Dear Doug,
I was running the
mri_glmfit-sim \ --glmdir GroupAnalysis_abs.glmdir \ --sim perm 1000 3 permcsd \ --sim-sign abs \ --cwpvalthresh .05
just to get an idea of how my FA Group Analysis data is looking and I got an error before to get the permcsd.sig.cluster.summary file from the mri_volcluster: mri_volcluster --in GroupAnalysis_abs.glmdir/contrast/sig.mgh --mask GroupAnalysis_abs.glmdir/mask.mgh --reg /usr/local/freesurfer/subjects/fsaverage/mri.2mm/reg.2mm.dat --no-fixmni --cwsig GroupAnalysis_abs.glmdir/contrast/permcsd.sig.cluster.mgh --sum GroupAnalysis_abs.glmdir/contrast/permcsd.sig.cluster.summary --ocn GroupAnalysis_abs.glmdir/contrast/permcsd.sig.ocn.mgh --cwpvalthresh .05 --seg fsaverage aparc+aseg.mgz --csdpdf GroupAnalysis_abs.glmdir/contrast/permcsd.pdf.dat --csd GroupAnalysis_abs.glmdir/csd/permcsd.j001-contrast.csd --vwsig GroupAnalysis_abs.glmdir/contrast/permcsd.sig.voxel.mgh ERROR: cannot find /media/2Tra/Antonella_2011/DTI_FS/GLM_CVS/fsaverage/mri/aparc+aseg.mgz
I replaced the fsaverage/mri/aparc+aseg.mgz with the cvs_avg35/mri/aseg.mgz and to get the cluster summary I was running:
mri_volcluster --in GroupAnalysis_abs.glmdir/contrast/sig.mgh --mask GroupAnalysis_abs.glmdir/mask.mgh --reg /usr/local/freesurfer/subjects/fsaverage/mri.2mm/reg.2mm.dat --no-fixmni --cwsig GroupAnalysis_abs.glmdir/contrast/permcsd.sig.cluster.mgh --sum GroupAnalysis_abs.glmdir/contrast/permcsd.sig.cluster.summary --ocn GroupAnalysis_abs.glmdir/contrast/permcsd.sig.ocn.mgh --cwpvalthresh .05 --seg cvs_avg35 aseg.mgz --csdpdf GroupAnalysis_abs.glmdir/contrast/permcsd.pdf.dat --csd GroupAnalysis_abs.glmdir/csd/permcsd.j001-contrast.csd --vwsig GroupAnalysis_abs.glmdir/contrast/permcsd.sig.voxel.mgh
but then I am getting the following error:
NOT Adjusting threshold for 1-tailed test Using defalt ctab /usr/local/freesurfer/FreeSurferColorLUT.txt INFO: loading mask volume: GroupAnalysis_abs.glmdir/mask.mgh INFO: MRIbinarize01(): nhits = 961260 Found 961260 voxels in mask ERROR: Input volume voxel dimensions do not match those in the registration file. If the input volume is in bshort/bfloat format, check that there is an accompanying bhdr file.
Can you please help me to solve this problem. I am not sure what to do.
Thank you very much! Antonella
PS. Just to remember you that I registered my FA data over the CVS space but I have for each subject from dt_recom the register.dat file so I suppose I need to keeo the fsaverage in the --reg COMMAND-LINE ARGUMENTS.