The 3dDeconvolve saves the full model design matrix. I was wondering if there was a way to also save a similar design matrix (with same formatting, including the # headers) that includes only the stim_times timeseries?
There is no way that I know of to get the header information,
however if you follow what an afni_proc.py script does, it at
least gets those columns. Consider this tcsh syntax to get
the ideals of interest and their sum:
set reg_cols = `1d_tool.py -infile X.nocensor.xmat.1D -show_indices_interest`
3dTstat -sum -prefix sum_ideal.1D X.nocensor.xmat.1D"[$reg_cols]"
1dcat X.nocensor.xmat.1D"[$reg_cols]" > X.stim.xmat.1D