<div dir="ltr">Many thanks, once again Nick! The full pipeline now works for me. </div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Apr 16, 2015 at 8:50 AM, Nick Oosterhof <span dir="ltr"><<a href="mailto:n.n.oosterhof@googlemail.com" target="_blank">n.n.oosterhof@googlemail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class=""><br>
On 12 Apr 2015, at 00:53, Mike E. Klein <<a href="mailto:michaeleklein@gmail.com">michaeleklein@gmail.com</a>> wrote:<br>
<br>
</span>> I’m trying to run a sample top-level analysis on surface-based searchlight data using Nick’s afni_surface_alphasim script. I’ve never tried this before, and am getting a ZeroDivisionError (bottom of output). […]<br>
<span class="">><br>
> Here is my bash input/output:<br>
</span>> […]<br>
<span class="">> ++ 3dbucket: AFNI version=AFNI_2011_12_21_1014 (May  1 2014) [64-bit]<br>
> ** FATAL ERROR: Dataset /Volumes/deep_thought/fmri/mvpa2/april_analyses/subs/06/suma_surfaces_deoblique/zeroed_soundVrest_NEW-minimal_ico32-128_mh_50vx.niml.dset differs in size from first one<br>
<br>
</span>Actually the issues turned out to be that a an anatomical surface file had to be specified. I’ve modified the code now to throw a more informative error message if that case is detected, and sent a PR [1] which will probably be merged soon. Also the padding actually had to go to 20483 to make it work, and increasing the smoothing bandwidth made it run significantly faster. So briefly, the options I added that work were:<br>
<br>
-N 20483 --surface_file intermediate_surf.asc -S 1.4<br>
<br>
where “-N” indicates the node padding (=(32^2*10+2)*2-1), and “-S” the smoothing bandwidth.<br>
<br>
Also, you may want to consider Threshold Free Cluster Enhancement [2]. Not in PyMVPA yet, but if you are happy to go into Matlab/Octave world, it is provided in CoSMoMVPA [3] through a monte carlo simulation [4]. There is an example script [5] illustrating its use (from "%% Run Threshold-Free Cluster Enhancement (TFCE)” onwards). Note: you would have to pad the nodes manually and convert to ASCII format, then use cosmo_surface_dataset to load each dataset and cosmo_stack to make it a single dataset.<br>
<br>
[1] <a href="https://github.com/PyMVPA/PyMVPA/pull/285" target="_blank">https://github.com/PyMVPA/PyMVPA/pull/285</a><br>
<span class="">[2] <a href="http://www.fmrib.ox.ac.uk/analysis/techrep/tr08ss1/tr08ss1.pdf" target="_blank">http://www.fmrib.ox.ac.uk/analysis/techrep/tr08ss1/tr08ss1.pdf</a><br>
[3] <a href="http://cosmomvpa.org" target="_blank">cosmomvpa.org</a><br>
</span>[4] <a href="https://github.com/CoSMoMVPA/CoSMoMVPA/blob/master/mvpa/cosmo_montecarlo_cluster_stat.m" target="_blank">https://github.com/CoSMoMVPA/CoSMoMVPA/blob/master/mvpa/cosmo_montecarlo_cluster_stat.m</a><br>
[5] <a href="http://cosmomvpa.org/contents_demo.html#demo-surface-tfce" target="_blank">http://cosmomvpa.org/contents_demo.html#demo-surface-tfce</a><br>
<br>
</blockquote></div><br></div>