I come up with two method,

one is using the 3dANOVA with option: -levels 5 , -DAFNI_FLOATIZE=YES and -contr 4 -1 -1 -1 -1, with every subject’s coef value from results of 3dDeconvolve

the other is using 3dttest++ with opiton: -setA, followed by every subject’s data corresponding to option -gltsym A -0.25B -0.25C -0.25D -0.25E of 3dDeconvolve.

the difference betwwen those two method.

I want to know which one is right and how to set the P value to get the cluster(ROI)

the following is part of my code

3dDeconvolve -input pb06.$subj.r01.scale+tlrc.HEAD

pb06.$subj.r02.scale+tlrc.HEAD pb06.$subj.r03.scale+tlrc.HEAD

…

-num_stimts 17

-stim_times 1 stimuli/AV1_A_halfo.txt ‘BLOCK(12,1)’

-stim_label 1 snake

-stim_times 2 stimuli/AV2_B_halfo.txt ‘BLOCK(12,1)’

-stim_label 2 mushroom

-stim_times 3 stimuli/AV3_C_halfo.txt ‘BLOCK(12,1)’

-stim_label 3 spider

-stim_times 4 stimuli/AV4_D_halfo.txt ‘BLOCK(12,1)’

-stim_label 4 animal

-stim_times 5 stimuli/AV5_E_halfo.txt ‘BLOCK(12,1)’

-stim_label 5 face

…

-gltsym ‘SYM: 0.2*A +0.2*B +0.2*C +0.2*D +0.2*E’
-glt_label 12 meanABCDE
-gltsym 'SYM: A -0.25*B -0.25

*C -0.25*D -0.25*E’

-glt_label 13 A-mean(BCDE)

…

-errts errts.${subj}.$cond

-bucket stats.$subj.$cond

tcsh -x stats.REML_cmd

3dANOVA -levels 5 -DAFNI_FLOATIZE=YES

-dset 1 ‘stats.S01.new_REML+tlrc[1]’

-dset 1 ‘stats.S02.new_REML+tlrc[1]’

-dset 1 ‘stats.S03.new_REML+tlrc[1]’

…

dset 5 ‘stats.S08.new_REML+tlrc[13]’

-dset 5 ‘stats.S09.new_REML+tlrc[13]’

-dset 5 ‘stats.S10.new_REML+tlrc[13]’

-ftr stim

-mean 1 A

-mean 2 B

-mean 3 C

-mean 4 D

-mean 5 E

-contr 4 －1 －1 －1 －1 A-mean(BCDE)

-contr 1 1 1 1 1 ABCD \

3dttest++ -prefix A-mean(BCDE)

-setA ‘stats.S01.new_REML+tlrc[52]’

‘stats.S02.new_REML+tlrc[52]’

‘stats.S03.new_REML+tlrc[52]’

‘stats.S04.new_REML+tlrc[52]’

‘stats.S05.new_REML+tlrc[52]’

‘stats.S06.new_REML+tlrc[52]’

‘stats.S07.new_REML+tlrc[52]’

‘stats.S08.new_REML+tlrc[52]’

‘stats.S09.new_REML+tlrc[52]’

‘stats.S10.new_REML+tlrc[52]’

-overwrite