[FieldTrip] plotting neighbours

[Tzvetan Popov]

Hi Victoria,
well if you work with the magnetometers yes. Maybe I misunderstood your previous posts. Yet this configuration from your post below cfg.channel = 'MEGGRAD’; results in time-frequency representation of power on the planar gradiometers only. In that case neither ft_megplanar nor ft_combineplanar are appropriate since the data stems from planar gradiometers. The tutorial you mentioned deals only with the magnetometers thus ft_megplanar and ft_combineplanar are used.

best
tzvetan
 
> Dear Tzvetan
> 
> Are you sure about that in the tutorial http://www.fieldtriptoolbox.org/tutorial/natmeg/preprocessing neuromag data is used and ft_combineplanar is used to combine the gradiometers. This is a paragraph from this tutorial:
> 
> In fact, we are now plotting the two different gradiometers together. You can see the channel locations being in pairs, one above the other. They are in reality, however, at the same location but oriented differently - radially and axially with respect to the surface of the helmet. They can thereby pick up both radial orientations of the magnetic fields. To use them properly for the purpose of plotting, we should therefor combine them first, adding their fields.
> 
> TFR_left_MEG_comb  = ft_combineplanar([],TFR_left_MEG);
> TFR_right_MEG_comb = ft_combineplanar([],TFR_right_MEG);
>  
> cfg = [];
> cfg.baseline     = [-0.5 -0.1];	
> cfg.baselinetype = 'absolute';
> cfg.xlim         = [0.4 0.8];   
> cfg.ylim         = [15 25];
> cfg.zlim         = [-4e-24 4e-24];	        
> cfg.marker       = 'on';
> cfg.layout       = 'neuromag306cmb.lay';
>  
> figure;
> ft_topoplotTFR(cfg, TFR_left_MEG_comb);
> print -dpng natmeg_freq13.png
> Thus, it should be possible to use ft_combineplanar for neuromag data.
> What do you think?
> 
> Cheers
> Victoria
> 
> 2016-01-16 18:00 GMT+00:00 Tzvetan Popov <tzvetan.popov at uni-konstanz.de>:
> 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
> 
> 
> _______________________________________________
> fieldtrip mailing list
> fieldtrip at donders.ru.nl
> http://mailman.science.ru.nl/mailman/listinfo/fieldtrip
> 
> _______________________________________________
> 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/20160118/22e5df87/attachment-0002.html>