<html><head></head><body><div style="font-family:Helvetica Neue, Helvetica, Arial, sans-serif;font-size:small;"><div>
<!--StartFragment-->
<p class="ydp3b17952cMsoNormal"><span style="font-size: 12pt; line-height: 107%; font-family: "Times New Roman", serif; color: rgb(38, 40, 42);">Dear Arno,</span></p>
<p class="ydp3b17952cMsoListParagraphCxSpFirst"><!--[if !supportLists]--><span style="font-size:10.0pt;mso-bidi-font-size:12.0pt;line-height:107%;font-family:"Helvetica Neue";mso-fareast-font-family:"Helvetica Neue";mso-bidi-font-family:"Helvetica Neue";color:#26282A">1)<span style="font-stretch: normal; font-size: 7pt; line-height: normal; font-family: "Times New Roman";"> </span></span><!--[endif]--><b><span style="font-size:12.0pt;line-height:107%;font-family:"Times New Roman",serif">On Phase Differences in the Original
vs the Delorme ICA Reconstruction: </span></b><span style="font-size:12.0pt;line-height:107%;font-family:"Times New Roman",serif">We can agree or disagree
about whether or not some small eye movement artifact was in the hand selection
that I did. But that misses the main
point here. That is the ICA
reconstruction alters each and every data point in the entire record including
all artifact free portions no matter what one selects. For example, the record is 6 minutes and 51
seconds = 411 seconds. The Mitsar sample
rate was 250 samples per second = 102,750 data samples. Phase difference for each frequency band for
each and every one of the 102,750 data samples has been altered by your own ICA
reconstruction in the EDF file that you emailed to me. Unless you were to sit next to me or if we
do a Team Viewer it is not possible for me to demonstrate this for all of the
data points and then create a power point for all of these data samples. However, I can show some exemplars, for
example, I have created two figures at 4 different time points (1 sec; 2:27 sec;
42 sec & 5:49 sec) that you can download.
You can extract each screen capture and expand them so that you can see that
the exact same time points were selected and the Hilbert transform JTFA for the
4 time points resulted in different phase differences in all channel
combinations with respect to O1 for all frequencies. The same is true no matter which channel is
selected to compute the phase differences in degrees. The same is true also if one computes averages of the instantaneous phase differences or if one uses the FFT. Here
is the download URL: </span></p>
<p class="ydp3b17952cMsoListParagraphCxSpMiddle"><span style="font-size:12.0pt;line-height:107%;font-family:"Times New Roman",serif"><a href="http://www.appliedneuroscience.com/Phase_Diff-Original_&_Delorme-Post-ICA-4_time_points.zip" rel="nofollow" target="_blank">http://www.appliedneuroscience.com/Phase_Diff-Original_&_Delorme-Post-ICA-4_time_points.zip</a></span></p>
<p class="ydp3b17952cMsoListParagraphCxSpMiddle"><br></p>
<p class="ydp3b17952cMsoListParagraphCxSpMiddle"><!--[if !supportLists]--><span class="ydp3b17952capple-converted-space"><span style="font-size:10.0pt;mso-bidi-font-size:12.0pt;line-height:107%;font-family:"Helvetica Neue";mso-fareast-font-family:"Helvetica Neue";mso-bidi-font-family:"Helvetica Neue";color:#26282A">2)<span style="font-stretch: normal; font-size: 7pt; line-height: normal; font-family: "Times New Roman";">
</span></span></span><!--[endif]--><span class="ydp3b17952capple-converted-space"><span style="font-size: 12pt; line-height: 107%; font-family: "Times New Roman", serif; color: rgb(38, 40, 42);"> </span><b>On the WinEEG ICA
Reconstruction: </b>I agree that having
access to ICA components themselves and the topography is critical in
understanding exactly what the WinEEG software did. Unfortunately, I personally do not have
access to the WinEEG software.
Clinician/Scientists in Australia use the WinEEG software and they were
the ones that expressed concern about phase difference distortion at a workshop
in Adelaide and gave me the original and the WinEEG ICA eye movement corrected
files in EDF format. They explained
that they removed only one ICA component for eye movement before they
reconstructed a new time series. At
first, I was impressed because the eye movements were absent in the
reconstructed time series. I then was
able to use JTFA (Hilbert transform) to compare the two edf files and
discovered that all of the phase differences for all channels for all
frequencies had been altered by the ICA reconstruction including artifact free
periods. I could demonstrate this by individual
time comparisons or averages of instantaneous phase differences or by the
FFT. A user of WinEEG explained that
they do not throw away the original raw digital data, however I was told that they
believe that the ICA reconstructed times series is artifact free and therefore
they compute means and standard deviations for their normative database using
the ICA reconstructed data and not the hand edited or artifact deleted original
data samples like other commercial companies do. Your ICA reconstructed time series is actually less different than the original phase difference in comparison to the WinEEG ICA. Nonetheless, both your ICA reconstruction and the WinEEG reconstructions are significantly different than the original recording.</span></p>
<p class="ydp3b17952cMsoListParagraphCxSpLast"><span class="ydp3b17952capple-converted-space"><span style="font-size:12.0pt;line-height:107%;font-family:"Times New Roman",serif"> </span></span></p>
<p class="ydp3b17952cMsoNormal">Best regards,</p>
<p class="ydp3b17952cMsoNormal"> </p>
<p class="ydp3b17952cMsoNormal">Robert</p>
<!--EndFragment-->
Cp ���<br></div><div><br></div><div><br></div><div id="yahoo_quoted_8604611870" class="yahoo_quoted"><div style="font-family:'Helvetica Neue', Helvetica, Arial, sans-serif;font-size:13px;color:#26282a;"><div>On Tuesday, June 20, 2017, 1:12:41 AM EDT, Arnaud Delorme <arno@ucsd.edu> wrote:</div><div><br></div><div><br></div><div><div id="yiv5249673184"><html><head></head><div>Dear Robert,<div class="yiv5249673184"><br clear="none" class="yiv5249673184"></div><div class="yiv5249673184">1) <b class="yiv5249673184">On my ICA decomposition analysis on your data.</b> You have selected a subset of the file where there is 1 minute and 41 second data of eye free data. I was only able to select 40 seconds in the same file, and I also showed that even in this short file, there was some residual eye movements. Jason and Stefan agreed with me. This is the reason why ICA components power spectrum over frontal channels (and frontal channels only) was affected below 10 Hz frequency band in my data analysis. So on my ICA decomposition, our disagreement comes from the interpretation. You feel that the power we remove at low frequency in frontal channel is not eye movement. In an attempt to convince you, I have picked up a clean region from your EDF dataset, and did some dipole localization at this latency. We see that in the clean data, the best dipolar fit (with 2 symmetrical dipoles) ends up near the eye balls with a residual variance of 6.9%. Hopefully this convinces you that your data is not free of eye movement artifacts. If you are willing to take a step further you might contemplate the idea that ICA can remove this residual spurious activity.</div><div class="yiv5249673184"><br clear="none" class="yiv5249673184"></div><div class="yiv5249673184">2) <b class="yiv5249673184">On the WinEEG ICA decomposition analysis.</b> It is critical for us to see the scalp topography (and if possible continuous activity) of the components the people at the Australia workshop selected. Without this, it is not possible for us to comment on the cleaned data. I agree with you that there was some phase distortion in alpha (visible directly in the raw data in the first email you sent) and that this should not be the case. However, without seing the ICA decomposition, it is not possible for us to conclude as to wether people selected the wrong ICA components or if the ICA decomposition implemented in this software is buggy (ICA is not a simple algorithm and it is sensitive to numerical imprecision and a lot of other parameters - a suboptimal implementation could easily explain the WinEEG results). Also, you seem to imply that the WinEEG people were running ICA on their data then throwing away the raw data (which is why their ICA biased neurofeedback database is useless for practical purposes). Is that correct? One should never throw away the raw data. If they did throw away the raw data, it is an indication that the WinEEG are not rigorous in their approach and therefore might not have implemented ICA in an optimal way. If it is not the case, one may easily reconstruct the database of measures with or without ICA decomposition (assuming ICA is done right which does not seem to be the case) then assess data measure distoritions (power, phase index, etc…) in a statistical fashion.</div><div class="yiv5249673184"><br clear="none" class="yiv5249673184"></div><div class="yiv5249673184">Best wishes,</div><div class="yiv5249673184"><br clear="none" class="yiv5249673184"></div><div class="yiv5249673184">Arno</div><div class="yiv5249673184"><br clear="none" class="yiv5249673184"></div><div class="yiv5249673184"><a rel="nofollow" shape="rect" class="yiv5249673184" target="_blank" href="http://sccn.ucsd.edu/~arno/download/clean_edf_file_analysis2.pdf">http://sccn.ucsd.edu/~arno/download/clean_edf_file_analysis2.pdf</a></div><div class="yiv5249673184"><br clear="none" class="yiv5249673184"><div><blockquote class="yiv5249673184" type="cite"><div class="yiv5249673184">On Jun 18, 2017, at 11:44 AM, Robert Thatcher <<a rel="nofollow" shape="rect" class="yiv5249673184" ymailto="mailto:rwthatcher2@yahoo.com" target="_blank" href="mailto:rwthatcher2@yahoo.com">rwthatcher2@yahoo.com</a>> wrote:</div><br clear="none" class="yiv5249673184Apple-interchange-newline"><div class="yiv5249673184"><span><Pre-ICA-Hand Artifact free selections.edf></span></div></blockquote></div><div class="yiv5249673184yqt0753035364" id="yiv5249673184yqtfd66224"><br clear="none" class="yiv5249673184"></div></div></div></html></div></div></div></div></div></body></html>