source coherence in one condition

Andrew Smart andrew.smart at NYU.EDU
Thu Apr 26 00:35:22 CEST 2007


Hi,
I would like to try to compute whether there is significant coherence across
subjects but within one condition and between two dipoles. This code
generates graphs but I am unclear about what they mean:

for i = 1 : length(the_sources) %this is coherence between pairs of BESA      
                                 %dipoles, imported *.tfc files
    cfg = [];
    cfg.keepindividual = 'yes';
    cfg.cohrefchannel= 'M100';   %I would like to test coherence between  
                                 % this dipole and the other dipoles in the 
                                 % BESA model on a group level
    cfg.channel = the_sources{i};
    cfg.parameter = 'cohspctrm';
    cfg.method = 'montecarlo';
    cfg.neighbours    = [];
    cfg.time          = 'all';
    cfg.freq        = 'all';
    cfg.avgovertime      = 'no';
    cfg.avgoverfreq      = 'no';
    cfg.avgoverchan      = 'no';
    cfg.computeprob      = 'yes';
    cfg.computecritval   = 'yes';
    cfg.statistic        = 'indepsamplesregrT' ; %only stat test that works
    cfg.correctm         = 'max';
    cfg.numrandomization = 500;
     cfg.design           = [ 1 2 3 4 5 6 7 8 9 10      %is this correct?
                              1 1 1 1 1 1 1 1 1 1 ] ;   %i have 10 subjects
           
     cfg.uvar = 1;                                       %one condition
     cfg.ivar = 1; 
     cfg.wvar = 1;                                                    
                                                                     

    Prepstat = freqstatistics(cfg, Prepgram{:});
    
    cfg = [];
    cfg.xparam       =  'time';
    cfg.yparam        = 'freq';
    cfg.zparam        = 'stat';    %what are the values in this field?
    cfg.zlim          = 'absmax';

    figure; singleplotTFR(cfg,Prepstat(:));
    title([the_sources{i} '_coherence']);
    colormap jet
end

Do you have any suggestions for improving this and making the results more
interpretable? (what do the results mean?) Are there other available
statistical tests in Fieldtrip, for example is 'indepsamplesZcoh' implemented?

Thanks,
andy

----------------------------------
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/



More information about the fieldtrip mailing list