[FieldTrip] Cluster-level permutation test problem

[Ana Laura Diez Martini]

Dear all,

I'm trying to run a cluster-level permutation test with an old script in
which I used the neighbourhood map of a Biosemi cap. The only thing I have
changed there is the neighbourhood that is now for a GSN-HydroCel-129 one,
so I just re wrote the neighbours for each channel the same way I did for
the Biosemi one. I use an old version of Fieldtrip here which was working
fine, but even with the newest versions and modifications, it doesn't work.
What do you think I could be doing wrong?

Thank you!!


*This is the output I get:*

Selecting and formatting the data.
selected 128 channels
selected 1000 time bins
selected 1 frequency bins
Calculating the neighbourhood structure of the channels.
Running the statistics engine.
Statistic-specific preprocessing (calculating critical values, initializing
draws from the randomization distribution, ...).
randomization 0 1000
Undefined function 'bwlabeln' for input arguments of type 'double'.

Error in findcluster (line 73)
  [labelmat(spatdimlev, :, :), num] = bwlabeln(reshape(onoff(spatdimlev, :,
:), nfreq, ntime), 4);

Error in clusterrandstatistics (line 851)
            [posclusterslabelmat,nposclusters] =
findcluster(onoff,data.channeighbstructmat,cfg.minnbchan);

Error in clusterrandanalysis (line 462)
[clusrand] = clusterrandstatistics(cfg, data);

Error in clusterrand_stats_ERP_CW (line 347)
  [clusrand] = clusterrandanalysis(cfg,GavgI,GavgC);


*And this is my script*

function [clusrand] = clusterrand_stats_ERP_WC()
S = load('4WC_1.avg','-MAT');
avgC41 = S.avgC;
avgI41 = S.avgI;

S = load('4WC_2.avg','-MAT');
avgC42 = S.avgC;
avgI42 = S.avgI;

etc...


cfg =[];
cfg.keepindividual = 'yes';

GavgC = timelockgrandaverage(cfg,avgC41,avgC42,.....);
GavgI = timelockgrandaverage(cfg,avgI41,avgI42.....);

save('grand_averages','GavgC','GavgI');

  cfg.latency = [0 1];
  cfg.channel = {'all'};
  cfg.avgovertime = 'no';

  cfg.statistic = 'depsamplesT';
  cfg.alpha = 0.05;
  cfg.clusterteststat = 'maxsum';

  cfg.makeclusters = 'yes';
  cfg.alphathresh = 0.05;
  cfg.minnbchan    = 0;

  cfg.elecfile = 'GSN-HydroCel-129.sfp';

  cfg.neighbours{1}.label = 'E1';
  cfg.neighbours{1}.neighblabel = {'E125', 'E8', 'E2', 'E122','E121'};
  cfg.neighbours{2}.label = 'E2';
  cfg.neighbours{2}.neighblabel = {'E1', 'E8', 'E14',
'E9','E3','E123','122'};
  cfg.neighbours{3}.label = 'E3';
  cfg.neighbours{3}.neighblabel = {'E2', 'E123', 'E124','E4','E10','E9'};

etc...


 cfg.nranddraws = 1000;

 [clusrand] = clusterrandanalysis(cfg,GavgI,GavgC);

save clusrand_WC clusrand;

end
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20140606/ea3f3d11/attachment-0001.html>