Hi Rick,

Thanks for the reply!

I just fixed the problem. It seems not to be related to the RAM.
In my script, for the output part, I have included the path in the output file name. When I removed those path and only included the file name. It worked.
Maybe 3dANOVA3 cannot tolerant path in the output setting?

The change is as follows:
%%%%%
-adiff 1 2 exper_vs_sham \
-fa stim_protocol \
-fb time \
-fab stim_time_interaction \
-bucket ANOVA

%%%%%

BTW, I have another question about the results. I'm wondering if the output F-statistics are normalized in the volume? In another way, is it z-transformed?
Or because that the number of voxels are too much, so we don't need to convert to z scores? Since I need to use another function to do the post analysis and it requires the input to be z-statistics. Sorry that I may be confused about these statistic terms...

Thanks and happy new year!

Tracy