<html xmlns:v="urn:schemas-microsoft-com:vml" 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=us-ascii"><meta name=Generator content="Microsoft Word 15 (filtered medium)"><style><!--
/* Font Definitions */
@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:0in;
        margin-bottom:.0001pt;
        font-size:12.0pt;
        font-family:"Times New Roman","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
        {mso-style-priority:99;
        mso-margin-top-alt:auto;
        margin-right:0in;
        mso-margin-bottom-alt:auto;
        margin-left:0in;
        font-size:12.0pt;
        font-family:"Times New Roman","serif";}
span.hoenzb
        {mso-style-name:hoenzb;}
span.EmailStyle18
        {mso-style-type:personal-reply;
        font-family:"Calibri","sans-serif";
        color:#1F497D;}
.MsoChpDefault
        {mso-style-type:export-only;}
@page WordSection1
        {size:8.5in 11.0in;
        margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
        {page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]--></head><body lang=EN-US link=blue vlink=purple><div class=WordSection1><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>I would like step in and add more comments which may be helpful (hopefully):<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><a name="_MailEndCompose"><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>The assumption of ICA is : The observed data is the sum of a set of inputs which have been mixed together in an unknown fashion and the aim of ICA is to discover both the inputs and how they were mixed. So, after ICA we have some sources which are temporally independent. In other words, they are independent at time t  McKeown, et al. (1998)<o:p></o:p></span></a></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>However and based on Clive Granger talk at 2003 Nobel Laureate in Economics “The basic "Granger Causality" definition is quite simple. Suppose that we have three terms, X<sub>t</sub>, Y<sub>t</sub>, and W<sub>t</sub>, and that we first attempt to forecast X<sub>t+1</sub> using past terms of Y<sub>t</sub> and W<sub>t</sub>. We then try to forecast X<sub>t+1</sub> using past terms of X<sub>t</sub>, Y<sub>t</sub>, and W<sub>t</sub>. If the second forecast is found to be more successful, according to standard cost functions, then the past of Y appears to contain information helping in forecasting X<sub>t+1</sub> that is not in past X<sub>t</sub> or W<sub>t. </sub>… Thus, Y<sub>t</sub> would "Granger cause" X<sub>t+1</sub> if (a) Y<sub>t</sub> occurs before X<sub>t+1</sub> ; and (b) it contains information useful in forecasting X<sub>t+1</sub> that is not found in a group of other appropriate variables.”  So, in Granger causality we try to relate time t+1 to t.<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>So, ICA and Granger causality are not contradicting each other and finding causality btw sources would not have anything to do with source space or channel space data. In my point of view, using ICA and source signal for Granger causality is good because you do not have to worry about the volume conductance problem. However, one can apply Granger causality in the channel space as well since the dipole localization has its own limitations. One clue code be transforming the channel space data to  current source density (CSD) format and then applying any causality/connectivity analysis you would like to study.<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>Best<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>Iman <o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><b><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>-------------------------------------------------------------<o:p></o:p></span></b></p><p class=MsoNormal><b><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>Iman M.Rezazadeh, Ph.D<o:p></o:p></span></b></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>Research Fellow<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>Semel Intitute, UCLA , Los Angeles<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>& Center for Mind and Brain, UC DAVIS, Davis<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><b><span style='font-size:11.0pt;font-family:"Calibri","sans-serif"'>From:</span></b><span style='font-size:11.0pt;font-family:"Calibri","sans-serif"'> eeglablist-bounces@sccn.ucsd.edu [mailto:eeglablist-bounces@sccn.ucsd.edu] <b>On Behalf Of </b>Makoto Miyakoshi<br><b>Sent:</b> Tuesday, February 18, 2014 3:54 PM<br><b>To:</b> mullen.tim@gmail.com<br><b>Cc:</b> eeglablist@sccn.ucsd.edu<br><b>Subject:</b> Re: [Eeglablist] Two step source connectivity analysis (as implemented in SIFT)<o:p></o:p></span></p><p class=MsoNormal><o:p> </o:p></p><div><p class=MsoNormal>Dear Tim,<o:p></o:p></p><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>Why don't you comment on the following question: If independent components are truly independent, how do causality analyses work?<o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>Dear Joe,<o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>Your inputs are too difficult for me to understand. In short, are you saying causality analysis works on independent components because they are not completely independent?<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><div><p class=MsoNormal style='margin-bottom:12.0pt'><o:p> </o:p></p><div><p class=MsoNormal>2014-02-18 15:46 GMT-08:00 Makoto Miyakoshi <<a href="mailto:mmiyakoshi@ucsd.edu" target="_blank">mmiyakoshi@ucsd.edu</a>>:<o:p></o:p></p><blockquote style='border:none;border-left:solid #CCCCCC 1.0pt;padding:0in 0in 0in 6.0pt;margin-left:4.8pt;margin-right:0in'><div><p class=MsoNormal>Dear Bethel,<o:p></o:p></p><div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>> say A=sunrise and B=ice-cream-sale, then the ICA in EEGLAB should find that A is maximally  temporaly independent from B.<o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div></div><div><p class=MsoNormal>ICA would find a correlation between sunrise and ice-cream-sale.<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><p class=MsoNormal>2014-02-10 4:57 GMT-08:00 Bethel Osuagwu <<a href="mailto:b.osuagwu.1@research.gla.ac.uk" target="_blank">b.osuagwu.1@research.gla.ac.uk</a>>:<o:p></o:p></p><div><div><p class=MsoNormal><o:p> </o:p></p><blockquote style='border:none;border-left:solid #CCCCCC 1.0pt;padding:0in 0in 0in 6.0pt;margin-left:4.8pt;margin-right:0in'><p class=MsoNormal>Hi<br>I am not an expert but I just want to give my own opinion!<br><br>I do not think that temporal independence of two variables (A and B) violets causality between them as implemented in SIFT. In fact if  say A=sunrise and B=ice-cream-sale, then the ICA in EEGLAB should find that A is maximally  temporaly independent from B. However we know there is causal flow from A to B.<br><br>This is what I think, but I wait to be corrected so that I can learn!<br><br>Thanks<br>Bethel<br>________________________________________<br>From: <a href="mailto:eeglablist-bounces@sccn.ucsd.edu" target="_blank">eeglablist-bounces@sccn.ucsd.edu</a> [<a href="mailto:eeglablist-bounces@sccn.ucsd.edu" target="_blank">eeglablist-bounces@sccn.ucsd.edu</a>] On Behalf Of IMALI THANUJA HETTIARACHCHI [<a href="mailto:ith@deakin.edu.au" target="_blank">ith@deakin.edu.au</a>]<br>Sent: 07 February 2014 01:27<br>To: <a href="mailto:mullen.tim@gmail.com" target="_blank">mullen.tim@gmail.com</a><br>Cc: <a href="mailto:eeglablist@sccn.ucsd.edu" target="_blank">eeglablist@sccn.ucsd.edu</a><br>Subject: [Eeglablist] Two step source connectivity analysis (as implemented     in SIFT)<o:p></o:p></p><div><div><p class=MsoNormal><br>Hi Tim and the list,<br><br>I am just in need of a clarification regarding the ICA source reconstruction and the subsequent MVAR –based effective connectivity analysis using the components, which is the basis of the SIFT toolbox. I was trying to use this approach in my work but was questioned on the validity using ICA and subsequent MVAR analysis by my colleagues.<br><br>“When using independent component analysis (ICA), we assume the mutual independence<br>of underlying sources, however when we try to estimate connectivity between EEG sources,<br>we implicitly assume that the sources may be  influenced by each other. This contradicts the<br>fundamental assumption of mutual independence between sources in ICA [Cheung et al., 2010, Chiang et al., 2012, Haufe et al., 2009 ]. “<br><br>So due to this reason different approaches such as MVARICA, CICAAR(convolution ICA+MVAR),  SCSA and state space-based methods have been proposed as ICA+MVAR based source connectivity analysis techniques.<br><br><br>·         So, how would you support the valid use of SIFT ( ICA+MVAR as a two-step procedure) for the source connectivity analysis?<br><br><br>·         If I argue that I do not assume independent sources but rely on the fact that ICA will decompose the EEG signals and output ‘maximally independent’ sources and then, I subsequently model for the dependency, will you agree with me? How valid would my argument be?<br><br>It would be really great to see different thoughts and opinions.<br><br>Kind regards<br><br>Imali<br><br><br>Dr. Imali Thanuja Hettiarachchi<br>Researcher<br>Centre for Intelligent Systems research<br>Deakin University, Geelong 3217, Australia.<br><br>Mobile : <a href="tel:%2B61430321972" target="_blank">+61430321972</a><o:p></o:p></p></div></div><p class=MsoNormal>Email: <a href="mailto:ith@deakin.edu.au" target="_blank">ith@deakin.edu.au</a><mailto:<a href="mailto:ith@deakin.edu.au" target="_blank">ith@deakin.edu.au</a>><br>Web :<a href="http://www.deakin.edu.au/cisr" target="_blank">www.deakin.edu.au/cisr</a><<a href="http://www.deakin.edu.au/cisr" target="_blank">http://www.deakin.edu.au/cisr</a>><br><br>[<a href="cid:image001.jpg@01CF23FF.F8259940">cid:image001.jpg@01CF23FF.F8259940</a>]<br><br><br><br><br><br><br>_______________________________________________<br>Eeglablist page: <a href="http://sccn.ucsd.edu/eeglab/eeglabmail.html" target="_blank">http://sccn.ucsd.edu/eeglab/eeglabmail.html</a><br>To unsubscribe, send an empty email to <a href="mailto:eeglablist-unsubscribe@sccn.ucsd.edu" target="_blank">eeglablist-unsubscribe@sccn.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" target="_blank">eeglablist-request@sccn.ucsd.edu</a><o:p></o:p></p></blockquote></div></div></div><p class=MsoNormal><span style='color:#888888'><br><br clear=all><span class=hoenzb><o:p></o:p></span></span></p><div><p class=MsoNormal><o:p> </o:p></p></div><p class=MsoNormal><span class=hoenzb><span style='color:#888888'>-- </span><o:p></o:p></span></p><div><p class=MsoNormal><span style='color:#888888'>Makoto Miyakoshi<br>Swartz Center for Computational Neuroscience<br>Institute for Neural Computation, University of California San Diego</span><o:p></o:p></p></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><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></body></html>