The MCZ tables were not built for fsaverage5. I think the best thing to do instead is a permutation test (since the MCZ is a little dodgy anyway:). Instructions are here:
https://surfer.nmr.mgh.harvard.edu/fswiki/FsTutorial/MultipleComparisonsV6.0...
You should run it from the development environment to get the up-to-date scripts.
It will take about 20 min to run using a single thread (maybe less with fsaverage5). If your computer has multiple cpus, you can specify that it run in parallel with --bg N where is the number of jobs (should not be more than number of cpus - 1). This will speed it up by a factor of N. You can also run it on launchpad using the secret --pbsubmit N flag. Don't pbsubmit mri_glmfit-sim, just specify --pbsubmit N and mri_glmfit-sim do submit all the jobs in the proper way.
On 3/22/18 5:18 PM, Tuominen, Lauri Johannes wrote:
Hello freesurfers!
My resting state data is in fsaverage5 space. mri_info gives dimensions of 3414 x 1 x 3 x 120 and voxelsize 1x1x1.
I do the group analysis and mri_glmfit creates a sig.mgh file that is has those same dimensions and results look good. Then I say:
mri_glmfit-sim --glmdir /autofs/space/bee_001/users/Lauri/test10.nosm.group/lh.surf.avg.rh.204 --3spaces --sim mc-z 10000 1.3 mc_1.3_3space --sim-sign abs —overwrite
And I get an error message saying dimension mismatch between input volume and seg. below is the tail of the message.
Thank you for your help! Lauri
number of voxels in search space = 9354 Done loading source values (nvtxs = 10242) overall max = 4.53945 at vertex 1977 overall min = -3.95722 at vertex 473 surface nvertices 10242 metric props tot surface area 66661.796875 group_avg_vtxarea_loaded 1 masked surface area 79140.789062 Computing voxel-wise significance CSDpvalMaxSigMap(): found 347/10242 above 0, max=0.485948 NOT Adjusting threshold for 1-tailed test thminadj = 1.3 Searching for Clusters ... thmin=1.300000 (1.300000), thmax=-1.000000 (-1), thsignid=0, minarea=0.000000 Found 46 clusters Max cluster size 1840.913696 Pruning by CW P-Value 0.05 Saving thresholded output to /autofs/space/bee_001/users/Lauri/test10.nosm.group/lh.surf.avg.rh.204/intercept/mc_1.3_3space.sig.masked.mgh Saving cluster numbers to /autofs/space/bee_001/users/Lauri/test10.nosm.group/lh.surf.avg.rh.204/intercept/mc_1.3_3space.sig.ocn.mgh Saving cluster pval /autofs/space/bee_001/users/Lauri/test10.nosm.group/lh.surf.avg.rh.204/intercept/mc_1.3_3space.sig.cluster.mgh Constructing output annotation Writing annotation /autofs/space/bee_001/users/Lauri/test10.nosm.group/lh.surf.avg.rh.204/intercept/mc_1.3_3space.sig.ocn.annot mri_segstats --seg /autofs/space/bee_001/users/Lauri/test10.nosm.group/lh.surf.avg.rh.204/intercept/mc_1.3_3space.sig.ocn.mgh --exclude 0 --i /autofs/space/bee_001/users/Lauri/test10.nosm.group/lh.surf.avg.rh.204.gm.nii.gz --avgwf /autofs/space/bee_001/users/Lauri/test10.nosm.group/lh.surf.avg.rh.204/intercept/mc_1.3_3space.y.ocn.dat --sum /tmp/mri_glmfit-sim.junk.8849
$Id: mri_segstats.c,v 1.121 2016/05/31 17:27:11 greve Exp $ cwd cmdline mri_segstats --seg /autofs/space/bee_001/users/Lauri/test10.nosm.group/lh.surf.avg.rh.204/intercept/mc_1.3_3space.sig.ocn.mgh --exclude 0 --i /autofs/space/bee_001/users/Lauri/test10.nosm.group/lh.surf.avg.rh.204.gm.nii.gz --avgwf /autofs/space/bee_001/users/Lauri/test10.nosm.group/lh.surf.avg.rh.204/intercept/mc_1.3_3space.y.ocn.dat --sum /tmp/mri_glmfit-sim.junk.8849 sysname Linux hostname bee.nmr.mgh.harvard.edu machine x86_64 user ljt22 UseRobust 0 Loading /autofs/space/bee_001/users/Lauri/test10.nosm.group/lh.surf.avg.rh.204/intercept/mc_1.3_3space.sig.ocn.mgh Loading /autofs/space/bee_001/users/Lauri/test10.nosm.group/lh.surf.avg.rh.204.gm.nii.gz ERROR: dimension mismatch between input volume and seg input 3414 1 3 seg 10242 1 1
Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer