<div dir="ltr"><div><div><div>Hello all,<br><br>I have not yet gotten a response to my question below, but in the meantime I have another question about the minimum norm estimate workflow--specifically, about the coordinate system for the skull-stripped anatomy in the step described at <a href="http://fieldtrip.fcdonders.nl/tutorial/minimumnormestimate#preprocessing_of_the_anatomical_mrisave_to_disk">http://fieldtrip.fcdonders.nl/tutorial/minimumnormestimate#preprocessing_of_the_anatomical_mrisave_to_disk</a>. I'm confused by the following bit of code:<br>
<br>% ensure that the skull-stripped anatomy is expressed in the same coordinate system as the anatomy<br>seg.transform = mri_tal.transform;<br><br></div>In my data, mri_tal.coordsys is 'spm' (I presume this is the result of re-aligning to Talairach in the previous step?) whereas seg.coordsys is 'ctf' (as a result of re-aligning to CTF several steps earlier). (But mri_tal also has a field mri_tal.transformorig, which seg does not have.) So should I really be using the same transform for both, as shown in the tutorial?<br>
<br></div>Apologies if this question is pretty basic; I'm just trying to pinpoint where the mis-alignment described in my message below occurred, so I want to make sure I understand each step of the workflow correctly<br>
<br></div>Best,<br>Steve<br><div><div><div><div><br><br>> Message: 1<br>> Date: Sat, 25 May 2013 08:11:18 -0500<br>> From: Stephen Politzer-Ahles <<a href="mailto:politzerahless@gmail.com">politzerahless@gmail.com</a>><br>
> To: <a href="mailto:fieldtrip@donders.ru.nl">fieldtrip@donders.ru.nl</a><br>> Subject: [FieldTrip] Sourcespace and volume conductor misaligned<br>> Message-ID:<br>>         <CAJT2k_9-hd_sM=hp4P-CUu+=aduSOMyZV7XPG0=<a href="mailto:fFk0ouR0wzA@mail.gmail.com">fFk0ouR0wzA@mail.gmail.com</a>><br>
> Content-Type: text/plain; charset="utf-8"<br>><br>> Hello all,<br>><br>> I am going through the workflow at<br>> <a href="http://fieldtrip.fcdonders.nl/tutorial/minimumnormestimate">http://fieldtrip.fcdonders.nl/tutorial/minimumnormestimate</a>. After making<br>
> the volume conduction model using ft_prepare_headmodel(), I noticed that<br>> although the volume conduction model and sourcespace have the same<br>> orientation and overall size/shape (after I converted the volume conduction<br>
> model to cm, which wasn't in the tutorial but my original model came out in<br>> mm), they don't quite line up, as you can see in this figure:<br>><br>> <a href="http://i.imgur.com/mGEtLOa.png">http://i.imgur.com/mGEtLOa.png</a><br>
><br>> I did interactively re-align the data to CTF (twice--in step 2 of<br>> "Preprocessing of the anatomical MRI" and in step 4 of "Source model")<br>> using fiducials, and to Talairach (step 5 of "Preprocessing of the<br>
> anatomical data"), so I'm not sure how it ended up this way. The code I've<br>> used at each step is basically the same as that in the tutorial.<br>><br>> Is there any way to line up my volume conduction model and sourcespace now,<br>
> without going back and re-running most of the workflow?<br>><br>> Best,<br>> Steve<br>><br>> --<br>> Stephen Politzer-Ahles<br>> University of Kansas<br>> Linguistics Department<br>> <a href="http://people.ku.edu/~sjpa/">http://people.ku.edu/~sjpa/</a><br>
<br><br><br>On Sat, May 25, 2013 at 1:56 PM, <<a href="mailto:fieldtrip-request@science.ru.nl">fieldtrip-request@science.ru.nl</a>> wrote:<br>><br>> Send fieldtrip mailing list submissions to<br>>         <a href="mailto:fieldtrip@science.ru.nl">fieldtrip@science.ru.nl</a><br>
><br>> To subscribe or unsubscribe via the World Wide Web, visit<br>>         <a href="http://mailman.science.ru.nl/mailman/listinfo/fieldtrip">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip</a><br>> or, via email, send a message with subject or body 'help' to<br>
>         <a href="mailto:fieldtrip-request@science.ru.nl">fieldtrip-request@science.ru.nl</a><br>><br>> You can reach the person managing the list at<br>>         <a href="mailto:fieldtrip-owner@science.ru.nl">fieldtrip-owner@science.ru.nl</a><br>
><br>> When replying, please edit your Subject line so it is more specific<br>> than "Re: Contents of fieldtrip digest..."<br>><br>><br>> Today's Topics:<br>><br>>    1. Sourcespace and volume conductor misaligned<br>
>       (Stephen Politzer-Ahles)<br>>    2. Re: fieldtrip Digest, Vol 30, Issue 31 (Johanna Zumer)<br>><br>><br>> ----------------------------------------------------------------------<br>><br>> Message: 1<br>
> Date: Sat, 25 May 2013 08:11:18 -0500<br>> From: Stephen Politzer-Ahles <<a href="mailto:politzerahless@gmail.com">politzerahless@gmail.com</a>><br>> To: <a href="mailto:fieldtrip@donders.ru.nl">fieldtrip@donders.ru.nl</a><br>
> Subject: [FieldTrip] Sourcespace and volume conductor misaligned<br>> Message-ID:<br>>         <CAJT2k_9-hd_sM=hp4P-CUu+=aduSOMyZV7XPG0=<a href="mailto:fFk0ouR0wzA@mail.gmail.com">fFk0ouR0wzA@mail.gmail.com</a>><br>
> Content-Type: text/plain; charset="utf-8"<br>><br>> Hello all,<br>><br>> I am going through the workflow at<br>> <a href="http://fieldtrip.fcdonders.nl/tutorial/minimumnormestimate">http://fieldtrip.fcdonders.nl/tutorial/minimumnormestimate</a>. After making<br>
> the volume conduction model using ft_prepare_headmodel(), I noticed that<br>> although the volume conduction model and sourcespace have the same<br>> orientation and overall size/shape (after I converted the volume conduction<br>
> model to cm, which wasn't in the tutorial but my original model came out in<br>> mm), they don't quite line up, as you can see in this figure:<br>><br>> <a href="http://i.imgur.com/mGEtLOa.png">http://i.imgur.com/mGEtLOa.png</a><br>
><br>> I did interactively re-align the data to CTF (twice--in step 2 of<br>> "Preprocessing of the anatomical MRI" and in step 4 of "Source model")<br>> using fiducials, and to Talairach (step 5 of "Preprocessing of the<br>
> anatomical data"), so I'm not sure how it ended up this way. The code I've<br>> used at each step is basically the same as that in the tutorial.<br>><br>> Is there any way to line up my volume conduction model and sourcespace now,<br>
> without going back and re-running most of the workflow?<br>><br>> Best,<br>> Steve<br>><br>> --<br>> Stephen Politzer-Ahles<br>> University of Kansas<br>> Linguistics Department<br>> <a href="http://people.ku.edu/~sjpa/">http://people.ku.edu/~sjpa/</a><br>
> -------------- next part --------------<br>> An HTML attachment was scrubbed...<br>> URL: <<a href="http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20130525/561df14f/attachment-0001.html">http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20130525/561df14f/attachment-0001.html</a>><br>
><br>> ------------------------------<br>><br>> Message: 2<br>> Date: Sat, 25 May 2013 20:54:59 +0200<br>> From: Johanna Zumer <<a href="mailto:johanna.zumer@donders.ru.nl">johanna.zumer@donders.ru.nl</a>><br>
> To: FieldTrip discussion list <<a href="mailto:fieldtrip@science.ru.nl">fieldtrip@science.ru.nl</a>><br>> Subject: Re: [FieldTrip] fieldtrip Digest, Vol 30, Issue 31<br>> Message-ID:<br>>         <CAL4kA6eTL_TwECF62W-OoWY=<a href="mailto:Y1mpBR6QoOLh%2Bug8AcXNbVxQ1A@mail.gmail.com">Y1mpBR6QoOLh+ug8AcXNbVxQ1A@mail.gmail.com</a>><br>
> Content-Type: text/plain; charset="utf-8"<br>><br>> Dear Jing,<br>><br>> 1.  Every MRI scanner is different, so I'm not familiar with that naming<br>> scheme, but it seems the 'GAD' at the end of the file name indicates that<br>
> was with the Gadolinium contrast agent, so I would suggest to not use that,<br>> but instead use the first (plain) T1.<br>><br>> 2. I'm not sure what to suggest.  I assume you can see the vitamin E marker<br>
> clearly when viewing the DICOM file?   Is the issue that some slices are<br>> missing in the conversion, or that the contrast on the image is not the<br>> same?  or right/left flipped?<br>><br>> Best,<br>
> Johanna<br>><br>><br>> 2013/5/23 WangJing <<a href="mailto:13681530640@139.com">13681530640@139.com</a>><br>><br>> > Dear Johanna,<br>> ><br>> >     Thank you for your suggestion.<br>
> >     1.I have two folder of T1 MRI, one is t1_vibe_tra_4,and another is<br>> > t1_vibe_traGAD-5. which one is best?<br>> >     2.when I convert MRI from DICOM into CTF format,the vitamin E marker<br>> > didn't be found. what can I do?<br>
> ><br>> > Best Regards,<br>> > Jing Wang<br>> ><br>> ><br>> ><br>> ><br>> > ------------------------------<br>> > ----The following is the content of the forwarded email----<br>
> > From?fieldtrip-request<<a href="mailto:fieldtrip-request@science.ru.nl">fieldtrip-request@science.ru.nl</a>><br>> > To?fieldtrip<<a href="mailto:fieldtrip@science.ru.nl">fieldtrip@science.ru.nl</a>><br>
> > Date?2013-05-22 15:33:13<br>> > Subject?fieldtrip Digest, Vol 30, Issue 31<br>> ><br>> > Send fieldtrip mailing list submissions to<br>> > <a href="mailto:fieldtrip@science.ru.nl">fieldtrip@science.ru.nl</a><br>
> ><br>> > To subscribe or unsubscribe via the World Wide Web, visit<br>> > <a href="http://mailman.science.ru.nl/mailman/listinfo/fieldtrip">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip</a><br>
> > or, via email, send a message with subject or body 'help' to<br>> > <a href="mailto:fieldtrip-request@science.ru.nl">fieldtrip-request@science.ru.nl</a><br>> ><br>> > You can reach the person managing the list at<br>
> > <a href="mailto:fieldtrip-owner@science.ru.nl">fieldtrip-owner@science.ru.nl</a><br>> ><br>> > When replying, please edit your Subject line so it is more specific<br>> > than "Re: Contents of fieldtrip digest..."<br>
> ><br>> ><br>> > Today's Topics:<br>> ><br>> >    1. Re: question About head model (Johanna Zumer)<br>> ><br>> ><br>> > ----------------------------------------------------------------------<br>
> ><br>> > Message: 1<br>> > Date: Wed, 22 May 2013 09:30:17 +0200<br>> > From: Johanna Zumer <<a href="mailto:johanna.zumer@donders.ru.nl">johanna.zumer@donders.ru.nl</a>><br>> > To: FieldTrip discussion list <<a href="mailto:fieldtrip@science.ru.nl">fieldtrip@science.ru.nl</a>><br>
> > Subject: Re: [FieldTrip] question About head model<br>> > Message-ID:<br>> > <CAL4kA6fyA18h3eFkBKu0V5+y4eA1rCZuaf8UGTZ1=<a href="mailto:HsAS_X0-w@mail.gmail.com">HsAS_X0-w@mail.gmail.com</a>><br>
> > Content-Type: text/plain; charset="utf-8"<br>> ><br>> > Dear Jing,<br>> ><br>> > It is best to use the T1 MRI rather than the other scans.   To open a DICOM<br>> > format file, you can try one of these options:<br>
> ><br>> > % read the DICOM files<br>> > mri = ft_read_mri('single_file_of_DICOM_series.IMA');<br>> >  % or use a graphical file selection[f, p] = uigetfile<br>> > <<a href="http://www.mathworks.com/access/helpdesk/help/techdoc/ref/uigetfile.html">http://www.mathworks.com/access/helpdesk/help/techdoc/ref/uigetfile.html</a><br>
> > >('*');<br>> > mri = ft_read_mri(fullfile<br>> > <<a href="http://www.mathworks.com/access/helpdesk/help/techdoc/ref/fullfile.html">http://www.mathworks.com/access/helpdesk/help/techdoc/ref/fullfile.html</a><br>
> > >(p,<br>> > f));<br>> ><br>> ><br>> > I copied the above from documentation on the FieldTrip wiki:<br>> ><br>> > <a href="http://fieldtrip.fcdonders.nl/faq/how_can_i_convert_an_anatomical_mri_from_dicom_into_ctf_format">http://fieldtrip.fcdonders.nl/faq/how_can_i_convert_an_anatomical_mri_from_dicom_into_ctf_format</a><br>
> ><br>> ><br>> > Then you can save the 'mri' into a *.mat file for future use with FieldTrip.<br>> ><br>> > Best regards,<br>> > Johanna<br>> ><br>> ><br>> > 2013/5/22 WangJing <<a href="mailto:13681530640@139.com">13681530640@139.com</a>><br>
> ><br>> > > Hi<br>> > > I am a freshman about field trip.<br>> > > I am trying to build the forward model. but my mri data is different from<br>> > > the data in the tutorial. So I don't know how to bulid the head model. I<br>
> > > use the date which is from the paper "Open Database of Epileptic EEG with<br>> > > MRI and Postoperational Assessment of Foci?a Real World Verification for<br>> > > the EEG Inverse Solutions",<a href="http://eeg.pl/epi">http://eeg.pl/epi</a><br>
> > . MRI recordings containing a<br>> ><br>> > > T1, T2 or fluid attenuated inversion recovery (FLAIR) weighted brain scans<br>> > > with morphologic substrate of the epilepsy (mostly cortical dysplasias,<br>
> > > dysplastic tumors etc.). In some cases there is also a scan with the<br>> ><br>> > > gadolinium (GAD) contrast. The images were collected by Siemens Sonata 1.5T<br>> > > scanner. The data is stored in the Digital Imaging and Communications in<br>
> > > Medicine (DICOM) format. The name of the scan?s folder indicates its<br>> > > weighting.The scans have the following resolutions: T1: 512 ? 512 pixels,<br>> > > pixel spacing 0.4687 mm, slice thickness 1.2 mm, T2: 256 ? 256 pixels,<br>
> > > pixel spacing 0.9375 mm,slice thickness 2.5 mm.<br>> ><br>> > > In the tutorial,the mri data stored in only one file, while the data which<br>> ><br>> > > I used have many files, I don't know how to process. who can help me,Thanks1<br>
> > >  Thanks very much!<br>> > > Best Reagards<br>> > > Jing Wang<br>> > ><br>> > ><br>> > ><br>> > ><br>> > ><br>> > > ------------------------------<br>
> > > ----The following is the content of the forwarded email----<br>> > > From?fieldtrip-request<<a href="mailto:fieldtrip-request@science.ru.nl">fieldtrip-request@science.ru.nl</a>><br>> > > To?fieldtrip<<a href="mailto:fieldtrip@science.ru.nl">fieldtrip@science.ru.nl</a>><br>
> > > Date?2013-05-21 23:59:18<br>> > > Subject?fieldtrip Digest, Vol 30, Issue 29<br>> > ><br>> > > Send fieldtrip mailing list submissions to<br>> > > <a href="mailto:fieldtrip@science.ru.nl">fieldtrip@science.ru.nl</a><br>
> > ><br>> > > To subscribe or unsubscribe via the World Wide Web, visit<br>> > > <a href="http://mailman.science.ru.nl/mailman/listinfo/fieldtrip">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip</a><br>
> > > or, via email, send a message with subject or body 'help' to<br>> > > <a href="mailto:fieldtrip-request@science.ru.nl">fieldtrip-request@science.ru.nl</a><br>> > ><br>> > > You can reach the person managing the list at<br>
> > > <a href="mailto:fieldtrip-owner@science.ru.nl">fieldtrip-owner@science.ru.nl</a><br>> > ><br>> > > When replying, please edit your Subject line so it is more specific<br>> > > than "Re: Contents of fieldtrip digest..."<br>
> > ><br>> > ><br>> > > Today's Topics:<br>> > ><br>> > >    1. Re: problem in connectivityanalysis (Gabriel Gonzalez Escamilla)<br>> > ><br>> > ><br>> > > ----------------------------------------------------------------------<br>
> > ><br>> > > Message: 1<br>> > > Date: Tue, 21 May 2013 17:52:21 +0200<br>> > > From: Gabriel Gonzalez Escamilla <<a href="mailto:ggonesc@upo.es">ggonesc@upo.es</a>><br>> > > To: FieldTrip discussion list <<a href="mailto:fieldtrip@science.ru.nl">fieldtrip@science.ru.nl</a>><br>
> > > Subject: Re: [FieldTrip] problem in connectivityanalysis<br>> > > Message-ID: <<a href="mailto:1348f6433edc26f9.519bb455@upo.es">1348f6433edc26f9.519bb455@upo.es</a>><br>> > > Content-Type: text/plain; charset="iso-8859-1"<br>
> > ><br>> > > Thank you so much for your answer J?rn,<br>> > ><br>> > ><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>> > ><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>> > ><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>> > ><br>> ><br>> > > So, both are apparently the same, but the data inside this matrices is different,<br>
> > ><br>> > ><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>> > ><br>> ><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>> > ><br>> ><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>> > ><br>> ><br>> > > I'm pretty sure I'm doing something wrong here, but I don't really have a clue,<br>> > ><br>> > ><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>> > ><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" <<a href="mailto:jm.horschig@donders.ru.nl">jm.horschig@donders.ru.nl</a>><br>> > > Fecha: Martes, 21 de Mayo de 2013, 9:03 am<br>> > > Asunto: Re: [FieldTrip] problem in connectivityanalysis<br>
> > > A: FieldTrip discussion list <<a href="mailto:fieldtrip@science.ru.nl">fieldtrip@science.ru.nl</a>><br>> > ><br>> > ><br>> > ><br>> > ><br>> > ><br>> > > > Dear Gabriel,<br>
> > >        ><br>> > ><br>> ><br>> > >        > 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 opposite, this is, to get the imaginary part of coherence for                   each epoch (repetition) on all                   my channels and not for each frequency, something like: chan_chan_rpt or chancbm_rpt,                 Is this possible?<br>
> > >                         ><br>> > ><br>> ><br>> > > >        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>
> > >        ><br>> > ><br>> ><br>> > >        > 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?       I see some terminology issues here. Let me try to explain in a       simple (aka least-mathematical I can think of) way:<br>
> > >        ><br>> > ><br>> ><br>> > > >        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>
> > ><br>> ><br>> > > >        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>
> > >        ><br>> > > >        I hope that somehow clarifies your questions.<br>> > >        ><br>> > > >        Best,<br>> > > >        J?rn<br>> > >        ><br>
> > > >        On 5/17/2013 5:45 PM, Gabriel Gonzalez Escamilla wrote:<br>> > >           > Dear Fieltrip experts,<br>> > >          ><br>> > ><br>> ><br>> > > >          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>
> > ><br>> ><br>> > > >          data.trials {1xNepochs}; % Nepochs = 7,         each epoch with 59 channels and 2000 samples<br>> > > >          data.label {59x1};<br>> > > >          data.time {1xNepochs};<br>
> > >          ><br>> > > >          then, calculate only the FFT of all sensors as:<br>> > > >          cfg=[];<br>> > > >          cfg.output='fourier';<br>> > > >          cfg.method='mtmfft';<br>
> > > >          cfg.taper='hanning';<br>> > > >          cfg.foilim=[8.5 9.5];<br>> > > >          cfg.tapsmofrq=0;<br>> > > >          cfg.trials='all';<br>
> > > >          cfg.keeptrials='yes';<br>> > > >          cfg.channel='all';<br>> > > >          fourier=ft_freqanalysis(cfg, data)<br>> > > >          as output I get:<br>
> > > >          dimord = 'rpttap_chan_freq'<br>> > > >          freq= [1x9]<br>> > > >          fourierspctrm=[7x59x9 double]<br>> > >          ><br>> > > >          Then compute the imaginary part of coherency as:<br>
> > >            > cfg=[];<br>> > >              > cfg.method ='cohe';<br>> > >              > cfg.complex='imag';<br>> > >                > cfg.channelcbm={'all'               'all'};<br>
> > >              > coher = ft_connectivityanalysis (cfg, fourier)<br>> > >          > as output get:<br>> > >            > dimord='chan_chan_freq';<br>> > >              > cohspctrm=[59x59x9 double];<br>
> > >              > dof=7;<br>> > >            ><br>> > >          > I have two main questions:<br>> > >          ><br>> > ><br>> ><br>> > >          > 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 opposite, this is, to               get the imaginary                 part of coherence for each epoch (repetition)                 on all my channels and not for each frequency, something like: chan_chan_rpt or chancbm_rpt,               Is this possible?<br>
> > >                  ><br>> > >                  ><br>> > >                  ><br>> > ><br>> ><br>> > >                  > 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?<br>
> > >                  ><br>> > >                  ><br>> > >                  ><br>> > >                  ><br>> > > >                  Many thanks in advanced,<br>
> > >                  > Gabriel.<br>> > >                  ><br>> > >                       ><br>> > >               ><br>> > ><br>> ><br>> > >        > _______________________________________________ > fieldtrip mailing list ><br>
> > > <a href="mailto:fieldtrip@donders.ru.nl">fieldtrip@donders.ru.nl</a><br>> > > <a href="http://mailman.science.ru.nl/mailman/listinfo/fieldtrip">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip</a><br>
> > >      ><br>> > ><br>> ><br>> > >      > --  > 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:<br>
> > > <a href="mailto:jm.horschig@donders.ru.nl">jm.horschig@donders.ru.nl</a> > Tel:    +31-(0)24-36-68493 > Web:<br>> > > <a href="http://www.ru.nl/donders">http://www.ru.nl/donders</a><br>> ><br>
> > >  >  > Visiting address: > Trigon, room 2.30 > Kapittelweg 29 > NL-6525 EN Nijmegen > The Netherlands     > _______________________________________________<br>> > > > fieldtrip mailing list<br>
> > > > <a href="mailto:fieldtrip@donders.ru.nl">fieldtrip@donders.ru.nl</a><br>> > > > <a href="http://mailman.science.ru.nl/mailman/listinfo/fieldtrip">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip</a><br>
> > ><br>> > ><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:<br>
> > > <a href="mailto:ggonesc@upo.es">ggonesc@upo.es</a><br /><a href="http://www.upo.es/neuroaging/es/">http://www.upo.es/neuroaging/es/</a></font><br>> > ><br>> > ><br>> > > -------------- next part --------------<br>
> > > An HTML attachment was scrubbed...<br>> > > URL: <<br>> > ><br>> > <a href="http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20130521/0886ac06/attachment.html">http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20130521/0886ac06/attachment.html</a><br>
> > > ><br>> > ><br>> > > ------------------------------<br>> > ><br>> > > _______________________________________________<br>> > > fieldtrip mailing list<br>> > > <a href="mailto:fieldtrip@donders.ru.nl">fieldtrip@donders.ru.nl</a><br>
> > > <a href="http://mailman.science.ru.nl/mailman/listinfo/fieldtrip">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip</a><br>> > ><br>> > > End of fieldtrip Digest, Vol 30, Issue 29<br>
> > > *****************************************<br>> > ><br>> > > _______________________________________________<br>> > > fieldtrip mailing list<br>> > > <a href="mailto:fieldtrip@donders.ru.nl">fieldtrip@donders.ru.nl</a><br>
> > > <a href="http://mailman.science.ru.nl/mailman/listinfo/fieldtrip">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip</a><br>> > ><br>> > -------------- next part --------------<br>> > An HTML attachment was scrubbed...<br>
> > URL: <<br>> > <a href="http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20130522/190aa246/attachment.html">http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20130522/190aa246/attachment.html</a><br>
> > ><br>> ><br>> > ------------------------------<br>> ><br>> > _______________________________________________<br>> > fieldtrip mailing list<br>> > <a href="mailto:fieldtrip@donders.ru.nl">fieldtrip@donders.ru.nl</a><br>
> > <a href="http://mailman.science.ru.nl/mailman/listinfo/fieldtrip">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip</a><br>> ><br>> > End of fieldtrip Digest, Vol 30, Issue 31<br>> > *****************************************<br>
> ><br>> > _______________________________________________<br>> > fieldtrip mailing list<br>> > <a href="mailto:fieldtrip@donders.ru.nl">fieldtrip@donders.ru.nl</a><br>> > <a href="http://mailman.science.ru.nl/mailman/listinfo/fieldtrip">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip</a><br>
> ><br>> -------------- next part --------------<br>> An HTML attachment was scrubbed...<br>> URL: <<a href="http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20130525/4259b000/attachment.html">http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20130525/4259b000/attachment.html</a>><br>
><br>> ------------------------------<br>><br>> _______________________________________________<br>> fieldtrip mailing list<br>> <a href="mailto:fieldtrip@donders.ru.nl">fieldtrip@donders.ru.nl</a><br>> <a href="http://mailman.science.ru.nl/mailman/listinfo/fieldtrip">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip</a><br>
><br>> End of fieldtrip Digest, Vol 30, Issue 36<br>> *****************************************<br><br><br><br><br>--<br>Stephen Politzer-Ahles<br>University of Kansas<br>Linguistics Department<br><a href="http://people.ku.edu/~sjpa/">http://people.ku.edu/~sjpa/</a><br>
<br><br></div></div></div></div></div>