design matrix only for stim_times

AFNI experts,

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?

Thank you,

Hi Chirag,

There is no way that I know of to get the header information,
however if you follow what an script does, it at
least gets those columns. Consider this tcsh syntax to get
the ideals of interest and their sum:

set reg_cols = ` -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
  • rick