I may be wrong but as far as i remember, the cluster analysis can only be done on one value of frequency, average or not.<div>There is no cfg.avgoverfreq option n your script.</div><div>But again i might be wrong.</div><div>
<br></div><div>Rodolphe <br><br><div class="gmail_quote">On Fri, Jan 14, 2011 at 11:58 AM, Jung, Ki-Young <span dir="ltr"><<a href="mailto:jungky@korea.ac.kr">jungky@korea.ac.kr</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Thanks! But I have still trouble with this message.<br>
<br>
??? Error using ==> statistics_montecarlo<br>
Too many input arguments.<br>
<div class="im"><br>
Error in ==> statistics_wrapper at 285<br>
    [stat, cfg] = statmethod(cfg, dat, cfg.design, 'issource',issource);<br>
<br>
Error in ==> ft_freqstatistics at 125<br>
[stat, cfg] = statistics_wrapper(cfg, varargin{:});<br>
<br>
</div>Best regards,<br>
<font color="#888888"><br>
KY<br>
</font><div class="im"><br>
On Fri, Jan 14, 2011 at 1:04 AM, Julian Keil <<a href="mailto:julian.keil@gmail.com">julian.keil@gmail.com</a>> wrote:<br>
> Hi,<br>
><br>
> did you check with 'which bwlabeln' if you have the respcetive function?<br>
> It's part of the image processing toolbox for matlab.<br>
><br>
> Good Luck.<br>
><br>
> Julian<br>
><br>
> Am 14.01.2011 um 09:58 schrieb Jung, Ki-Young:<br>
><br>
</div><div><div></div><div class="h5">>> Dear all,<br>
>><br>
>> I'd like to do statistical comparison using cluster based permuation<br>
>> between two groups on ERP data.<br>
>> But I am receiving following error messages<br>
>> ??? Undefined function or method 'bwlabeln' for input arguments of<br>
>> type 'double'.<br>
>><br>
>> Error in ==> findcluster at 89<br>
>>  [labelmat(spatdimlev, :, :), num] = bwlabeln(reshape(onoff(spatdimlev, :, :),<br>
>>  nfreq, ntime), 4);<br>
>><br>
>> Error in ==> clusterstat at 194<br>
>>     posclusobs = findcluster(reshape(postailobs,<br>
>>     [cfg.dim,1]),channeighbstructmat,cfg.minnbchan);<br>
>><br>
>> Error in ==> statistics_montecarlo at 321<br>
>>  [stat, cfg] = clusterstat(cfg, statrand, statobs,'issource',issource);<br>
>><br>
>> Error in ==> statistics_wrapper at 285<br>
>>   [stat, cfg] = statmethod(cfg, dat, cfg.design, 'issource',issource);<br>
>><br>
>> Error in ==> ft_freqstatistics at 125<br>
>> [stat, cfg] = statistics_wrapper(cfg, varargin{:});<br>
>><br>
>> My configuration for ft_freqstatistics are<br>
>> cfg = [];<br>
>> cfg.channel          = 'all';<br>
>> cfg.latency          = 'all';<br>
>> cfg.frequency        = [30 50];<br>
>> cfg.method           = 'montecarlo';<br>
>> cfg.statistic        = 'indepsamplesT';<br>
>> cfg.correctm         = 'cluster';<br>
>> cfg.clusteralpha     = 0.025;<br>
>> cfg.clusterstatistic = 'maxsum';<br>
>> cfg.minnbchan        = 2;<br>
>> cfg.tail             = 0;<br>
>> cfg.clustertail      = 0;<br>
>> cfg.alpha            = 0.025;<br>
>> cfg.numrandomization = 500;<br>
>> cfg.neighbours       = [ ];<br>
>> design = [ones(1, 16), ones(1, 16)*2];<br>
>> cfg.design           = design;<br>
>> cfg.ivar             = 1;<br>
>> [stat] = ft_freqstatistics(cfg, mig_grand_avg, con_grand_avg);<br>
>><br>
>><br>
>> Data structure of con_grand_avg are<br>
>><br>
>> con_grand_avg =<br>
>>       label: {27x1 cell}<br>
>>        freq: [30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50]<br>
>>        time: [1x21 double]<br>
>>      dimord: 'subj_chan_freq_time'<br>
>>   powspctrm: [4-D double]<br>
>>         cfg: [1x1 struct]<br>
>><br>
>><br>
>> Any comment helps me.<br>
>><br>
>><br>
>> Best regards,<br>
>><br>
>> Ki-Young Jung<br>
>><br>
>> _______________________________________________<br>
>> fieldtrip mailing list<br>
>> <a href="mailto:fieldtrip@donders.ru.nl">fieldtrip@donders.ru.nl</a><br>
>> <a href="http://mailman.science.ru.nl/mailman/listinfo/fieldtrip" target="_blank">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip</a><br>
><br>
><br>
> _______________________________________________<br>
> fieldtrip mailing list<br>
> <a href="mailto:fieldtrip@donders.ru.nl">fieldtrip@donders.ru.nl</a><br>
> <a href="http://mailman.science.ru.nl/mailman/listinfo/fieldtrip" target="_blank">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip</a><br>
><br>
><br>
><br>
<br>
_______________________________________________<br>
fieldtrip mailing list<br>
<a href="mailto:fieldtrip@donders.ru.nl">fieldtrip@donders.ru.nl</a><br>
<a href="http://mailman.science.ru.nl/mailman/listinfo/fieldtrip" target="_blank">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip</a><br>
</div></div></blockquote></div><br></div>