<div dir="ltr"><div><div><div></div>Thanks for the links to the Nima's paper and your pipeline Makoto, these were extremely helpful. I have started using the PREP pipeline and I like it so far. I assume for subsequent ICA analysis that interpolated channels must be excluded for proper ICA component extraction, correct? I was wondering about the case where no channels are excluded. Since the data have been re-referenced to an average reference the data dimensionality must be reduced by 1 to keep the dataset at full rank. I have found in the EEGLAB wiki that with the runica function, using the flag 'pca' followed by number of principal components to keep (n-1) is a good way to reduce data dimension, but I have also seen suggestions to just remove a data channel and not do a reconstruction from a principal component subspace. What are the pros and cons of each? I'd guess that for the dimension reduction when the 'pca' flag is thrown, runica would project the data in the PC subspace consisting of the eigenvectors with the maximum eigenvalues up to the number of eigenvectors specified in the function argument, but I'm not really sure if this is true.<br><br></div>Thanks!<br></div>Michael<br></div><br><div class="gmail_quote"><div dir="ltr">On Tue, Aug 25, 2015 at 6:00 AM Andreas Widmann <<a href="mailto:widmann@uni-leipzig.de">widmann@uni-leipzig.de</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear Makoto, Michael, and list,<br>
<br>
> Standard filtering policy does not go well with variations of (multivariate) Granger Causality analysis. Many says 1-Hz high-pass filter we do is no good for analyzing averaged ERP. High-pass filter should be done before CleanLine.<br>
Nima introduced a clever solution for the Cleanline vs. highpass filter issue in the prep-pipeline/paper: The low frequency signal components are preserved during highpass filtering (i.e., the difference between filtered and unfiltered data) and re-added to the cleaned data after applying Cleanline.<br>
<br>
Best,<br>
Andreas<br>
<br>
> > 0. Start with 128-channel EEG from EGI HydroCel Geodesic Sensor Net, with flags indicating events and epochs around events predefined (-0.1s to 0.9s around the events).<br>
><br>
> Use -1 to 2 sec to epoch. This is because you want to use 3-Hz 3-cycle wavelet later. It does not matter even if this causes overlap of the trials; it always happen as long as you use sliding windows.<br>
><br>
> Also, you may want to use continuous data as long as it is cleaned.<br>
><br>
> > 1. Keep all valid (even if lots of line noise) channels in dataset, only remove channels which clearly have no physiological information.<br>
><br>
> I agree, unless you discard all the channels.<br>
><br>
> > 2. Epoch the data, reject bad epochs by eye (significant movement or rare, not-stereotyped artifacts)<br>
><br>
> No, make every effort not to reject epochs due to eye activity (saccade, blink, etc) This is because ICA is best capable of identifying them. Save your trials.<br>
><br>
> > 3. Run ICA on epoched data.<br>
><br>
> Make sure that you clean your data before running ICA.<br>
><br>
> > 4. Reject components related to blinks, stereotypical muscle activity, and 60 Hz noise, and EKG artifacts.<br>
><br>
> No, do not reject anything. You'll create STUDY later, in doing which ICs with 15% residual variance and those with outside-brain dipoles are kicked out. Usually this will reject around 70% of your ICs (depending on your data).<br>
><br>
> > 5. Reconstruct sensor-level signals from remaining ICA weights<br>
><br>
> It's a forward projection. Fine.<br>
><br>
> > 6. Re-reference all cleaned EEG channels to the average reference<br>
><br>
> Nima emphasized that when you use average reference you should choose only clean channels to include. See his paper for details. Jason also said average reference should be done before ICA, but it will cause rank reduction (and does not change the results except for zero-centering the scalp maps) so I would recommend you do it after ICA.<br>
><br>
> > 7. Interpolate (using spherical splines from neighboring electrodes) any channels that were rejected in step 1.<br>
><br>
> No it's not necessary if you stick to ICA-centric analysis to the end.<br>
><br>
> > 8. Calculate ERPs, time-frequency analysis, source localization<br>
><br>
> Use STUDY for the group-level analysis.<br>
><br>
> Additional info can be found here.<br>
> <a href="http://sccn.ucsd.edu/wiki/Makoto's_preprocessing_pipeline" rel="noreferrer" target="_blank">http://sccn.ucsd.edu/wiki/Makoto's_preprocessing_pipeline</a><br>
><br>
> Makoto<br>
><br>
> On Tue, Aug 18, 2015 at 11:54 AM, Michael Boyle <<a href="mailto:mrboyle@live.unc.edu" target="_blank">mrboyle@live.unc.edu</a>> wrote:<br>
> Dear EEGLABers,<br>
><br>
> I would like to know whether or not the particular pre-processing steps listed are recommended for ERP analysis, time-frequency analysis (wavelet amplitude spectrograms or DFT spectrograms, ITPC), and source localization analysis.<br>
><br>
> In particular, I am uncertain about the ordering of steps such as interpolating bad channels before or after re-referencing or ICA, or if particular steps work for some analysis strategies but not others (e.g. if a preprocessing step would work fine for time-frequency analysis but not be appropriate for source localization). Any feedback is much appreciated!<br>
><br>
> 0. Start with 128-channel EEG from EGI HydroCel Geodesic Sensor Net, with flags indicating events and epochs around events predefined (-0.1s to 0.9s around the events).<br>
> 1. Keep all valid (even if lots of line noise) channels in dataset, only remove channels which clearly have no physiological information.<br>
> 2. Epoch the data, reject bad epochs by eye (significant movement or rare, not-stereotyped artifacts)<br>
> 3. Run ICA on epoched data.<br>
> 4. Reject components related to blinks, stereotypical muscle activity, and 60 Hz noise, and EKG artifacts.<br>
> 5. Reconstruct sensor-level signals from remaining ICA weights<br>
> 6. Re-reference all cleaned EEG channels to the average reference<br>
> 7. Interpolate (using spherical splines from neighboring electrodes) any channels that were rejected in step 1.<br>
> 8. Calculate ERPs, time-frequency analysis, source localization<br>
><br>
> Best and many thanks,<br>
> Michael<br>
><br>
> _______________________________________________<br>
> Eeglablist page: <a href="http://sccn.ucsd.edu/eeglab/eeglabmail.html" rel="noreferrer" target="_blank">http://sccn.ucsd.edu/eeglab/eeglabmail.html</a><br>
> To unsubscribe, send an empty email to <a href="mailto:eeglablist-unsubscribe@sccn.ucsd.edu" target="_blank">eeglablist-unsubscribe@sccn.ucsd.edu</a><br>
> For digest mode, send an email with the subject "set digest mime" to <a href="mailto:eeglablist-request@sccn.ucsd.edu" target="_blank">eeglablist-request@sccn.ucsd.edu</a><br>
><br>
><br>
><br>
> --<br>
> Makoto Miyakoshi<br>
> Swartz Center for Computational Neuroscience<br>
> Institute for Neural Computation, University of California San Diego<br>
> _______________________________________________<br>
> Eeglablist page: <a href="http://sccn.ucsd.edu/eeglab/eeglabmail.html" rel="noreferrer" target="_blank">http://sccn.ucsd.edu/eeglab/eeglabmail.html</a><br>
> To unsubscribe, send an empty email to <a href="mailto:eeglablist-unsubscribe@sccn.ucsd.edu" target="_blank">eeglablist-unsubscribe@sccn.ucsd.edu</a><br>
> For digest mode, send an email with the subject "set digest mime" to <a href="mailto:eeglablist-request@sccn.ucsd.edu" target="_blank">eeglablist-request@sccn.ucsd.edu</a><br>
<br>
</blockquote></div>