Thanks for the feedback!
To debug, we collapsed correct by contrasting the two values (errorCorrect-correct): 'Ext_t*ageatvisit+reward+meanmotion' as well as trying the suggested 'reward+Ext_t+ageatvisit+meanmotion' with MVM .
But both still fail the model test.
The datatable is on
pastebin, if that is useful.
Ultimately we would like to model with missing values. Not all subjects came back for every visit. (50 for visit 3, 142 for 2, 168 for visit 1).
Here is our 3dMVM call without visit and with the input file as the difference score (contrast errorcorrect - correct for each reward type)
3dMVM -overwrite -prefix ../ncanda_scripts/ncanda_3dmvm_outputs/3dmvm_output_extxagexcor_continuous.nii.gz \
-mask gm_50mask+tlrc \
-model 'Ext_t+ageatvisit+reward+meanmotion' \
-qVars "Ext_t,ageatvisit,meanmotion" \
-ranEff '~1' \
-jobs 8 \
-SS_type 3 \
-dataTable @datatable_3dlme_resp_diffscore2018524.txt
And it failed the model test:
Quote
Warning messages:
1: In process.MVM.opts(lop, verb = lop$verb) : NAs introduced by coercion
2: In process.MVM.opts(lop, verb = lop$verb) : NAs introduced by coercion
...
~~~~~~~~~~~~~~~~~~~ Model test failed! ~~~~~~~~~~~~~~~~~~~
Possible reasons:
....
** Error:
Quitting due to model test failure...