Name: PIPEMAT-RS: Development and Validation of a Standardized MATLAB Pipeline for Resting-State EEG Preprocessing (Video) | JoVE Journal
Uploaded: 2025-06-06T15:55:00.000+00:00
Duration: 6 min 51 s
Description: Learn PIPEMAT-RS for standardized resting-state EEG preprocessing using MATLAB and ICA-based artifact removal. This pipeline improves signal quality and supports reproducible neurophysiological research across clinical studies.

ai_agent_IP

ai_agent

maintenance

instantTrialDurationMinutes

isBiopharma

enableNewMyJoveSectionDesign

educationFeatures

seoMetaApi

ai_copilot_institution_blacklist

ai_chat_institution_blacklist

freeAccessWindow

ssoRoleCapture

scheduledMaintenanceBanner

joveLabsFeatures

disable_autoplay_for_bots

extended_trial_duration

5 minutes

Specifically, LMM was supported by a postdoctoral research grant #21/05897-5, S&#227;o Paulo Research Foundation (FAPESP). AC was supported by a scientific initiation grant #21/12790-2, S&#227;o Paulo Research Foundation (FAPESP). SPB was supported by a postdoctoral research grant #20/08512-4, S&#227;o Paulo Research Foundation (FAPESP). FF and LRB are supported by research grant #17/12943-8, S&#227;o Paulo Research Foundation (FAPESP). Individually, FF received support from NHI 2020 R01 AT, Project #1R01AT009491-01A1.

PIPEMAT-RS was developed as part of our ongoing research project, approved by the Ethics Committee of the Clinics Hospital, University of S&#227;o Paulo Medical School (CAAE: 86832518.7.0000.0068). This project aims to enhance the quality and reliability of resting-state EEG data for various neurological and psychiatric studies. The data used in this study are derived from a cohort study, the protocol of which was previously published by Simis et al. (2021)31. The main study is still in progress, with data collection ongoing for the remaining clinical groups. The script for the pipeline is provided in Supplementary File 1.
1. File format conversion
<ol>
	<li>Prepare the data directory.
	<ol>
		<li>Open MATLAB and navigate to the directory containing raw EEG data files.</li>
		<li>Define the path using pathname = &#39;C:\\Path\\To\\RawData\\&#39;; and list subjects with name_subj_raw = {&#39;subject1&#39;, &#39;subject2&#39;, ...};.</li>
		<li>If the dataset is already in .mat format, skip this step and proceed to EEG Montage.</li>
	</ol>
	</li>
	<li>Convert files.
	<ol>
		<li>Load each raw data file using the appropriate MATLAB function for the file type. For example, use load for .mat files or file-specific EEGLAB functions such as pop_biosig (for .edf, .bdf) or pop_fileio (for various proprietary formats). 
		NOTE: Users should adapt this step based on the format of their raw EEG files and refer to EEGLAB&#39;s documentation for additional file support.</li>
		<li>Save each file in .mat format, ensuring data integrity and retaining all relevant metadata.</li>
	</ol>
	</li>
</ol>
2. EEG montage
<ol>
	<li>Load EEG data.
	<ol>
		<li>Open MATLAB and navigate to the directory with the .mat files.</li>
		<li>Use the pop_loadset function from the EEGLAB toolbox to load each subject&#39;s EEG dataset (e.g., EEG = pop_loadset(&#39;filename.set&#39;);).</li>
	</ol>
	</li>
	<li>Apply electrode locations.
	<ol>
		<li>Ensure the appropriate electrode location file (e.g., EEGSystem_128.loc) corresponding to the EEG cap configuration is available.</li>
		<li>Use the pop_chanedit function to apply electrode locations. For example: 
		EEG = pop_chanedit(EEG, &#39;load&#39;, {&#39;EEGSystem_128.loc&#39;, &#39;filetype&#39;, &#39;loc&#39;});</li>
		<li>Visually confirm that the electrode positions were loaded correctly by reviewing the channel layout using EEGLAB&#39;s GUI.</li>
	</ol>
	</li>
	<li>Save the dataset.
	<ol>
		<li>Save the updated EEG dataset using the pop_saveset function (e.g., pop_saveset(EEG, &#39;filename&#39;, &#39;subject1_loc.set&#39;)).</li>
	</ol>
	</li>
</ol>
3. Downsampling and filtering
<ol>
	<li>Downsample the EEG Data.
	<ol>
		<li>Reduce the sampling rate of the EEG data to 250 Hz using the pop_resample function (optional). 
		NOTE: The pop_resample function automatically applies an appropriate low-pass anti-aliasing filter before downsampling. This step ensures that high-frequency components are attenuated prior to resampling, preventing aliasing artifacts and preserving the quality of the EEG signal.</li>
	</ol>
	</li>
	<li>Apply band-pass filtering.
	<ol>
		<li>Apply a band-pass filter from 1 Hz to 50 Hz to channels 1-64 using the pop_basicfilter function with a Butterworth design and filter order of 2.</li>
		<li>Use zero-phase forward and reverse filtering to prevent phase distortion.</li>
	</ol>
	</li>
	<li>Apply notch filtering (Optional).
	<ol>
		<li>Apply a notch filter at 50 or 60 Hz to channels 1-64 using pop_basicfilter with the &#39;PMnotch&#39; design and filter order of 180. 
		NOTE: This step is optional and should be used when analyzing oscillatory activity above 30 Hz (e.g., gamma band). In such cases, researchers may increase the upper limit of the band-pass filter (e.g., up to 80 or 100 Hz), making the notch filter crucial for eliminating powerline interference. The command is provided in the pipeline for convenience and can be easily adjusted depending on the local electrical frequency (e.g., 50 Hz or 60 Hz).</li>
	</ol>
	</li>
</ol>
4. Artifact rejection and rereferencing
<ol>
	<li>Perform automatic artifact rejection.
	<ol>
		<li>Use the clean_rawdata function with parameters (5, [-1], 0.7, -1, 8, 0.85) to automatically detect and remove flatline channels, noisy segments, and low-frequency drifts. These parameters correspond respectively to: 
		FlatlineCriterion (5): removes channels flat for more than 5 s. 
		ChannelCriterion ([-1]): disables channel correlation-based rejection. 
		LineNoiseCriterion (0.7): removes channels with excessive line noise above this threshold. 
		Highpass (&#8722;1): disables additional high-pass filtering during cleaning. 
		BurstCriterion (8): removes data bursts that exceed 8 standard deviations. 
		WindowCriterion (0.85): requires 85% of data in a window to be clean to retain it. 
		NOTE: Users may adjust these values depending on dataset characteristics. For instance, reducing the BurstCriterion will result in more aggressive artifact rejection, while lowering the WindowCriterion will increase the tolerance for noisy data windows.</li>
	</ol>
	</li>
	<li>Perform manual artifact inspection (optional).
	<ol>
		<li>If necessary, visually inspect the EEG data using eegplot from EEGLAB to identify residual artifacts not detected by automatic methods.</li>
		<li>Manually mark and remove remaining artifacts to ensure high data quality. 
		NOTE: This step is optional and primarily recommended when dealing with datasets that present non-standard noise patterns or when greater precision is needed beyond automated procedures.</li>
	</ol>
	</li>
	<li>Rereference the data.
	<ol>
		<li>Use the pop_reref function to rereference EEG signals to the average of all electrodes, ensuring that the original reference electrode is retained in the dataset. 
		NOTE: In this pipeline, average re-referencing is performed using the available channels without reconstructing the original reference electrode. However, users aiming for enhanced rank consistency and cross-subject comparability may consider temporarily reinserting the original reference electrode (e.g., Cz) as a zero-valued time series before re-referencing. This allows the reference to contribute to the average calculation. Afterward, the added channel can be removed to avoid rank deficiency before performing ICA. This process is exemplified in the RELAX_average_rereference.m function of the RELAX pipeline, which also includes interpolation of all channels prior to referencing. While PIPEMAT-RS maintains a simplified version for clarity and accessibility, users are encouraged to adapt this step as needed depending on the goals of their analysis.</li>
	</ol>
	</li>
</ol>
5. Independent component analysis (ICA)
<ol>
	<li>Load preprocessed EEG data.
	<ol>
		<li>Load the file with the suffix _loc_filt_cleanraw_reref.set using the pop_loadset function.</li>
	</ol>
	</li>
	<li>Perform ICA.
	<ol>
		<li>Run ICA using the pop_runica function with the runica algorithm to decompose the data into independent components.The &#39;runica&#39; algorithm, based on the Infomax approach, maximizes the statistical independence of the components and is widely used in EEG preprocessing. 
		NOTE: For users seeking faster ICA computation, alternative algorithms such as picard may be used by specifying icatype, picard in the function call. These alternatives offer comparable performance and may reduce computation time.</li>
	</ol>
	</li>
	<li>Save the ICA dataset.
	<ol>
		<li>Save the dataset with the suffix _loc_filt_cleanraw_reref_ICA using pop_saveset to indicate that ICA has been performed.</li>
	</ol>
	</li>
</ol>
6. ICLabel classification
<ol>
	<li>Load the ICA dataset.
	<ol>
		<li>Import the dataset with ICA components using the pop_loadset function.</li>
	</ol>
	</li>
	<li>Apply ICLabel classification.
	<ol>
		<li>Run the pop_iclabel function with the default model to classify independent components into brain and artifact categories.</li>
		<li>Identify components with a probability greater than 0.7 of representing brain activity.</li>
	</ol>
	</li>
	<li>Remove artifact components.
	<ol>
		<li>Use the pop_subcomp function to remove all components with ICLabel brain probability below 0.7.</li>
		<li>Retain only components with a brain activity probability above 0.7. The threshold of 0.7 balances the need to preserve genuine neural signals while ensuring the effective removal of artifacts.</li>
	</ol>
	</li>
	<li>Save the cleaned dataset.
	<ol>
		<li>Save the cleaned dataset with the suffix _loc_filt_cleanraw_reref_ICA_Cleaned using pop_saveset to indicate the removal of artifact components.</li>
	</ol>
	</li>
</ol>
7. Data normalization
<ol>
	<li>Load cleaned EEG data.
	<ol>
		<li>Import the dataset with artifact-free ICA components using the pop_loadset function.</li>
	</ol>
	</li>
	<li>Apply Z-Score normalization.
	<ol>
		<li>Normalize each EEG channel using Z-score transformation by computing (x - mean(x)) / std(x) for each channel individually. 
		NOTE: Z-score normalization is optional and should be applied based on the analytical goals of the study. While it reduces inter-subject variability and standardizes data for group-level comparisons, it may also attenuate topographical distinctions and eliminate meaningful amplitude differences across channels. Researchers interested in preserving raw signal magnitude or spatial distribution are advised to skip this step.</li>
	</ol>
	</li>
	<li>Save the final dataset.
	<ol>
		<li>Save the fully preprocessed and normalized dataset with the suffix _loc_filt_cleanraw_reref_ICA_Normalized using pop_saveset to indicate the completion of preprocessing.</li>
	</ol>
	</li>
</ol>

<ol><li>Niedermeyer, E., da Silva, F. L. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aNiedermeyer%2C+E+%22Electroencephalography%3A+basic+principles%2C+clinical+applications%2C+and+related+fields%22">Electroencephalography: basic principles, clinical applications, and related fields</a>. , Lippincott Williams and Wilkins. Philadelphia, PA. (2004).</li>
<li>Delorme, A., Makeig, S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((EEGLAB%3A+an+open-source+toolbox+for+analysis+of+single-trial+EEG+dynamics+including+independent+component+analysis%5BTitle%5D)+AND+%22J+Neurosci+Methods%22%5BJournal%5D)">EEGLAB: an open-source toolbox for analysis of single-trial EEG dynamics including independent component analysis.</a> J Neurosci Methods. 134 (1), 9-21 (2004).</li>
<li>Cohen, M. X. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aCohen%2C+MX+%22Analyzing+neural+time+series+data%3A+theory+and+practice%22">Analyzing neural time series data: theory and practice</a>. , MIT Press. (2014).</li>
<li>Buzsáki, G., Anastassiou, C. A., Koch, C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+origin+of+extracellular+fields+and+currents-EEG%2C+ECoG%2C+LFP+and+spikes%5BTitle%5D)+AND+%22Nat+Rev+Neurosci%22%5BJournal%5D)">The origin of extracellular fields and currents-EEG, ECoG, LFP and spikes.</a> Nat Rev Neurosci. 13 (6), 407-420 (2012).</li>
<li>Nunez, P. L., Srinivasan, R. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aNunez%2C+PL+%22Electric+fields+of+the+brain%3A+the+neurophysics+of+EEG%22">Electric fields of the brain: the neurophysics of EEG</a>. , Oxford University Press. (2006).</li>
<li>Gratton, G. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Brain+reflections%3A+a+circuit-based+framework+for+understanding+information+processing+and+cognitive+control%5BTitle%5D)+AND+%22Psychophysiology%22%5BJournal%5D)">Brain reflections: a circuit-based framework for understanding information processing and cognitive control.</a> Psychophysiology. 55 (3), e13038(2018).</li>
<li>Bigdely-Shamlo, N., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+PREP+pipeline%3A+standardized+preprocessing+for+large-scale+EEG+analysis%5BTitle%5D)+AND+%22Front+Neuroinform%22%5BJournal%5D)">The PREP pipeline: standardized preprocessing for large-scale EEG analysis.</a> Front Neuroinform. 9, 16(2015).</li>
<li>Pernet, C. R., Chauveau, N., Gaspar, C., Rousselet, G. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((LIMO+EEG%3A+a+toolbox+for+hierarchical+linear+modeling+of+electroencephalographic+data%5BTitle%5D)+AND+%22Comput+Intell+Neurosci%22%5BJournal%5D)">LIMO EEG: a toolbox for hierarchical linear modeling of electroencephalographic data.</a> Comput Intell Neurosci. 2011, 831409(2011).</li>
<li>Onton, J., Makeig, S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Independent+component+analysis+for+EEG+artifacts%5BTitle%5D)+AND+%22Neuroimage%22%5BJournal%5D)">Independent component analysis for EEG artifacts.</a> Neuroimage. 15 (3), 789-801 (2006).</li>
<li>Delorme, A., Miyakoshi, M., Jung, T. P., Makeig, S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Grand+average+ERP-image+plotting+and+statistics%3A+a+method+for+comparing+variability+in+event-related+single-trial+EEG+activities+across+subjects+and+conditions%5BTitle%5D)+AND+%22J+Neurosci+Methods%22%5BJournal%5D)">Grand average ERP-image plotting and statistics: a method for comparing variability in event-related single-trial EEG activities across subjects and conditions.</a> J Neurosci Methods. 210 (2), 132-145 (2012).</li>
<li>Debener, S., Thorne, J., Schneider, T. R., Viola, F. C. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aDebener%2C+S+%22Using+ICA+for+the+analysis+of+multi-channel+EEG+data.+Simultaneous+EEG+and+FMRI%22">Using ICA for the analysis of multi-channel EEG data. Simultaneous EEG and FMRI</a>. , Springer. New York, NY. (2010).</li>
<li>Luck, S. J. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aLuck%2C+SJ+%22An+introduction+to+the+event-related+potential+technique%22">An introduction to the event-related potential technique</a>. , MIT Press. (2014).</li>
<li>Cohen, M. X. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Where+does+EEG+come+from+and+what+does+it+mean%5BTitle%5D)+AND+%22Trends+Neurosci%22%5BJournal%5D)">Where does EEG come from and what does it mean.</a> Trends Neurosci. 40 (4), 208-218 (2017).</li>
<li>Keil, A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Committee+report%3A+publication+guidelines+and+recommendations+for+studies+using+electroencephalography+and+magnetoencephalography%5BTitle%5D)+AND+%22Psychophysiology%22%5BJournal%5D)">Committee report: publication guidelines and recommendations for studies using electroencephalography and magnetoencephalography.</a> Psychophysiology. 51 (1), 1-21 (2014).</li>
<li>Gramfort, A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((MEG+and+EEG+data+analysis+with+MNE-Python%5BTitle%5D)+AND+%22Front+Neurosci%22%5BJournal%5D)">MEG and EEG data analysis with MNE-Python.</a> Front Neurosci. 7, 7(2013).</li>
<li>Pernet, C., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Issues+and+recommendations+from+the+OHBM+COBIDAS+MEEG+committee+for+reproducible+EEG+and+MEG+research%5BTitle%5D)+AND+%22Nat+Neurosci%22%5BJournal%5D)">Issues and recommendations from the OHBM COBIDAS MEEG committee for reproducible EEG and MEG research.</a> Nat Neurosci. 23 (12), 1473-1483 (2020).</li>
<li>Oostenveld, R., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((FieldTrip%3A+open-source+software+for+advanced+analysis+of+MEG%2C+EEG%2C+and+invasive+electrophysiological+data%5BTitle%5D)+AND+%22Comput+Intell+Neurosci%22%5BJournal%5D)">FieldTrip: open-source software for advanced analysis of MEG, EEG, and invasive electrophysiological data.</a> Comput Intell Neurosci. 2011, 156869(2011).</li>
<li>Bailey, N. W., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Introducing+RELAX%3A+an+automated+preprocessing+pipeline+for+cleaning+EEG+data-Part+1%3A+algorithm+and+application+to+oscillations%5BTitle%5D)+AND+%22Clin+Neurophysiol%22%5BJournal%5D)">Introducing RELAX: an automated preprocessing pipeline for cleaning EEG data-Part 1: algorithm and application to oscillations.</a> Clin Neurophysiol. 149, 178-201 (2023).</li>
<li>Pedroni, A., Bahreini, A., Langer, N. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Automagic%3A+standardized+preprocessing+of+big+EEG+data%5BTitle%5D)+AND+%22Neuroimage%22%5BJournal%5D)">Automagic: standardized preprocessing of big EEG data.</a> Neuroimage. 200, 460-473 (2019).</li>
<li>da Cruz, J. R., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((An+automatic+preprocessing+pipeline+for+EEG+analysis+%28APP%29+based+on+robust+statistics%5BTitle%5D)+AND+%22Clin+Neurophysiol%22%5BJournal%5D)">An automatic preprocessing pipeline for EEG analysis (APP) based on robust statistics.</a> Clin Neurophysiol. 129 (7), 1427-1437 (2018).</li>
<li>Bailey, N. W., et al. EEG is better when cleaning effectively targets artifacts. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aBailey%2C+NW+%22bioRxiv%22">bioRxiv</a>. , 2024-2106 (2024).</li>
<li>Global EEG devices market report. , Grand View Research. Available at: 
 <a target="_blank" href="https://www.grandviewresearch.com">https://www.grandviewresearch.com</a> (2021).</li>
<li>Makoto's preprocessing pipeline. , Available at: 
 <a target="_blank" href="https://sccn.ucsd.edu/wiki/Makoto's_preprocessing_pipeline">https://sccn.ucsd.edu/wiki/Makoto's_preprocessing_pipeline</a> (2015).</li>
<li>Winkler, I., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((On+the+influence+of+high-pass+filtering+on+ICA-based+artifact+reduction+in+EEG-ERP%5BTitle%5D)+AND+%22Annu+Int+Conf+IEEE+Eng+Med+Biol+Soc%22%5BJournal%5D)">On the influence of high-pass filtering on ICA-based artifact reduction in EEG-ERP.</a> Annu Int Conf IEEE Eng Med Biol Soc. 2015, 4101-4105 (2015).</li>
<li>Gramfort, A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Time-frequency+mixed-norm+estimates%3A+sparse+M%2FEEG+imaging+with+non-stationary+source+activations%5BTitle%5D)+AND+%22Neuroimage%22%5BJournal%5D)">Time-frequency mixed-norm estimates: sparse M/EEG imaging with non-stationary source activations.</a> Neuroimage. 70, 410-422 (2013).</li>
<li>Delorme, A., Sejnowski, T., Makeig, S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Enhanced+detection+of+artifacts+in+EEG+data+using+higher-order+statistics+and+independent+component+analysis%5BTitle%5D)+AND+%22Neuroimage%22%5BJournal%5D)">Enhanced detection of artifacts in EEG data using higher-order statistics and independent component analysis.</a> Neuroimage. 34 (4), 1443-1449 (2007).</li>
<li>Marques, L. M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Resting-state+EEG+as+biomarker+of+maladaptive+motor+function+and+depressive+profile+in+stroke+patients%5BTitle%5D)+AND+%22Clin+EEG+Neurosci%22%5BJournal%5D)">Resting-state EEG as biomarker of maladaptive motor function and depressive profile in stroke patients.</a> Clin EEG Neurosci. 55 (4), 496-507 (2024).</li>
<li>Lacerda, G. J., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+neural+signature+for+brain+compensation+in+stroke+with+EEG+and+TMS%3A+insights+from+the+DEFINE+cohort+study%5BTitle%5D)+AND+%22Neurophysiol+Clin%22%5BJournal%5D)">A neural signature for brain compensation in stroke with EEG and TMS: insights from the DEFINE cohort study.</a> Neurophysiol Clin. 54 (5), 102985(2024).</li>
<li>Camargo, L., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Adaptive+and+compensatory+neural+signatures+in+fibromyalgia%3A+an+analysis+of+resting-state+and+stimulus-evoked+EEG+oscillations%5BTitle%5D)+AND+%22Biomedicines%22%5BJournal%5D)">Adaptive and compensatory neural signatures in fibromyalgia: an analysis of resting-state and stimulus-evoked EEG oscillations.</a> Biomedicines. 12 (7), 1428(2024).</li>
<li>Barbosa, S. P., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Resting-state+electroencephalography+delta+and+theta+bands+as+compensatory+oscillations+in+chronic+neuropathic+pain%3A+a+secondary+data+analysis%5BTitle%5D)+AND+%22Brain+Netw+Modul%22%5BJournal%5D)">Resting-state electroencephalography delta and theta bands as compensatory oscillations in chronic neuropathic pain: a secondary data analysis.</a> Brain Netw Modul. 3 (2), 52-60 (2024).</li>
<li>Simis, M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Deficit+of+inhibition+as+a+marker+of+neuroplasticity+%28DEFINE+study%29+in+rehabilitation%3A+a+longitudinal+cohort+study+protocol%5BTitle%5D)+AND+%22Front+Neurol%22%5BJournal%5D)">Deficit of inhibition as a marker of neuroplasticity (DEFINE study) in rehabilitation: a longitudinal cohort study protocol.</a> Front Neurol. 12, 695406(2021).</li>
<li>Pion-Tonachini, L., Kreutz-Delgado, K., Makeig, S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((ICLabel%3A+an+automated+electroencephalographic+independent+component+classifier%2C+dataset%2C+and+website%5BTitle%5D)+AND+%22Neuroimage%22%5BJournal%5D)">ICLabel: an automated electroencephalographic independent component classifier, dataset, and website.</a> Neuroimage. 198, 181-197 (2019).</li>
<li>Jaramillo-Jimenez, A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Reproducible+python+workflow+for+multi-site+resting-state+EEG+analysis%5BTitle%5D)+AND+%22Alzheimer%27s+Dement%22%5BJournal%5D)">Reproducible python workflow for multi-site resting-state EEG analysis.</a> Alzheimer's Dement. 19, e076353(2023).</li>
<li>Levin, A. R., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((BEAPP%3A+the+batch+electroencephalography+automated+processing+platform%5BTitle%5D)+AND+%22Front+Neurosci%22%5BJournal%5D)">BEAPP: the batch electroencephalography automated processing platform.</a> Front Neurosci. 12, 513(2018).</li>
</ol>

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

The PIPEMAT-RS pipeline was developed to provide a standardized, efficient method for preprocessing resting-state EEG data. Critical steps in this protocol include artifact rejection and Independent Component Analysis (ICA), both of which significantly enhance the signal-to-noise ratio and ensure the extraction of meaningful neural signals. The combination of automatic artifact rejection using the clean_rawdata function and manual inspection via eegplot ensures comprehensive artifact management, balancing efficiency and ...

The electroencephalogram (EEG) is a critical tool in both neuroscience research and clinical practice, offering valuable insights into brain electrical activity. Resting-state EEG data, in particular, are extensively utilized to investigate various neurological and psychiatric conditions, such as stroke, fibromyalgia, and chronic neuropathic pain. Despite its widespread application, the analysis of EEG data necessitates rigorous preprocessing to eliminate artifacts and noise, ensuring the integrity and reliability of the results1,2,3. Over the years, EEG has been invaluable in providing non-invasive and high-temporal resolution measurements of brain activity, making it indispensable for studying the dynamic processes of the brain4,5. However, the raw EEG signals are often contaminated by a variety of noise sources that can obscure the neural signals of interest, thereby complicating the interpretation of the data6.
EEG signals are highly susceptible to contamination from sources such as muscle activity, eye movements, and electrical interference7. These artifacts can obscure neural signals, making preprocessing a critical step in EEG analysis. This process typically involves multiple stages-file format conversion, downsampling, filtering, artifact rejection, independent component analysis (ICA), and exclusion of noise-related components8,9. Each step contributes to improving signal quality and ensuring that the data accurately reflect neural activity10. The complexity of these processes necessitates the use of sophisticated algorithms and software tools to ensure that the preprocessing is both efficient and effective11.
One of the primary challenges in EEG data preprocessing is the variability in the data, which can differ significantly between subjects and recording sessions10,12. This variability arises due to differences in the physiological and anatomical characteristics of the subjects, as well as variations in the experimental conditions and equipment used. Additionally, the lack of standardized preprocessing protocols can lead to inconsistencies in data analysis and interpretation13,14. While several pipelines and preprocessing scripts are available, many of them are either tailored to specific datasets or lack comprehensive documentation, making them less accessible to the broader research community14. Platforms like GitHub host numerous EEG preprocessing scripts and plugins, facilitating the sharing and collaborative improvement of these tools15. However, the fragmented nature of these resources underscores the need for a robust, efficient, and standardized editable preprocessing pipeline that can be broadly applied across different datasets and research contexts16. Such a pipeline would not only enhance the reproducibility of EEG studies but also provide a foundation upon which the research community can build and adapt for specific needs17.
In addition to EEGLAB, which forms the basis of PIPEMAT-RS due to its extensive user base and integration with MATLAB, other EEG processing toolboxes are also widely used in the field. For instance, MNE-Python is an open-source Python-based platform that offers advanced functionalities for source localization, sensor-space analyses, and integration with machine learning workflows. FieldTrip, a MATLAB-based toolbox like EEGLAB, is known for its flexible scripting environment and comprehensive support for time-frequency and connectivity analyses. While these platforms provide powerful alternatives, they typically require steeper learning curves and more advanced programming skills. PIPEMAT-RS was designed to bridge accessibility and transparency, especially for researchers seeking an editable and modular solution for resting-state EEG preprocessing using EEGLAB. Nonetheless, the conceptual structure of PIPEMAT-RS can be adapted to other toolboxes as needed, fostering interoperability and customization across different software environments.
Several automated EEG preprocessing pipelines have been developed in recent years, such as RELAX18, Automagic19, APP20, and PREP7. These pipelines incorporate advanced algorithms for artifact removal and data cleaning, often employing ICA-based strategies with limited transparency or customization for non-expert users. For instance, RELAX introduces a multi-step ICA cleaning method that improves artifact rejection but relies on complex configurations that may not be easily modifiable21. In contrast, PIPEMAT-RS is designed to provide a flexible, transparent, and educationally accessible alternative, enabling researchers to understand and manually adjust each preprocessing step. This transparency allows PIPEMAT-RS to be used as both a research and teaching tool while maintaining rigor and reproducibility. Moreover, PIPEMAT-RS emphasizes documentation, modular execution, and adaptability across datasets-features that are often underreported or absent in other pipelines.
Economically, the development and implementation of a standardized editable preprocessing pipeline for EEG data can have significant impacts. The global EEG devices market was valued at approximately USD 1.2 billion in 2020 and is expected to grow at a compound annual growth rate (CAGR) of 7.5% from 2021 to 202822. Efficient preprocessing pipelines, particularly those that are editable and modular, can reduce the time and cost associated with data analysis-an important consideration in large-scale clinical trials and public health research. By minimizing human error and automating repetitive tasks, such pipelines can improve the accuracy of EEG data analysis and lead to more reliable research outcomes13. This, in turn, can accelerate the development of new diagnostic tools and therapeutic interventions, potentially reducing healthcare costs and improving patient outcomes1.
This paper aims to present and validate a standardized editable Preprocessing Integrated Pipeline for Resting-State EEG using MATLAB (PIPEMAT-RS). Designed to be both robust and efficient, this editable pipeline is versatile and can be widely applied to various datasets and research scenarios. The seven comprehensive steps include file format conversion, downsampling, filtering, artifact rejection, independent component analysis (ICA), and noise related components exclusion, all integrated into a cohesive workflow23,24. The development of PIPEMAT-RS was guided by the need to address common challenges in EEG preprocessing and to provide a tool that researchers can use to streamline their data analysis processes25. By integrating with commonly adopted practices in the field and leveraging the capabilities of MATLAB, PIPEMAT-RS aims to offer a reliable, editable, and user-friendly framework to support EEG data preprocessing26.
To validate the effectiveness of PIPEMAT-RS, we have applied it in several studies. In one study, PIPEMAT-RS was utilized to investigate EEG biomarkers in stroke patients, identifying significant correlations between maladaptive motor function and depressive profiles27. Another study applied PIPEMAT-RS to analyze neural signatures of brain compensation in stroke patients using EEG and TMS, revealing critical insights into the brain&#39;s adaptive mechanisms post-stroke28. Further validation comes from a study that explored the neural adaptations and compensations in fibromyalgia patients29. Additionally, another study demonstrated how delta and theta bands in resting-state EEG serve as compensatory mechanisms in chronic neuropathic pain30.
These studies collectively highlight PIPEMAT-RS&#39;s utility in processing and analyzing EEG data, establishing its efficacy across various conditions and applications. By providing a standardized approach, this pipeline not only enhances the reproducibility of EEG studies but also facilitates the broader adoption of consistent preprocessing methods in the research community. Through this paper, we aim to contribute a valuable resource to the field, promoting more accurate and reliable EEG data analysis16,17. In doing so, we hope to support the advancement of research in neurophysiology and the development of novel therapeutic approaches for neurological and psychiatric conditions.
The foundational preprocessing steps used in PIPEMAT-RS are rooted in robust, widely validated methods that have become essential in EEG research for ensuring data quality and reproducibility2,8,14. Although these methodologies are well-established, they remain crucial due to their effectiveness in addressing common EEG artifacts and noise, which continues to support reliable analyses across studies. However, few studies have recently provided a fully documented, standardized pipeline tailored for resting-state EEG. By integrating these proven techniques with recent advances, such as ICLabel for automated artifact classification, PIPEMAT-RS delivers a contemporary, structured approach that enhances replicability and accessibility for researchers, fulfilling an ongoing need for rigorously documented preprocessing workflows in the field.
While PIPEMAT-RS pipeline follows established preprocessing steps commonly employed in EEG studies, its documentation in the form of a scientific paper fills an important gap in the literature. Many EEG studies rely on similar preprocessing methods but often lack a comprehensive, step-by-step description that allows for full reproducibility and ease of adoption by other researchers. By detailing each component of PIPEMAT-RS and presenting a standardized sequence for resting-state EEG analysis, this work aims to enhance transparency, minimize user variability, and provide a replicable framework applicable across studies. Rather than introducing novel methods, this contribution is intended to serve as an instructional and practical guide that helps standardize preprocessing practices and improve accessibility for a wide range of EEG researchers, particularly those working with resting-state data.
The preprocessing of EEG data is a critical step that significantly influences the quality and reliability of subsequent analyses. The Preprocessing Integrated Pipeline for Resting-State EEG using MATLAB (PIPEMAT-RS) was developed to address common challenges in EEG data preprocessing by providing a comprehensive, standardized approach applicable to various datasets. This structured, editable pipeline consists of seven key steps: i) File Format Conversion; ii) EEG Montage; iii) Downsampling and Filtering; iv) Artifact Rejection and Rereferencing; v) Independent Component Analysis (ICA); vi) ICLabel Classification; vii) Data Normalization. Each step is meticulously designed to enhance signal quality and facilitate the extraction of meaningful information from raw EEG data.
A critical feature of PIPEMAT-RS is that each preprocessing step concludes with saving the dataset under a distinct filename. This structure ensures that, for each processed file, a derivative version corresponding to each specific step is generated. This systematic approach allows users to access data at any point within the seven-step preprocessing pipeline. The saved files include identifiers that clearly indicate the stage of preprocessing, making it easier to track and manage data throughout the workflow.
PIPEMAT-RS begins with file format conversion, where raw EEG data files are transformed into a MATLAB-compatible format. This conversion is essential for the seamless handling of data in MATLAB&#39;s robust computational environment. Following this, the EEG montage step assigns precise spatial locations to each electrode on the scalp based on standardized electrode placement systems. This spatial information is crucial for the accurate interpretation of EEG data in neurophysiological studies.
Next, the pipeline applies downsampling and filtering to reduce the data&#39;s sampling rate and eliminate noise while retaining the frequency components relevant to neurophysiological research. Downsampling decreases the computational load and storage requirements without compromising data quality. Filtering removes specific frequency components known to be associated with artifacts, thereby preserving the essential neural signals.
Subsequently, artifact rejection and rereferencing are performed. This step identifies and removes segments of data contaminated by noise sources such as eye blinks, muscle activity, and electrode movement. Automatic artifact rejection algorithms are utilized to minimize manual intervention, reducing the potential for human error and increasing preprocessing efficiency. Rereferencing the data to the average of all electrodes helps to mitigate the influence of any single electrode and provides a stable reference for the EEG signals.
The Independent Component Analysis (ICA) step further refines the data by separating mixed signals into their independent sources. This decomposition facilitates the identification and removal of artifacts such as eye blinks, muscle artifacts, and line noise. ICA is a powerful technique that ensures the remaining data accurately reflects underlying neural activity.
Finally, ICLabel classification is applied to automatically identify and remove components classified as artifacts. ICLabel assigns probabilistic labels to independent components, categorizing them as brain activity, muscle activity, eye blinks, heartbeats, line noise, or channel noise. Components with a high probability of representing brain activity are retained, while those identified as artifacts are removed. This automated approach significantly reduces manual effort and ensures consistent and objective classification across datasets.
Each of these preprocessing steps is essential for improving the quality of EEG data. By standardizing the preprocessing workflow, PIPEMAT-RS minimizes variability introduced by different preprocessing methods, facilitating the comparison of results across studies. The pipeline&#39;s implementation in MATLAB, a widely-used platform in the neuroscience community, ensures accessibility and ease of integration into existing research workflows.
The development of PIPEMAT-RS was guided by commonly adopted practices in EEG data analysis, incorporating methods and techniques validated in previous research. Each step was carefully designed and rigorously tested to ensure its effectiveness in enhancing data quality (Figure 1). This comprehensive approach not only improves the reliability of EEG data but also supports the identification of neurophysiological markers critical for understanding various neurological and psychiatric conditions.

<table><tbody><tr><td>MATLAB&amp;nbsp;</td><td>MathWorks</td><td>R2020a or newer</td><td>Required for executing the PIPEMAT-RS script.</td></tr><tr><td>EEG Recording System</td><td>Brain Products</td><td>https://www.brainproducts.com/</td><td>EEG acquisition system used for data collection. Other brands&amp;nbsp; include BioSemi / ANT Neuro (or equivalent)&amp;nbsp;</td></tr><tr><td>EEGLAB Toolbox</td><td>Swartz Center for Computational Neuroscience</td><td>N/A</td><td>Open-source MATLAB toolbox for EEG analysis.</td></tr><tr><td>Electrodes (Ag/AgCl)</td><td>Brain Products&amp;nbsp;</td><td>https://www.brainproducts.com/solutions/r-net/</td><td>Used in EEG data acquisition. Other brands&amp;nbsp; include BioSemi / ANT Neuro (or equivalent)&amp;nbsp;</td></tr><tr><td>ICLabel Plugin</td><td>Swartz Center for Computational Neuroscience</td><td>N/A</td><td>Automated artifact classification tool for EEG.</td></tr><tr><td>PIPEMAT-RS Script</td><td>N/A</td><td>N/A</td><td>Custom MATLAB script for standardized EEG preprocessing.</td></tr><tr><td>Signal Amplifier</td><td>Brain Products&amp;nbsp;</td><td>https://www.brainproducts.com/solutions/#amplifiers</td><td>Amplifies EEG signals for processing. Other brands&amp;nbsp; include BioSemi / ANT Neuro (or equivalent)&amp;nbsp;</td></tr><tr><td>Standard 64-Channel EEG Cap</td><td>Brain Products&amp;nbsp;</td><td>https://www.brainproducts.com/solutions/#electrodes-caps</td><td>Electrode cap for EEG recording. Other brands&amp;nbsp; include BioSemi / ANT Neuro (or equivalent)&amp;nbsp;</td></tr></tbody></table>

pipemat-rs: development and validation of a standardized matlab pipeline for resting-state eeg preprocessing

PIPEMAT-RS: Development and Validation of a Standardized MATLAB Pipeline for Resting-State EEG Preprocessing

Electroencephalography (EEG) is a crucial tool in neuroscience research and clinical applications, but raw EEG data often contain noise and artifacts that compromise signal quality. To address this, we developed PIPEMAT-RS, a standardized MATLAB-based preprocessing pipeline for resting-state EEG data. PIPEMAT-RS follows a structured seven-step workflow: file format conversion, EEG montage configuration, downsampling, filtering, artifact rejection, independent component analysis (ICA), and ICLabel classification for automated artifact removal. This protocol enhances EEG data quality by minimizing human intervention while maintaining high accuracy in artifact rejection. It was validated using multiple datasets, demonstrating its robustness in improving signal integrity. PIPEMAT-RS provides a systematic approach that facilitates reproducibility and reliability in EEG studies, aligning with commonly adopted practices in the field and offering a clearly documented structure that can complement existing pipelines. By standardizing EEG preprocessing, PIPEMAT-RS facilitates neurophysiological research and clinical applications, allowing for more accurate interpretations of resting-state brain activity and its associations with neurological and psychiatric conditions.

The validation of the EEG preprocessing pipeline is essential to ensure the reliability and effectiveness of the methods employed for extracting meaningful neurophysiological data. This process involves empirical evidence from published studies, statistical assessments, and comparisons with established benchmarks in the field.
The robustness and utility of PIPEMAT-RS have been demonstrated through its successful application in several published studies spanning a range of neurological and psyc...

This protocol presents PIPEMAT-RS, a standardized MATLAB-based preprocessing pipeline for resting-state EEG data. It ensures artifact removal, improves signal quality, and enhances data reproducibility across studies. The pipeline automates key preprocessing steps, including filtering, independent component analysis (ICA), and artifact classification, facilitating consistent and reliable EEG analysis for neurophysiological research.

Electroencephalography (EEG) is a crucial tool in neuroscience research and clinical applications, but raw EEG data often contain noise and artifacts that ...

pipemat-rs-development-validation-standardized-matlab-pipeline-for

Research

JoVE Journal

Neuroscience

Assessment of the EEG mu rhythm provides a unique methodology for examining brain activity and when combined with behaviorally based assays, can be a powerful tool for elucidating aspects of social cognition, such as imitation, in clinical populations.

EEG Mu Rhythm in Typical and Atypical Development

This video introduces the preparation, recording, and source analysis procedures of high-resolution EEG on sedated rats with a particular preclinical model of focal epilepsy under noninvasive conditions.

Brain Source Imaging in Preclinical Rat Models of Focal Epilepsy using High-Resolution EEG Recordings

This article describes the protocol underlying electroencephalography (EEG) microstate analysis and omega complexity analysis, which are two reference-free EEG measures and highly valuable to explore the neural mechanisms of brain disorders.

Microstate and Omega Complexity Analyses of the Resting-state Electroencephalography

This protocol describes a method for obtaining stable resting-state functional magnetic resonance imaging (rs-fMRI) data from a rat using low dose isoflurane in combination with low dose dexmedetomidine.

Acquisition of Resting-State Functional Magnetic Resonance Imaging Data in the Rat

A standardized pipeline is presented for examining cerebellum grey matter morphometry. The pipeline combines high-resolution, state-of-the-art approaches for optimized and automated cerebellum parcellation and voxel-based registration of the cerebellum for volumetric quantification.

A Standardized Pipeline for Examining Human Cerebellar Grey Matter Morphometry using Structural Magnetic Resonance Imaging

Here, we describe a protocol combining adeno-associated virus injection with cranial window implantation for simultaneous imaging of microglial dynamics and neuronal activity in awake mice.

Simultaneous Imaging of Microglial Dynamics and Neuronal Activity in Awake Mice

The BrainBeats toolbox is an open-source EEGLAB plugin designed to jointly analyze EEG and cardiovascular (ECG/PPG) signals. It includes heartbeat-evoked potentials (HEP) assessment, feature-based analysis, and heart artifact extraction from EEG signals. The protocol will aid in studying brain-heart interplay through two lenses (HEP and features), enhancing reproducibility and accessibility.

BrainBeats as an Open-Source EEGLAB Plugin to Jointly Analyze EEG and Cardiovascular Signals

Presented here is a video analysis pipeline that overcomes challenges in behavioral monitoring within MRI environments, allowing for the detection of uninstructed behavioral responses to external cues. This analysis will facilitate a more comprehensive understanding of evoked internal state changes and brain-wide activity.

Profiling Maternal Behavior Responses During Whole-Brain Imaging

A novel protocol for personalized transcranial magnetic stimulation in combination with electroencephalography (TMS–EEG) mapping of cortical networks, informed by multiple magnetic resonance imaging (MRI) data modalities.

Neuroimaging-Guided TMS&#8211;EEG for Real-Time Cortical Network Mapping

This study presents a standardized and reproducible protocol for implementing the Spatial-Temporal-Frequency EEG analysis tool (STFEEG) for motor imagery EEG decoding, incorporating configurable spatial-temporal-frequency segmentation, Common Spatial Patterns (CSP)-based feature extraction, multiple classification algorithms, and visualization capabilities.

STFEEG-Tool: A Spatial-Temporal-Frequency EEG Analysis Tool for Motor Imagery Brain-Computer Interfaces

This protocol demonstrates how physics-based neural simulations can be used to interpret electrophysiological biomarkers of neurotherapeutics and uncover their effect on neural circuits, providing a mechanistically grounded approach for neurotherapeutic development.

A Protocol For Uncovering Neural Mechanisms Of Neurotherapeutic Effects On Electroencephalography Using The Human Neocortical Neurosolver

Learn PIPEMAT-RS for standardized resting-state EEG preprocessing using MATLAB and ICA-based artifact removal. This pipeline improves signal quality and supports reproducible neurophysiological research across clinical studies.

PIPEMAT-RS: Development and Validation of a Standardized MATLAB Pipeline for Resting-State EEG Preprocessing

Explore More Videos