Hi,
I am looking at resting state connectivity with an amygdala ROI longitudinally and I want to see if amygdala connectivity changes with age (inverse,invcage). I am also including mean absolute motion displacement as a nuisance covariate (motion). I tried to set up the model as below, but 3dlme gets hung up and does not continue running. I think I am doing something incorrect with the nuisance covariate (motion)? I've already centered age prior to including it in the model. Do you have any suggestions on what I am doing wrong?
Thank you for your time!
Best, Maria
data=/data/Maria
conditions="Ramyg Lamyg"
for condition in $conditions; do
3dLME -prefix $data/stats/20150206_3dlme_invcage_mot/3dLME_${condition}_invcage -jobs 8 \
-model 'invcage+motion' \
-qVars 'invcage, motion' \
-ranEff '~1' \
-dataTable \
Subj invcage motion InputFile
10152 0.000136066 0.035470598 /data/Maria/10152_20100514/10152_20100514_${condition}_REML_r_Z+tlrc.[0] \
10152 -0.004708593 0.027585698 /data/Maria/10152_20111123/10152_20111123_${condition}_REML_r_Z+tlrc.[0] \
10153 0.008060631 0.025299051 /data/Maria/10153_20100604/10153_20100604_${condition}_REML_r_Z+tlrc.[0] \
10153 0.002051341 0.035811259 /data/Maria/10153_20111128/10153_20111128_${condition}_REML_r_Z+tlrc.[0] \
10156 -0.000361161 0.02155011 /data/Maria/10156_20100313/10156_20100313_${condition}_REML_r_Z+tlrc.[0] \
10156 -0.005135271 0.030034241 /data/Maria/10156_20110923/10156_20110923_${condition}_REML_r_Z+tlrc.[0] \
10173 -0.001777112 0.027198092 /data/Maria/10173_20100313/10173_20100313_${condition}_REML_r_Z+tlrc.[0] \
10173 -0.006075057 0.03254388 /data/Maria/10173_20110820/10173_20110820_${condition}_REML_r_Z+tlrc.[0] \
10181 0.005629966 0.12263516 /data/Maria/10181_20090604/10181_20090604_${condition}_REML_r_Z+tlrc.[0] \
10181 7.01E-05 0.051321493 /data/Maria/10181_20101124/10181_20101124_${condition}_REML_r_Z+tlrc.[0] \