Hi everyone!
I'm looking to set up a 3lme with two covariates of interest that are continuous variables. Below is my model and the error message that I'm getting. Any ideas what I could be doing wrong? The script runs fine if I only use one covariate of interest under my gltcode. Thanks!!
3dLME -prefix LME_TaskxSexxRunxPercHDDxCravProv_covaryAgeDaysLastDrink.nii.gz -jobs 30 \
-mask /home/mr2649/YSC1/bioimagesuite35/images/MNI_T1_2mm_mask.nii.gz \
-model 'Age+DaysDrk+Task*Sex*Time*percHDD*CravP' \
-SS_type 3 \
-qVars 'Age,DaysDrk,percHDD,CravP' \
-qVarCenters '35.9,17.26,11.35,0' \
-ranEff '~1+CravP' \ \
-resid LME_TaskxSexxRunxPercHDDxCravProv_covaryAgeDaysLastDrink_resid.nii.gz \
-num_glt 2 \
-gltLabel 1 'Male_stress_HDD_CravP' -gltCode 1 'Gender : 1*Male Task : 1*stress percHDD : CravP : ' \
-gltLabel 2 'Male_drug_HDD_CravP' -gltCode 2 'Gender : 1*Male Task : 1*drug percHDD : CravP :' \
-dataTable \
Subj Sex Age percHDD DaysDrk CravP Time Task InputFile \
pa0080 Male 47 34 85 -0.14 early stress ../../../Common_paXXX_stress_orig.nii.gz \
pa0086 Male 27 12 31 -2.14 early stress ../../../Common_paXXX_stress_orig.nii.gz \
...
Sex = male or female
Time = early, mid or late
Task = stress, drug or neutral
These are just sample subjects and contrasts.
** Error:
Incorrect variable name in GLF coding: 1*stress
Warning messages:
1: In if (is.na(code[[n]][QVpos + 2])) { :
the condition has length > 1 and only the first element will be used
2: In if (QVpos == 1) outList[[1]][[n]] <- NA else outList[[1]][[n]] <- glfConstr(code[[n]][-(QVpos:(QVpos + :
the condition has length > 1 and only the first element will be used
3: In QVpos:(QVpos + 2) :
numerical expression has 2 elements: only the first used
4: In QVpos:(QVpos + 2) :
numerical expression has 2 elements: only the first used
