AFNI Message Board

Dear AFNI users-

We are very pleased to announce that the new AFNI Message Board framework is up! Please join us at:

https://discuss.afni.nimh.nih.gov

Existing user accounts have been migrated, so returning users can login by requesting a password reset. New users can create accounts, as well, through a standard account creation process. Please note that these setup emails might initially go to spam folders (esp. for NIH users!), so please check those locations in the beginning.

The current Message Board discussion threads have been migrated to the new framework. The current Message Board will remain visible, but read-only, for a little while.

Sincerely, AFNI HQ

History of AFNI updates  

Previous Message Next Message
Message List New Topic Print View
|
cnweis
February 04, 2020 10:45AM
Hello,

I'm struggling to get my 3dLME to run. I keep getting stuck at the model test, with the suggestion that there's inappropriate coding in -model, -qVars, or -gltCode. I've combed through the message board and tried various things other people have had problems with but I can't seem to come to any resolution. Any insights as to what I may be doing wrong?

*I know the 3dLME command looks a little crazy with all the posthoc tests, but we're interested in the specific results of the various contrasts as they pertain to the 3 covariates (PCL_2wk, ETV_2wk, and PEDQ_2WK_TOTAL) and not just "controlling for their effects".

Thanks in advance!!


AFNI version: 
Precompiled binary linux_openmp_64: Jul  3 2019 (Version AFNI_19.2.01 'Claudius')


3dLME -prefix unc_clock_ETC_PEDQ -jobs 24 \
        -model "Valence*Predict+PCL_2wk+ETV_2wk+PEDQ_2WK_TOTAL+Age+Gender" \
        -qVars "Age,ETV_2wk,PCL_2wk,PEDQ_2WK_TOTAL" \
       -qVarsCenters "32.47, 9.01, 23.77, 28.74" \
        -ranEff "~1+Age" \
        -SS_type 3 \
        -num_glt 28 \
        -gltLabel 1 'neg-neu' -gltCode 1 'Valence : 1*neg -1*neu' \
        -gltLabel 2 'u-p' -gltCode 2 'Predict : 1*u -1*p' \
        -gltLabel 3 'Uneg-Uneu' -gltCode 3 'Valence : 1*neg -1*neu Predict : 1*u' \
        -gltLabel 4 'Pneg-Pneu' -gltCode 4 'Valence : 1*neg -1*neu Predict : 1*p' \
        -gltLabel 5 'Uneg-Pneg' -gltCode 5 'Predict : 1*u -1*p Valence : 1*neg' \
        -gltLabel 6 'Uneu-Pneu' -gltCode 6 'Predict : 1*u -1*p Valence : 1*neu' \
        -gltLabel 7 'Uneg-Pneu' -gltCode 7 'Predict : 1*u -1*p Valence : 1*neg -1*neu' \
        -gltLabel 8 'neg-neu_PCL_2wk' -gltCode 8 'Valence : 1*neg -1*neu PCL_2wk :' \
        -gltLabel 9 'u-p_PCL_2wk' -gltCode 9 'Predict : 1*u -1*p PCL_2wk :' \
        -gltLabel 10 'Uneg-Uneu_PCL_2wk' -gltCode 10 'Valence : 1*neg -1*neu Predict : 1*u PCL_2wk :' \
        -gltLabel 11 'Pneg-Pneu_PCL_2wk' -gltCode 11 'Valence : 1*neg -1*neu Predict : 1*p PCL_2wk :' \
        -gltLabel 12 'Uneg-Pneg_PCL_2wk' -gltCode 12 'Predict : 1*u -1*p Valence : 1*neg PCL_2wk :' \
        -gltLabel 13 'Uneu-Pneu_PCL_2wk' -gltCode 13 'Predict : 1*u -1*p Valence : 1*neu PCL_2wk :' \
        -gltLabel 14 'Uneg-Pneu_PCL_2wk' -gltCode 14 'Predict : 1*u -1*p Valence : 1*neg -1*neu PCL_2wk :' \
        -gltLabel 15 'neg-neu_ETV_2wk' -gltCode 15 'Valence : 1*neg -1*neu ETV_2wk :' \
        -gltLabel 16 'u-p_ETV_2wk' -gltCode 16 'Predict : 1*u -1*p ETV_2wk :' \
        -gltLabel 17 'Uneg-Uneu_ETV_2wk' -gltCode 17 'Valence : 1*neg -1*neu Predict : 1*u ETV_2wk :' \
        -gltLabel 18 'Pneg-Pneu_ETV_2wk' -gltCode 18 'Valence : 1*neg -1*neu Predict : 1*p ETV_2wk :' \
        -gltLabel 19 'Uneg-Pneg_ETV_2wk' -gltCode 19 'Predict : 1*u -1*p Valence : 1*neg ETV_2wk :' \
        -gltLabel 20 'Uneu-Pneu_ETV_2wk' -gltCode 20 'Predict : 1*u -1*p Valence : 1*neu ETV_2wk :' \
        -gltLabel 21 'Uneg-Pneu_ETV_2wk' -gltCode 21 'Predict : 1*u -1*p Valence : 1*neg -1*neu ETV_2wk :' \
        -gltLabel 22 'neg-neu_PEDQ_2WK_TOTAL' -gltCode 22 'Valence : 1*neg -1*neu PEDQ_2WK_TOTAL :' \
        -gltLabel 23 'u-p_PEDQ_2WK_TOTAL' -gltCode 23 'Predict : 1*u -1*p PEDQ_2WK_TOTAL :' \
        -gltLabel 24 'Uneg-Uneu_PEDQ_2WK_TOTAL' -gltCode 24 'Valence : 1*neg -1*neu Predict : 1*u PEDQ_2WK_TOTAL :' \
        -gltLabel 25 'Pneg-Pneu_PEDQ_2WK_TOTAL' -gltCode 25 'Valence : 1*neg -1*neu Predict : 1*p PEDQ_2WK_TOTAL :' \
        -gltLabel 26 'Uneg-Pneg_PEDQ_2WK_TOTAL' -gltCode 26 'Predict : 1*u -1*p Valence : 1*neg PEDQ_2WK_TOTAL :' \
        -gltLabel 27 'Uneu-Pneu_PEDQ_2WK_TOTAL' -gltCode 27 'Predict : 1*u -1*p Valence : 1*neu PEDQ_2WK_TOTAL :' \
        -gltLabel 28 'Uneg-Pneu_PEDQ_2WK_TOTAL' -gltCode 28 'Predict : 1*u -1*p Valence : 1*neg -1*neu PEDQ_2WK_TOTAL :' \
       -num_glf 2 \
       -glfLabel 1 'cond_valence' -glfCode 1 'Valence : 1*neg & 1*neu' \
       -glfLabel 2 'cond_pred' -glfCode 2 'Predict : 1*u & 1*p' \
        -dataTable @LME_table_unc_trl_ETV_PEDQ_PCL.txt



Here's a preview of my dataTable "LME_table_unc_trl_ETV_PEDQ_PCL.txt", in total there are 65 subjects with complete data. I've included 3 subjects below. 

Subj    Age     Gender  ETV_2wk PCL_2wk PEDQ_2WK_TOTAL  Valence Predict InputFile \
s270    25      0       4       41      51      neg     p       /raid-06/LS/Data/istar/270_day1/unc_fullN_SSWarp_trl_results/stats.270+tlrc[61] \
s270    25      0       4       42      52      neu     p       /raid-06/LS/Data/istar/270_day1/unc_fullN_SSWarp_trl_results/stats.270+tlrc[64] \
s270    25      0       4       43      53      neg     u       /raid-06/LS/Data/istar/270_day1/unc_fullN_SSWarp_trl_results/stats.270+tlrc[67] \
s270    25      0       4       44      54      neu     u       /raid-06/LS/Data/istar/270_day1/unc_fullN_SSWarp_trl_results/stats.270+tlrc[70] \
s283    24.6    1       0       36      40      neg     p       /raid-06/LS/Data/istar/283_day1/unc_fullN_SSWarp_trl_results/stats.283+tlrc[61] \
s283    24.6    1       0       37      41      neu     p       /raid-06/LS/Data/istar/283_day1/unc_fullN_SSWarp_trl_results/stats.283+tlrc[64] \
s283    24.6    1       0       38      42      neg     u       /raid-06/LS/Data/istar/283_day1/unc_fullN_SSWarp_trl_results/stats.283+tlrc[67] \
s283    24.6    1       0       39      43      neu     u       /raid-06/LS/Data/istar/283_day1/unc_fullN_SSWarp_trl_results/stats.283+tlrc[70] \
s294    39.4    0       10      40      45      neg     p       /raid-06/LS/Data/istar/294_day1/unc_fullN_SSWarp_trl_results/stats.294+tlrc[61] \
s294    39.4    0       10      40      45      neu     p       /raid-06/LS/Data/istar/294_day1/unc_fullN_SSWarp_trl_results/stats.294+tlrc[64] \
s294    39.4    0       10      40      45      neg     u       /raid-06/LS/Data/istar/294_day1/unc_fullN_SSWarp_trl_results/stats.294+tlrc[67] \
s294    39.4    0       10      40      45      neu     u       /raid-06/LS/Data/istar/294_day1/unc_fullN_SSWarp_trl_results/stats.294+tlrc[70] \
.
.
.
Reply Quote
RSS
Subject Author Posted

3dLME model test failure

cnweis February 04, 2020 10:45AM

Re: 3dLME model test failure

gang February 04, 2020 05:13PM