I tried this flag, but only outputted a set of subbricks (2). I don't know if that represents only the mean for 1 cell only. I would like to obtain 4 sets of subbricks just to compare each cell. This is how I added the flag for xmeans
3dANOVA2 -type 1 \
-mask ${maskdir}/automask_all.nii.gz \
-alevels 2 \
-blevels 2 \
-dset 1 1 ${dir_kids}/01253B/01253B_stats.nii.gz'[22]' \
-dset 1 1 ${dir_kids}/01254B/01254B_stats.nii.gz'[22]' \
......
-dset 2 1 ${dir_kids}/01178B/01178B_stats.nii.gz'[22]' \
-dset 2 1 ${dir_kids}/01299B/01299B_stats.nii.gz'[22]' \
......
-dset 1 2 ${dir_adults}/00753B/00753B_stats.nii.gz'[22]' \
-dset 1 2 ${dir_adults}/00769B/00769B_stats.nii.gz'[22]' \
........
-dset 2 2 ${dir_adults}/00328B/00328B_stats.nii.gz'[22]' \
-dset 2 2 ${dir_adults}/00365B/00365B_stats.nii.gz'[22]' \
.......
-xmean 1 2 means \
-amean 1 TD \
-amean 2 ASD \
-bmean 1 Kids \
-bmean 2 Adults \
-fa Main_effect_A \
-fb Main_effect_B \
-fab Interaction_AxB \
-bucket ${output}/ANOVA_results.nii.gz