[FieldTrip] Fwd: question about ft_statistics errors

Julian Keil julian.keil at gmail.com
Fri Jan 14 10:04:19 CET 2011


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





More information about the fieldtrip mailing list