[Freesurfer] adding covariates to a longitudinal glmfit model