[FieldTrip] plotting neighbours

[Tzvetan Popov]

Dear Victoria,
you should remove the ft_combineplanar step from your analysis pipeline. This function takes input from ft_megplanar. The latter is applied only to magnetometers or axial gradiometers. 

best
tzvetan
 

> Hi !
> 
> I am working with a neuromag system and want to do cluster based permutation tests with time frequency data. However, before running the stats i would like to check the clusters which are formed by ft_prepare_neighbours. Thus i use ft_neighbourplot. When looking at the plot and clicking on the channels i would expect to see something like :
> 
> Selected channel MEG1032+1033, which has 7 neighbours: MEG0622+0623, MEG1012+1013, MEG1022+1023, MEG1042+1043, MEG1112+1113, MEG1122+1123, MEG1242+1243
> Selected channel MEG0742+0743, which has 7 neighbours: MEG0432+0433, MEG0712+0713, MEG0722+0723, MEG0732+0733, MEG1822+1823, MEG1832+1833, MEG2242+2243
> 
> however what i see is:
> 
> Selected channel MEG0741, which has 0 neighbours: 
> Selected channel MEG0731, which has 0 neighbours: 
> Selected channel MEG0742+0743, which has 7 neighbours: MEG0432+0433, MEG0712+0713, MEG0722+0723, MEG0732+0733, MEG1822+1823, MEG1832+1833, MEG2242+2243
> Selected channel MEG0711, which has 0 neighbours: 
> Selected channel MEG0741, which has 0 neighbours: 
> 
> It is absolutely random when the correct output occurs. When clicking on a particular channel it would randomly show me the following results:
> 
> Selected channel MEG1032+1033, which has 7 neighbours: MEG0622+0623, MEG1012+1013, MEG1022+1023, MEG1042+1043, MEG1112+1113, MEG1122+1123, MEG1242+1243
> 
> or
> 
> Selected channel MEG1031, which has 0 neighbours: 
> 
> Do you have any suggestions why this happen? Are the channels not grouped correctly?
> 
> Below is the code starting from the time-frequency analysis of the data.
> 
> Thank you very much
> Victoria
> 
> cfg = [];
> cfg.output = 'pow';
> cfg.channel = 'MEGGRAD';
> cfg.method = 'mtmconvol';
> cfg.taper = 'hanning';
> cfg.foi = [2:2:30];
> %cfg.keeptrials = 'yes';   for stats keep trials
> cfg.t_ftimwin    = ones(length(cfg.foi),1).*0.5;   % length of time window = 0.5 sec
> cfg.toi          = [-1.3:0.05:1];   
> cfg.keeptrials   = 'yes';
> FreqCon = ft_freqanalysis(cfg, Con_clean);
> FreqIncon = ft_freqanalysis(cfg, Incon_clean);
> 
> %combine planar
> cfg = [];
> cfg.method ='sum';
> ComCon = ft_combineplanar(cfg,FreqCon);
> ComIncon = ft_combineplanar(cfg, FreqIncon);
> 
> % plot neighbours
> cfg=[];
> cfg.method        = 'triangulation';
> cfg.layout           = 'neuromag306cmb.lay';
> ft_neighbourplot(cfg,ComCon)
> _______________________________________________
> fieldtrip mailing list
> fieldtrip at donders.ru.nl
> http://mailman.science.ru.nl/mailman/listinfo/fieldtrip

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20160116/dd53a3b9/attachment-0002.html>