<div dir="ltr">Dear Frederik,<div><br></div><div>It's such an honor to receive a kind explanation from the first author of the critical paper. Thank you.</div><div><br></div><div>I read your comment but did not see the answer. I read discussion section of your paper and found the following description.</div><div><br></div><div><div>> However, since the mixing of sources is not identical for the different channels, the inclusion of past values of another (driver) channel can improve the prediction of another (target) channel’s current state.</div></div><div><br></div><div>Yeah, this makes sense... I confused sources with channels. After applying MVAR on channel signals, of course the result is across-channel connectivity... who wants to know it.</div><div><br></div><div>However, this means that preprocessing methods to obtain source EEG activity is critical. The author of SIFT, Tim Mullen, participated to write this paper, which claims that it does better job than AMICA.</div><div><a href="https://www.hindawi.com/journals/cin/2016/9754813/">https://www.hindawi.com/journals/cin/2016/9754813/</a><br></div><div>My point is that if the selection of this kind of preprocessing methods and algorithms can determine the basic quality of the connectivity ananlysis, then the EEG connectivity analysis itself if very much dependent on its source-decomposing preprocessing. I kind of knew it, but your publication, together with Clemens's, proved it. This may promote the use of ICA actually.</div><div><br></div><div>Anyway, thank you very much Frederik. I understood it now.</div><div><br></div><div>Makoto</div><div> </div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Dec 7, 2016 at 3:22 AM, Frederik Van de Steen <span dir="ltr"><<a href="mailto:Frederik.VandeSteen@ugent.be" target="_blank">Frederik.VandeSteen@ugent.be</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">




<div dir="ltr" style="font-size:12pt;color:#000000;background-color:#ffffff;font-family:Calibri,Arial,Helvetica,sans-serif">

<p align="LEFT">Dear Makoto,</p>
<p><br>
<br>
</p>
<p>We investigated this issue in our recent paper: </p>
<p>Van de Steen, F., Faes, L., Karahan, E., Songsiri, J., Valdes-Sosa, P. A., & Marinazzo, D. (2016). Critical comments on EEG sensor space dynamical connectivity analysis.
<em>Brain Topography</em>, 1-12.</p>
<p align="JUSTIFY">This work was conducted independently from and simultaneously as
<font color="#000000"><font face="helvetica, sans-serif, ariel"><span style="background:#ffffee">Brunner</span></font></font> et al (2016). Both papers show how spurious connectivity on sensor space analysis can, in many cases, occur (not always though).
</p>
<p>I'll try my very best to make it intuitively:</p>
<p>first of all, you should be aware that non negative (d)DTF, (r)PDC, time domain granger causality etc. all imply at least one non-negative off diagonal coefficient of the coefficient matrix (<b>B</b>) of the MVAR model.</p>
<p>To clarify: time series are modelled with a multivariate autoregressive model (MVAR) in which current values of time series are explained by a linear combination of past values of that time series itself and past values of other time series:
<span style="font-weight:normal">X</span>(t)= <b>B</b>1*<span style="font-weight:normal">X</span>(t-1) +
<b>B</b><span style="font-weight:normal">2*X</span>(t-2) + etc+ error where X(t) is a column vector contain the current values of all time series (X(t) = [X1(t) X2(t) ...]T). Diagonal elements of
<b>B</b> relate past values (t-1, t-2 etc) of a time series to the current value (t) of the time series itself (e.g X1(t-1) to X1(t) )). Off diagonal elements relate past values of another time series to the current value (e.g X2(t-1) to X1(t)))</p>
<p>secondly, two independent time series ≠ two independent time series. what i mean is that there are different ways in which time series can be independent of one another. e.g two white noise processes vs. two independent autoregressive time series. In the
 former case the time series cannot be modelled with past values  of another time series nor can it be modelled with past values of that time series itself (all elements in
<b>B</b><span style="font-weight:normal"> are zero</span>). In the latter case, the situation is different. There you cannot model one time series with past values of the other but you can model it with past values of the time series itself (i.e. only non
 zero elements on diagonal of <b>B</b>). Now when you mix up (linear superposition) these time series (by volume conduction), the past contains information of both original time series.
</p>
<p><a name="m_5072483075482788299_ms-rterangepaste-start"></a>since the mixing of time series is not identical (i.e. each EEG channel is a unique linear combination of sources), the inclusion of past values of another mixed time series can improve the prediction of the other mixed
 time series. This translates, in terms of the MVAR model, in a (some) non-negative off diagonal elements in
<b>B</b> and thus also non negative DTF, PDC, etc. When mixing up white noise, no spurious connectivity will occur because the past wasn't informative in the first place.
</p>
<p>You need to be aware that even though volume conduction implies instantaneously mixing of time series, the continuation of time does not unmix the past...</p>
<p>So basically these measures are not insensitive to volume conduction as claimed by some.
</p>
<p>Hope it is a bit clear. Please do not hesitate to ask more questions</p>
<p><br>
<br>
</p>
<p>Kind regards,</p>
<p><br>
<br>
</p>
<p>Frederik </p>
<p><br>
<br>
</p>
<div style="color:rgb(33,33,33)">
<hr style="display:inline-block;width:98%">
<div id="m_5072483075482788299divRplyFwdMsg" dir="ltr"><font style="font-size:11pt" color="#000000" face="Calibri, sans-serif"><b>Van:</b> Makoto Miyakoshi <<a href="mailto:mmiyakoshi@ucsd.edu" target="_blank">mmiyakoshi@ucsd.edu</a>><br>
<b>Verzonden:</b> woensdag 7 december 2016 03:46<br>
<b>Aan:</b> Scott Makeig<br>
<b>CC:</b> <a href="mailto:eeglablist@sccn.ucsd.edu" target="_blank">eeglablist@sccn.ucsd.edu</a><br>
<b>Onderwerp:</b> Re: [Eeglablist] paper published</font>
<div> </div>
</div><div><div class="h5">
<div>
<div dir="ltr">Dear Scott and list,
<div><br>
</div>
<div>During the EEGLAB Workshop 2016, I asked this question to Tim Mullen during his lecture, without knowing this ongoing debate. Intuitively, it does not make sense to me why just linear mixing affects connectivity calculation, if dDTF or rPDC can suppress
 the spurious connections... Tim mentioned that there is good reason for this, but did not explain it during the lecture due to limited time. Can anyone give me an intuitive explanation why it is bad?</div>
<div><br>
</div>
<div>Makoto</div>
<div><br>
</div>
<div><br>
</div>
<div class="gmail_extra"><br>
<div class="gmail_quote">On Wed, Nov 30, 2016 at 6:02 PM, Scott Makeig <span dir="ltr">
<<a href="mailto:smakeig@ucsd.edu" target="_blank">smakeig@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 dir="ltr">Some of us noticed a paper published recently that claimed that effective connectivity measures between scalp EEG channels suffer no ill effects of volume conduction -- and immediately questioned its conclusions!  Clement Brunner of Graz mounted
 an effort to publish a rebuttal in the same journal, which has now appeared:
<div><br>
</div>
<div><font color="#000000"><font style="font-family:helvetica,sans-serif,ariel;text-align:justify;background-color:rgb(255,255,238)">C. Brunner, M. Billinger, M. Seeber, T.R. Mullen, <font>S Makeig</font>, </font><font style="font-family:helvetica,sans-serif,ariel;text-align:justify;background-color:rgb(255,255,238)">Volume
 conduction influences scalp-based connectivity estimates<a href="https://sccn.ucsd.edu/~scott/pdf/brunner16.pdf" style="text-decoration:none" target="_blank"> </a>(a rebuttal). <em>Frontiers in Computational Neuroscience</em>, </font></font></div>
<div><font color="#000000"><font style="font-family:helvetica,sans-serif,ariel;text-align:justify;background-color:rgb(255,255,238)">doi:10.3389/fncom.2016.00121, 22 November 2016.</font></font><br clear="all">
<div><br>
</div>
<div>We have learned that another group is publishing a separate rebuttal ...<span class="m_5072483075482788299HOEnZb"><font color="#888888"><br>
</font></span></div>
<span class="m_5072483075482788299HOEnZb"><font color="#888888">
<div><br>
</div>
<div>Scott Makeig</div>
<div><br>
</div>
-- <br>
<div class="m_5072483075482788299m_5057240027199685336gmail_signature">Scott Makeig, Research Scientist and Director, Swartz Center for Computational Neuroscience, Institute for Neural Computation, University of California San Diego, La Jolla CA 92093-0961,
<a href="http://sccn.ucsd.edu/~scott" target="_blank">http://sccn.ucsd.edu/~scott</a></div>
</font></span></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/ee<wbr>glabmail.html</a><br>
To unsubscribe, send an empty email to <a href="mailto:eeglablist-unsubscribe@sccn.ucsd.edu" target="_blank">
eeglablist-unsubscribe@sccn.uc<wbr>sd.edu</a><br>
For digest mode, send an email with the subject "set digest mime" to <a href="mailto:eeglablist-request@sccn.ucsd.edu" target="_blank">
eeglablist-request@sccn.ucsd.e<wbr>du</a><br>
</blockquote>
</div>
<br>
<br clear="all">
<div><br>
</div>
-- <br>
<div class="m_5072483075482788299gmail_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>
</div>
</div></div></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><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>