Dear AFNI users,

I have run into errors for running 3dMVM,

./3dMVM_runscript.sh
Loading required package: afex
Loading required package: lme4
Loading required package: Matrix
Loading required package: emmeans
************
Welcome to afex. For support visit: [afex.singmann.science]
- Functions for ANOVAs: aov_car(), aov_ez(), and aov_4()
- Methods for calculating p-values with mixed(): 'KR', 'S', 'LRT', and 'PB'
- 'afex_aov' and 'mixed' objects can be passed to emmeans() for follow-up tests
- Get and set global package options with: afex_options()
- Set orthogonal sum-to-zero contrasts globally: set_sum_contrasts()
- For example analyses see: browseVignettes("afex")
************

Attaching package: ‘afex’

The following object is masked from ‘package:lme4’:

lmer

Loading required package: phia
Loading required package: car
Error in `$<-.data.frame`(`*tmp*`, Subj, value = integer(0)) :
replacement has 0 rows, data has 532
Calls: $<- -> $<-.data.frame
Execution halted

The command lines i used are


3dMVM -prefix TBI_PVSinteraction \
-bsVars 'TBI*PVS+Gender+Age' \
-qVars 'Age' \
-num_glt 7 \
-gltLabel 1 TBI_PVSpresent -gltCode 1 'TBI : 1*TBI PVS : 1*PVS_present' \
-gltLabel 2 TBI_PVSabsent -gltCode 2 'TBI : 1*TBI PVS : 1*PVS_absent' \
-gltLabel 3 nonTBI_PVSpresent -gltCode 3 'TBI : 1*nonTBI PVS : 1*PVS_present' \
-gltLabel 4 nonTBI_PVSabsent -gltCode 4 'TBI : 1*nonTBI PVS : 1*PVS_absent' \
-gltLabel 5 TBI_by_PVS -gltCode 5 'TBI : 1*TBI -1*nonTBI PVS : 1*PVS_present -1*PVS_absent' \
-gltLabel 6 Gender_TBI_interaction -gltCode 6 'TBI : 1*TBI -1*nonTBI Gender : 1*M -1*F' \
-gltLabel 7 Gender_PVS_interaction -gltCode 7 'PVS : 1*PVS_present -1*PVS_absent Gender : 1*M -1*F' \
-dataTable \
ID TBI PVS Gender Age InputFile \
NCNC0062 TBI PVS_present F 44 /data/NEURO2/ping/SWI/QSM/NCNC0062/nonlinear_t2map_inQSM_MNI_sgima2.1233226_5mm_smooth.nii.gz .
.
.
NCNC9074 TBI PVS_present M 42 /data/NEURO2/ping/SWI/QSM/NCNC9074/nonlinear_t2map_inQSM_MNI_sgima2.1233226_5mm_smooth.nii.gz \