On Jul 09, 2016, at 09:00 AM, freesurfer-request@nmr.mgh.harvard.edu wrote:
Send Freesurfer mailing list submissions to freesurfer@nmr.mgh.harvard.edu
To subscribe or unsubscribe via the World Wide Web, visit https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer or, via email, send a message with subject or body 'help' to freesurfer-request@nmr.mgh.harvard.edu
You can reach the person managing the list at freesurfer-owner@nmr.mgh.harvard.edu
When replying, please edit your Subject line so it is more specific than "Re: Contents of Freesurfer digest..."
Today's Topics:
1. Longitudinal LME contrast vectors (K) 2. Re: Longitudinal LME contrast vectors (Martin Reuter) 3. Question about NULL event with Optseq2 (Feng Zhou) 4. Cerebral volume without ventricles (Kaiming Yin) 5. Re: FSFAST error: ERROR: fast_selxavg3() failed\n (Douglas Greve)
----------------------------------------------------------------------
Message: 1 Date: Fri, 8 Jul 2016 15:02:58 -0700 From: K kam.ntnu@gmail.com Subject: [Freesurfer] Longitudinal LME contrast vectors To: freesurfer@nmr.mgh.harvard.edu Message-ID: CAH54DzgqKWSCwm0gu_ecbB3jrQmLc=Y6Nun5j4DujKwaxn78UA@mail.gmail.com Content-Type: text/plain; charset="utf-8"
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!