<html xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"><head><meta http-equiv=Content-Type content="text/html; charset=utf-8"><meta name=Generator content="Microsoft Word 15 (filtered medium)"><style><!--
/* Font Definitions */
@font-face
        {font-family:Helvetica;
        panose-1:2 11 6 4 2 2 2 2 2 4;}
@font-face
        {font-family:"Cambria Math";
        panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
        {font-family:Calibri;
        panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
        {margin:0cm;
        margin-bottom:.0001pt;
        font-size:11.0pt;
        font-family:"Calibri",sans-serif;}
a:link, span.MsoHyperlink
        {mso-style-priority:99;
        color:blue;
        text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
        {mso-style-priority:99;
        color:purple;
        text-decoration:underline;}
p.msonormal0, li.msonormal0, div.msonormal0
        {mso-style-name:msonormal;
        mso-margin-top-alt:auto;
        margin-right:0cm;
        mso-margin-bottom-alt:auto;
        margin-left:0cm;
        font-size:11.0pt;
        font-family:"Calibri",sans-serif;}
p.gmail-m-1834111884996008056ydp82770b41msonormal, li.gmail-m-1834111884996008056ydp82770b41msonormal, div.gmail-m-1834111884996008056ydp82770b41msonormal
        {mso-style-name:gmail-m_-1834111884996008056ydp82770b41msonormal;
        mso-margin-top-alt:auto;
        margin-right:0cm;
        mso-margin-bottom-alt:auto;
        margin-left:0cm;
        font-size:11.0pt;
        font-family:"Calibri",sans-serif;}
span.E-mailStijl19
        {mso-style-type:personal;
        color:#1F497D;}
span.E-mailStijl20
        {mso-style-type:personal-compose;
        font-family:"Calibri",sans-serif;
        color:windowtext;}
.MsoChpDefault
        {mso-style-type:export-only;
        font-family:"Calibri",sans-serif;
        mso-fareast-language:EN-US;}
@page WordSection1
        {size:612.0pt 792.0pt;
        margin:70.85pt 70.85pt 70.85pt 70.85pt;}
div.WordSection1
        {page:WordSection1;}
--></style></head><body lang=NL link=blue vlink=purple><div class=WordSection1><p class=MsoNormal><span lang=EN-US style='font-size:14.0pt;color:#1F497D;mso-fareast-language:EN-US'>Dear Makoto<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:14.0pt;color:#1F497D;mso-fareast-language:EN-US'><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:14.0pt;color:#1F497D;mso-fareast-language:EN-US'>These are wise words born from scientific humbleness and relentless search for a deeper truth. I look forward to hear your further views on ICA as I am to very worried about phase contortion in reconstruction of signals from a limited set of components. I am aware that most researchers on the EEGLab list are working with 256 or more high density channel EEG recordings but we must also  be aware that in the clinical scene 19 channels are conventionally used and reconstruction after rejection of 4-5 components is of course quite a different story than rejecting one component in a 256 ch recording. New techniques will at the end of the day pop up in the clinical theater and I think it is wise that controversies are straightened out or at least be documented and discussed so that the necessary caveats can be drawn. <o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:14.0pt;color:#1F497D;mso-fareast-language:EN-US'><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:14.0pt;color:#1F497D;mso-fareast-language:EN-US'>That being said I wonder a bit about your remark on EEG models. Are you not a bit too pessimistic? As a clinician and neurologist, I had the impression that after all this time and basic scientific work we do have quite clear ideas about the generating mechanisms in pyramidal cells embedded in cortical columns, synaptic channel dynamics and its equations, dendritic and axonal conduction and its correlation with the electrical activity recorded at the scalp.  <o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:14.0pt;color:#1F497D;mso-fareast-language:EN-US'><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:14.0pt;color:#1F497D;mso-fareast-language:EN-US'>What is missing? (or what am I missing?) <o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:14.0pt;color:#1F497D;mso-fareast-language:EN-US'><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:14.0pt;color:#1F497D;mso-fareast-language:EN-US'>Sincerely<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:14.0pt;color:#1F497D;mso-fareast-language:EN-US'><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:14.0pt;color:#1F497D;mso-fareast-language:EN-US'>Georges<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:14.0pt;color:#1F497D;mso-fareast-language:EN-US'><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:14.0pt;color:#1F497D;mso-fareast-language:EN-US'><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:14.0pt;color:#1F497D;mso-fareast-language:EN-US'>PS My latest mail to this list (dated start of July) concerning the transit from EEG (or MEG) channel recordings to network parameters in a graph theory model  was not published. Was something wrong with it?<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:14.0pt;color:#1F497D;mso-fareast-language:EN-US'><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-US style='color:#1F497D;mso-fareast-language:EN-US'><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-US style='color:#1F497D;mso-fareast-language:EN-US'><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-US style='color:#1F497D;mso-fareast-language:EN-US'><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-US style='color:#1F497D;mso-fareast-language:EN-US'><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-US style='color:#1F497D;mso-fareast-language:EN-US'><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-US style='color:#1F497D;mso-fareast-language:EN-US'><o:p> </o:p></span></p><p class=MsoNormal><b><span lang=EN-US>From:</span></b><span lang=EN-US> Makoto Miyakoshi [mailto:mmiyakoshi@ucsd.edu] <br><b>Sent:</b> Wednesday, July 26, 2017 3:01 AM<br><b>To:</b> Rob</span>ert Thatcher <rwthatcher2@yahoo.com><br><b>Cc:</b> EEGLAB List <eeglablist@sccn.ucsd.edu>; Georges Otte <georges.otte@telenet.be><br><b>Subject:</b> Re: How phase is calculated in linear decomposition<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><div><div><p class=MsoNormal>Dear Robert,<o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>Sorry for belated response.<o:p></o:p></p></div><div><p class=MsoNormal>I updated my wiki page to explain how phase is synthesized and split when components are rejected.<o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal><a href="https://sccn.ucsd.edu/wiki/How_phase_is_calculated_in_linear_decomposition">https://sccn.ucsd.edu/wiki/How_phase_is_calculated_in_linear_decomposition</a><o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>As I explained above, you are calling the difference between 18 Bitcoins and 19 Bitcoins 'distortion'.<o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><p class=MsoNormal>> In addition to mathematics, there is empirical proof of “distortion” or “adulteration” that has been repeated over and over again on this forum and accepted to be true by Arnaud and Stephan and Georges and Robert Larson and Ramesh but apparently not yet by you.<o:p></o:p></p><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>I don't think Arno and Stefan accepted 'the proof' because they will write a rebuttal paper with me. It's a good opportunity for us anyways to refute the Montefucso-Siegmunt paper.<o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>By the way Robert, I like your claim that people tend to be trapped by their own dogma and make errors in circular logic. That's a good point. I believe that this criticism actually applies to the entire EEG research field, since we do not have a <i>ground truth</i> in generative model of EEG signals; and do we not tend to compensate the lack of the ground truth by using fancy signal processings? All computational EEG analysts need a good introspection on this issue. We also need to know the difference between engineering and science, models and truth, etc; Good engineering is not necessarily a good science!<o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>I also enjoy and sympathize your criticism against dogmatic people. However, your criticism against ICA is invalid and unsatisfactory to me. Let me invite you to the next thread to share my real criticisms against ICA. <o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>Makoto<o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p><div><p class=MsoNormal>On Fri, Jul 7, 2017 at 11:07 AM, Robert Thatcher <<a href="mailto:rwthatcher2@yahoo.com" target="_blank">rwthatcher2@yahoo.com</a>> wrote:<o:p></o:p></p><blockquote style='border:none;border-left:solid #CCCCCC 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-right:0cm'><div><div><div><p class=gmail-m-1834111884996008056ydp82770b41msonormal><span style='font-size:12.0pt;font-family:"Times New Roman",serif'>Dear Makoto,</span><span style='font-size:12.0pt;font-family:"Helvetica",sans-serif'><o:p></o:p></span></p><p class=gmail-m-1834111884996008056ydp82770b41msonormal><span style='font-size:12.0pt;font-family:"Times New Roman",serif'>      It is correct that a linear decomposition where one decomposes a 19 channel time series into 19 ICA components then, like the inverse Fourier transform, one can recover the phase information if one uses all 19 ICA components.   However, one cannot recover the phase information in 19 channels if one uses 18 or 17 or 16 or any number less than 19 components to reconstruct a new 19 channel time series.  </span><span style='font-size:12.0pt;font-family:"Helvetica",sans-serif'><o:p></o:p></span></p><p class=gmail-m-1834111884996008056ydp82770b41msonormal><span style='font-size:12.0pt;font-family:"Times New Roman",serif'>The mathematics proving “distortion” by ICA reconstruction with a small number of components is explained by Taken’s Theorem (& other theorems in differential geometry, e.g. Shannon's) where a lower dimensional time series does not preserve phase when used to reconstruct a higher dimensional time series.  Therefore mathematically (including the inverse Fourier transform) there is always a loss of information when using a smaller number of components to reconstruct a larger number of channels. Also, commonsense tells one that it is not possible to create something out of nothing, e.g., if I give you 18 Bitcoins will you give me 19 Bitcoins? (credit to Georges for asking this question).</span><span style='font-size:12.0pt;font-family:"Helvetica",sans-serif'><o:p></o:p></span></p><p class=gmail-m-1834111884996008056ydp82770b41msonormal><span style='font-size:12.0pt;font-family:"Times New Roman",serif'> </span><span style='font-size:12.0pt;font-family:"Helvetica",sans-serif'><o:p></o:p></span></p><p class=gmail-m-1834111884996008056ydp82770b41msonormal><span style='font-size:12.0pt;font-family:"Times New Roman",serif'>In addition to mathematics, there is empirical proof of “distortion” or “adulteration” that has been repeated over and over again on this forum and accepted to be true by Arnaud and Stephan and Georges and Robert Larson and Ramesh but apparently not yet by you.   The empirical proof can be demonstrated by yourself if you compute the phase differences between all channel combinations in the original EEG time series vs the ICA reconstructed time series using a smaller number of ICA components.  Anyone can demonstrate the phase difference “distortion” for themselves using the Hilbert transform for instantaneous phase differences (i.e., each time point) or the FFT for an average.  </span><span style='font-size:12.0pt;font-family:"Helvetica",sans-serif'><o:p></o:p></span></p><p class=gmail-m-1834111884996008056ydp82770b41msonormal><span style='font-size:12.0pt;font-family:"Times New Roman",serif'>Here is a URL of a You Tube Video that provides additional empirical proof that ICA reconstruction distorts or adulterates the physiologically based phase differences in the original EEG recording.</span><span style='font-size:12.0pt;font-family:"Helvetica",sans-serif'><o:p></o:p></span></p><p class=gmail-m-1834111884996008056ydp82770b41msonormal><span style='font-size:12.0pt;font-family:"Times New Roman",serif'><a href="https://youtu.be/Q36ojib5OZE" target="_blank">https://youtu.be/Q36ojib5OZE</a></span><span style='font-size:12.0pt;font-family:"Helvetica",sans-serif'><o:p></o:p></span></p><p class=gmail-m-1834111884996008056ydp82770b41msonormal><span style='font-size:12.0pt;font-family:"Times New Roman",serif'> </span><span style='font-size:12.0pt;font-family:"Helvetica",sans-serif'><o:p></o:p></span></p><p class=gmail-m-1834111884996008056ydp82770b41msonormal><span style='font-size:12.0pt;font-family:"Times New Roman",serif'>Here is a url to four power points of two different time points that compares the phase difference before and after Arnaud’s ICA reconstruction.  </span><span style='font-size:12.0pt;font-family:"Helvetica",sans-serif'><o:p></o:p></span></p><p class=gmail-m-1834111884996008056ydp82770b41msonormal><span style='font-size:12.0pt;font-family:"Times New Roman",serif'><a href="http://www.appliedneuroscience.com/Phase_Differences-between-Original_&_Delorme-Post-ICA.pptx" target="_blank">http://www.appliedneuroscience.com/Phase_Differences-between-Original_&_Delorme-Post-ICA.pptx</a></span><span style='font-size:12.0pt;font-family:"Helvetica",sans-serif'><o:p></o:p></span></p><p class=gmail-m-1834111884996008056ydp82770b41msonormal><span style='font-size:12.0pt;font-family:"Times New Roman",serif'>Please download and expand the power points so that you can verify for yourself that the phase differences in the original recording have been distorted or altered.  I would be happy to issue a temporary license to Neuroguide so that you can launch two Neuroguides and import and compare the original time series and Arnaud’s ICA reconstructed time series to verify for yourself that phase differences for all channel combinations for each and every time  point have been distorted by ICA reconstruction.</span><span style='font-size:12.0pt;font-family:"Helvetica",sans-serif'><o:p></o:p></span></p><p class=gmail-m-1834111884996008056ydp82770b41msonormal><span style='font-size:12.0pt;font-family:"Times New Roman",serif'> </span><span style='font-size:12.0pt;font-family:"Helvetica",sans-serif'><o:p></o:p></span></p><p class=gmail-m-1834111884996008056ydp82770b41msonormal><span style='font-size:12.0pt;font-family:"Times New Roman",serif'>I offer this same opportunity to all of the members of this forum.  Download, install and launch NG-2.9.4 and paste the key A to <a href="mailto:qeeg@appliedneuroscience.com" target="_blank">qeeg@appliedneuroscience.com</a>.   Here is the download url:</span><span style='font-size:12.0pt;font-family:"Helvetica",sans-serif'><o:p></o:p></span></p><p class=gmail-m-1834111884996008056ydp82770b41msonormal><span style='font-size:12.0pt;font-family:"Times New Roman",serif'><a href="http://www.appliedneuroscience.com/Download_NeuroGuide.htm" target="_blank">http://www.appliedneuroscience.com/Download_NeuroGuide.htm</a></span><span style='font-size:12.0pt;font-family:"Helvetica",sans-serif'><o:p></o:p></span></p><p class=gmail-m-1834111884996008056ydp82770b41msonormal><span style='font-size:12.0pt;font-family:"Times New Roman",serif'>I will mail step by step mouse instructions so that anyone interested can verify for themselves that ICA reconstruction distorts the phase differences for all channel combinations and all time points.</span><span style='font-size:12.0pt;font-family:"Helvetica",sans-serif'><o:p></o:p></span></p><p class=gmail-m-1834111884996008056ydp82770b41msonormal><span style='font-size:12.0pt;font-family:"Times New Roman",serif'>This is important because if the primary level of measurement is distorted then all subsequent analyses are invalid as a matter of verifiable fact.</span><span style='font-size:12.0pt;font-family:"Helvetica",sans-serif'><o:p></o:p></span></p><p class=gmail-m-1834111884996008056ydp82770b41msonormal><span style='font-size:12.0pt;font-family:"Times New Roman",serif'>Best wishes,</span><span style='font-size:12.0pt;font-family:"Helvetica",sans-serif'><o:p></o:p></span></p><p class=gmail-m-1834111884996008056ydp82770b41msonormal><span style='font-size:12.0pt;font-family:"Helvetica",sans-serif'><o:p> </o:p></span></p><p class=gmail-m-1834111884996008056ydp82770b41msonormal><span style='font-size:12.0pt;font-family:"Times New Roman",serif'>Robert</span><span style='font-size:12.0pt;font-family:"Helvetica",sans-serif'><o:p></o:p></span></p><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Helvetica",sans-serif'><o:p> </o:p></span></p></div><div><div><div><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Helvetica",sans-serif'><o:p> </o:p></span></p></div><div><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Helvetica",sans-serif'><o:p> </o:p></span></p></div><div id="gmail-m_-1834111884996008056yahoo_quoted_0378714856"><div><div><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Helvetica",sans-serif;color:#26282A'>On Thursday, July 6, 2017, 2:31:56 PM PDT, Makoto Miyakoshi <<a href="mailto:mmiyakoshi@ucsd.edu" target="_blank">mmiyakoshi@ucsd.edu</a>> wrote:<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Helvetica",sans-serif;color:#26282A'><o:p> </o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Helvetica",sans-serif;color:#26282A'><o:p> </o:p></span></p></div><div><div id="gmail-m_-1834111884996008056yiv5626669134"><div><div><div><div><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Helvetica",sans-serif;color:#26282A'>Dear Robert and list,<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Helvetica",sans-serif;color:#26282A'><o:p> </o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Helvetica",sans-serif;color:#26282A'>If I remember correctly, one of the critical issue was how phase changes in 'channel data' after removing independent components. Let me focus on this point only.<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Helvetica",sans-serif;color:#26282A'><o:p> </o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Helvetica",sans-serif;color:#26282A'>1. To show my rebuttal and mathematical proof, I prepared a wiki page. In this page. I described how a phase is split/synthesized before and after linear decomposition. Note that this applies to any linear decomposition and not specific to ICA. In doing so, made it nice and educational to show my respect to Arno's effort. Hopefully this mathematical clarification helps.<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Helvetica",sans-serif;color:#26282A'><o:p> </o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Helvetica",sans-serif;color:#26282A'>          <a href="https://sccn.ucsd.edu/wiki/How_phase_is_calculated_in_linear_decomposition" target="_blank">https://sccn.ucsd.edu/wiki/How_phase_is_calculated_in_linear_decomposition</a><o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Helvetica",sans-serif;color:#26282A'><o:p> </o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Helvetica",sans-serif;color:#26282A'>Robert, if you do not agree,<b> please show the right process in a form of math</b>.<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Helvetica",sans-serif;color:#26282A'><o:p> </o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Helvetica",sans-serif;color:#26282A'>2. Apart from this scientific/mathematical dispute, I insist that one should not use the word 'distortion' in this context because it includes subjective judgement.<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Helvetica",sans-serif;color:#26282A'><o:p> </o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Helvetica",sans-serif;color:#26282A'>Makoto<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Helvetica",sans-serif;color:#26282A'><o:p> </o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Helvetica",sans-serif;color:#26282A'><o:p> </o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Helvetica",sans-serif;color:#26282A'><o:p> </o:p></span></p></div></div><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Helvetica",sans-serif;color:#26282A'>-- <o:p></o:p></span></p><div><div><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Helvetica",sans-serif;color:#26282A'>Makoto Miyakoshi<br>Swartz Center for Computational Neuroscience<br>Institute for Neural Computation, University of California San Diego<o:p></o:p></span></p></div></div></div></div></div></div></div></div></div></div></div></div></blockquote></div><p class=MsoNormal><br><br clear=all><o:p></o:p></p><div><p class=MsoNormal><o:p> </o:p></p></div><p class=MsoNormal>-- <o:p></o:p></p><div><div><p class=MsoNormal>Makoto Miyakoshi<br>Swartz Center for Computational Neuroscience<br>Institute for Neural Computation, University of California San Diego<o:p></o:p></p></div></div></div></div></div></body></html>