Hi Gang,

Thanks for your help last year. I'm finally at the point of running this analysis, and I've run into problems. I tried to follow your suggestion above about how to set up the model statement and data table, but the model test failed. I have the necessary R packages, and I believe the data table is set up correctly, as it is able to read the input files (see below; note I added a subject column because without it, I got the error "The content under -dataTable is not rectangular!"), but I suspect I haven't set up the rest of the script correctly. I wasn't sure if I needed a -gltCode line, and if so, how to write it. I'm interested in the effect of a behavior (cToM_Negative_A) on a specific contrast, PNM-CNM; the input files are the sub-brick with that contrast. Does the contrast need to be indicated somewhere in the model statement?

3dLMEr -prefix test_LME_cToMneg -jobs 24 \
-model 'cToM_Negative_A+(1|dyad)+(0+relationship|dyad)' \
-qVars 'cToM_Negative_A' \
-dataTable \
Subj dyad cToM_Negative_A cToM_Overall_A cToM_Positive_A relationship InputFile \
1 1 0 3 3 0 REDCAT256_stats_dm_REML+tlrc'[PNM-CNM_G#0_Coef]' \
2 2 0 4 4 0 REDCAT181_stats_dm_REML+tlrc'[PNM-CNM_G#0_Coef]' \
3 3 0 4 4 0 REDCAT184_stats_dm_REML+tlrc'[PNM-CNM_G#0_Coef]' \
4 4 0 1 1 0 REDCAT174_stats_dm_REML+tlrc'[PNM-CNM_G#0_Coef]' \
5 5 0 2 2 0 REDCAT270_stats_dm_REML+tlrc'[PNM-CNM_G#0_Coef]' \
6 8 0 3 3 0 REDCMNT166_stats_dm_REML+tlrc'[PNM-CNM_G#0_Coef]' \
...