Hi Gang,
I obtain the same results both with single-subject analysis and with group analysis. Since I have many regressors, I usually run GLM analysis and contrast computation separately. I report below the 3dDeconvolve scripts, so that you can see if there are some mistakes:
Analysis with BLOCK (…,1)
set sub=XX
3dDeconvolve -input {$sub}_all+orig -concat /data/project/concat.txt \
-mask mask_all+orig -polort 2 -basis_normall 1 \
-num_stimts 46 \
-stim_times 1 '1D: 15' 'BLOCK(5.5,1)' -stim_label 1 s1 \
-stim_times 2 '1D: 45' 'BLOCK(5.5,1)' -stim_label 2 s2 \
-stim_times 3 '1D: 75' 'BLOCK(14.4,1)' -stim_label 3 s3 \
-stim_times 4 '1D: 105' 'BLOCK(11.1,1)' -stim_label 4 s4 \
-stim_times 5 '1D: 135' 'BLOCK(5.1,1)' -stim_label 5 s5 \
-stim_times 6 '1D: 165' 'BLOCK(14.8,1)' -stim_label 6 s6 \
-stim_times 7 '1D: 195' 'BLOCK(7.9,1)' -stim_label 7 s7 \
-stim_times 8 '1D: 225' 'BLOCK(4.9,1)' -stim_label 8 s8 \
-stim_times 9 '1D: 285' 'BLOCK(3.7,1)' -stim_label 9 s9 \
..............
-stim_times 38 '1D: 1245' 'BLOCK(2.5,1)' -stim_label 38 s38 \
-stim_times 39 '1D: 1275' 'BLOCK(10.8,1)' -stim_label 39 s39 \
-stim_times 40 '1D: 1305' 'BLOCK(30,1)' -stim_label 40 s40 \
-stim_file 41 {$sub}_mov.out'[1]' -stim_base 41 -stim_file 42 {$sub}_mov.out'[2]' -stim_base 42 \
-stim_file 43 {$sub}_mov.out'[3]' -stim_base 43 -stim_file 44 {$sub}_mov.out'[4]' -stim_base 44 \
-stim_file 45 {$sub}_mov.out'[5]' -stim_base 45 -stim_file 46 {$sub}_mov.out'[6]' -stim_base 46 \
-iresp 1 {$sub}_BLOCK_iresp_s1 \
-iresp 2 {$sub}_BLOCK_iresp_s2 \
-iresp 3 {$sub}_BLOCK_iresp_s3 \
-iresp 4 {$sub}_BLOCK_iresp_s4 \
-iresp 5 {$sub}_BLOCK_iresp_s5 \
-iresp 6 {$sub}_BLOCK_iresp_s6 \
-iresp 7 {$sub}_BLOCK_iresp_s7 \
-iresp 8 {$sub}_BLOCK_iresp_s8 \
-iresp 9 {$sub}_BLOCK_iresp_s9 \
..............
-iresp 38 {$sub}_BLOCK_iresp_s38 \
-iresp 39 {$sub}_BLOCK_iresp_s39 \
-iresp 40 {$sub}_BLOCK_iresp_s40 \
-fout -tout \
-full_first -xsave -nobout\
-xjpeg {$sub}_x_BLOCK.jpeg -x1D {$sub}_x_BLOCK.1D -cbucket {$sub}_cBLOCK \
-fitts {$sub}_fittsBLOCK -bucket {$sub}_resultsBLOCK
3dDeconvolve -input {$sub}_all+orig -concat /data4/mangina/concat.txt -mask mask_all+orig -polort 2 -xrestore {$sub}_resultsBLOCK.xsave \
-num_glt 5 \
-gltsym 'SYM: +s22 +s24 +s28 +s35 +s36 +s37' -glt_label 1 mix \
-gltsym 'SYM: +s1 +s3 +s4 +s5 +s8 +s9 +s10 +s12 +s13 +s14 +s15 +s16' -glt_label 2 dif1 \
-gltsym 'SYM: +s2 +s6 +s7 +s11 +s17 +s19 +s23 +s24 +s31 +s32 +s33 +s35 +s36 +s38' -glt_label 3 dif2 \
-gltsym 'SYM: +s18 +s20 +s21 +s22 +s27 +s28 +s37 +s39' -glt_label 4 dif3 \
-gltsym 'SYM: +s25 +s26 +s29 +s30 +s34 +s40' -glt_label 5 dif4 \
-fout -tout -bucket {$sub}_resultsBLOCK_glt3
Analysis with BLOCK ()
set sub=XX
3dDeconvolve -input {$sub}_all+orig -concat /data/project/concat.txt \
-mask mask_all+orig -polort 2 -basis_normall 1 \
-num_stimts 46 \
-stim_times 1 '1D: 15' 'BLOCK(5.5)' -stim_label 1 s1 \
-stim_times 2 '1D: 45' 'BLOCK(5.5)' -stim_label 2 s2 \
-stim_times 3 '1D: 75' 'BLOCK(14.4)' -stim_label 3 s3 \
-stim_times 4 '1D: 105' 'BLOCK(11.1)' -stim_label 4 s4 \
-stim_times 5 '1D: 135' 'BLOCK(5.1)' -stim_label 5 s5 \
-stim_times 6 '1D: 165' 'BLOCK(14.8)' -stim_label 6 s6 \
-stim_times 7 '1D: 195' 'BLOCK(7.9)' -stim_label 7 s7 \
-stim_times 8 '1D: 225' 'BLOCK(4.9)' -stim_label 8 s8 \
-stim_times 9 '1D: 285' 'BLOCK(3.7)' -stim_label 9 s9 \
.........
-stim_file 41 {$sub}_mov.out'[1]' -stim_base 41 -stim_file 42 {$sub}_mov.out'[2]' -stim_base 42 \
-stim_file 43 {$sub}_mov.out'[3]' -stim_base 43 -stim_file 44 {$sub}_mov.out'[4]' -stim_base 44 \
-stim_file 45 {$sub}_mov.out'[5]' -stim_base 45 -stim_file 46 {$sub}_mov.out'[6]' -stim_base 46 \
-iresp 1 {$sub}_BLOCKw_iresp_s1 \
-iresp 2 {$sub}_BLOCKw_iresp_s2 \
-iresp 3 {$sub}_BLOCKw_iresp_s3 \
-iresp 4 {$sub}_BLOCKw_iresp_s4 \
-iresp 5 {$sub}_BLOCKw_iresp_s5 \
-iresp 6 {$sub}_BLOCKw_iresp_s6 \
-iresp 7 {$sub}_BLOCKw_iresp_s7 \
-iresp 8 {$sub}_BLOCKw_iresp_s8 \
-iresp 9 {$sub}_BLOCKw_iresp_s9 \
.......
-iresp 38 {$sub}_BLOCKw_iresp_s38 \
-iresp 39 {$sub}_BLOCKw_iresp_s39 \
-iresp 40 {$sub}_BLOCKw_iresp_s40 \
-fout -tout \
-full_first -xsave -nobout\
-xjpeg {$sub}_x_BLOCKw.jpeg -x1D {$sub}_x_BLOCKw.1D -cbucket {$sub}_cBLOCKw \
-fitts {$sub}_fittsBLOCKw -bucket {$sub}_resultsBLOCKw
3dDeconvolve -input {$sub}_all+orig -concat /data4/mangina/concat.txt -mask mask_all+orig -polort 2 -xrestore {$sub}_resultsBLOCKw.xsave \
-num_glt 5 \
-gltsym 'SYM: +s22 +s24 +s28 +s35 +s36 +s37' -glt_label 1 mix \
-gltsym 'SYM: +s1 +s3 +s4 +s5 +s8 +s9 +s10 +s12 +s13 +s14 +s15 +s16' -glt_label 2 dif1 \
-gltsym 'SYM: +s2 +s6 +s7 +s11 +s17 +s19 +s23 +s24 +s31 +s32 +s33 +s35 +s36 +s38' -glt_label 3 dif2 \
-gltsym 'SYM: +s18 +s20 +s21 +s22 +s27 +s28 +s37 +s39' -glt_label 4 dif3 \
-gltsym 'SYM: +s25 +s26 +s29 +s30 +s34 +s40' -glt_label 5 dif4 \
-fout -tout -bucket {$sub}_resultsBLOCKw_glt3
The coefficients and statistics contained in
XX_resultsBLOCK and XX_resultsBLOCKw;
XX_BLOCK_iresp_sx and XX_BLOCKw_iresp_sx;
XX_resultsBLOCK_glt3 and XX_resultsBLOCKw_glt3
are identical.
The ANOVA scripts for group analysis are the following:
3dANOVA3 -type 4 -alevels 2 -blevels 2 -clevels 12 \
-dset 1 1 1 AT_resultsBLOCK_glt2beta+tlrc'[0]' \
-dset 2 1 1 AT_resultsBLOCK_glt2beta+tlrc'[1]' \
-dset 1 2 1 AT_resultsBLOCK_glt2beta+tlrc'[2]' \
-dset 2 2 1 AT_resultsBLOCK_glt2beta+tlrc'[3]' \
-dset 1 1 2 AV_resultsBLOCK_glt2beta+tlrc'[0]' \
-dset 2 1 2 AV_resultsBLOCK_glt2beta+tlrc'[1]' \
-dset 1 2 2 AV_resultsBLOCK_glt2beta+tlrc'[2]' \
-dset 2 2 2 AV_resultsBLOCK_glt2beta+tlrc'[3]' \
........................................
-dset 1 1 12 ST_resultsBLOCK_glt2beta+tlrc'[0]' \
-dset 2 1 12 ST_resultsBLOCK_glt2beta+tlrc'[1]' \
-dset 1 2 12 ST_resultsBLOCK_glt2beta+tlrc'[2]' \
-dset 2 2 12 ST_resultsBLOCK_glt2beta+tlrc'[3]' \
-diskspace \
-fa main_ability -fb main_difficulty -fab interact_abdi \
-amean 1 math -amean 2 read \
-bmean 1 easy -bmean 2 diff \
-adiff 1 2 math-read \
-bdiff 1 2 easy-diff \
-abmean 1 1 matheasy -abmean 1 2 mathdiff -abmean 2 1 readeasy -abmean 2 2 readdiff \
-aBdiff 1 2 : 1 math-read_easy -aBdiff 1 2 : 2 math-read_diff \
-Abdiff 1 : 1 2 easy-diff_math -Abdiff 2 : 1 2 easy-diff_read \
-bucket group12_3danova3_2ab2di_glt
3dANOVA3 -type 4 -alevels 2 -blevels 2 -clevels 12 \
-dset 1 1 1 AT_resultsBLOCKw_glt2beta+tlrc'[0]' \
-dset 2 1 1 AT_resultsBLOCKw_glt2beta+tlrc'[1]' \
-dset 1 2 1 AT_resultsBLOCKw_glt2beta+tlrc'[2]' \
-dset 2 2 1 AT_resultsBLOCKw_glt2beta+tlrc'[3]' \
-dset 1 1 2 AV_resultsBLOCKw_glt2beta+tlrc'[0]' \
-dset 2 1 2 AV_resultsBLOCKw_glt2beta+tlrc'[1]' \
-dset 1 2 2 AV_resultsBLOCKw_glt2beta+tlrc'[2]' \
-dset 2 2 2 AV_resultsBLOCKw_glt2beta+tlrc'[3]' \
..........................................
-dset 1 1 12 ST_resultsBLOCKw_glt2beta+tlrc'[0]' \
-dset 2 1 12 ST_resultsBLOCKw_glt2beta+tlrc'[1]' \
-dset 1 2 12 ST_resultsBLOCKw_glt2beta+tlrc'[2]' \
-dset 2 2 12 ST_resultsBLOCKw_glt2beta+tlrc'[3]' \
-diskspace \
-fa main_ability -fb main_difficulty -fab interact_abdi \
-amean 1 math -amean 2 read \
-bmean 1 easy -bmean 2 diff \
-adiff 1 2 math-read \
-bdiff 1 2 easy-diff \
-abmean 1 1 matheasy -abmean 1 2 mathdiff -abmean 2 1 readeasy -abmean 2 2 readdiff \
-aBdiff 1 2 : 1 math-read_easy -aBdiff 1 2 : 2 math-read_diff \
-Abdiff 1 : 1 2 easy-diff_math -Abdiff 2 : 1 2 easy-diff_read \
-bucket group12_3danova3_2ab2di_gltw
Also in this case the coefficients and statistics in group12_3danova3_2ab2di_glt and group12_3danova3_2ab2di_gltw are identical.
Thank you very much for your help,
Silvia
