[FieldTrip] plotting neighbours

victoria schroeder vic.schroeder2 at gmail.com
Sat Jan 16 18:26:53 CET 2016


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)
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20160116/3204dfe4/attachment-0001.html>


More information about the fieldtrip mailing list