Hello,
I was trying to use 3dlme to analyze our longitudinal data but the model failed to run. In brief, we have one within-subject variable (ScanPoint, 6 levels with missing data), one within-subject nuisance variable (Scanner, 14 levels with missing data), and one between-subject variable of interest (Group, 3 levels). We are interested in the main effects of Group and ScanPoint. I would like to model random intercept for each subject and random slopes for ScanPoint and Scanner. Here is part of my code:
3dLME -prefix compare_volume.nii \
-jobs 10 \
-model "Group*ScanPoint+Sex+Age+TIV" \
-qVars "Age,TIV" \
-ranEff "~1+Scanner ~1+ScanPoint" \
-SS_type 3 \
-num_glt 3 \
-gltLabel 1 "Conv+Nonconv>Controls" -gltCode 1 "Group : 1*conv +1*nonconv -1*control" \
-gltLabel 2 "Conv>Nonconv+Controls" -gltCode 2 "Group : 1*conv -1*nonconv -1*control" \
-gltLabel 3 "Conv>Nonconv" -gltCode 3 "Group : 1*conv -1*nonconv" \
-num_glf 2 \
-glfLabel 1 'StartPoint_maineffect' -glfCode 1 'Group : 1*conv -1*nonconv & 1*conv -1*control ScanPoint : 1*p1' \
-glfLabel 2 'EndPoint_maineffect' -glfCode 2 'Group : 1*conv -1*nonconv & 1*conv -1*control ScanPoint : 1*p6' \
-dataTable \
Subj Group Sex Age TIV Scanner ScanPoint InputFile \
sub1 nonconv 2 19 1515.41000000000 s1 p1 smwp101_S0301_00_T1.nii \
sub1 nonconv 2 19 1515.41000000000 s1 p5 smwp101_S0301_08_T1.nii \
sub1 nonconv 2 19 1515.41000000000 s2 p6 smwp101_S0301_OSA1_T1.nii \
sub2 control 1 13 1530.95000000000 s1 p1 smwp101_S0310_00_T1.nii \
sub2 control 1 13 1530.95000000000 s1 p2 smwp101_S0310_02_T1.nii \
sub2 control 1 13 1530.95000000000 s1 p3 smwp101_S0310_04_T1.nii \
sub2 control 1 13 1530.95000000000 s1 p4 smwp101_S0310_06_T1.nii \
sub2 control 1 13 1530.95000000000 s1 p5 smwp101_S0310_08_T1.nii \
sub3 conv 2 16 1486.07000000000 s2 p1 smwp101_S0333_00_T1.nii \
sub3 conv 2 16 1486.07000000000 s2 p2 smwp101_S0333_02_T1.nii \
sub3 conv 2 16 1486.07000000000 s2 p3 smwp101_S0333_04_T1.nii \
sub3 conv 2 16 1486.07000000000 s2 p4 smwp101_S0333_06_T1.nii \
...
When the code was running, it stuck at here with no response:
"If the program hangs here for more than, for example, half an hour,
kill the process because the model specification or the GLT coding
is likely inappropriate."
I am wondering if this is because something is wrong with my code? If so, how can I change it?
Thank you!
Best,
Hengyi
