Hello freesurfers,
I would like to use mri_glmfit with --mask but I get this error: dimension mismatch 1 between y and mask. I created a binary mask from 22 aparc labels (using mri_binarize) and I would like to run a surface-based analysis only on those regions. The command I use is:
mri_glmfit --y lh.g2v0.thickness.15.mgh --fsgd g2v0.fsgd dods --C lh.group.diff.mtx --surf fsaverage lh --mask lh22mask.mgz --glmdir lh.group_diff_g2v0.glmdir
The input file was generated with this command: mris_preproc --fsgd g2v0.fsgd --cache-in thickness.fwhm15.fsaverage --target fsaverage --hemi lh --out lh.g2v0.thickness.15.mgh
From the error, it seems that I am not using the right input file. What
kind of input file should be used with a mask?
I also thought to use --label but I am not sure how to combine my 22 labels in one single label. mris_label2annot can combine them but the output is an annotation file not a label.
As a note, we will want to do FDR correction, not MC.
Thanks! Mihaela