<div dir="ltr">Dear Alexandra,<div><br></div><div>> (Minor question: At this point I do not baseline correct each epoch, is this recommended?)<br></div><div class="gmail_extra"><br></div><div class="gmail_extra">Yes, I recommend it.</div><div class="gmail_extra"><a href="https://sccn.ucsd.edu/wiki/Makoto's_preprocessing_pipeline#Epoch_data_to_-1_to_2_sec">https://sccn.ucsd.edu/wiki/Makoto's_preprocessing_pipeline#Epoch_data_to_-1_to_2_sec</a><br></div><div class="gmail_extra"><br></div><div class="gmail_extra">> Here is where I am losing track. For this example, there are 116 non-interpolated channels and the data set was re-referenced to common reference before ICA.<br>So I guess this might be behind the rank reduction, although I am not sure how runica knows about this (where is such information retrieved?)</div><div class="gmail_extra"><br></div><div class="gmail_extra">It computes rank in pop_runica() line 523 getrank().<br><br>> What precisely is used as the input data now? Any first 115 channels or only the non-interpolated channels?</div><div class="gmail_extra"><br></div><div class="gmail_extra">It's not the channels, but the rank of data that matters. See this.</div><div class="gmail_extra"><a href="https://sccn.ucsd.edu/wiki/Makoto's_preprocessing_pipeline#Re-reference_the_data_to_average">https://sccn.ucsd.edu/wiki/Makoto's_preprocessing_pipeline#Re-reference_the_data_to_average</a><br></div><div class="gmail_extra"><br></div><div class="gmail_extra">> What happens when the resulting (clean) components are projected back onto channel space if some of the channels are not actually included in the data set used to calculate the components?<br></div><div class="gmail_extra"><br></div><div class="gmail_extra">All of your scalp channels will be there, so no worries.</div><div class="gmail_extra">It's just that you have less number of *independent* sources than the number of channels. Even if you have 1 independent component, it can project to 100 channels. In your case, you have 99 components projecting to 100 channels (i.e. one less).</div><div class="gmail_extra"><br></div><div class="gmail_extra">Makoto</div><div class="gmail_extra"><br></div><div class="gmail_extra"><br></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Aug 16, 2016 at 9:18 AM, Alexandra Yvonne Vossen <span dir="ltr"><<a href="mailto:a.vossen.1@research.gla.ac.uk" target="_blank">a.vossen.1@research.gla.ac.uk</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><br>
Dear all,<br>
<br>
I realise this is a popular question but I am really uncertain what happens to my data in my preprocessing pipeline. I am quite new to ICA/PCA decomposition.<br>
<br>
First, the raw data are run through the PREP pipeline (in EEGLAB 13.5.4b), including the (spherical) interpolation of noisy channels.<br>
<br>
Then after rejection of very noisy epochs I run ICA (for identification of artefact components) with EEG = pop_runica(EEG,'extended',1) %this is somebody else's code that I am re-using.<br>
(Minor question: At this point I do not baseline correct each epoch, is this recommended?)<br>
<br>
I get the info: "Data rank (115) is smaller than the number of channels (128).  Input data size [115,525723] = 115 channels, 525723 frames/nAfter PCA dimension reduction,<br>
finding 115 ICA components using extended ICA."<br>
<br>
Here is where I am losing track. For this example, there are 116 non-interpolated channels and the data set was re-referenced to common reference before ICA.<br>
So I guess this might be behind the rank reduction, although I am not sure how runica knows about this (where is such information retrieved?)<br>
<br>
What precisely is used as the input data now? Any first 115 channels or only the non-interpolated channels?<br>
What happens when the resulting (clean) components are projected back onto channel space if some of the channels are not actually included in the data set used to calculate the components?<br>
<br>
Many thanks in advance,<br>
<br>
<br>
Alexandra Vossen<br>
PhD student<br>
School of Psychology<br>
College of Science & Engineering<br>
University of Glasgow<br>
______________________________<wbr>_________________<br>
Eeglablist page: <a href="http://sccn.ucsd.edu/eeglab/eeglabmail.html" rel="noreferrer" target="_blank">http://sccn.ucsd.edu/eeglab/<wbr>eeglabmail.html</a><br>
To unsubscribe, send an empty email to <a href="mailto:eeglablist-unsubscribe@sccn.ucsd.edu">eeglablist-unsubscribe@sccn.<wbr>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.<wbr>edu</a><br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr">Makoto Miyakoshi<br>Swartz Center for Computational Neuroscience<br>Institute for Neural Computation, University of California San Diego<br></div></div>
</div></div>