I have several sessions from 30 subjects, and I want to test if the HDR changes from session 0 to 1 (S00 vs S01). So far, I'm trying to use the idea from Gang Chen's 2015 paper on this and am using 3dLME. I have missing data and only one group. When looking at the data, the HDR is decreased at most points along the curve in S01 compared to S00 (20sec total at 0.5sec TR, giving 39 time points). To test for this, I'm using the command:
3dLME -prefix lmeS00vsS01_3 -jobs 8 \
-model '0+Time+Session' \
-ranEff '~1' \
-num_glf 1 \
-qVars Order \
-qVarCenters 0 \
-corSt 'Order : AR1' \
-SS_type 3 \
-num_glt 1 \
-gltLabel 1 'S00minusS01' -gltCode 1 'Session : 1*S00 -1*S01' \
-dataTable Subj Session Time Order InputFile \
P0204 S00 t00 0 iresp_stim01.P0204_S00_masked+tlrc[00] \
P0204 S00 t01 1 iresp_stim01.P0204_S00_masked+tlrc[01] \
P0204 S00 t02 2 iresp_stim01.P0204_S00_masked+tlrc[02] \
...
Is this the appropriate way to test for this difference?
Thanks,
Trey
