REVIEW 4 major objections 6 minor 24 references
Implementation of tools for lessening the influence of artifacts in EEG signal analysis
T0 review · 4 major / 6 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read This paper provides ready-to-run scripts that reduce ocular artifacts in EEG by removing or attenuating independent components correlated with EOG signals.
desk verdict A useful open-source Fieldtrip implementation of EOG-artifact ICA methods that honestly demonstrates rather than validates itself, with the effectiveness claim resting on visual inspection and a tautological metric. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the artifact membership function (MSF), a logical vector of the same length as the data in which artifact samples are set to 1 and all others to 0. The paper provides a function that builds the MSF automatically using thresholded z-scores on the EOG channels after a 2-15 Hz band-pass filter, with a cut-off of 0.5 and 0.1 seconds of padding; it also uses the sum of absolute correlation coefficients between ICA components and the vertical and horizontal EOG channels to select which components to remove. The MSF makes the pipeline modular, because any artifact detector that outputs such a marker vector can be substituted.
What would settle it
Take an EEG recording with simultaneous video-based eye tracking or expert-annotated blink and saccade intervals, build the membership function with the provided defaults, and compare its marked samples against the annotations; if many real artifacts go unmarked or many non-artifact segments are marked, the partial-removal methods will either leave artifacts in or attenuate brain responses. A quantitative version would add simulated ocular components with known timing into clean EEG and check whether the partial-removal methods attenuate exactly those components.
Extended reading notes
Core claim
The paper's central claim is that its provided functions, built around an artifact membership function and EOG-correlation measures, successfully reduce ocular artifact influence in EEG channels. In the complete-removal variant, ICA components whose summed absolute correlation with the vertical and horizontal EOG signals is an outlier are discarded and the channels are reconstructed from the rest. In the two partial-removal variants, the same component selection is applied but only samples marked by the artifact membership function are zeroed, so non-artifact brain activity outside those windows is meant to be preserved; the third variant recomputes the ICA unmixing matrix on artifact-free samples before the partial removal. The paper shows example results on frontal channels, where large artifact peaks disappear and the summed correlation between independent components and EOG channels drops after processing.
Load-bearing premise
The automatic artifact detector that builds the membership function is assumed to mark exactly the eye-artifact samples, so that zeroing them removes ocular noise without also deleting genuine brain activity.
Editorial extensions
If this is right
- Researchers can apply the functions to their own EEG recordings and keep contaminated trials in the analysis instead of discarding them.
- Because the artifact membership function is a simple binary marker vector, alternative artifact detectors can be plugged in without changing the component selection or reconstruction code.
- In the two partial-removal methods, only samples identified as artifacts are altered, so genuine brain activity in clean segments should survive the cleaning.
- The complete-removal method offers a baseline that removes whole components and may discard some non-artifact signal, which the partial methods are designed to avoid.
Reading between the lines
- An implication the paper leaves implicit is that the default artifact-detection parameters are shown on one subject and one stimulus; testing them across many recordings would show whether they generalize or need per-recording tuning.
- A testable extension would replace the sum-of-absolute-correlation selection with component topographies or equivalent dipole locations, which often separate ocular sources more cleanly than raw correlation.
- The same membership-function machinery could be pointed at other artifact types with reference channels, such as muscle or cardiac artifacts, though the paper does not demonstrate this.
- Allowing the membership function to carry continuous weights instead of zeros and ones would enable graded attenuation, a possibility the authors mention but do not implement.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript describes and releases Matlab/Fieldtrip implementations of three methods, originally from [14], for reducing ocular artifacts in EEG: (1) complete removal of ICA components most correlated with EOG channels; (2) partial removal of those components only at samples flagged by an artifact membership function (MSF); and (3) partial removal using an unmixing matrix estimated from artifact-free data. It also provides a function to build the MSF using ft_artifact_zvalue, sample scripts (Boxes 1-3), and example results for one subject (001). The central claim is that these implementations lessen the influence of ocular artifacts and that the code and data allow exact replication.
Significance. The paper addresses a real need identified in [13]: lack of open-source artifact-removal implementations. Its strengths are the detailed function-by-function description, the release of data and code at a persistent handle, and the attempt to keep the MSF format simple and replaceable. If the methods were quantitatively validated, the contribution would be a useful reference implementation for practitioners. However, the current evidence for effectiveness is limited to visual inspection of two channels in one subject and a correlation metric that is partly self-referential; the paper therefore does not yet substantiate its central claim.
major comments (4)
- [Section 5, Figure 2] The only quantitative evidence is the reduction in the sum of absolute correlations between ICs and EOG before and after processing. Because the algorithm selects components precisely to minimize this quantity, the decrease is expected by construction and does not independently show that ocular artifacts in EEG channels are reduced or that brain signals are preserved. I recommend adding validation against expert-labeled artifact intervals, a downstream EEG analysis endpoint, or comparison with an established method such as ADJUST, MARA, or ASR, with quantitative metrics such as artifact attenuation, signal-to-error ratio, or preservation of known evoked responses.
- [Section 5, Method 1] The text states that in some areas not influenced by artifacts there is a decrease in amplitude with respect to the original shape of the signal. This indicates possible over-removal, but it is not quantified or compared across methods. Since the paper claims lessening artifact influence, it must report the trade-off between artifact energy removed and non-artifact signal distortion, for example by computing RMS error or spectral error on artifact-free segments, or by showing preservation of known neural responses.
- [Section 4.5] The automatic MSF created by fta_create_msf_fieldtrip uses ft_artifact_zvalue with a cut-off of 0.5, a 2-15 Hz band-pass, and padding of 0.1. These parameters are asserted without validation. Methods 2 and 3 depend entirely on this MSF, so the paper should report the sensitivity of the results to these choices and, ideally, compare the MSF against expert annotation or another artifact detector, such as those in [20-22], in terms of detection accuracy.
- [Section 5] All examples are from a single subject (001) and a single stimulus ('truck'). No across-subject statistics, error bars, or reproducibility checks are provided. For a software/tool paper, a small multi-subject demonstration with quantitative summary measures is necessary to support the claim that the implementations generally lessen artifact influence rather than merely working on one selected example.
minor comments (6)
- [Abstract] The phrase 'describes and implementation of scripts of code' should be corrected to 'describes an implementation of scripts of code' or 'describes implementations of scripts of code'.
- [Section 2.1] The units '10 KQ' should be '10 kΩ', and the same correction should be applied throughout the equipment description.
- [Section 4.1, step 4] The term 'modulus of correlation coefficients' should be replaced with 'absolute value' or 'magnitude' for clarity and consistency with standard statistical terminology.
- [Boxes 1-3] The example code contains typos and inconsistencies with the function signatures described in the text, such as 'clean.data' instead of 'cleaned_data', 'cfg.hpfilttype -' instead of 'cfg.hpfilttype =', and 'Corretion' instead of 'Correction'. These should be fixed so that the examples run exactly as printed.
- [Sections 4.1-4.9] Function names are rendered inconsistently, for example 'fta_select_compjrom_coeff_correlation' versus 'fta_select_compjrom_coeffcorrelation' and 'fta_channels_jrom_ica_reject_components' versus 'fta_channels_from_ica_reject_components'. The text should use exactly the names that appear in the released code to avoid confusion for readers attempting to replicate the pipeline.
- [Figure captions] The captions for Figures 1 and 3 use 'Fpl', but the standard and presumably intended label is 'Fp1'.
Circularity Check
Quantitative artifact-reduction evidence (Figure 2) is tautological: the plotted correlation sum is the algorithm's own selection criterion, so its decrease is guaranteed by construction.
-
fitted input called prediction
[Section 5, Method 1, Figure 2; selection rule in Section 4.1, steps 4–5]
"The reduction of the correlation coefficients between the independent components and EOG data after using the method is illustrated in Figure 2; specifically, Figure 2b shows a significant reduction of the correlation coefficients with respect to Figure 2a, after the application of the procedure."
The selection rule in Section 4.1 uses the same quantity as the reported outcome metric: 'the sum of the correlation coefficients is ordered from highest to lowest to show these components can be selected using a different criterion. These components are the ones selected for elimination.' Removing the components with the largest summed |correlation| with EOG guarantees that the post-removal sum in Figure 2b is smaller than the pre-removal sum in Figure 2a. The reported 'significant reduction' is therefore true by construction; it is the method's objective function renamed as a result, and it provides no independent evidence that ocular influence in reconstructed EEG channels is reduced or that non-artifact brain activity is preserved.
full rationale
The paper is primarily a software-implementation and reproducibility contribution, not a claim of a newly derived scientific law. The central deliverable—Matlab/Fieldtrip code, data, and examples—is self-contained and externally reproducible, so the availability claim is not circular. Reference [14] is a self-citation, but it is cited as the method being implemented rather than as a proof of efficacy; no load-bearing derivation rests on that citation. The circular element is concentrated in the quantitative support of Figure 2: the before/after correlation drop is the method's own selection objective, so it reduces by construction. Figures 1, 3, and 4 are single-subject visual illustrations with no independent artifact annotation, downstream-task metric, or comparison to a gold standard, and the automatically built MSF is not validated against expert labeling. Those are validation gaps and correctness risks rather than circular derivations. Overall, the artifact-reduction claim is partially supported by a tautological metric, giving a score of 6.
Assumptions & free parameters
free parameters (5)
- MSF z-value cut-off =
0.5
- MSF band-pass filter =
Butterworth order 3, 2-15 Hz
- artifact padding =
0.1 s
- trial padding =
0 s
- IC outlier selection criterion
assumptions (4)
- domain assumption EEG channels are a linear mixture of statistically independent sources (ICA model), so unmixing separates ocular artifacts into distinct components.
- domain assumption Independent components with high absolute correlation with VEOG/HEOG are ocular artifact components.
- domain assumption The artifact membership function produced by ft_artifact_zvalue with the selected parameters correctly identifies ocular artifact samples.
- domain assumption An ICA unmixing matrix estimated from artifact-free samples is valid for the entire recording.
Cite this review
Pith. "Pith review of Implementation of tools for lessening the influence of artifacts in EEG signal analysis." pith.science (2026). https://pith.science/paper/BKUTV55O
@misc{pith2026241113931,
author = {Pith},
title = {Pith review of: Implementation of tools for lessening the influence of artifacts in EEG signal analysis},
year = {2026},
howpublished = {\url{https://pith.science/paper/BKUTV55O}},
note = {Machine review of arXiv:2411.13931}
}
read the original abstract
This manuscript describes and implementation of scripts of code aimed at reducing the influence of artifacts, specifically focused on ocular artifacts, in the measurement and processing of electroencephalogram (EEG) signals. This process is of importance because it benefits the analysis and study of long trial samples when the appearance of ocular artifacts cannot be avoided by simply discarding trials. The implementations provided to the reader illustrate, with slight modifications, previously proposed methods aimed at the partial or complete elimination of EEG channels or components are those that resemble the electro-oculogram (EOG) signals in which artifacts are detected. In addition to the description of each of the provided functions, examples of utilization and illustrative figures will be included to show the expected results and processing pipeline.
Reference graph
Works this paper leans on
-
[14]
Preprocessing for Lessening the Influence of Eye Artifacts in EEG Analysis
Villena, A.; Tardón, L.J.; Barbancho, I.; Barbancho, A.M.; Brattico, E.; Haumann, N.T. Preprocessing for Lessening the Influence of Eye Artifacts in EEG Analysis. Appl. Sci. 2019, 9,1757. [CrossRef]
work page 2019
-
[13]
Review of Challenges Associated with the EEG Artifact Removal Methods
Mumtaz, W; Rasheed, S.; Irfan, A. Review of Challenges Associated with the EEG Artifact Removal Methods. Biomed. Signal Process. Control 2021, 68,102741. [CrossRef]
work page 2021
-
[1]
Verleger, R.; Gasser, T.; Mocks, J. Correction of EOG Artifacts in Event-Related Potentials of the EEG: Aspects of Reliability and Validity. Psychophysiology 1982,19,472-480. [CrossRef] [PubMed]
work page 1982
-
[2]
The Removal of the Eye-Movement Artifact from the EEG by Regression Analysis in the Frequency Domain
Woestenburg, J.; Verbaten, M.; Slangen, J. The Removal of the Eye-Movement Artifact from the EEG by Regression Analysis in the Frequency Domain. Biol. Psychol. 1983,16,127-147. [CrossRef] [PubMed]
work page 1983
-
[3]
Adaptive AR Modeling of Nonstationary Time Series by Means of Kalman Filtering
Arnold, M.; Milner, X.H.R.; Witte, H.; Bauer, R.; Braun, C. Adaptive AR Modeling of Nonstationary Time Series by Means of Kalman Filtering. IEEE Trans. Biomed. Eng. 1998,45,553-562. [CrossRef] [PubMed]
work page 1998
-
[4]
Automatic Removal of Eye Movement and Blink Artifacts from EEG Data using Blind Component Separation
Joyce, C.; Gorodnitsky, L; Kutas, M. Automatic Removal of Eye Movement and Blink Artifacts from EEG Data using Blind Component Separation. Psychophysiology 2004, 41,313-325. [CrossRef] [PubMed]
work page 2004
-
[5]
Removal of Artifacts from EEG Signals: A Review
Jiang, X.; Bian, G.B.; Tian, Z. Removal of Artifacts from EEG Signals: A Review. Sensors 2019,19, 987. [CrossRef] [PubMed]
work page 2019
-
[6]
Isolation of Multiple Electrocardiogram Artifacts using Independent Vector Analysis
Uddin, Z.; Altaf, M.; Ahmad, A.; Qamar, A.; Orakzai, F. Isolation of Multiple Electrocardiogram Artifacts using Independent Vector Analysis. Peer] Comput. Sci. 2023, 9, ell89. [CrossRef] [PubMed]
work page 2023
Show all 24 references
-
[7]
Taxonomy on EEG Artifacts Removal Methods, Issues, and Healthcare Applications
Shukla, P.K.; Roy, V; Shukla, P.K.; Gupta, A.K.; Goel, V; Shukla, S. Taxonomy on EEG Artifacts Removal Methods, Issues, and Healthcare Applications. J. Organ. End User Comput. 2021,33,19-46. [CrossRef] Appl. Sci. 2024,14, 971 14 of 14
2021
-
[8]
Removal of Movement-Induced EEG Artifacts: Current State of the Art and Guidelines
Gorjan, D.; Gramann, K.; Pauw, K.D.; Marusic, U. Removal of Movement-Induced EEG Artifacts: Current State of the Art and Guidelines. J. Neural Eng. 2022, 19, 011004. [CrossRef] [PubMed]
2022
-
[9]
Automatic Detection and Classification of EEG Artifacts using Fuzzy Kernel SVM and Wavelet ICA (WICA)
Yasoda, K.; Ponmagal, R.S.; Bhuvaneshwari, K.S.; Venkatachalam4, K. Automatic Detection and Classification of EEG Artifacts using Fuzzy Kernel SVM and Wavelet ICA (WICA). Soft Comput. 2020, 24,16011-16019. [CrossRef]
2020
-
[10]
Automatic and Tunable Algorithm for EEG Artifact Removal using Wavelet Decomposition with Applications in Predictive Modeling during Auditory Tasks
Bajaj, N.; Requena-Carrión, J.; Bellotti, E; Berta, R.; De-Gloria, A. Automatic and Tunable Algorithm for EEG Artifact Removal using Wavelet Decomposition with Applications in Predictive Modeling during Auditory Tasks. Biomed. Signal Process. Control 2020, 55,101624. [CrossRef]
2020
-
[11]
Noorbasha, S.K.; Sudha, G.F. Removal of EOG Artifacts and Separation of Different Cerebral Activity Components from Single Channel EEG—An Efficient Approach Combining SSA-ICA with Wavelet Thresholding for BCI Applications. Biomed. Signal Process. Control 2021, 63,102168. [CrossRef]
2021
-
[12]
Adaptive Single-Channel EEG Artifact Removal With Applications to Clinical Monitoring
Dora, M.; Holcman, D. Adaptive Single-Channel EEG Artifact Removal With Applications to Clinical Monitoring. IEEE Trans. Neural. Syst. Rehabil. Eng. 2022,30, 286-295. [CrossRef] [PubMed]
2022
-
[15]
E-Prime; Psychology Software Tools: Pittsburgh, PA, USA, 2022
Psychology Software Tools. E-Prime; Psychology Software Tools: Pittsburgh, PA, USA, 2022
2022
-
[16]
FieldTrip: Open Source Software for Advanced Analysis of MEG, EEG, and Invasive Electrophysiological Data
Oostenveld, R.; Fries, E; Maris, E.; Schoffelen, J.M. FieldTrip: Open Source Software for Advanced Analysis of MEG, EEG, and Invasive Electrophysiological Data. Comput. Intell. Neurosci. 2011,2011,1-9. [CrossRef] [PubMed]
2011
-
[17]
Fieldtrip Matlab Software Toolbox for MEG, EEG and iEEG Analysis; Donders Institute for Brain, Cognition and Behaviour: Nijmegen, The Netherlands, 2024
Donders Institute for Brain, Cognition and Behaviour. Fieldtrip Matlab Software Toolbox for MEG, EEG and iEEG Analysis; Donders Institute for Brain, Cognition and Behaviour: Nijmegen, The Netherlands, 2024
2024
-
[18]
MATLAB; version: 9.13.0 (R2022b); MathWorks Inc.: Natick, MA, USA, 2022
The MathWorks Inc. MATLAB; version: 9.13.0 (R2022b); MathWorks Inc.: Natick, MA, USA, 2022
2022
-
[19]
An Improved Algorithm for Automatic Detection of Saccades in Eye Movement Data and for Calculating Saccade Parameters
Behrens, E; MacKeben, M.; Schróder-Preikschat, W. An Improved Algorithm for Automatic Detection of Saccades in Eye Movement Data and for Calculating Saccade Parameters. Behav. Res. Methods 2010, 42, 701-708. [CrossRef] [PubMed]
2010
-
[20]
An Algorithm Separating Saccadic from Nonsaccadic Eye Movements Automatically by use of the Acceleration Signal
Behrens, E; Weiss, L.R. An Algorithm Separating Saccadic from Nonsaccadic Eye Movements Automatically by use of the Acceleration Signal. Vis. Res. 1992, 32,889-893. [CrossRef] [PubMed]
1992
-
[21]
A Probabilistic Real-Time Algorithm for Detecting Blinks, Saccades, and Fixations from EOG Data
Toivanen, M.; Pettersson, K.; Lukander, K. A Probabilistic Real-Time Algorithm for Detecting Blinks, Saccades, and Fixations from EOG Data. J. Eye Mov. Res. 2015, 8,1-14. [CrossRef]
2015
-
[22]
Feature Extraction from Regenerated EEG - A Better Approach for ICA Based Eye Blink Artifact Detection
Rashida, M.; Habib, M.A. Feature Extraction from Regenerated EEG - A Better Approach for ICA Based Eye Blink Artifact Detection. Front. Biomed. Technol. 2023,11,569-576. Available online: , (accessed on 3 December 2023). https://fbt.tums.ac.ir/index.php/fbt/article/view/569
2023
-
[23]
Optimizing the ICA-based Removal of Ocular EEG Artifacts from Free Viewing Experiments
Dimigen, O. Optimizing the ICA-based Removal of Ocular EEG Artifacts from Free Viewing Experiments. NeuroImage 2020, 207,116117. [CrossRef] [PubMed]
2020
-
[24]
EEG Artifact Removal by Bayesian Deep Learning & ICA
Lee, S.S.; Lee, K.; Kang, G. EEG Artifact Removal by Bayesian Deep Learning & ICA. In Proceedings of the 2020 42nd Annual International Conference of the IEEE Engineering in Medicine Biology Society (EMBC), Montreal, QC, Canada, 20-24 July 2020; pp. 932-935. [CrossRef] Disclai...
2020
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.