<div dir="ltr">Hey Julian,<div><br></div><div>thanks for the answer.<br><div><br></div><div>We are using some sort of an Intel iCore i7 with 16 GB of RAM as well as 40 GB of swap and Fedora 16.</div><div><br></div><div>We do interpolate the data on an MRI. In fact I'm not sure if that is the source of the problem. We downsampled the data and it did not change anything. The problem seems to be that there is a number generated that is too big for MATLAB to process it with the zeros(x) instruction.</div><div><br></div><div>Around 1-2 years ago I did nearly the same thing in fieldtrip but including an LCMV-Beamformer, the resulting data structures where much bigger and it worked without a problem. </div></div><div><br></div><div>Thanks, Martin</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Feb 23, 2015 at 5:14 PM, Julian Keil <span dir="ltr"><<a href="mailto:julian.keil@gmail.com" target="_blank">julian.keil@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word">Dear Martin,<div><br></div><div>what kind of machine are you using?</div><div>Did you interpolate your data to an MRI?</div><div>What is your grid resolution?</div><div><br></div><div>You have quite a high number of grid points that you want to compare.</div><div>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.</div><div><br></div><div>Good luck,</div><div><br></div><div>Julian</div><div><br><div>
<span style="border-collapse:separate;color:rgb(0,0,0);font-family:Helvetica;font-style:normal;font-variant:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:-webkit-auto;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;font-size:medium"><div><div style="word-wrap:break-word;font-size:12px">********************</div><div style="word-wrap:break-word"><b>Dr. Julian Keil</b></div><div style="word-wrap:break-word"><b><br></b></div>AG Multisensorische Integration<br>Psychiatrische Universitätsklinik<br>der Charité im St. Hedwig-Krankenhaus<br>Große Hamburger Straße 5-11, Raum E 307<br>10115 Berlin</div><div><br>Telefon: <a href="tel:%2B49-30-2311-1879" value="+493023111879" target="_blank">+49-30-2311-1879</a><br>Fax: <a href="tel:%2B49-30-2311-2209" value="+493023112209" target="_blank">+49-30-2311-2209</a> </div><div><a href="http://psy-ccm.charite.de/forschung/bildgebung/ag_multisensorische_integration" target="_blank">http://psy-ccm.charite.de/forschung/bildgebung/ag_multisensorische_integration</a></div></span>
</div>
<br><div><div>Am 23.02.2015 um 17:00 schrieb Martin Winkels:</div><br><blockquote type="cite"><div><div class="h5"><div dir="ltr"><div><div><div>Dear Fieldtrippers,<br><br></div>we encountered a problem during our DICS Beamformer-Statistics. <br><br></div>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). <br></div><div><br></div><div>The code we used is as follows:<br></div><div><br> cfg = [];<br> <br> cfg.method = 'montecarlo';<br> %cfg.statistic = 'depsamplesT';<br> cfg.statistic = 'ft_statfun_indepsamplesT';<br> cfg.correctm = 'cluster';<br> cfg.clusteralpha = 0.05; %ft default 0.05; <br> cfg.clusterstatistic = 'maxsum';<br> cfg.minnbchan = 2;<br> cfg.tail = 0;<br> cfg.clustertail = 0;<br> cfg.alpha = 0.025; %ft hat hier 0,025<br> <br> cfg.parameter = 'pow';<br> cfg.dim = grandavgA.dim;<br> <br> cfg.numrandomization = 1; % number of draws from the permutation distribution<br> <br> design(1,:) = [1 1 1 2 2 2];<br> design(2,:) = [1 1 1 1 1 1];<br> <br> cfg.design = design;<br> cfg.ivar = 1;<br> <br> stat = ft_sourcestatistics(cfg, grandavgA, grandavgB);<br> <br><br><br></div><div>The input data structure (grandavgA, grandavgB) is as follows:<br><br>grandavgA = <br><br> pow: [3x116380 double]<br> dim: [46 55 46]<br> inside: [116380x1 logical]<br> pos: [116380x3 double]<br> cfg: [1x1 struct]<br><br>grandavgB = <br><br> pow: [3x116380 double]<br> dim: [46 55 46]<br> inside: [116380x1 logical]<br> pos: [116380x3 double]<br> cfg: [1x1 struct]<br><br><br></div><div>Fieldtrip version: current (02/23/2015)<br><br><br></div><div>Thanks, Martin<br><br><p class="MsoNormal"><span style="font-size:10pt;font-family:"Courier New"">--
</span></p><p class="MsoNormal"><span style="font-size:10pt;font-family:"Courier New"">M.Sc. Martin Winkels</span></p><p class="MsoNormal"><span style="font-size:10pt;font-family:"Courier New"">Universitätsklinikum Münster</span></p><p class="MsoNormal"><span style="font-size:10pt;font-family:"Courier New"">Institut für Biomagnetismus & Biosignalanalyse</span></p><p class="MsoNormal"><span style="font-size:10pt;font-family:"Courier New"">Malmedyweg 15</span></p><p class="MsoNormal"><span style="font-size:10pt;font-family:"Courier New"">48149 Münster</span></p><p class="MsoNormal"><span style="font-size:10pt;font-family:"Courier New"">GERMANY</span></p><div><span style="font-size:10pt;font-family:"Courier New""> </span><br></div><p class="MsoNormal"><span style="font-size:10pt;font-family:"Courier New"">Telefon: <a href="tel:%2B49%20251%2083%2056%20846" value="+492518356846" target="_blank">+49 251 83 56 846</a></span></p>
<span style="font-size:10pt;font-family:"Courier New"">Web: <a href="http://biomag.uni-muenster.de/" target="_blank">http://biomag.uni-muenster.de</a></span><br></div><div><br></div></div></div></div>
_______________________________________________<br>fieldtrip mailing list<br><a href="mailto:fieldtrip@donders.ru.nl" target="_blank">fieldtrip@donders.ru.nl</a><br><a href="http://mailman.science.ru.nl/mailman/listinfo/fieldtrip" target="_blank">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip</a></blockquote></div><br></div></div><br>_______________________________________________<br>
fieldtrip mailing list<br>
<a href="mailto:fieldtrip@donders.ru.nl">fieldtrip@donders.ru.nl</a><br>
<a href="http://mailman.science.ru.nl/mailman/listinfo/fieldtrip" target="_blank">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip</a><br></blockquote></div><br></div>