Hello,
I use 3dttest to calculate group difference activation maps between two groups of subjects, where each t-test command compares one specific sub-brik of individual subjects' statistical maps, as follows:
# using 3dttest on Rew vs Neu anticipation group differences-
if ( -e COArvnttestglt+tlrc.BRIK ) then
rm -rf COArvnttestglt+tlrc.*
endif
3dttest -session . -prefix COArvnttestglt \
-set1 \
../cr083101/ttestgltreg+tlrc'[46]' \
../la081001/ttestgltreg+tlrc'[46]' \
../sw020802/ttestgltreg+tlrc'[46]' \
../ks082802/ttestgltreg+tlrc'[46]' \
../jg091002/ttestgltreg+tlrc'[46]' \
../jb062802/ttestgltreg+tlrc'[46]' \
../mw082802/ttestgltreg+tlrc'[46]' \
../mrm071202/ttestgltreg+tlrc'[46]' \
../mn051302/ttestgltreg+tlrc'[46]' \
../yt042202/ttestgltreg+tlrc'[46]' \
../vw101901/ttestgltreg+tlrc'[46]' \
../tc040502/ttestgltreg+tlrc'[46]' \
-set2 \
../COA-DBD/ar073101/ttestgltreg+tlrc'[46]' \
../COA-DBD/mb011802/ttestgltreg+tlrc'[46]' \
../COA-DBD/me030802/ttestgltreg+tlrc'[46]' \
../COA-DBD/ra092602/ttestgltreg+tlrc'[46]' \
../COA-DBD/tp111702/ttestgltreg+tlrc'[46]' \
../COA-DBD/md031203/ttestgltreg+tlrc'[46]' \
../COA-DBD/mk051403/ttestgltreg+tlrc'[46]' \
../COA-DBD/ad052803/ttestgltreg+tlrc'[46]' \
../COA-DBD/nw070903/ttestgltreg+tlrc'[46]' \
../COA-DBD/dy082703/ttestgltreg+tlrc'[46]' \
# using 3dttest on Pun vs Neu anticipation group differences-
if ( -e COApvnttestglt+tlrc.BRIK ) then
rm -rf COApvnttestglt+tlrc.*
endif
3dttest -session . -prefix COApvnttestglt -set1 ../cr083101/ttestgltreg+tlrc'[48]' \
../la081001/ttestgltreg+tlrc'[48]' \
../sw020802/ttestgltreg+tlrc'[48]' \
../ks082802/ttestgltreg+tlrc'[48]' \
../jg091002/ttestgltreg+tlrc'[48]' \
../jb062802/ttestgltreg+tlrc'[48]' \
../mw082802/ttestgltreg+tlrc'[48]' \
../mrm071202/ttestgltreg+tlrc'[48]' \
../mn051302/ttestgltreg+tlrc'[48]' \
../yt042202/ttestgltreg+tlrc'[48]' \
../vw101901/ttestgltreg+tlrc'[48]' \
../tc040502/ttestgltreg+tlrc'[48]' \
-set2 \
../COA-DBD/ar073101/ttestgltreg+tlrc'[48]' \
../COA-DBD/mb011802/ttestgltreg+tlrc'[48]' \
../COA-DBD/me030802/ttestgltreg+tlrc'[48]' \
../COA-DBD/ra092602/ttestgltreg+tlrc'[48]' \...
etc., where sub-brick 46 is the "coef" of the event (contrast) of interest which is contained in each subject's "ttestrgltreg" statistical bucket dataset.
Based on a sample script I have been working from, I essentially repeat the same command lines for each of the other sub-bricks to be compared between sets, e.g. 48 and 50. This results in multiple "difference map" files in the directory for AFNI-- eg "rewvsneuttestglt", followed by "punvsneuttestglt" etc.
Considering how all the event-related contrasts in which I am interested in comparing between groups are originally found in the same bucket dataset, is there a way to redo my script so that the individual difference maps for each sub-brik can be generated as separate sub-briks within the same "group difference" bucket dataset? It is cluttering to my files and confusing to keep having to switch functional dataset files in afni viewer. I would prefer to have a single file with multiple group comparison maps that could be accessed with the old pull down menus under "define function."
Perhaps there is some obvious after-the-fact command I could use at the end of my ttests in the script to mesh the individual file maps all together.
Please advise.
Jim B