[FieldTrip] Fwd: question about ft_statistics errors

[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