[FieldTrip] Fwd: question about ft_statistics errors

[Jung, Ki-Young]

Thanks! But I have still trouble with this message.

??? Error using ==> statistics_montecarlo
Too many input arguments.

Error in ==> statistics_wrapper at 285
    [stat, cfg] = statmethod(cfg, dat, cfg.design, 'issource',issource);

Error in ==> ft_freqstatistics at 125
[stat, cfg] = statistics_wrapper(cfg, varargin{:});

Best regards,

KY

On Fri, Jan 14, 2011 at 1:04 AM, Julian Keil <julian.keil at gmail.com> wrote:
> Hi,
>
> did you check with 'which bwlabeln' if you have the respcetive function?
> It's part of the image processing toolbox for matlab.
>
> Good Luck.
>
> Julian
>
> Am 14.01.2011 um 09:58 schrieb Jung, Ki-Young:
>
>> Dear all,
>>
>> I'd like to do statistical comparison using cluster based permuation
>> between two groups on ERP data.
>> But I am receiving following error messages
>> ??? Undefined function or method 'bwlabeln' for input arguments of
>> type 'double'.
>>
>> Error in ==> findcluster at 89
>>  [labelmat(spatdimlev, :, :), num] = bwlabeln(reshape(onoff(spatdimlev, :, :),
>>  nfreq, ntime), 4);
>>
>> Error in ==> clusterstat at 194
>>     posclusobs = findcluster(reshape(postailobs,
>>     [cfg.dim,1]),channeighbstructmat,cfg.minnbchan);
>>
>> Error in ==> statistics_montecarlo at 321
>>  [stat, cfg] = clusterstat(cfg, statrand, statobs,'issource',issource);
>>
>> Error in ==> statistics_wrapper at 285
>>   [stat, cfg] = statmethod(cfg, dat, cfg.design, 'issource',issource);
>>
>> Error in ==> ft_freqstatistics at 125
>> [stat, cfg] = statistics_wrapper(cfg, varargin{:});
>>
>> My configuration for ft_freqstatistics are
>> cfg = [];
>> cfg.channel          = 'all';
>> cfg.latency          = 'all';
>> cfg.frequency        = [30 50];
>> cfg.method           = 'montecarlo';
>> cfg.statistic        = 'indepsamplesT';
>> cfg.correctm         = 'cluster';
>> cfg.clusteralpha     = 0.025;
>> cfg.clusterstatistic = 'maxsum';
>> cfg.minnbchan        = 2;
>> cfg.tail             = 0;
>> cfg.clustertail      = 0;
>> cfg.alpha            = 0.025;
>> cfg.numrandomization = 500;
>> cfg.neighbours       = [ ];
>> design = [ones(1, 16), ones(1, 16)*2];
>> cfg.design           = design;
>> cfg.ivar             = 1;
>> [stat] = ft_freqstatistics(cfg, mig_grand_avg, con_grand_avg);
>>
>>
>> Data structure of con_grand_avg are
>>
>> con_grand_avg =
>>       label: {27x1 cell}
>>        freq: [30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50]
>>        time: [1x21 double]
>>      dimord: 'subj_chan_freq_time'
>>   powspctrm: [4-D double]
>>         cfg: [1x1 struct]
>>
>>
>> Any comment helps me.
>>
>>
>> Best regards,
>>
>> Ki-Young Jung
>>
>> _______________________________________________
>> fieldtrip mailing list
>> fieldtrip at donders.ru.nl
>> http://mailman.science.ru.nl/mailman/listinfo/fieldtrip
>
>
> _______________________________________________
> fieldtrip mailing list
> fieldtrip at donders.ru.nl
> http://mailman.science.ru.nl/mailman/listinfo/fieldtrip
>
>
>