Hello! We are trying to run a longitudinal morphometry analysis with a
control group and clinical group scanned twice. We'd like to look at group
x time interactions and would like help figuring out the best way to do so.
As the time variable we use centered age -- the deviation from the mean age
across all subjects and timepoints. To do so, we use the following
lme model:
intercept + centeredage.controls + centeredage.clinical + group(0,1)
+ sex(0,1)
Contrast vector: [0 1 -1 0 0]
where the centeredage column contains only the centered age for that group
and zeros for the opposite group. The thinking behind this is to
investigate which cortical regions show opposite effects of time
between the groups.
However, the FreeSurfer lme software requests a specified "time
column". Obviously, we cannot add a centered age for all subjects, nor does
it seem to make sense to use one of the two group-specific
centeredage columns as the time column.
Our alternative model for investigating interatction is:
intercept + centeredage + group(0,1) + sex(0,1) + groupxcenteredage
where centeredage is the time columns and is a common regressor for both
groups and thus, groupxcenteredage becomes 0 for control subjects.
Contrast vector: [0 0 0 0 1]
This would test whether there is an added effect of centered age in the
clinical group.
Is there a better way to do an group x time interaction analysis here, or
is there a solution for the time column? Thanks!