<div dir="ltr"><div class="gmail_default" style="font-family:courier new,monospace">I first noticed the problem with phase distortion more than a decade ago when Marco Congedo implemented ICA artifact removal in software he developed. I pointed out this phase problem to him and he was unable to solve it.</div><div class="gmail_default" style="font-family:courier new,monospace"><br></div><div class="gmail_default" style="font-family:courier new,monospace">I told Arnaud about it in 2004 and we have been debating this ever since. I am glad to see you all solving the question empirically.</div><div class="gmail_default" style="font-family:courier new,monospace"><br></div><div class="gmail_default" style="font-family:courier new,monospace">I had an interesting talk with Phil Zeman about this last year. He said that phase distortion may be more of a problem with fewer channels. Removing 2 components from a 19 channel recording is 10 percent of the information. Removing 2 components from a 256 channel recording is less than 0.5 percent of the information. Another good experiment would be to do ICA on a high density recording- 128 or 256 channels. Remove the artifact components, reconstruct the data set and then pick out the 19 10-20 channels. Compare the phase distortion to the same channels with all components in it. Do the same experiment but pick out the 19 channels before ICA and then reconstruct this data with 17 or 18 components. It would also be interesting to compare the artifact components themselves from 256 channels of data and only the 19 channels of 10-20 system.</div><div class="gmail_default" style="font-family:courier new,monospace"><br></div><div class="gmail_default" style="font-family:courier new,monospace">Robert Lawson</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Jun 13, 2017 at 11:15 PM, Arnaud Delorme <span dir="ltr"><<a href="mailto:arno@ucsd.edu" target="_blank">arno@ucsd.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word">Dear Robert,<div><br></div><div>Below is the EDF file with the components removed.</div><div>There was 19 channels, so I got 19 components, of which I removed 2 (components 1 and 2). 1 was blink. 2 was lateral eye movements.</div><div>I have rerun the decomposition keeping all the data. This way you can align it with the dataset you sent me.</div><div><br></div><div><a href="https://sccn.ucsd.edu/~arno/download/Australia-Post-ICA_delorme.edf" target="_blank">https://sccn.ucsd.edu/~arno/<wbr>download/Australia-Post-ICA_<wbr>delorme.edf</a></div><div><br></div><div>Best wishes,</div><div><br></div><div>Arno</div><div><br><div><blockquote type="cite"><div>On Jun 13, 2017, at 4:11 PM, Robert Thatcher <<a href="mailto:rwthatcher2@yahoo.com" target="_blank">rwthatcher2@yahoo.com</a>> wrote:</div><br class="m_-4019525169749825605Apple-interchange-newline"><div>
<div><div style="font-family:Helvetica Neue,Helvetica,Arial,sans-serif;font-size:small"><div><div>Arnaud,</div></div><div>   Thank you. Please attach the raw digital data produced after ICA reconstruction.  Images are not sufficient and also what is of concern is the phase differences between all of the channels and this cannot be evaluated without having access to the raw digital data following ICA reconstruction..</div><div><br></div><div>I look forward to receiving an EDF formatted data file with the post ICA reconstruction.  Also, how many and which ICA components did you delete to re-create the 19 channels, was it 18 or 17 or 16, etc?</div><div><br></div><div>Thank you,</div><div><br></div><div>Robert</div><div><br></div><div><br></div><div id="m_-4019525169749825605yahoo_quoted_8109086840" class="m_-4019525169749825605yahoo_quoted"><div style="font-family:'Helvetica Neue',Helvetica,Arial,sans-serif;font-size:13px;color:#26282a"><div>On Tuesday, June 13, 2017, 6:47:15 PM EDT, Arnaud Delorme <<a href="mailto:arno@ucsd.edu" target="_blank">arno@ucsd.edu</a>> wrote:</div><div><br></div><div><br></div><div><div id="m_-4019525169749825605yiv9205346509"><div>Hi Robert,<div class="m_-4019525169749825605yiv9205346509"><br clear="none" class="m_-4019525169749825605yiv9205346509"></div><div class="m_-4019525169749825605yiv9205346509">I have imported your dataset in EEGLAB and I have looked at the data. The ICA components I get are actually beautiful.</div><div class="m_-4019525169749825605yiv9205346509"><br clear="none" class="m_-4019525169749825605yiv9205346509"></div><div class="m_-4019525169749825605yiv9205346509">I have taken a couple of screen captures. The first 2 components are eye blinks and lateral eye movements. </div><div class="m_-4019525169749825605yiv9205346509"><br clear="none" class="m_-4019525169749825605yiv9205346509"></div><div class="m_-4019525169749825605yiv9205346509">Then I used one of the EEGLAB functions to compare before (in black) and after (red) ICA rejection. There are small captures of alpha oscillation before and after removing ICA. You can see that the phase of alpha is preserved (there is no need for spectral decomposition to see that the peaks are almost perfectly aligned which means that the phase is preserved).</div><div class="m_-4019525169749825605yiv9205346509"><br clear="none" class="m_-4019525169749825605yiv9205346509"></div><div class="m_-4019525169749825605yiv9205346509">I attached the code that you can use to get these results in EEGLAB (you will need to install the BIOSIG extension to import EDF files). I have only used EEGLAB default (except manual rejection of 2 portions of data totaling 3 seconds or so but you can try without that and it will most likely return very similar results). Imported the data, ran ICA (Infomax as in the paper) using default parameters and plotted results using EEGLAB menus. I did not tweak the data in any way.</div><div class="m_-4019525169749825605yiv9205346509"><br clear="none" class="m_-4019525169749825605yiv9205346509"></div><div class="m_-4019525169749825605yiv9205346509">Let me know if anything is unclear.</div><div class="m_-4019525169749825605yiv9205346509"><br clear="none" class="m_-4019525169749825605yiv9205346509"></div><div class="m_-4019525169749825605yiv9205346509">Best wishes,</div><div class="m_-4019525169749825605yiv9205346509"><br clear="none" class="m_-4019525169749825605yiv9205346509"></div><div class="m_-4019525169749825605yiv9205346509">Arno</div><div class="m_-4019525169749825605yiv9205346509"><span id="m_-4019525169749825605cid:oL1SLJDJIWU3q4zGrQEV"><Screen Shot 2017-06-13 at 3.31.50 PM.png></span><span id="m_-4019525169749825605cid:iUQ249Z1Ai7BLDwqCKOS"><Screen Shot 2017-06-13 at 3.31.28 PM.png></span><span id="m_-4019525169749825605cid:rU8Ek3YB98h0AwdX65PX"><Screen Shot 2017-06-13 at 3.30.56 PM.png></span><span id="m_-4019525169749825605cid:JK5dKNjSUpJ1nrVJ6ZSI"><Screen Shot 2017-06-13 at 3.30.04 PM.png></span><span id="m_-4019525169749825605cid:XmDwYIsMyU3MExAhL4te"><Screen Shot 2017-06-13 at 3.29.02 PM.png></span><span id="m_-4019525169749825605cid:kt6JQYpxskLc4fPCiaNS"><Screen Shot 2017-06-13 at 3.28.26 PM.png></span></div><div class="m_-4019525169749825605yiv9205346509"><br clear="none" class="m_-4019525169749825605yiv9205346509"></div><div class="m_-4019525169749825605yiv9205346509">% import the data</div><div class="m_-4019525169749825605yiv9205346509">EEG = pop_biosig('/Users/arno/<wbr>Downloads/Australia-Pre-ICA.<wbr>edf’);</div><div class="m_-4019525169749825605yiv9205346509"><br clear="none" class="m_-4019525169749825605yiv9205346509"></div><div class="m_-4019525169749825605yiv9205346509">% remove some bad portion of data</div><div class="m_-4019525169749825605yiv9205346509">EEG = eeg_eegrej( EEG, [11 922;21141 21618]);</div><div class="m_-4019525169749825605yiv9205346509"><br clear="none" class="m_-4019525169749825605yiv9205346509"></div><div class="m_-4019525169749825605yiv9205346509">% look up channel locations<br clear="none" class="m_-4019525169749825605yiv9205346509">EEG=pop_chanedit(EEG, 'lookup'<wbr>,'/data/matlab/eeglab/plugins/<wbr>dipfit2.3/standard_BESA/<wbr>standard-10-5-cap385.elp','<wbr>rplurchanloc',1,'changefield',<wbr>{1 'labels' 'FP1'},'<wbr>changefield',{2 'labels' 'FP2'<wbr>},'changefield',{3 'labels' '<wbr>F7'},'changefield',{4 'labels'<wbr> 'F3'},'changefield',{5 '<wbr>labels' 'Fz'},'changefield',{<wbr>6 'labels' 'F4'},'changefield'<wbr>,{7 'labels' 'F8'},'<wbr>changefield',{8 'labels' 'T3'}<wbr>,'changefield',{9 'labels' '<wbr>C3'},'changefield',{10 '<wbr>labels' 'Cz'},'changefield',{<wbr>11 'labels' 'C4'},'<wbr>changefield',{12 'labels' 'T4'<wbr>},'changefield',{13 'labels' '<wbr>T5'},'changefield',{14 '<wbr>labels' 'P3'},'changefield',{<wbr>15 'labels' 'Pz'},'<wbr>changefield',{16 'labels' 'P4'<wbr>},'changefield',{17 'labels' '<wbr>T6'},'changefield',{18 '<wbr>labels' 'O1'},'changefield',{<wbr>19 'labels' 'O2'},'lookup','/<wbr>data/matlab/eeglab/plugins/<wbr>dipfit2.3/standard_BESA/<wbr>standard-10-5-cap385.elp','<wbr>rplurchanloc',1);<br clear="none" class="m_-4019525169749825605yiv9205346509"><br clear="none" class="m_-4019525169749825605yiv9205346509"></div><div class="m_-4019525169749825605yiv9205346509">% run ICA (default)</div><div class="m_-4019525169749825605yiv9205346509">EEG = pop_runica(EEG, 'extended',1,'<wbr>interupt','on’);</div><div class="m_-4019525169749825605yiv9205346509"><br clear="none" class="m_-4019525169749825605yiv9205346509"></div><div class="m_-4019525169749825605yiv9205346509">% plot properties of the first two components</div><div class="m_-4019525169749825605yiv9205346509">pop_prop( EEG, 0, [1  2], NaN, {'freqrange' [2 50] });<br clear="none" class="m_-4019525169749825605yiv9205346509"><br clear="none" class="m_-4019525169749825605yiv9205346509">% Use menu Tools > Remove components to compare the data before and after ICA<br clear="none" class="m_-4019525169749825605yiv9205346509"></div><div class="m_-4019525169749825605yiv9205346509yqt0588457362" id="m_-4019525169749825605yiv9205346509yqt78581"><div class="m_-4019525169749825605yiv9205346509"><br clear="none" class="m_-4019525169749825605yiv9205346509"><div><blockquote class="m_-4019525169749825605yiv9205346509" type="cite"><div class="m_-4019525169749825605yiv9205346509">On Jun 13, 2017, at 2:32 PM, Robert Thatcher <<a rel="nofollow" shape="rect" class="m_-4019525169749825605yiv9205346509" href="mailto:rwthatcher2@yahoo.com" target="_blank">rwthatcher2@yahoo.com</a>> wrote:</div><br clear="none" class="m_-4019525169749825605yiv9205346509Apple-interchange-newline"><div class="m_-4019525169749825605yiv9205346509">
</div></blockquote></div></div></div></div><div class="m_-4019525169749825605yiv9205346509yqt0588457362" id="m_-4019525169749825605yiv9205346509yqt37625"><div><div class="m_-4019525169749825605yiv9205346509"><div class="m_-4019525169749825605yiv9205346509" style="font-family:Helvetica Neue,Helvetica,Arial,sans-serif;font-size:small"><div class="m_-4019525169749825605yiv9205346509"><div class="m_-4019525169749825605yiv9205346509">Hi Arno,</div></div><div class="m_-4019525169749825605yiv9205346509">   I agree that the study by Montefusco-Siegmund  et al does not give a lot of detail and they used simulated EEG and simulated eye movement and a very large time series and they used phase unwrapping which is not necessary when one uses JTFA like the Hilbert transform.  They produced reliable phase changes due to the ICA reconstruction, albeit small in magnitude.  The analyses by Georges and colleaques used real human EEG and standard ICA software and more realistic simulations with much larger changes in phase differences due to ICA reconstruction.  Large changes in phase differences by ICA is best demonstrated by regular type EEG data samples and with more realistic mathematical simulations.  The EEG samples from the Australian workshop are better examples of ICA phase distortion and these examples are commonplace and WinEEG is commercial software used world wide.   Georges and colleaques as well as myself and colleaques can easily reproduce large phase distortion and will do so in a future publication.  I expect that we will share our analyses on this forum because of the rare expertise and interest by member of this forum.</div><div class="m_-4019525169749825605yiv9205346509"><br clear="none" class="m_-4019525169749825605yiv9205346509"></div><div class="m_-4019525169749825605yiv9205346509">It would be helpful for all of us together to mathematically and empirically test the changes in phase differences due to ICA reconstruction to allegedly remove artifact.   I think that at the end of the day we will find agreement and this will strengthen the proper uses of ICA.   </div><div class="m_-4019525169749825605yiv9205346509"><br clear="none" class="m_-4019525169749825605yiv9205346509"></div><div class="m_-4019525169749825605yiv9205346509">An example of important advantages of ICA are the last few postings on this forum about 3D source vectors and the best way to use the cross-spectrum and the imaginary numbers to derive reliable estimates of network hub/node coupling.  ICA feature detection is excellent, the problem is limited to the mathematically "ill-posed" reconstruction when used for the purposes of artifact rejection.</div><div class="m_-4019525169749825605yiv9205346509"><br clear="none" class="m_-4019525169749825605yiv9205346509"></div><div class="m_-4019525169749825605yiv9205346509">Best wishes,</div><div class="m_-4019525169749825605yiv9205346509"><br clear="none" class="m_-4019525169749825605yiv9205346509"></div><div class="m_-4019525169749825605yiv9205346509">Robert</div><div class="m_-4019525169749825605yiv9205346509"><br clear="none" class="m_-4019525169749825605yiv9205346509"></div><div class="m_-4019525169749825605yiv9205346509"><br clear="none" class="m_-4019525169749825605yiv9205346509"></div><div class="m_-4019525169749825605yiv9205346509ydpae948b6yahoo_quoted" id="m_-4019525169749825605yiv9205346509ydpae948b6yahoo_quoted_7835141380"><div class="m_-4019525169749825605yiv9205346509" style="font-family:'Helvetica Neue',Helvetica,Arial,sans-serif;font-size:13px;color:#26282a"><div class="m_-4019525169749825605yiv9205346509">On Tuesday, June 13, 2017, 3:15:51 PM EDT, Arnaud Delorme <<a rel="nofollow" shape="rect" class="m_-4019525169749825605yiv9205346509" href="mailto:arno@ucsd.edu" target="_blank">arno@ucsd.edu</a>> wrote:</div><div class="m_-4019525169749825605yiv9205346509"><br clear="none" class="m_-4019525169749825605yiv9205346509"></div><div class="m_-4019525169749825605yiv9205346509"><br clear="none" class="m_-4019525169749825605yiv9205346509"></div><div class="m_-4019525169749825605yiv9205346509"><div class="m_-4019525169749825605yiv9205346509" id="m_-4019525169749825605yiv9205346509ydpae948b6yiv3412171908"><div class="m_-4019525169749825605yiv9205346509">Hi Robert,<div class="m_-4019525169749825605yiv9205346509ydpae948b6yiv3412171908"><br clear="none" class="m_-4019525169749825605yiv9205346509ydpae948b6yiv3412171908"></div><div class="m_-4019525169749825605yiv9205346509ydpae948b6yiv3412171908">Thank you for sharing the manuscript.</div><div class="m_-4019525169749825605yiv9205346509ydpae948b6yiv3412171908"><br clear="none" class="m_-4019525169749825605yiv9205346509ydpae948b6yiv3412171908"></div><div class="m_-4019525169749825605yiv9205346509ydpae948b6yiv3412171908">Regarding the actual EEG data (figure 1), component activities seem quite noisy for blink components and the scalp topographies also appear patchy which is not typical of ICA blink or eye movement components which are usually quite smooth. It might be that the ICA decomposition was suboptimal (preprocessing of the data is not detailed and this might be the reason why - although sometimes ICA also fails to return meaningful components for reasons which are not well understood). The units on this figure are wrong by the way (ICA component activity or scalp topography are not in microvolt). So, the poor quality of this decomposition might be responsible for the phase distortion you are mentioning.</div><div class="m_-4019525169749825605yiv9205346509ydpae948b6yiv3412171908"><br clear="none" class="m_-4019525169749825605yiv9205346509ydpae948b6yiv3412171908"></div><div class="m_-4019525169749825605yiv9205346509ydpae948b6yiv3412171908">Regarding the simulated data, projection to scalp channels are usually done using a leadfield matrix that models the conductivity of the different brain layers. In this paper, it seems to be a simple linear combination (and the coefficients and the positions of the sources are not indicated which makes it hard to reproduce these results). The phase difference shown in Figure 5 (on the order of 0.001 radian for most of them) even if significant, do not appear convincing - is a change in phase of 0.001 radian really that critical. Changing the ICA stop threshold might actually dramatically affect these results.</div><div class="m_-4019525169749825605yiv9205346509ydpae948b6yiv3412171908"><br clear="none" class="m_-4019525169749825605yiv9205346509ydpae948b6yiv3412171908"></div><div class="m_-4019525169749825605yiv9205346509ydpae948b6yiv3412171908">I will take some time to download your dataset and have a look at it.</div><div class="m_-4019525169749825605yiv9205346509ydpae948b6yiv3412171908">Best wishes,</div><div class="m_-4019525169749825605yiv9205346509ydpae948b6yiv3412171908"><br clear="none" class="m_-4019525169749825605yiv9205346509ydpae948b6yiv3412171908"></div><div class="m_-4019525169749825605yiv9205346509ydpae948b6yiv3412171908">Arno</div><div class="m_-4019525169749825605yiv9205346509ydpae948b6yiv3412171908"><br clear="none" class="m_-4019525169749825605yiv9205346509ydpae948b6yiv3412171908"><div class="m_-4019525169749825605yiv9205346509"><blockquote class="m_-4019525169749825605yiv9205346509ydpae948b6yiv3412171908" type="cite"><div class="m_-4019525169749825605yiv9205346509ydpae948b6yiv3412171908yqt2544582234" id="m_-4019525169749825605yiv9205346509ydpae948b6yiv3412171908yqtfd67336"><div class="m_-4019525169749825605yiv9205346509ydpae948b6yiv3412171908">On Jun 13, 2017, at 11:34 AM, Robert Thatcher <<a rel="nofollow" shape="rect" class="m_-4019525169749825605yiv9205346509ydpae948b6yiv3412171908" href="mailto:rwthatcher2@yahoo.com" target="_blank">rwthatcher2@yahoo.com</a>> wrote:</div><br clear="none" class="m_-4019525169749825605yiv9205346509ydpae948b6yiv3412171908Apple-interchange-newline"><div class="m_-4019525169749825605yiv9205346509ydpae948b6yiv3412171908"><div class="m_-4019525169749825605yiv9205346509ydpae948b6yiv3412171908" style="font-family:'Helvetica Neue',Helvetica,Arial,sans-serif;font-style:normal;font-weight:normal;letter-spacing:normal;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">ICA is a wonderful tool for feature extraction but it has limitations when used to allegedly remove artifact.   The brain is not a large balloon filled with saline – there are hubs and clusters of neurons connected in networks (Brodmann areas, etc).  The phase differences between hubs are due to differences in the synaptic rise times and synaptic integration times and differences in conduction velocity, etc.   These differences are vital and critical to understanding brain function (see Nunez, “Electrical Fields of the Brain”, Oxford Univ. Press, 1981 or Walter Freeman and many others).</div><br clear="none" class="m_-4019525169749825605yiv9205346509ydpae948b6yiv3412171908Apple-interchange-newline"></div></div></blockquote></div><div class="m_-4019525169749825605yiv9205346509ydpae948b6yiv3412171908yqt2544582234" id="m_-4019525169749825605yiv9205346509ydpae948b6yiv3412171908yqtfd42203"><br clear="none" class="m_-4019525169749825605yiv9205346509ydpae948b6yiv3412171908"></div></div></div></div></div></div></div></div></div><br clear="none" class="m_-4019525169749825605yiv9205346509"></div></div></div></div></div></div></div></div><span id="m_-4019525169749825605cid:JK5dKNjSUpJ1nrVJ6ZSI"><Screen Shot 2017-06-13 at 3.30.04 PM.png></span><span id="m_-4019525169749825605cid:rU8Ek3YB98h0AwdX65PX"><Screen Shot 2017-06-13 at 3.30.56 PM.png></span><span id="m_-4019525169749825605cid:kt6JQYpxskLc4fPCiaNS"><Screen Shot 2017-06-13 at 3.28.26 PM.png></span><span id="m_-4019525169749825605cid:XmDwYIsMyU3MExAhL4te"><Screen Shot 2017-06-13 at 3.29.02 PM.png></span><span id="m_-4019525169749825605cid:iUQ249Z1Ai7BLDwqCKOS"><Screen Shot 2017-06-13 at 3.31.28 PM.png></span><span id="m_-4019525169749825605cid:oL1SLJDJIWU3q4zGrQEV"><Screen Shot 2017-06-13 at 3.31.50 PM.png></span></div></blockquote></div><br></div></div><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></div>