The next operation would be to actually output smallest FDR q-values for each t-/F-stat volume.
Consider adding -regress_compute_fitts to the AP command, which will save RAM in 3dDeconvolve by computing the fitts dataset via 3dcalc rather than via 3dDeconvolve (which uses the same RAM as for the linear regression).