Hi Gang,
Thank you, that makes sense. I have run the following 3dDeconvolve command to generate the coefficients and have pasted the output of one of the conditions. There are 24 sub-bricks associated with one condition and I was wondering which one represents the "intercept" that corresponds to the condition effect when the modulatory variable is adjusted at its center value? In addition, within the 3dDeconvolve command, does the gltlsym flag automatically use sub-brick 0 to create the glt? As an example, I am currently using coefficients generated from glt_label 9 as input to the population level and assumed it was the condition effect when the modulatory variable is adjusted at its center value. Is this correct?
-- At sub-brick #50 'hour_targets#0_Coef' datum type is float: -10.4489 to 12.8519
-- At sub-brick #51 'hour_targets#0_Tstat' datum type is float: -7.96445 to 10.5607
statcode = fitt; statpar = 1374
-- At sub-brick #52 'hour_targets#1_Coef' datum type is float: -14.9663 to 11.7733
-- At sub-brick #53 'hour_targets#1_Tstat' datum type is float: -11.5405 to 14.3313
statcode = fitt; statpar = 1374
-- At sub-brick #54 'hour_targets#2_Coef' datum type is float: -11.499 to 13.5667
-- At sub-brick #55 'hour_targets#2_Tstat' datum type is float: -9.02144 to 24.0982
statcode = fitt; statpar = 1374
-- At sub-brick #56 'hour_targets#3_Coef' datum type is float: -17.0224 to 10.1095
-- At sub-brick #57 'hour_targets#3_Tstat' datum type is float: -10.5102 to 26.6022
statcode = fitt; statpar = 1374
-- At sub-brick #58 'hour_targets#4_Coef' datum type is float: -21.7766 to 13.1005
-- At sub-brick #59 'hour_targets#4_Tstat' datum type is float: -10.0235 to 18.7194
statcode = fitt; statpar = 1374
-- At sub-brick #60 'hour_targets#5_Coef' datum type is float: -14.5087 to 12.9383
-- At sub-brick #61 'hour_targets#5_Tstat' datum type is float: -7.91893 to 9.23477
statcode = fitt; statpar = 1374
-- At sub-brick #62 'hour_targets#6_Coef' datum type is float: -8.64132 to 9.26615
-- At sub-brick #63 'hour_targets#6_Tstat' datum type is float: -5.35049 to 5.07596
statcode = fitt; statpar = 1374
-- At sub-brick #64 'hour_targets#7_Coef' datum type is float: -6.8888 to 7.93592
-- At sub-brick #65 'hour_targets#7_Tstat' datum type is float: -4.55072 to 4.17768
statcode = fitt; statpar = 1374
-- At sub-brick #66 'hour_targets#8_Coef' datum type is float: -0.026589 to 0.0286437
-- At sub-brick #67 'hour_targets#8_Tstat' datum type is float: -4.11187 to 5.41255
statcode = fitt; statpar = 1374
-- At sub-brick #68 'hour_targets#9_Coef' datum type is float: -0.0281863 to 0.0315514
-- At sub-brick #69 'hour_targets#9_Tstat' datum type is float: -3.79148 to 4.3071
statcode = fitt; statpar = 1374
-- At sub-brick #70 'hour_targets#10_Coef' datum type is float: -0.0288239 to 0.0401153
-- At sub-brick #71 'hour_targets#10_Tstat' datum type is float: -4.09525 to 4.27468
statcode = fitt; statpar = 1374
-- At sub-brick #72 'hour_targets#11_Coef' datum type is float: -0.0301019 to 0.0427169
-- At sub-brick #73 'hour_targets#11_Tstat' datum type is float: -3.97246 to 5.50767
statcode = fitt; statpar = 1374
-- At sub-brick #74 'hour_targets#12_Coef' datum type is float: -0.0275101 to 0.0282634
-- At sub-brick #75 'hour_targets#12_Tstat' datum type is float: -4.53677 to 5.98939
statcode = fitt; statpar = 1374
-- At sub-brick #76 'hour_targets#13_Coef' datum type is float: -0.0256773 to 0.0365262
-- At sub-brick #77 'hour_targets#13_Tstat' datum type is float: -3.90988 to 4.54588
statcode = fitt; statpar = 1374
-- At sub-brick #78 'hour_targets#14_Coef' datum type is float: -0.0301191 to 0.0580204
-- At sub-brick #79 'hour_targets#14_Tstat' datum type is float: -4.39122 to 4.33509
statcode = fitt; statpar = 1374
-- At sub-brick #80 'hour_targets#15_Coef' datum type is float: -0.0282869 to 0.0264222
-- At sub-brick #81 'hour_targets#15_Tstat' datum type is float: -3.80869 to 4.21798
statcode = fitt; statpar = 1374
-- At sub-brick #82 'hour_targets#16_Coef' datum type is float: -13.7903 to 10.0837
-- At sub-brick #83 'hour_targets#16_Tstat' datum type is float: -4.12466 to 4.50353
statcode = fitt; statpar = 1374
-- At sub-brick #84 'hour_targets#17_Coef' datum type is float: -10.6915 to 10.858
-- At sub-brick #85 'hour_targets#17_Tstat' datum type is float: -4.01089 to 4.65436
statcode = fitt; statpar = 1374
-- At sub-brick #86 'hour_targets#18_Coef' datum type is float: -15.2205 to 10.6691
-- At sub-brick #87 'hour_targets#18_Tstat' datum type is float: -4.0562 to 4.32744
statcode = fitt; statpar = 1374
-- At sub-brick #88 'hour_targets#19_Coef' datum type is float: -11.9866 to 11.6637
-- At sub-brick #89 'hour_targets#19_Tstat' datum type is float: -3.89083 to 4.91938
statcode = fitt; statpar = 1374
-- At sub-brick #90 'hour_targets#20_Coef' datum type is float: -10.585 to 15.5468
-- At sub-brick #91 'hour_targets#20_Tstat' datum type is float: -3.97106 to 4.45029
statcode = fitt; statpar = 1374
-- At sub-brick #92 'hour_targets#21_Coef' datum type is float: -12.8017 to 10.3374
-- At sub-brick #93 'hour_targets#21_Tstat' datum type is float: -3.96983 to 4.52196
statcode = fitt; statpar = 1374
-- At sub-brick #94 'hour_targets#22_Coef' datum type is float: -7.64403 to 9.91451
-- At sub-brick #95 'hour_targets#22_Tstat' datum type is float: -3.68258 to 3.85035
statcode = fitt; statpar = 1374
-- At sub-brick #96 'hour_targets#23_Coef' datum type is float: -11.5217 to 11.9084
-- At sub-brick #97 'hour_targets#23_Tstat' datum type is float: -4.00909 to 3.8033
3dDeconvolve -input pb01.$subject.r01.scale+orig \
pb01.$subject.r02.scale+orig \
pb01.$subject.r03.scale+orig \
pb01.$subject.r04.scale+orig \
pb01.$subject.r05.scale+orig \
pb01.$subject.r06.scale+orig \
-polort A \
-CENSORTR *:0-2 \
-x1D_stop \
-num_stimts 17 \
-mask /usr/local/mridata/Consolidation_2019/atlases/Kirby/dil+orig \
-stim_file 1 /usr/local/mridata/Consolidation_2019/${subject}/${subject}.nmi/dfile_rall.1D'[0]' -stim_base 1 -stim_label 1 roll \
-stim_file 2 /usr/local/mridata/Consolidation_2019/${subject}/${subject}.nmi/dfile_rall.1D'[1]' -stim_base 2 -stim_label 2 pitch \
-stim_file 3 /usr/local/mridata/Consolidation_2019/${subject}/${subject}.nmi/dfile_rall.1D'[2]' -stim_base 3 -stim_label 3 yaw \
-stim_file 4 /usr/local/mridata/Consolidation_2019/${subject}/${subject}.nmi/dfile_rall.1D'[3]' -stim_base 4 -stim_label 4 dS \
-stim_file 5 /usr/local/mridata/Consolidation_2019/${subject}/${subject}.nmi/dfile_rall.1D'[4]' -stim_base 5 -stim_label 5 dL \
-stim_file 6 /usr/local/mridata/Consolidation_2019/${subject}/${subject}.nmi/dfile_rall.1D'[5]' -stim_base 6 -stim_label 6 dP \
-stim_file 7 /usr/local/mridata/Consolidation_2019/${subject}/${subject}.nmi/motion.deriv.1D'[0]' -stim_base 7 -stim_label 7 rolldx \
-stim_file 8 /usr/local/mridata/Consolidation_2019/${subject}/${subject}.nmi/motion.deriv.1D'[1]' -stim_base 8 -stim_label 8 pitchdx \
-stim_file 9 /usr/local/mridata/Consolidation_2019/${subject}/${subject}.nmi/motion.deriv.1D'[2]' -stim_base 9 -stim_label 9 yawdx \
-stim_file 10 /usr/local/mridata/Consolidation_2019/${subject}/${subject}.nmi/motion.deriv.1D'[3]' -stim_base 10 -stim_label 10 dSdx \
-stim_file 11 /usr/local/mridata/Consolidation_2019/${subject}/${subject}.nmi/motion.deriv.1D'[4]' -stim_base 11 -stim_label 11 dLdx \
-stim_file 12 /usr/local/mridata/Consolidation_2019/${subject}/${subject}.nmi/motion.deriv.1D'[5]' -stim_base 12 -stim_label 12 dPdx \
-stim_times_AM2 13 ../timing_files/'1dMarry_foil_demean_rt_strength_AM2.1D' 'TENT(0,14,8)' -stim_label 13 all_foils \
-stim_times_AM2 14 ../timing_files/'1dMarry_hour_demean_rt_strength_AM2.1D' 'TENT(0,14,8)' -stim_label 14 hour_targets \
-stim_times_AM2 15 ../timing_files/'1dMarry_day_demean_rt_strength_AM2.1D' 'TENT(0,14,8)' -stim_label 15 day_targets \
-stim_times_AM2 16 ../timing_files/'1dMarry_week_demean_rt_strength_AM2.1D' 'TENT(0,14,8)' -stim_label 16 week_targets \
-stim_times_AM2 17 ../timing_files/'1dMarry_month_demean_rt_strength_AM2.1D' 'TENT(0,14,8)' -stim_label 17 month_targets \
-num_glt 18 \
-gltsym 'SYM: +.25*hour_targets +.25*day_targets +.25*week_targets +.25*month_targets' \
-glt_label 1 "all_targets" \
-gltsym 'SYM: 0.5*hour_targets -0.167*month_targets -0.167*day_targets -0.167*week_targets' \
-glt_label 2 "recent-all_remote" \
-gltsym 'SYM: +.25*hour_targets +.25*day_targets +.25*week_targets +.25*month_targets -all_foils' \
-glt_label 3 "targ-foils" \
-gltsym 'SYM: +hour_targets' \
-glt_label 4 "hour_targets" \
-gltsym 'SYM: +day_targets' \
-glt_label 5 "day_targets" \
-gltsym 'SYM: +week_targets' \
-glt_label 6 "week_targets" \
-gltsym 'SYM: +month_targets' \
-glt_label 7 "month_targets" \
-gltsym 'SYM: +all_foils' \
-glt_label 8 "all_foils" \
-gltsym 'SYM: +hour_targets[3..4]' \
-glt_label 9 "peak_hour_targets" \
-gltsym 'SYM: +day_targets[3..4]' \
-glt_label 10 "peak_day_targets" \
-gltsym 'SYM: +week_targets[3..4]' \
-glt_label 11 "peak_week_targets" \
-gltsym 'SYM: +month_targets[3..4]' \
-glt_label 12 "peak_month_targets" \
-gltsym 'SYM: +all_foils[3..4]' \
-glt_label 13 "peak_foils" \
-gltsym 'SYM: +.25*hour_targets[3..4] +.25*day_targets[3..4] +.25*week_targets[3..4] +.25*month_targets[3..4]' \
-glt_label 14 "peak_all_targets" \
-gltsym 'SYM: +0.5*hour_targets[3..4] -.167*day_targets[3..4] -.167*week_targets[3..4] -.167*month_targets[3..4] ' \
-glt_label 15 "peak_recent-all_remote" \
-gltsym 'SYM: +.25*hour_targets[3..4] +.25*day_targets[3..4] +.25*week_targets[3..4] +.25*month_targets[3..4] -all_foils[3..4]' \
-glt_label 16 "peak_targ-foils" \
-gltsym 'SYM: +hour_targets[3..4] -day_targets[3..4] \ +hour_targets[3..4] -week_targets[3..4] \ +hour_targets[3..4] -month_targets[3..4]' \
-glt_label 17 "time_ANOVA_peak" \
-gltsym 'SYM: +hour_targets -day_targets \ +hour_targets -week_targets \ +hour_targets -month_targets' \
-glt_label 18 "time_ANOVA" \
-jobs 16 \
-errts ../AM2_rt_conf_tr2/errts \
-fout -tout -x1D ../AM2_rt_conf_tr2/X.xmat.1D -xjpeg ../AM2_rt_conf_tr2/X.jpg \
-bucket ../AM2_rt_conf_tr2/${subject}_decon_results_all_targets_AM \
-cbucket ../AM2_rt_conf_tr2/${subject}_decon_results_all_targets_cbucket_AM \
Thank you.
Edited 1 time(s). Last edit at 12/22/2020 07:35PM by Catherine Tallman.
