<META HTTP-EQUIV="Content-Type" CONTENT="text/html; charset=iso-8859-1">
<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 name=Generator content="Microsoft Word 12 (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;}
@font-face
{font-family:Tahoma;
panose-1:2 11 6 4 3 5 4 4 2 4;}
@font-face
{font-family:Verdana;
panose-1:2 11 6 4 3 5 4 4 2 4;}
@font-face
{font-family:Consolas;
panose-1:2 11 6 9 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman","serif";
color:black;}
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;}
pre
{mso-style-priority:99;
mso-style-link:"HTML Preformatted Char";
margin:0cm;
margin-bottom:.0001pt;
font-size:10.0pt;
font-family:"Courier New";
color:black;}
span.HTMLPreformattedChar
{mso-style-name:"HTML Preformatted Char";
mso-style-priority:99;
mso-style-link:"HTML Preformatted";
font-family:Consolas;
color:black;}
span.EmailStyle19
{mso-style-type:personal-reply;
font-family:"Verdana","sans-serif";
color:#1F497D;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:70.85pt 70.85pt 70.85pt 70.85pt;}
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 bgcolor=white lang=EN-US link=blue vlink=purple><div class=WordSection1><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Verdana","sans-serif";color:#1F497D'>Hi Ramin and Jörn,<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Verdana","sans-serif";color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Verdana","sans-serif";color:#1F497D'>In multiple-subject study, it is straightforward to statistically test differences in imaginary coherence using Fieldtrip. This is analogous to the statistical testing of regular coherence, which I outlined in a reply to Ana.<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Verdana","sans-serif";color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Verdana","sans-serif";color:#1F497D'>For a single-subject study, you would have to write your own statfun. This requires some skill, but it is conceptually straightforward. Importantly, for any test statistic that you would like to use for comparing imaginary coherence, you do not(!) have to know it’s parametric reference distribution (thanks to the beauty of permutation inference (-;)<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Verdana","sans-serif";color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Verdana","sans-serif";color:#1F497D'>Best,<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Verdana","sans-serif";color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Verdana","sans-serif";color:#1F497D'>Eric Maris<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Verdana","sans-serif";color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Verdana","sans-serif";color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Verdana","sans-serif";color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Verdana","sans-serif";color:#1F497D'><o:p> </o:p></span></p><div style='border:none;border-left:solid blue 1.5pt;padding:0cm 0cm 0cm 4.0pt'><div><div style='border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm'><p class=MsoNormal><b><span style='font-size:10.0pt;font-family:"Tahoma","sans-serif";color:windowtext'>From:</span></b><span style='font-size:10.0pt;font-family:"Tahoma","sans-serif";color:windowtext'> "Jörn M. Horschig" [mailto:jm.horschig@donders.ru.nl] <br><b>Sent:</b> vrijdag 5 april 2013 16:35<br><b>To:</b> FieldTrip discussion list<br><b>Subject:</b> Re: [FieldTrip] Imaginary part of coherence, 'indepsamplesZcoh'<o:p></o:p></span></p></div></div><p class=MsoNormal><o:p> </o:p></p><div><p class=MsoNormal>Hi Ramin,<br><br>hope you are doing fine. For the imaginary part of coherence you need to call ft_connectivityanalysis with cfg.output = 'imag'.<br>Afaik, z-scoring for the imaginary part is different from doing it for the plain coherence, so you should not use indepsamplesZcoh. You could have a look at Guido Nolte's paper on imaginary coherence and look it up in the methods section how he is computing it - similar but different from what I know how the z-scoring for coherence is computed. Thus, I am not sure how to best do cluster-based statistics on that (that means, I know how you could do it, but not what the best way to do it would be). Of course, you write an own statistics function, something like 'indepsamplesZicoh'... but that would be tons of work.<br><br>Good luck! I hope anyone out there has other/better ideas :)<br><br>Best,<br>Jörn<br><br>On 4/5/2013 9:49 AM, Ramin Azodi wrote:<o:p></o:p></p></div><blockquote style='margin-top:5.0pt;margin-bottom:5.0pt'><div><div><div><div><div><p class=MsoNormal>Hello,<br>In case that I want just to calculate imaginary part of coherence instead of coherence with 'indepsamplesZcoh', what kind of operation should I do for that? can anybody lead me through the right way ?<o:p></o:p></p></div><p class=MsoNormal style='margin-bottom:12.0pt'>You can find my code as following,<br><br>cfg = [];<br>cfg.method = 'mtmfft';<br>cfg.output = 'fourier';<br>cfg.taper = 'hanning';<br>cfg.foilim = [12 20];<br>freq1 = ft_freqanalysis(cfg, data1);<br>cfg = [];<br>cfg.method = 'mtmfft';<br>cfg.output = 'fourier';<br>cfg.complex = 'imag';<br>cfg.taper = 'hanning';<br>cfg.foilim = [12 20];<br>freq2 = ft_freqanalysis(cfg, data2);<br>cfg = [];<br>cfg.method = 'montecarlo';<br>cfg.statistic = 'indepsamplesZcoh';<br>cfg.parameter = 'fourierspctrm';<br>cfg.computecritval = 'yes';<br>cfg.numrandomization = 500;<br>cfg.design = [ones(size(freq1.fourierspctrm,1),1);<br>2*ones(size(freq2.fourierspctrm,1),1)]';<br>cfg.label = freq.label;<br>stat = ft_freqstatistics(cfg, freq,freq1);<o:p></o:p></p></div><p class=MsoNormal>Thanks in advanced for your consideration.<o:p></o:p></p></div><p class=MsoNormal>Best Regards,<o:p></o:p></p></div><p class=MsoNormal>Ramin Azodi<o:p></o:p></p><div><div><div><div><p class=MsoNormal style='margin-bottom:12.0pt'><br><br><br><br><o:p></o:p></p></div></div></div></div></div><p class=MsoNormal><br><br><br><o:p></o:p></p><pre>_______________________________________________<o:p></o:p></pre><pre>fieldtrip mailing list<o:p></o:p></pre><pre><a href="mailto:fieldtrip@donders.ru.nl">fieldtrip@donders.ru.nl</a><o:p></o:p></pre><pre><a href="http://mailman.science.ru.nl/mailman/listinfo/fieldtrip">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip</a><o:p></o:p></pre></blockquote><p class=MsoNormal><br><br><br><o:p></o:p></p><pre>-- <o:p></o:p></pre><pre>Jörn M. Horschig<o:p></o:p></pre><pre>PhD Student<o:p></o:p></pre><pre>Donders Institute for Brain, Cognition and Behaviour <o:p></o:p></pre><pre>Centre for Cognitive Neuroimaging<o:p></o:p></pre><pre>Radboud University Nijmegen <o:p></o:p></pre><pre>Neuronal Oscillations Group<o:p></o:p></pre><pre>FieldTrip Development Team<o:p></o:p></pre><pre><o:p> </o:p></pre><pre>P.O. Box 9101<o:p></o:p></pre><pre>NL-6500 HB Nijmegen<o:p></o:p></pre><pre>The Netherlands<o:p></o:p></pre><pre><o:p> </o:p></pre><pre>Contact:<o:p></o:p></pre><pre>E-Mail: <a href="mailto:jm.horschig@donders.ru.nl">jm.horschig@donders.ru.nl</a><o:p></o:p></pre><pre>Tel: +31-(0)24-36-68493<o:p></o:p></pre><pre>Web: <a href="http://www.ru.nl/donders">http://www.ru.nl/donders</a><o:p></o:p></pre><pre><o:p> </o:p></pre><pre>Visiting address:<o:p></o:p></pre><pre>Trigon, room 2.30<o:p></o:p></pre><pre>Kapittelweg 29<o:p></o:p></pre><pre>NL-6525 EN Nijmegen<o:p></o:p></pre><pre>The Netherlands<o:p></o:p></pre></div></div></body></html>