AW: [FIELDTRIP] BUG in ft_topoplotER (cfg.xlim) version 20100127

jan-mathijs schoffelen jan.schoffelen at DONDERS.RU.NL
Tue Feb 2 16:28:23 CET 2010


Dear Hanneke,

Just to follow up: it seems there was a discrepancy between the part  
of code you copied (and adjusted) into one of your previous e-mails,  
and the code of the latest ft_topoplotER.

> 444 % make dat structure with one value for each channel
> 445 dat = getsubfield(data, cfg.zparam);
> 446 if ~isempty(cfg.yparam),
> 447   dat = dat(:, ymin:ymax, xmin:xmax);
> 448   dat = nanmean(nanmean(dat, 2), 3);
> 449 elseif ~isempty(cfg.component) | ~isempty(cfg.xlim);%this part  
> Hanneke added
> 450  dat = dat(:, xmin:xmax);
> 451  dat = nanmean(dat, 2);
> 452 end
> 453
> 454 dat = dat(:);


There seems to be an 'else' missing after line 449.

I guess your problem will be solved when you upgrade to the latest  
fieldtrip version.
Sorry for the inconvenience,

Jan-Mathijs



On Feb 2, 2010, at 12:43 PM, Hanneke van Dijk wrote:

> Maybe like this then sorry for the large amount of mails...
>
> On Tue, Feb 2, 2010 at 12:40 PM, Hanneke Van Dijk <Hanneke.vanDijk at med.uni-duesseldorf.de 
> > wrote:
> Hi Arjen,
>
> This time I am plotting cortico-muscular coherence obtained from  
> ft_freqanalysis_mtmfft with fourier output and then using  
> ft_connectivityanalysis. I wanted to plot the coherence on helmet- 
> level only for 8-12 Hz. I attached the figures; When I do the  
> average over the .cohspctrm field myself (or add the line I reported  
> before) I get the left plot (figure 2) and if I say cfg.xlim I get  
> the right plot (figure 1).
>
> Groetjes Hanneke
>
>
>
> --------------------------------------------------
> Institut für Klinische Neurowissenschaften und Medizinische  
> Psychologie
> Gebäude-Nr.: 23.02
> Ebene: 03 Zimmer-Nr.: 47
> Tel.:  0211-81-13074
> Mail : hanneke.vandijk at med.uni-duesseldorf.de
> http://www.uniklinik-duesseldorf.de/deutsch/unternehmen/institute/KlinNeurowiss/Team/HannekevanDijk/page.html
>
>
>
> -----Ursprüngliche Nachricht-----
> Von: FieldTrip discussion list im Auftrag von a.stolk at fcdonders.ru.nl
> Gesendet: Di 02.02.2010 12:03
> An: FIELDTRIP at NIC.SURFNET.NL
> Betreff: Re: [FIELDTRIP] BUG in ft_topoplotER (cfg.xlim) version  
> 20100127
>
> Hi Hanneke,
>
> Thank you for your report. Could you specify some more details to  
> help locating the problem?
>
> Are you trying to plot event-related fields or components (I'm  
> asking this because of our earlier mail exchanges)? For the latter,  
> I could recommend the ft_topoplotIC function. This one has maxmin  
> values for the z-axis by default, but is more intuitive.
>
> An example:
>
> cfg = [];
> cfg.layout = 'NM122_ver.lay';
> cfg.component = 1;
> ft_topoplotIC(cfg,comp);
>
>
> Best,
> Arjen
>
>
>
>
>
>
>
>
>
>
> ----- Original Message -----
> From: "Hanneke Van Dijk" <Hanneke.vanDijk at MED.UNI-DUESSELDORF.DE>
> To: FIELDTRIP at NIC.SURFNET.NL
> Sent: Tuesday, February 2, 2010 11:16:40 AM GMT +01:00 Amsterdam /  
> Berlin / Bern / Rome / Stockholm / Vienna
> Subject: [FIELDTRIP] BUG in ft_topoplotER (cfg.xlim) version 20100127
>
> Dear all,
>
> I wanted to report a bug in ft_topoplotER (Fieldtrip version  
> 20100127), regarding the use of cfg.xlim. It doesn't seem to use  
> cfg.xlim when you choose one, but just to use the 'maxmin'. When I  
> average the data over xlim myself I get a completely different  
> topoplot then when I use cfg.xlim for the same range. I wouldn't  
> have known if I didn't need to average over xlim for other reasons.
>
> It seems to happen between lines 445 and 452. The data are selected  
> for cfg.ylim and cfg.component but not for cfg.xlim. This changes  
> the outlook on my data completely, and I think this could be  
> important for all users of this version. Therefore I send it to the  
> entire mailing list.
>
> I have changed it now to:
> 444 % make dat structure with one value for each channel
> 445 dat = getsubfield(data, cfg.zparam);
> 446 if ~isempty(cfg.yparam),
> 447   dat = dat(:, ymin:ymax, xmin:xmax);
> 448   dat = nanmean(nanmean(dat, 2), 3);
> 449 elseif ~isempty(cfg.component) | ~isempty(cfg.xlim);%this part  
> Hanneke added
> 450  dat = dat(:, xmin:xmax);
> 451  dat = nanmean(dat, 2);
> 452 end
> 453
> 454 dat = dat(:);
>
> This helps for me and it is very important for me as well, my data  
> looks meaningfull now!
>
> Best and good luck!
>
> Hanneke
>
>
>
>
> --------------------------------------------------
> Institut für Klinische Neurowissenschaften und Medizinische  
> Psychologie
> Gebäude-Nr.: 23.02
> Ebene: 03 Zimmer-Nr.: 47
> Tel.:  0211-81-13074
> Mail : hanneke.vandijk at med.uni-duesseldorf.de
> http://www.uniklinik-duesseldorf.de/deutsch/unternehmen/institute/KlinNeurowiss/Team/HannekevanDijk/page.html
>
>
> ----------------------------------
> The aim of this list is to facilitate the discussion between users  
> of the FieldTrip  toolbox, to share experiences and to discuss new  
> ideas for MEG and EEG analysis. See also http://listserv.surfnet.nl/archives/fieldtrip.html 
>  and http://www.ru.nl/neuroimaging/fieldtrip.
>
> ----------------------------------
> The aim of this list is to facilitate the discussion between users  
> of the FieldTrip  toolbox, to share experiences and to discuss new  
> ideas for MEG and EEG analysis. See also http://listserv.surfnet.nl/archives/fieldtrip.html 
>  and http://www.ru.nl/neuroimaging/fieldtrip.
>
> figure;
> cfgp=[];
> cfgp.cohrefchannel = {'MEG063+064'};
> cfgp.xparam = 'freq';
> cfgp.zparam = 'cohspctrm';
> cfgp.layout = 'NM122combined_test.lay';
> cfgp.showlabels = 'yes';
> cfgp.xlim = [8 12];
> %cfgp.zlim = [0 1];
> ft_topoplotER(cfgp,fd_ccc_off);colorbar
>
>
> ----------------------------------
> The aim of this list is to facilitate the discussion between users  
> of the FieldTrip  toolbox, to share experiences and to discuss new  
> ideas for MEG and EEG analysis. See also http://listserv.surfnet.nl/archives/fieldtrip.html 
>  and http://www.ru.nl/neuroimaging/fieldtrip.
>
> ----------------------------------
>
> The aim of this list is to facilitate the discussion between users  
> of the FieldTrip toolbox, to share experiences and to discuss new  
> ideas for MEG and EEG analysis.
>
> http://listserv.surfnet.nl/archives/fieldtrip.html
>
> http://www.ru.nl/fcdonders/fieldtrip/
>
> <cfgxlimplots.jpeg>


----------------------------------
The aim of this list is to facilitate the discussion between users of the FieldTrip  toolbox, to share experiences and to discuss new ideas for MEG and EEG analysis. See also http://listserv.surfnet.nl/archives/fieldtrip.html and http://www.ru.nl/neuroimaging/fieldtrip.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20100202/4357248b/attachment-0001.html>


More information about the fieldtrip mailing list