I have a question about running mri_glmfit with 2 different surface modalities as the independent variable and a single value per participant as the dependent variable. I'm trying to examine the unique and shared contribution of the 2 surfaces to predict a cognitive measure.

E.g. Cognition ~ ?h.thickness + ?h.PET

I tried mri_glmfit with the "--table" option, giving it a column of scores per participant (in the aseg2stats table format), and then using "--pvr lh.thickness-stack.fsaverage.mgh" and a second "--pvr lh.PET-stack.fsaverage.mgh", using contrasts matrices as [0 0 1], [0 1 0], and [0 1 1].

ERROR: mri_reshape: number of elements cannot change
  nv1 = 163842, nv1 = 1

It continues to run, but all surface maps are full of zeroes. Any suggestions? Is there a better way to do this?

Thanks,

Adam