<div dir="ltr">I am using this code:<div><br></div><div><div>cfg = [];</div><div>cfg.method = 'mtmfft';</div><div>cfg.output = 'pow';</div><div>cfg.taper     = 'hanning';</div><div>%cfg.toi       = 1.6 : 0.02 : 6.0;</div><div>cfg.foi       = 2:60;</div><div>%cfg.t_ftimwin = ones(size(cfg.foi)) * 1.2;</div><div>cfg.keeptrials = 'yes';</div><div><br></div><div>freq1          = ft_freqanalysis(cfg, data1);</div><div>freq2          = ft_freqanalysis(cfg, data2);</div><div><br></div><div>%layout = ft_prepare_layout(cfg, ftadata1)</div><div><br></div><div>cfg = [];</div><div>cfg.channel          = 'all';</div><div>%cfg.latency          = 'all';</div><div>cfg.frequency        = 'all';</div><div>cfg.method           = 'montecarlo';</div><div>cfg.statistic        = 'ft_statfun_indepsamplesT';</div><div>cfg.correctm         = 'cluster';</div><div>cfg.clusteralpha     = 0.05;</div><div>cfg.clusterstatistic = 'maxsum';</div><div>cfg.minnbchan        = 2;</div><div>cfg.tail             = 0;</div><div>cfg.clustertail      = 0;</div><div>cfg.alpha            = 0.025;</div><div>cfg.numrandomization = 1000;</div><div>% prepare_neighbours determines what sensors may form clusters</div><div>cfg_neighb.method    = 'distance';</div><div>cfg.neighbours       = ft_prepare_neighbours(cfg_neighb, ftadata1);</div><div>design = zeros(1,size(freq1,1) + size(freq2.powspctrm,1));</div><div>design(1,1:size(freq1.powspctrm,1)) = 1;</div><div>design(1,(size(freq1.powspctrm,1)+1):(size(freq1.powspctrm,1)+size(freq2.powspctrm,1))) = 2;</div><div>cfg.design           = design;</div><div>cfg.ivar             = 1;</div><div>[stat] = ft_freqstatistics(cfg, freq1, freq2);</div><div><br></div><div><br></div><div>cfg = []; </div><div>freq11 = ft_freqdescriptives(cfg, freq1);</div><div>freq22  = ft_freqdescriptives(cfg, freq2);</div><div>stat.raweffect = freq11.powspctrm - freq22.powspctrm;</div><div><br></div><div>cfg = [];</div><div>cfg.alpha  = 0.05;</div><div>cfg.parameter = 'raweffect';</div><div>%cfg.zlim   = [-1e-27 1e-27];</div><div>cfg.layout = layout ;</div><div>ft_clusterplot(cfg, stat);</div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Jul 17, 2015 at 3:52 PM, Jose Rebola <span dir="ltr"><<a href="mailto:jrebola@gmail.com" target="_blank">jrebola@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">How can I use the outputs of mtmfft ft_freqanalysis to compare two spectra, by means of non-parametric randomization procedures?<div><br></div><div>I have the two outputs with the keeptrials option on</div><div>Does anybody have sample code for this or an example?</div><div>I only can find examples for TFR. I have adapted the code but not successfully.</div><div>I can't interpret the output of ft_freqstatistics so well..._</div><div><br></div><div>thank you all so much</div><div><br></div><div>José</div></div>
</blockquote></div><br></div>