<div dir="ltr"><div><div><div><div><div><div>Hi all,<br><br></div>I'm having an issue running freqstatistics using the indepsamplesZcoh and a cluster correction. The indepsamplesZcoh portion seems to run just fine but I get an error when I try to perform the cluster statistics. The first issue is error1: The error describes a dimension mismatch at line 297 of ft_statistics_montecarlo: statrand(:,:,i) = getfield(dum, 'stat');<br>
</div>Apparently the output (dum) from indepsamplesZcoh collapses my 100 frequencies and 100 timepoints into 10000 points (nchancmb/2 x 10000) but the dimensions of statrand are nchancmb/2 x 100 x Nrand which obviously creates a dimension mismatch. I can get past that part if I set avgovertime = 'yes' which takes me to error 2: "To RESHAPE the number of elements must not change" refers to line 202 in clusterstat.m, posclusobs = findcluster(reshape(postailobs, [cfg.dim,1]),channeighbstructmat,cfg.minnbchan);<br>
</div>Here, cfg.dim is nchan x 100 freq x 1 (avgtime) which obviously doesn't match postailobs dimensions of nchancmb/2 x 100.<br><br></div>I thought perhaps the ndepsamplesZcoh function is not supposed to be used in conjunction with the cluster statistics but the Maris 2007 paper appears to use the two together. Any help with this would be much appreciated.<br>
<br></div>Below is the code that that I'm using:<br><br>cfg = [];<br>%cfg.avgovertime = 'yes';<br>cfg.correctm =           'cluster';<br>cfg.clusteralpha =       0.05;    <br>cfg.clusterstatistic =   'maxsum';   <br>
cfg.clusterthreshold =   'parametric';  <br>cfg.minnbchan =          2;         <br>cfg.correcttail =        'alpha';<br>cfg.method    = 'montecarlo';<br>cfg.statistic = 'indepsamplesZcoh';<br>
cfg.parameter = 'fourierspctrm';<br>cfg.computecritval = 'yes';<br>cfg.neighbours =          neighbours;<br>cfg.numrandomization = 99;<br>cfg.alpha =              0.05;       <br>cfg.tail =               0;          <br>
design = ones(1, 90); design(91:180)=2;<br>cfg.design =             design;<br>cfg.label     = freqL.label;<br>stat = ft_freqstatistics(cfg, freqL,freqR);<br><br><br></div>Thanks,<br><div><div><div><div>-Tony<br></div></div>
</div></div></div>