Hello,

I am trying to run a correlational analysis with glm_fit on a masked area with the code:

foreach hemi (lh)

  foreach smoothness (10)

    foreach meas (volume thickness)

        mri_glmfit \

        --y {$hemi}.{$meas}.{$study}.{$smoothness}.mgh \

        --fsgd FSGD/{$study}.fsgd \

        --C Contrasts/Scores_Slope.mtx  \

        --surf fsaverage {$hemi}  \

        --mask lhYeo2011_7NetworksConfidence_N1000.mgz  \

        --glmdir {$hemi}.{$meas}.{$study}.{$smoothness}.glmdir

    end

  end

end