Hi Gang,

It turns out the error messages derived from both -wsE2 option and wrong path for R package.
Thank you so much.

I have one question regarding output.
I expected I could get the same number of t-stats image as the no. of gltLabel (in my case, 12).
However, just one output file was produced after performing 3dMVM script.
Here is the output which was printed in terminal.

++ Smallest FDR q [0 (Intercept) F] = 6.99541e-07
++ Smallest FDR q [1 c01 F] = 7.63635e-08
++ Smallest FDR q [2 c02 F] = 4.77009e-13
++ Smallest FDR q [3 c01:c02 F] = 3.48735e-05
++ Smallest FDR q [5 .... t] = 6.70835e-07
++ Smallest FDR q [7 ... t] = 3.05506e-08
++ Smallest FDR q [9 ... t] = 0.00627071
++ Smallest FDR q [11 _vs_ t] = 1.56036e-08
++ Smallest FDR q [13 _vs_ t] = 1.54795e-06
++ Smallest FDR q [15 _vs_ t] = 8.3346e-06
++ Smallest FDR q [17 _vs_ t] = 0.00029295
++ Smallest FDR q [19 _vs_ t] = 1.13479e-06
++ Smallest FDR q [21 2way_interaction t] = 0.00208255
*+ WARNING: Smallest FDR q [23 2way_interaction t] = 0.832452 ==> few true single voxel detections
++ Smallest FDR q [25 2way_interaction t] = 0.00813095
++ Smallest FDR q [27 2way_interaction t] = 0.00293417
*+ WARNING: varying brick factors, writing NIfTI as float
In summary.Anova.mlm(object$Anova, multivariate = FALSE) :
HF eps > 1 treated as 1

Thank you for your great help.

Jeong