Sara, Given enough channels (i.e., deg of freedom), ICA separates out each scalp/neck muscle activity -- Note there may be information of interest in their independent activities as well! Therefore, EEG montages including low electrodes need relatively high density there, since the spatial source (muscle) density is high around the insertions of the neck muscles into the skull (this is where the EMG activity appears on the scalp, with a polarity reversal in the direction of the muscle and centered at the muscle end). This is quite likely why in your experience muscle activity "usually ... spreads on half of the components," as you say.<div>
<div><br></div><div>Scott Makeig<br><br><div class="gmail_quote">On Thu, Oct 6, 2011 at 2:50 AM, Sara Graziadio <span dir="ltr"><<a href="mailto:sara.graziadio@newcastle.ac.uk">sara.graziadio@newcastle.ac.uk</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">Hello,<br>
Thanks for your suggestion.<br>
<br>
As I was planning to do also a PSD analysis on the data I guess that to remove the mean is not the best method if it works as a non-selective high pass filter, am I right?<br>
<br>
I am applying the PCA before applying the ICA to reduce the number of components. How the data rank would be modified in this case?<br>
I have to admit that it never happened to me that the muscle artefact is put in a single source with the ICA. Usually it spreads on half of the components, is this only my experience?<br>
<br>
Thanks again<br>
<br>
Best wishes<br>
<br>
Sara<br>
<br>
<br>
>-----Original Message-----<br>
>From: David Groppe [mailto:<a href="mailto:david.m.groppe@gmail.com">david.m.groppe@gmail.com</a>]<br>
>Sent: 05 October 2011 23:10<br>
>To: Sara Graziadio<br>
>Cc: <a href="mailto:eeglablist@sccn.ucsd.edu">eeglablist@sccn.ucsd.edu</a><br>
>Subject: Re: [Eeglablist] filters, ICA and erp<br>
><br>
>Hi Sara,<br>
> I found that a good way to improve the performance of ICA for ERP<br>
>analysis is to<br>
>1) Epoch your data into one or two second chunks time locked to the<br>
>event of interest<br>
>2) Remove the mean of each epoch at each channel<br>
>3) Run ICA to remove artifacts<br>
>4) Use a standard pre-event time window to baseline your data<br>
>5) Compute your ERPs<br>
><br>
>Removing the mean of each epoch acts as a crude high-pass filter.<br>
>It's not nearly as selective as a "true" high pass filter but it<br>
>doesn't distort the ERP waveforms as much either. Moreover we've<br>
>found that the procedure described above massively improves the<br>
>reliability of ICA when compared to standard ERP prestimulus<br>
>baselines:<br>
><br>
>Groppe, D.M., Makeig, S., & Kutas, M. (2009) Identifying reliable<br>
>independent components via split-half comparisons. NeuroImage, 45<br>
>pp.1199-1211.<br>
><br>
>Hope this helps,<br>
> -David<br>
<div><div></div><div class="h5">><br>
><br>
><br>
>On Wed, Oct 5, 2011 at 10:46 AM, Sara Graziadio<br>
><<a href="mailto:sara.graziadio@newcastle.ac.uk">sara.graziadio@newcastle.ac.uk</a>> wrote:<br>
>> Hello,<br>
>> I would like just a suggestion about some data cleaning/analysis I am doing. I<br>
>am doing an ERP analysis and I want to clean my data first with the ICA. In<br>
>theory, though, I should not use an high-pass cutoff higher than 0.1 Hz to not<br>
>reduce the erp amplitude. On the other side the ICA does not work well if the<br>
>high-pass cutoff is lower than 0.5 Hz...what is then the best method to apply?<br>
>Has anybody tested how robust the ica is with a 0.1Hz filter?<br>
>> I have also another question: I am doing the analysis on 94 electrodes<br>
>referenced to Fz. I planned to average reference the data but actually there is<br>
>quite a large spread of noise on all the electrodes with this method (muscular<br>
>artefacts for example from the temporal electrodes). But actually almost all<br>
>the papers are using the average reference so I was surprised, am I the only<br>
>one having this problem of noise? Would not be better just to keep the Fz<br>
>reference and then perhaps to average the erps for every different cortical<br>
>area and do the analysis on these averaged erps?<br>
>><br>
>> Thank you very much<br>
>><br>
>> Best wishes<br>
>><br>
>> Sara Graziadio<br>
>> Research Associate<br>
>> Newcastle University<br>
>><br>
>><br>
>><br>
>> _______________________________________________<br>
>> Eeglablist page: <a href="http://sccn.ucsd.edu/eeglab/eeglabmail.html" target="_blank">http://sccn.ucsd.edu/eeglab/eeglabmail.html</a><br>
>> To unsubscribe, send an empty email to eeglablist-<br>
><a href="mailto:unsubscribe@sccn.ucsd.edu">unsubscribe@sccn.ucsd.edu</a><br>
>> For digest mode, send an email with the subject "set digest mime" to<br>
><a href="mailto:eeglablist-request@sccn.ucsd.edu">eeglablist-request@sccn.ucsd.edu</a><br>
>><br>
><br>
><br>
><br>
</div></div>>--<br>
>David Groppe, Ph.D.<br>
>Postdoctoral Researcher<br>
>North Shore LIJ Health System<br>
>New Hyde Park, New York<br>
><a href="http://www.cogsci.ucsd.edu/~dgroppe/" target="_blank">http://www.cogsci.ucsd.edu/~dgroppe/</a><br>
<div><div></div><div class="h5"><br>
_______________________________________________<br>
Eeglablist page: <a href="http://sccn.ucsd.edu/eeglab/eeglabmail.html" 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">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">eeglablist-request@sccn.ucsd.edu</a><br>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br>Scott Makeig, Research Scientist and Director, Swartz Center for Computational Neuroscience, Institute for Neural Computation; Prof. of Neurosciences (Adj.), University of California San Diego, La Jolla CA 92093-0559, <a href="http://sccn.ucsd.edu/%7Escott" target="_blank">http://sccn.ucsd.edu/~scott</a><br>
</div></div>