[FieldTrip] ft_clusterstat OUT OF MEMORY - DICS

[Julian Keil]

Dear Martin,

what kind of machine are you using?
Did you interpolate your data to an MRI?
What is your grid resolution?

You have quite a high number of grid points that you want to compare.
So in case you run out of memory, I'd suggest not interpolating to an MRI (in case you have done this) but to stay on the grid-point level for your stats. Otherwise, you could use a less dense grid which obviously results in smaller data structures.

Good luck,

Julian

********************
Dr. Julian Keil

AG Multisensorische Integration
Psychiatrische Universitätsklinik
der Charité im St. Hedwig-Krankenhaus
Große Hamburger Straße 5-11, Raum E 307
10115 Berlin

Telefon: +49-30-2311-1879
Fax:        +49-30-2311-2209
http://psy-ccm.charite.de/forschung/bildgebung/ag_multisensorische_integration

Am 23.02.2015 um 17:00 schrieb Martin Winkels:

> Dear Fieldtrippers,
> 
> we encountered a problem during our DICS Beamformer-Statistics.
> 
> After calculating a beamformer (DICS), normalisation and building grandaverages across subjects (here exemplarily 3 subjects) we try to calculate cluster based permutation statistic (in this study: between groups - one condition).
> 
> The code we used is as follows:
> 
>         cfg = [];
> 
>         cfg.method           = 'montecarlo';
>         %cfg.statistic        = 'depsamplesT';
>         cfg.statistic        = 'ft_statfun_indepsamplesT';
>         cfg.correctm         = 'cluster';
>         cfg.clusteralpha     = 0.05; %ft default 0.05;
>         cfg.clusterstatistic = 'maxsum';
>         cfg.minnbchan        = 2;
>         cfg.tail             = 0;
>         cfg.clustertail      = 0;
>         cfg.alpha            = 0.025; %ft hat hier 0,025
> 
>         cfg.parameter   = 'pow';
>         cfg.dim         =  grandavgA.dim;
> 
>         cfg.numrandomization = 1;      % number of draws from the permutation distribution
> 
>         design(1,:) = [1 1 1 2 2 2];
>         design(2,:) = [1 1 1 1 1 1];
> 
>         cfg.design = design;
>        cfg.ivar   = 1;
> 
>         stat = ft_sourcestatistics(cfg, grandavgA, grandavgB);
> 
> 
> 
> The input data structure (grandavgA, grandavgB) is as follows:
> 
> grandavgA =
> 
>        pow: [3x116380 double]
>        dim: [46 55 46]
>     inside: [116380x1 logical]
>        pos: [116380x3 double]
>        cfg: [1x1 struct]
> 
> grandavgB =
> 
>        pow: [3x116380 double]
>        dim: [46 55 46]
>     inside: [116380x1 logical]
>        pos: [116380x3 double]
>        cfg: [1x1 struct]
> 
> 
> Fieldtrip version: current (02/23/2015)
> 
> 
> Thanks, Martin
> 
> --
> 
> M.Sc. Martin Winkels
> 
> Universitätsklinikum Münster
> 
> Institut für Biomagnetismus & Biosignalanalyse
> 
> Malmedyweg 15
> 
> 48149 Münster
> 
> GERMANY
> 
> 
> Telefon: +49 251 83 56 846
> 
> Web: http://biomag.uni-muenster.de
> 
> _______________________________________________
> fieldtrip mailing list
> fieldtrip at donders.ru.nl
> http://mailman.science.ru.nl/mailman/listinfo/fieldtrip

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20150223/0e072034/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 495 bytes
Desc: Message signed with OpenPGP using GPGMail
URL: <http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20150223/0e072034/attachment-0001.sig>