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