Dear Afni experts,
I set up a 3dMVM script to find effects with 2 (med) x 2 (cue) x 3 (TENT-basis), all within-ss factors (attached at the bottom). However, the script runs a bit, but throws an error (below). It would be greatly appreciated if there is any problem with the model that I have.
Thank you very much in advance!
Error:
[1] "Great, test run passed at voxel (20, 36, 30)!"
[1] "Start to compute 61 slices along Z axis. You can monitor the progress"
[1] "and estimate the total run time as shown below."
[1] "03/30/16 17:35:50.294"
Z slice 1 done: 03/30/16 17:35:50.425
Z slice 2 done: 03/30/16 17:35:50.515
Error in if (p_wsmvt < uvfm$anova[ii, "Pr(>F)"]) out[ii] <- qf(p_wsmvt, :
argument is of length zero
Calls: aperm -> apply -> FUN
In addition: Warning message:
In summary.Anova.mlm(object$Anova, multivariate = FALSE) :
HF eps > 1 treated as 1
Execution halted
My script:
3dMVM -prefix ../ANOVA_results/3dMVM_VvsN_CueRetTENT-2-4TRs_x_CueType_x_med \
-jobs 1 \
-wsVars "med*cueType*Time" \
-wsE2 \
-mask ../ANOVA_results/mask_group_N=22_fromNative+tlrc. \
-num_glt 8 \
-gltLabel 1 'D-P' -gltCode 1 'med : 1*D -1*P' \
-gltLabel 2 'Valid-Neutral' -gltCode 2 'cueType : 1*V -1*N' \
-gltLabel 3 'Valid-Neutral_4s' -gltCode 3 'cueType : 1*V -1*N Time : 1*t2' \
-gltLabel 4 'Valid-Neutral_6s' -gltCode 4 'cueType : 1*V -1*N Time : 1*t3' \
-gltLabel 5 'Valid-Neutral_8s' -gltCode 5 'cueType : 1*V -1*N Time : 1*t4' \
-gltLabel 6 'CueType_x_med_4s' -gltCode 6 'cueType : 1*V -1*N med : 1*D -1*P Time : 1*t2' \
-gltLabel 7 'CueType_x_med_6s' -gltCode 7 'cueType : 1*V -1*N med : 1*D -1*P Time : 1*t3' \
-gltLabel 8 'CueType_x_med_8s' -gltCode 8 'cueType : 1*V -1*N med : 1*D -1*P Time : 1*t4' \
-num_glf 2 \
-glfLabel 1 'CueRetTime_x_CueType' -glfCode 1 'Time : 1*t2 -1*t3 & 1*t2 -1*t4 cueType : 1*V -1*N' \
-glfLabel 2 'CueRetTime_x_med' -glfCode 2 'Time : 1*t2 -1*t3 & 1*t2 -1*t4 med : 1*D -1*P' \
-dataTable \
Subj med cueType Time InputFile \
...