Hello. I am having an issue with 3dMEMA with a covariate that I have been unable to solve. Can anyone see what I am doing wrong? The message is:
Have lop$centerType 1
Have lop$centerVal[ii] 28 for bmi
Warning message:
In cbind(lop$xMat, lop$covData) :
number of rows of result is not a multiple of vector length (arg 1)
Warning in DOF.check.MEMA(lop) : Too few subjects for options selected.
Error:
Execution halted
Here is the command:
3dMEMA -prefix x_3dMEMA_CAF_H_OLD_BMI_COV_eqV \
-jobs 4 \
-mask x_OLD_Overall_Mask+tlrc \
-set OlderAdults \
DRU 'DRU_2runs_H_P_bucket+tlrc[18]' 'DRU_2runs_H_P_bucket+tlrc[19]' \
WL 'WL_2runs_H_P_bucket+tlrc[18]' 'WL_2runs_H_P_bucket+tlrc[19]' \
AP 'AP_2runs_H_P_bucket+tlrc[18]' 'AP_2runs_H_P_bucket+tlrc[19]' \
JTR 'JTR_2runs_H_P_bucket+tlrc[18]' 'JTR_2runs_H_P_bucket+tlrc[19]' \
EPO 'EPO_2runs_H_P_bucket+tlrc[18]' 'EPO_2runs_H_P_bucket+tlrc[19]' \
CS 'CS_2runs_H_P_bucket+tlrc[18]' 'CS_2runs_H_P_bucket+tlrc[19]' \
KJ 'KJ_2runs_H_P_bucket+tlrc[18]' 'KJ_2runs_H_P_bucket+tlrc[19]' \
NS 'NS_2runs_H_P_bucket+tlrc[18]' 'NS_2runs_H_P_bucket+tlrc[19]' \
LF 'LF_2runs_H_P_bucket+tlrc[18]' 'LF_2runs_H_P_bucket+tlrc[19]' \
CEG 'CEG_2runs_H_P_bucket+tlrc[18]' 'CEG_2runs_H_P_bucket+tlrc[19]' \
ESH 'ESH_2runs_H_P_bucket+tlrc[18]' 'ESH_2runs_H_P_bucket+tlrc[19]' \
LD 'LD_2runs_H_P_bucket+tlrc[18]' 'LD_2runs_H_P_bucket+tlrc[19]' \
TK 'TK_2runs_H_P_bucket+tlrc[18]' 'TK_2runs_H_P_bucket+tlrc[19]' \
JD 'JD_2runs_H_P_bucket+tlrc[18]' 'JD_2runs_H_P_bucket+tlrc[19]' \
MRJ 'MRJ_2runs_H_P_bucket+tlrc[18]' 'MRJ_2runs_H_P_bucket+tlrc[19]' \
TKA 'TKA_2runs_H_P_bucket+tlrc[18]' 'TKA_2runs_H_P_bucket+tlrc[19]' \
TMA 'TMA_2runs_H_P_bucket+tlrc[18]' 'TMA_2runs_H_P_bucket+tlrc[19]' \
KB 'KB_2runs_H_P_bucket+tlrc[18]' 'KB_2runs_H_P_bucket+tlrc[19]' \
DBA 'DBA_2runs_H_P_bucket+tlrc[18]' 'DBA_2runs_H_P_bucket+tlrc[19]' \
MAR 'MAR_2runs_H_P_bucket+tlrc[18]' 'MAR_2runs_H_P_bucket+tlrc[19]' \
DJO 'DJO_2runs_H_P_bucket+tlrc[18]' 'DJO_2runs_H_P_bucket+tlrc[19]' \
TDO 'TDO_2runs_H_P_bucket+tlrc[18]' 'TDO_2runs_H_P_bucket+tlrc[19]' \
-missing_data 0 \
-HKtest \
-covariates x_OLD_BMI_COVAR_FILE \
-covariates_center bmi=28
And here is the covariate file:
subj bmi
TKA 19.37
MAR 20.04
DJO 20.64
KB 22.59
JD 22.77
TK 24.33
TMA 25.49
DBA 26.05
MRJ 26.99
LD 28.08
TDO 26.9
LF 25.68
AP 25.82
JTR 27.15
CS 30.47
EPO 31.45
NS 32.26
ESH 34.59
WL 35.48
DRU 35.96
KJ 40.17
CEG 27.54