3dDeconvolve GLT output

AFNI Message Board
1

I’m running 3dDeconvolve (through afni_proc.py) with eight conditions (2x4), the 3dDeconvolve command is attached below.

I’m used to getting three subricks in the stats.subj_REML file for each GLT, namely #0_Coef, #0_Tstat and _Fstat, however this time, for all the GLTs that have all 8 conditions (but not for the other GLTs), I’m getting five output subbricks: #0_Coef, #0_Tstat, #1_Coef, #1_Tstat and _Fstat.
The #0 and #1 subbricks show completely different results from each other.

What are the #1 subbricks, and how should they be interpreted?


3dDeconvolve -input pb04.$subj.r*.scale+tlrc.HEAD                          \
     -censor censor_${subj}_combined_2.1D                                   \
     -polort 2                                                              \
     -local_times                                                           \
     -num_stimts 14                                                         \
     -stim_times 1 stimuli/SL_01_ling_frst.txt 'GAM'                        \
     -stim_label 1 ling_frst                                                \
     -stim_times 2 stimuli/SL_02_nonl_frst.txt 'GAM'                        \
     -stim_label 2 nonl_frst                                                \
     -stim_times 3 stimuli/SL_03_ling_rand.txt 'GAM'                        \
     -stim_label 3 ling_rand                                                \
     -stim_times 4 stimuli/SL_04_ling_pred.txt 'GAM'                        \
     -stim_label 4 ling_pred                                                \
     -stim_times 5 stimuli/SL_05_ling_unpd.txt 'GAM'                        \
     -stim_label 5 ling_unpd                                                \
     -stim_times 6 stimuli/SL_06_nonl_rand.txt 'GAM'                        \
     -stim_label 6 nonl_rand                                                \
     -stim_times 7 stimuli/SL_07_nonl_pred.txt 'GAM'                        \
     -stim_label 7 nonl_pred                                                \
     -stim_times 8 stimuli/SL_08_nonl_unpd.txt 'GAM'                        \
     -stim_label 8 nonl_unpd                                                \
     -stim_file 9 motion_demean.1D'[0]' -stim_base 9 -stim_label 9 roll     \
     -stim_file 10 motion_demean.1D'[1]' -stim_base 10 -stim_label 10 pitch \
     -stim_file 11 motion_demean.1D'[2]' -stim_base 11 -stim_label 11 yaw   \
     -stim_file 12 motion_demean.1D'[3]' -stim_base 12 -stim_label 12 dS    \
     -stim_file 13 motion_demean.1D'[4]' -stim_base 13 -stim_label 13 dL    \
     -stim_file 14 motion_demean.1D'[5]' -stim_base 14 -stim_label 14 dP    \
     -num_glt 20                                                            \
     -gltsym 'SYM: +ling_frst +nonl_frst +ling_rand +ling_pred +ling_unpd   \
     +nonl_rand +nonl_pred +nonl_unpd'                                      \
     -glt_label 1 all_cond                                                  \
     -gltsym 'SYM: +ling_frst +ling_rand +ling_pred +ling_unpd'             \
     -glt_label 2 ling_all                                                  \
     -gltsym 'SYM: +nonl_frst +nonl_rand +nonl_pred +nonl_unpd'             \
     -glt_label 3 nonl_all                                                  \
     -gltsym 'SYM: +ling_rand +nonl_rand'                                   \
     -glt_label 4 rand                                                      \
     -gltsym 'SYM: +ling_pred +ling_unpd +nonl_pred +nonl_unpd'             \
     -glt_label 5 strc                                                      \
     -gltsym 'SYM: +ling_pred +nonl_pred'                         \
     -glt_label 6 pred   									               \
     -gltsym 'SYM: +ling_unpd +nonl_unpd'                         \
     -glt_label 7 unpd                                                      \
     -gltsym 'SYM: +ling_rand +ling_pred +ling_unpd'                        \
     -glt_label 8 ling                                                      \
     -gltsym 'SYM: +nonl_rand +nonl_pred +nonl_unpd'                        \
     -glt_label 9 nonl                                                      \
     -gltsym 'SYM: +ling_pred +ling_unpd -ling_rand +nonl_pred +nonl_unpd   \
     -nonl_rand'                                                            \
     -glt_label 10 strc-rand                                                \
     -gltsym 'SYM: +ling_pred +ling_unpd -ling_rand'                        \
     -glt_label 11 ling_strc-rand                                           \
     -gltsym 'SYM: +ling_pred -ling_unpd'                                   \
     -glt_label 12 ling_pred-unpd                                           \
     -gltsym 'SYM: +ling_pred -ling_unpd -ling_rand'                        \
     -glt_label 13 ling_pred-unpd+rand                                      \
     -gltsym 'SYM: +nonl_pred +nonl_unpd -nonl_rand'                        \
     -glt_label 14 nonl_strc-rand                                           \
     -gltsym 'SYM: +nonl_pred -nonl_unpd'                                   \
     -glt_label 15 nonl_pred-unpd                                           \
     -gltsym 'SYM: +nonl_pred -nonl_unpd -nonl_rand'                        \
     -glt_label 16 nonl_pred-unpd+rand                                      \
     -gltsym 'SYM: +ling_frst -nonl_frst +ling_rand +ling_pred +ling_unpd   \
     -nonl_rand -nonl_pred -nonl_unpd'                                      \
     -glt_label 17 ling-nonl_all                                            \
     -gltsym 'SYM: +ling_rand +ling_pred +ling_unpd -nonl_rand -nonl_pred   \
     -nonl_unpd'                                                            \
     -glt_label 18 ling-nonl                                                \
     -gltsym 'SYM: +ling_pred +ling_unpd -ling_rand -nonl_pred -nonl_unpd   \
     +nonl_rand'                                                            \
     -glt_label 19 stimXstrc                                                \
     -gltsym 'SYM: +ling_pred -ling_unpd -nonl_pred +nonl_unpd'             \
     -glt_label 20 stimXpred                                                \
     -jobs 4                                                                \
     -fout -tout -x1D X.xmat.1D -xjpeg X.jpg                                \
     -x1D_uncensored X.nocensor.xmat.1D                                     \
     -fitts fitts.$subj                                                     \
     -errts errts.${subj}                                                   \
     -bucket stats.$subj
2

Hi Henry,

Would you check to see that perhaps the nonl GLT has 3 output volumes, while strc-rand has 5? Those should be sequential.

Thanks,

  • rick
3

The nonl GLT has three, the strc-rand has 5.

This is the list of subbricks.


Number of values stored at each pixel = 95
  -- At sub-brick #0 'Full_Fstat' datum type is float:            0 to       284.502
     statcode = fift;  statpar = 8 1447
  -- At sub-brick #1 'ling_frst#0_Coef' datum type is float:     -1.81978 to         2.096
  -- At sub-brick #2 'ling_frst#0_Tstat' datum type is float:     -11.3147 to       16.4216
     statcode = fitt;  statpar = 1447
  -- At sub-brick #3 'ling_frst_Fstat' datum type is float:            0 to        269.67
     statcode = fift;  statpar = 1 1447
  -- At sub-brick #4 'nonl_frst#0_Coef' datum type is float:     -2.26268 to       1.79856
  -- At sub-brick #5 'nonl_frst#0_Tstat' datum type is float:     -12.6531 to        13.426
     statcode = fitt;  statpar = 1447
  -- At sub-brick #6 'nonl_frst_Fstat' datum type is float:            0 to       180.257
     statcode = fift;  statpar = 1 1447
  -- At sub-brick #7 'ling_rand#0_Coef' datum type is float:      -1.5234 to       1.45699
  -- At sub-brick #8 'ling_rand#0_Tstat' datum type is float:     -10.8983 to       27.4944
     statcode = fitt;  statpar = 1447
  -- At sub-brick #9 'ling_rand_Fstat' datum type is float:            0 to       755.942
     statcode = fift;  statpar = 1 1447
  -- At sub-brick #10 'ling_pred#0_Coef' datum type is float:     -2.12394 to        2.5055
  -- At sub-brick #11 'ling_pred#0_Tstat' datum type is float:     -8.82936 to       19.6803
     statcode = fitt;  statpar = 1447
  -- At sub-brick #12 'ling_pred_Fstat' datum type is float:            0 to       387.313
     statcode = fift;  statpar = 1 1447
  -- At sub-brick #13 'ling_unpd#0_Coef' datum type is float:     -4.08426 to       2.32614
  -- At sub-brick #14 'ling_unpd#0_Tstat' datum type is float:     -6.63447 to       17.2248
     statcode = fitt;  statpar = 1447
  -- At sub-brick #15 'ling_unpd_Fstat' datum type is float:            0 to       296.694
     statcode = fift;  statpar = 1 1447
  -- At sub-brick #16 'nonl_rand#0_Coef' datum type is float:     -2.13064 to       1.77059
  -- At sub-brick #17 'nonl_rand#0_Tstat' datum type is float:     -13.7465 to       22.7572
     statcode = fitt;  statpar = 1447
  -- At sub-brick #18 'nonl_rand_Fstat' datum type is float:            0 to       517.889
     statcode = fift;  statpar = 1 1447
  -- At sub-brick #19 'nonl_pred#0_Coef' datum type is float:     -1.93336 to       2.14809
  -- At sub-brick #20 'nonl_pred#0_Tstat' datum type is float:     -13.4094 to       16.5902
     statcode = fitt;  statpar = 1447
  -- At sub-brick #21 'nonl_pred_Fstat' datum type is float:            0 to       275.235
     statcode = fift;  statpar = 1 1447
  -- At sub-brick #22 'nonl_unpd#0_Coef' datum type is float:     -4.15211 to       4.70986
  -- At sub-brick #23 'nonl_unpd#0_Tstat' datum type is float:     -7.02976 to        15.377
     statcode = fitt;  statpar = 1447
  -- At sub-brick #24 'nonl_unpd_Fstat' datum type is float:            0 to       236.451
     statcode = fift;  statpar = 1 1447
  -- At sub-brick #25 'all_cond#0_Coef' datum type is float:     -6.95755 to       6.26203
  -- At sub-brick #26 'all_cond#0_Tstat' datum type is float:     -19.7129 to       35.1732
     statcode = fitt;  statpar = 1447
  -- At sub-brick #27 'all_cond#1_Coef' datum type is float:     -3.30265 to       4.93838
  -- At sub-brick #28 'all_cond#1_Tstat' datum type is float:     -16.3772 to       30.3611
     statcode = fitt;  statpar = 1447
  -- At sub-brick #29 'all_cond_Fstat' datum type is float:            0 to       946.723
     statcode = fift;  statpar = 2 1447
  -- At sub-brick #30 'ling_all#0_Coef' datum type is float:     -6.19311 to       5.20055
  -- At sub-brick #31 'ling_all#0_Tstat' datum type is float:     -15.8389 to       35.2692
     statcode = fitt;  statpar = 1447
  -- At sub-brick #32 'ling_all_Fstat' datum type is float:            0 to       1243.92
     statcode = fift;  statpar = 1 1447
  -- At sub-brick #33 'nonl_all#0_Coef' datum type is float:     -4.56772 to       6.05644
  -- At sub-brick #34 'nonl_all#0_Tstat' datum type is float:     -19.2577 to       29.9529
     statcode = fitt;  statpar = 1447
  -- At sub-brick #35 'nonl_all_Fstat' datum type is float:            0 to       897.174
     statcode = fift;  statpar = 1 1447
  -- At sub-brick #36 'rand#0_Coef' datum type is float:     -2.81011 to       2.95263
  -- At sub-brick #37 'rand#0_Tstat' datum type is float:     -17.3873 to       33.3413
     statcode = fitt;  statpar = 1447
  -- At sub-brick #38 'rand_Fstat' datum type is float:            0 to       1111.64
     statcode = fift;  statpar = 1 1447
  -- At sub-brick #39 'strc#0_Coef' datum type is float:     -6.74291 to       5.75183
  -- At sub-brick #40 'strc#0_Tstat' datum type is float:     -15.1926 to       32.1429
     statcode = fitt;  statpar = 1447
  -- At sub-brick #41 'strc_Fstat' datum type is float:            0 to       1033.17
     statcode = fift;  statpar = 1 1447
  -- At sub-brick #42 'pred#0_Coef' datum type is float:      -4.0049 to       3.36988
  -- At sub-brick #43 'pred#0_Tstat' datum type is float:     -15.0482 to       23.4258
     statcode = fitt;  statpar = 1447
  -- At sub-brick #44 'pred_Fstat' datum type is float:            0 to       548.769
     statcode = fift;  statpar = 1 1447
  -- At sub-brick #45 'unpd#0_Coef' datum type is float:     -7.55168 to       6.94643
  -- At sub-brick #46 'unpd#0_Tstat' datum type is float:     -8.93151 to       20.2873
     statcode = fitt;  statpar = 1447
  -- At sub-brick #47 'unpd_Fstat' datum type is float:            0 to       411.575
     statcode = fift;  statpar = 1 1447
  -- At sub-brick #48 'ling#0_Coef' datum type is float:     -5.25008 to       4.04901
  -- At sub-brick #49 'ling#0_Tstat' datum type is float:     -13.6882 to       34.8114
     statcode = fitt;  statpar = 1447
  -- At sub-brick #50 'ling_Fstat' datum type is float:            0 to       1211.83
     statcode = fift;  statpar = 1 1447
  -- At sub-brick #51 'nonl#0_Coef' datum type is float:     -3.30265 to       4.93838
  -- At sub-brick #52 'nonl#0_Tstat' datum type is float:     -16.3772 to       30.3611
     statcode = fitt;  statpar = 1447
  -- At sub-brick #53 'nonl_Fstat' datum type is float:            0 to       921.796
     statcode = fift;  statpar = 1 1447
  -- At sub-brick #54 'strc-rand#0_Coef' datum type is float:     -6.12877 to        4.3347
  -- At sub-brick #55 'strc-rand#0_Tstat' datum type is float:     -12.1645 to       22.5615
     statcode = fitt;  statpar = 1447
  -- At sub-brick #56 'strc-rand#1_Coef' datum type is float:     -1.77059 to       2.13064
  -- At sub-brick #57 'strc-rand#1_Tstat' datum type is float:     -22.7572 to       13.7465
     statcode = fitt;  statpar = 1447
  -- At sub-brick #58 'strc-rand_Fstat' datum type is float:            0 to       494.799
     statcode = fift;  statpar = 2 1447
  -- At sub-brick #59 'ling_strc-rand#0_Coef' datum type is float:     -4.02179 to       2.13153
  -- At sub-brick #60 'ling_strc-rand#0_Tstat' datum type is float:     -5.75157 to        12.578
     statcode = fitt;  statpar = 1447
  -- At sub-brick #61 'ling_strc-rand_Fstat' datum type is float:            0 to       158.206
     statcode = fift;  statpar = 1 1447
  -- At sub-brick #62 'ling_pred-unpd#0_Coef' datum type is float:     -4.39767 to       6.58146
  -- At sub-brick #63 'ling_pred-unpd#0_Tstat' datum type is float:     -5.50436 to       5.04205
     statcode = fitt;  statpar = 1447
  -- At sub-brick #64 'ling_pred-unpd_Fstat' datum type is float:            0 to        30.298
     statcode = fift;  statpar = 1 1447
  -- At sub-brick #65 'ling_pred-unpd+rand#0_Coef' datum type is float:     -4.18745 to       5.94839
  -- At sub-brick #66 'ling_pred-unpd+rand#0_Tstat' datum type is float:     -13.8196 to       5.43438
     statcode = fitt;  statpar = 1447
  -- At sub-brick #67 'ling_pred-unpd+rand_Fstat' datum type is float:            0 to       190.982
     statcode = fift;  statpar = 1 1447
  -- At sub-brick #68 'nonl_strc-rand#0_Coef' datum type is float:     -3.57638 to       3.47854
  -- At sub-brick #69 'nonl_strc-rand#0_Tstat' datum type is float:     -7.52991 to       11.2293
     statcode = fitt;  statpar = 1447
  -- At sub-brick #70 'nonl_strc-rand_Fstat' datum type is float:            0 to       126.098
     statcode = fift;  statpar = 1 1447
  -- At sub-brick #71 'nonl_pred-unpd#0_Coef' datum type is float:     -6.63696 to       5.95475
  -- At sub-brick #72 'nonl_pred-unpd#0_Tstat' datum type is float:     -6.02616 to       5.39157
     statcode = fitt;  statpar = 1447
  -- At sub-brick #73 'nonl_pred-unpd_Fstat' datum type is float:            0 to       36.3146
     statcode = fift;  statpar = 1 1447
  -- At sub-brick #74 'nonl_pred-unpd+rand#0_Coef' datum type is float:     -6.84601 to       5.82035
  -- At sub-brick #75 'nonl_pred-unpd+rand#0_Tstat' datum type is float:     -11.3323 to        5.4652
     statcode = fitt;  statpar = 1447
  -- At sub-brick #76 'nonl_pred-unpd+rand_Fstat' datum type is float:            0 to       128.421
     statcode = fift;  statpar = 1 1447
  -- At sub-brick #77 'ling-nonl_all#0_Coef' datum type is float:     -5.42867 to       4.13906
  -- At sub-brick #78 'ling-nonl_all#0_Tstat' datum type is float:     -8.68324 to       29.4352
     statcode = fitt;  statpar = 1447
  -- At sub-brick #79 'ling-nonl_all#1_Coef' datum type is float:     -4.93838 to       3.30265
  -- At sub-brick #80 'ling-nonl_all#1_Tstat' datum type is float:     -30.3611 to       16.3772
     statcode = fitt;  statpar = 1447
  -- At sub-brick #81 'ling-nonl_all_Fstat' datum type is float:            0 to       772.022
     statcode = fift;  statpar = 2 1447
  -- At sub-brick #82 'ling-nonl#0_Coef' datum type is float:      -6.9395 to       2.50938
  -- At sub-brick #83 'ling-nonl#0_Tstat' datum type is float:     -12.3136 to       20.0322
     statcode = fitt;  statpar = 1447
  -- At sub-brick #84 'ling-nonl#1_Coef' datum type is float:     -4.70986 to       4.15211
  -- At sub-brick #85 'ling-nonl#1_Tstat' datum type is float:      -15.377 to       7.02976
     statcode = fitt;  statpar = 1447
  -- At sub-brick #86 'ling-nonl_Fstat' datum type is float:            0 to        217.73
     statcode = fift;  statpar = 2 1447
  -- At sub-brick #87 'stimXstrc#0_Coef' datum type is float:     -3.70064 to       3.07961
  -- At sub-brick #88 'stimXstrc#0_Tstat' datum type is float:     -13.1413 to       5.84163
     statcode = fitt;  statpar = 1447
  -- At sub-brick #89 'stimXstrc#1_Coef' datum type is float:     -2.13064 to       1.77059
  -- At sub-brick #90 'stimXstrc#1_Tstat' datum type is float:     -13.7465 to       22.7572
     statcode = fitt;  statpar = 1447
  -- At sub-brick #91 'stimXstrc_Fstat' datum type is float:            0 to       325.205
     statcode = fift;  statpar = 2 1447
  -- At sub-brick #92 'stimXpred#0_Coef' datum type is float:     -4.94181 to       4.14377
  -- At sub-brick #93 'stimXpred#0_Tstat' datum type is float:     -3.98362 to       3.95166
     statcode = fitt;  statpar = 1447
  -- At sub-brick #94 'stimXpred_Fstat' datum type is float:            0 to       15.8692
     statcode = fift;  statpar = 1 1447
4

Do you have the output.proc.SUBJECT text file, that contains all of the text output from running the proc script? If you have that, would you mail it to me, along with the proc script itself? Click on my name for my address.

Thanks,

  • rick