[FieldTrip] between-subject cluster-stats wont average over specified frequency

George McKenzie Opie george.opie at adelaide.edu.au
Wed Oct 26 22:51:57 CEST 2016


Wow, cant believe I missed that! Thanks Julian, much appreciated.

George Opie
ARC Research Associate
Discipline of Physiology
School of Medicine
The University of Adelaide, AUSTRALIA 5005
Ph    : +61 8 8313 4157
Fax   : +61 8 8303 5384
e-mail: george.opie at adelaide.edu.au<mailto:george.opie at adelaide.edu.au>
________________________________
From: fieldtrip-bounces at science.ru.nl <fieldtrip-bounces at science.ru.nl> on behalf of Julian Keil <julian.keil at gmail.com>
Sent: Thursday, 27 October 2016 7:20:40 AM
To: FieldTrip discussion list
Subject: Re: [FieldTrip] between-subject cluster-stats wont average over specified frequency

Dear George,

You have a typo in the cfg at cfg.frequency.

Hope this helps,

Julian

Am Mittwoch, 26. Oktober 2016 schrieb George McKenzie Opie :

Hi all,


Im trying to run some between-subject cluster-based analyses on some time-frequency data, but am having some issues getting the analysis to average over a specified frequency range. For some reason this only happens with between-subject comparisons and not within-subject. My cfg structure is shown below. D1 and D2 are grandaverage data from two groups calculated using ft_freqgrandaverage. When I call ft_freqstatistics with this cfg, it averages over all frequencies (5 - 45 Hz), instead of the specified frequency range (31 - 45 Hz). Any help would be very much appreciated.


cfg = [];
cfg.channel     = {'all'};
cfg.minnbchan        = 2;
cfg.clusteralpha = 0.01;
cfg.clusterstatistic = 'maxsum';
cfg.alpha       = 0.05;
cfg.latency     = [0.025, 0.220];
cfg.avgoverfreq = 'yes';
cfg.frequnecy = [31 45];
cfg.avgovertime = 'yes';
cfg.avgoverchan = 'no';
cfg.statistic   = 'indepsamplesT';
cfg.numrandomization = 2000;
cfg.correctm    = 'cluster';
cfg.method      = 'montecarlo';
cfg.tail             = 0;
cfg.clustertail      = 0;
cfg.neighbours  = neighbours;
cfg.parameter   = 'powspctrm';

design = zeros(1,size(D1.powspctrm,1) + size(D2.powspctrm,1));
design(1,1:size(D1.powspctrm,1)) = 1;
design(1,(size(D1.powspctrm,1)+1):(size(D1.powspctrm,1) + size(D2.powspctrm,1)))= 2;

cfg.design = design;
cfg.ivar  = 1;

[stat] = ft_freqstatistics(cfg, D1, D2);



kind regards,

George Opie
ARC Research Associate
Discipline of Physiology
School of Medicine
The University of Adelaide, AUSTRALIA 5005
Ph    : +61 8 8313 4157
Fax   : +61 8 8303 5384
e-mail: george.opie at adelaide.edu.au<javascript:_e(%7B%7D,'cvml','george.opie at adelaide.edu.au');>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20161026/16906edd/attachment-0001.html>


More information about the fieldtrip mailing list