Well, unforunately, my traditional contrast setup:
3dDeconvolve \
-input normalized+orig \
-concat ../vectors/runs.1D -nfirst 0 -num_stimts 20 \
-mask 3dCL_mask+orig \
-num_glt 21 \
-polort 1 \
-stim_file 1 ../vectors/'cue_noresp_neu.1D' -stim_label 1 'antnr_neu' \
-stim_file 2 ../vectors/'cue_noresp_mys.1D' -stim_label 2 'antnr_mys' \
-stim_file 3 ../vectors/'cue_noresp_gain.1D' -stim_label 3 'antnr_rew' \
-stim_file 4 ../vectors/'cue_resp_neu.1D' -stim_label 4 'antres_neu' \
-stim_file 5 ../vectors/'cue_resp_mys.1D' -stim_label 5 'antres_mys' \
-stim_file 6 ../vectors/'cue_resp_gain.1D' -stim_label 6 'antres_rew' \
-stim_file 7 ../vectors/'fbk_noresp_neu.1D' -stim_label 7 'fbknr_neu' \
-stim_file 8 ../vectors/'fbk_noresp_mys_nogain.1D' -stim_label 8 'fbknr_mnw' \
-stim_file 9 ../vectors/'fbk_noresp_mys_gain.1D' -stim_label 9 'fbknr_mw' \
-stim_file 10 ../vectors/'fbk_noresp_gain.1D' -stim_label 10 'fbknr_rew' \
-stim_file 11 ../vectors/'fbk_resp_neu.1D' -stim_label 11 'fbkres_neu' \
-stim_file 12 ../vectors/'fbk_resp_mys_nogain.1D' -stim_label 12 'fbkres_mnw' \
-stim_file 13 ../vectors/'fbk_resp_mys_gain.1D' -stim_label 13 'fbkres_mw' \
-stim_file 14 ../vectors/'fbk_resp_gain.1D' -stim_label 14 'fbkres_rew' \
-stim_file 15 '3dmotion.1D[1]' -stim_label 15 'roll' \
-stim_file 16 '3dmotion.1D[2]' -stim_label 16 'pitch' \
-stim_file 17 '3dmotion.1D[3]' -stim_label 17 'yaw' \
-stim_file 18 '3dmotion.1D[4]' -stim_label 18 'dS' \
-stim_file 19 '3dmotion.1D[5]' -stim_label 19 'dL' \
-stim_file 20 '3dmotion.1D[6]' -stim_label 20 'dP' \
-tout -bucket ttestreg_cue \
-glt 1 ../vectors/rvn_noresp.mat -glt_label 1 'ant_rvn_nr' \
-glt 1 ../vectors/mvn_noresp.mat -glt_label 2 'ant_mvn_nr' \
-glt 1 ../vectors/rmvn_noresp.mat -glt_label 3 'ant_rmvn_nr' \
-glt 1 ../vectors/mvrn_noresp.mat -glt_label 4 'ant_mvrn_nr' \
-glt 1 ../vectors/rvn_resp.mat -glt_label 5 'ant_rvn_res' \
-glt 1 ../vectors/mvn_resp.mat -glt_label 6 'ant_mvn_res' \
-glt 1 ../vectors/rmvn_resp.mat -glt_label 7 'ant_rmvn_r' \
-glt 1 ../vectors/mvrn_resp.mat -glt_label 8 'ant_mvrn_res' \
-glt 1 ../vectors/mvr_all.mat -glt_label 9 'ant_mvr_all' \
-glt 1 ../vectors/n_rspvnorsp.mat -glt_label 10 'neuant_resvnr' \
-glt 1 ../vectors/m_rspvnorsp.mat -glt_label 11 'mysant_resvnr' \
-glt 1 ../vectors/r_rspvnorsp.mat -glt_label 12 'rewant_resvnr' \
-glt 1 ../vectors/resvnr_rvn.mat -glt_label 13 'HO_resvnr_rvn' \
-glt 1 ../vectors/resvnr_rmvn.mat -glt_label 14 'HO_resvnr_rmvn' \
-glt 1 ../vectors/fbk_mysvneu.mat -glt_label 15 'fbk_mysvneu' \
-glt 1 ../vectors/fbk_rewvneu.mat -glt_label 16 'fbk_rewvneu' \
-glt 1 ../vectors/fbk_mysvrew.mat -glt_label 17 'fbk_mysvrew' \
-glt 1 ../vectors/fbk_winvnowin_resp.mat -glt_label 18 'fbk_winvnw_resp' \
-glt 1 ../vectors/fbk_winvnowin_noresp.mat -glt_label 19 'fbk_winvnw_nr' \
-glt 1 ../vectors/fbk_winvnowin_allmys.mat -glt_label 20 'fbk_winvnw_all' \
-glt 1 ../vectors/fbk_gainvnongain.mat -glt_label 21 'fbk_gainvnongain' \
-iresp 1 antnr_neu \
-iresp 2 antnr_mys \
-iresp 3 antnr_rew \
-iresp 4 antres_neu \
-iresp 5 antres_mys \
-iresp 6 antres_rew \
-iresp 7 fbknr_neu \
-iresp 8 fbknr_mnw \
-iresp 9 fbknr_mw \
-iresp 10 fbknr_rew \
-iresp 11 fbkres_neu \
-iresp 12 fbkres_mnw \
-iresp 13 fbkres_mw \
-iresp 14 fbkres_rew
does NOT give me what I need with -iresp, and instead, I get this error message:
++ WARNING: Only 1 time point for output dataset antnr_neu (etc for all iresps)
And datasets that have only a unitary value in each voxel, and are NOT 3d+time datasets
This error led me to conclude that I had to start specifying a -maxlag in order to actually get plots of a full HRF in the output.
I chose maxlag 7 since I use a 2 s TR, and the 0-8 = 16 s for a full hemodynamic response. Again, the script above is what I normally do to calculate AUC linear contrants, only now I would like to obtain an -iresp dump lasting across the HRF for each of the stim_files.
Jim