Thank you so much for your answer Jörn,<br><br>Now I got it, so, I will let the 'fourier' as input to ft_connectivityanalysis, because for some indices I will need the auto-spectrum, and if I fully undertood what you said, it will be the same, as they always take the phase information from the fourier transform. And fieltrip will automatically compute the csd if the input is the frequency data.<br><br>Nevertheless, it give me some error when I tried to use the output from 'powandcsd' of freqanalysis, always gives me a matrix dimensions error, and when I use the output 'fourier' from connectivityanalysis, it computes almost all indices to ok, raging from  -1 to 1 (I'm saw that they are signed, so this is ok), the problem are the ImC ranging from -inf to inf, and the PLV values ranging from 0 to inf, but when I use the 'fourier' the values are. Which for me is really weird.<br><br><br>I have checked internally, and even when I always use the 'fourier' as input, when it internally calculates cross-spectrum for the ImC and PLV, the data.crsspctrm is different, e.g: <br>data.dimord = 'rpt_rpt_chan_chan_freq'<br>ans size(data.crsspctrm) = [1 7 59 59]; <br>So, both are apparently the same, but the data inside this matrices is different,<br><br>for example, when calculating the ImC the first 7 columns for the first two rows are: <br>val(:,:,1,1)<br> [0.9323  0.3679  0.8629  1.6436  0.6404  1.8782  0.1139] <br>val(:,:,1,2)<br> [0.8509-0.1555i  0.2507-0.1066i  0.8173-0.0290i  1.4328-0.2427i  0.5457-0.0500i  1.6505-0.0161i  0.0993+0.0016i]<br><br>while when calculating the PLV are:<br>val(:,:,1,1)<br> [1  1  1  1  1  1  1] <br>val(:,:,1,2)<br> [0.9694-0.2133i  0.7011-0.3564i  0.9951-0.0104i  0.9014-0.3481i  0.9781-0.0819i  0.9982-0.0046i  0.9507-0.0483i]<br><br><br>I'm pretty sure I'm doing something wrong here, but I don't really have a clue, <br><br>If I use the fourier transform to calculate the 'csd' in ft_connectivity analysis and use this csd as input to all other ft_connectivityanalysis indices all of them are set to NaN or 1, not ranging between any numbers.<br><br><br><br><br>On the other hand, I have seen searching on the net, that fieltrip is capable of create surrogate data, by shifting the phase information on every trial for only one electrode of each pair while using the original samples of each trials on the second channel, but I couldn't find on the manual how to do such a thing. Is there any way to compute this kind of surrogates for every sensor? or was it something that some one suggest to somebody?<br><br><br>Many thanks in advanced,<br>Gabriel.<br><br><br><br><br>----- Mensaje original -----<br>De: "Jörn M. Horschig" <jm.horschig@donders.ru.nl><br>Fecha: Martes, 21 de Mayo de 2013, 9:03 am<br>Asunto: Re: [FieldTrip] problem in connectivityanalysis<br>A: FieldTrip discussion list <fieldtrip@science.ru.nl><br><br><span><p>        <meta content="text/html; charset=ISO-8859-1" http-equiv="Content-Type">      </p><p>     </p><div class="moz-cite-prefix"><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>Dear Gabriel,<br>       <font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font><br>       <blockquote type="cite"><font size="3"><font size="4"><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>A) As you             see, this makes an average over the epochs of each channel,             and not over the samples (which in this case are the 9             frequencies of the FFT), And I would like to have the op</font><font size="3"><font size="4">posite, </font><font size="3"><font size="4">t<font size="4">his is, to</font> get the </font><font size="3"><font size="4">imaginary part of coherence for                   each epoch (repeti<font size="4">t</font>ion) on all                   my channels and not for each frequency, som<font size="4">ething like: chan_chan_rpt or chancbm_rpt</font>,                 </font><font size="3"><font size="4">Is this possible?</font><br>                 </font></font></font></font></font></blockquote>       <font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font><br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>       No. You should think of coherence more as a measure for the       consistency of the phase relation between channels across trials,       so coherence is just not defined per trial,neither would it make       sense to compute coherence for a single trial.<br>       <font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font><br>       <blockquote type="cite"><font size="3"><font size="3"><font size="3"><font size="3"><font size="3"><font size="4"><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>B)                     Another question I have is: for the Imaginary part                     of Coherence is not necessary to compute the                     cross-spectrum with 'powandcsd', only the FFT, but                     for the PLV, PPC or WPLI is that way too?</font></font></font></font></font></font></blockquote>       I see some terminology issues here. Let me try to explain in a       simple (aka least-mathematical I can think of) way:<br>       <font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font><br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>       For all these connectivity measures, we need the phase information       from the Fourier transform (FFT) to compute the cross-spectral       density matrix (CSD). The CSD can be regarded as the equivalent in       the frequency domain to the covariance matrix in the time domain,       thus it is a measures how a certain channels activity is       co-modulated with another channels activity (for a particular       frequency) - makes pretty much sense to make use of this when       computing connectivity, right? :)<br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>       As you might know, a Fourier transform returns complex numbers,       where the imaginary part of this number contains phase       information. When calling ft_freqanalysis, you can decide if it       shall only return the (squared) real part of the frequency       spectrum (cfg.output='pow'), or instead the full CSD (cfg.output =       'powandcsd'). Alternatively, you can also let ft_freqanalysis       return the raw Fourier coefficients (cfg.output='fourier'). From       the Fourier coefficients, however, you can easily obtain the CSD       or the power spectrum basically by multiplication the Fourier       matrix with itself (transposed). All of these measures quantify       the phase relation between channels across trials, and as phase       information is coded in the imaginary part of the fourier output,       thus for connectivityanalysis you need the CSD which can be       obtained either by 'powandcsd' or by 'fourier'. <br>       <font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font><br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>       I hope that somehow clarifies your questions.<br>       <font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font><br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>       Best,<br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>       Jörn<br>       <font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font><br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>       On 5/17/2013 5:45 PM, Gabriel Gonzalez Escamilla wrote:<br>     </div>     <blockquote cite="mid:9f89709b303dae2.51966cc4@upo.es" type="cite"><font size="3"><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>Dear Fieltrip experts,<br>         <font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font><br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>         I have my continuous data imported to fieldtrip with matlab, I'm         looking for performing connectivity analysis between pairs of         sensors, I have succed as:<br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>         data.trials {1xNepochs}; % <font size="3">N</font>epochs = 7,         each epoch with 59 channels and 2000 samples<br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>         data.label {59x1};<br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>         data.time {1xNepochs};<br>         <font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font><br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>         then, calculate only the FFT of all sensors as:<br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>         cfg=[];<br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>         cfg.output='fourier';<br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>         cfg.method='mtmfft';<br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>         cfg.taper='hanning';<br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>         cfg.foilim=[8.5 9.5];<br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>         cfg.tapsmofrq=0;<br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>         cfg.trials='all';<br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>         cfg.keeptrials='yes';<br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>         cfg.channel='all';<br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>         fourier=ft_freqanalysis(cfg, data)<br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>         as output I get:<br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>         dimord = 'rpttap_chan_freq'<br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>         freq= [1x9]<br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>         fourierspctrm=[7x59x9 double]<br>         <font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font><br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>         Then compute the ima<font size="3">ginary part of coherency as:<br>           <font size="3"><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>cfg=[];<br>             <font size="3"><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>cfg.method ='cohe'</font>;<br>             <font size="3"><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>cfg.complex='imag';<br>               <font size="3"><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>cfg.channelc<font size="3">bm={'all'</font></font>               'all'};<br>             </font></font></font><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>coher = ft_connect<font size="3">iv<font size="3">ityanalysis </font></font>(cfg, fourier)<br>         <font size="3"><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>as output get: <br>           <font size="3"><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>dimord='chan_chan_freq';<br>             <font size="3"><font size="3"><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>cohspctrm=[59x59x9 double]</font></font>;<br>             <font size="3"><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>dof=7;</font><br>           </font></font><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font><br>         <font size="4"><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>I have two main questions:</font><br>         <font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font><br>         <font size="4"><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>A) As you see, this makes an average over the           epochs of each channel, and not over the samples (which in           this case are the 9 frequencies of the FFT), And I would like           to have the op</font><font size="3"><font size="4">posite, </font><font size="3"><font size="4">t<font size="4">his is, to</font>               get the </font><font size="3"><font size="4">imaginary                 part of coherence for each epoch (repeti<font size="4">t</font>ion)                 on all my channels and not for each frequency, som<font size="4">ething like: chan_chan_rpt or chancbm_rpt</font>,               </font><font size="3"><font size="4">Is this possible?</font><br>                 <font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font><br>                 <font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font><br>                 <font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font><br>                 <font size="4"><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>B) Another question I have is: for the                   Imaginary part of Coherence is not necessary to                   compute the cross-spectrum with 'powandcsd', only the                   FFT, but for the PLV, PPC or WPLI is that way too?</font><br>                 <font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font><br>                 <font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font><br>                 <font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font><br>                 <font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font><br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>                 Many thanks in advanced,<br>                 <font size="3"><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>Gabriel.</font><br>                 <font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font><br>               </font></font></font></font></font>       <font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font><br>       <fieldset class="mimeAttachmentHeader"></fieldset>       <font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font><br>       <pre wrap=""><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>_______________________________________________ > fieldtrip mailing list > <a href="javascript:main.compose('new', 't=fieldtrip@donders.ru.nl')">fieldtrip@donders.ru.nl</a> <a class="moz-txt-link-freetext" href="http://mailman.science.ru.nl/mailman/listinfo/fieldtrip" target="1">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip</a></pre>     </blockquote>     <br>     <font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font><br>     <pre class="moz-signature" cols="72"><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">> </font>--  > Jörn M. Horschig > PhD Student > Donders Institute for Brain, Cognition and Behaviour  > Centre for Cognitive Neuroimaging > Radboud University Nijmegen  > Neuronal Oscillations Group > FieldTrip Development Team >  > P.O. Box 9101 > NL-6500 HB Nijmegen > The Netherlands >  > Contact: > E-Mail: <a href="javascript:main.compose('new', 't=jm.horschig@donders.ru.nl')">jm.horschig@donders.ru.nl</a> > Tel:    +31-(0)24-36-68493 > Web: <a class="moz-txt-link-freetext" href="http://www.ru.nl/donders" target="1">http://www.ru.nl/donders</a> >  > Visiting address: > Trigon, room 2.30 > Kapittelweg 29 > NL-6525 EN Nijmegen > The Netherlands</pre>    </span> > _______________________________________________<br>> fieldtrip mailing list<br>> fieldtrip@donders.ru.nl<br>> http://mailman.science.ru.nl/mailman/listinfo/fieldtrip<br><br><font size="3">--------------------------<br />PhD. student Gabriel González-Escamilla<br />Laboratory of Functional Neuroscience<br />Department of Physiology, Anatomy, and Cell Biology<br />University Pablo de Olavide<br />Ctra. de Utrera, Km.1<br />41013 - Seville<br />- Spain -<br /><br />Email: ggonesc@upo.es<br />http://www.upo.es/neuroaging/es/</font><br>