Hi everybody,
I'm using 3dDeconvolve to clean up some EPI data for resting state analyses. I'm using example 9 from the output of afni_proc.py as my guide and am running 3dDeconvolve indicated below. The problem I'm having is that my errts timeseries is entirely zero. Can any one see any problem with my command that might give rise to this issue?
subjectNumber=330_A
prefixPrefix=${subjectNumber}.restingstate
nvals=$( 3dnvals ${prefixPrefix}.preprocessed.MNI+tlrc.HEAD ) ## 255
tr=$( 3dinfo -tr ${prefixPrefix}.preprocessed.MNI+tlrc.HEAD ) ## 2 sec
echo "*** Found $nvals with a TR of ${tr} sec in the preprocessed dataset"
echo "*** Generating bandpass regressors"
# create bandpass regressors (instead of using 3dBandpass, say)
1dBport -nodata $nvals $tr -band 0.01 0.1 -invert -nozero > bandpass_rall.1D
JOBS="-jobs 4"
3dDeconvolve $JOBS \
-polort $polort \
-float \
-fout -tout \
-censor motion_${subjectNumber}_censor.1D \
-ortvec bandpass_rall.1D bandpass \
-input ${prefixPrefix}.preprocessed.MNI+tlrc.HEAD \
-bucket ${prefixPrefix}.dec \
-errts ${prefixPrefix}.errts \
-fitts ${prefixPrefix}.fitts \
-mask ${prefixPrefix}.preprocessed.mask.MNI+tlrc.HEAD \
-xjpeg ${prefixPrefix}.matrix \
-x1D X.xmat.1D \
-x1D_uncensored X.nocensor.xmat.1D \
-num_stimts 12 \
-stim_file 1 motion_demean.1D'[0]' -stim_base 1 -stim_label 1 roll_01 \
-stim_file 2 motion_demean.1D'[1]' -stim_base 2 -stim_label 2 pitch_01 \
-stim_file 3 motion_demean.1D'[2]' -stim_base 3 -stim_label 3 yaw_01 \
-stim_file 4 motion_demean.1D'[3]' -stim_base 4 -stim_label 4 dS_01 \
-stim_file 5 motion_demean.1D'[4]' -stim_base 5 -stim_label 5 dL_01 \
-stim_file 6 motion_demean.1D'[5]' -stim_base 6 -stim_label 6 dP_01 \
-stim_file 7 motion_deriv.1D'[0]' -stim_base 7 -stim_label 7 roll_02 \
-stim_file 8 motion_deriv.1D'[1]' -stim_base 8 -stim_label 8 pitch_02 \
-stim_file 9 motion_deriv.1D'[2]' -stim_base 9 -stim_label 9 yaw_02 \
-stim_file 10 motion_deriv.1D'[3]' -stim_base 10 -stim_label 10 dS_02 \
-stim_file 11 motion_deriv.1D'[4]' -stim_base 11 -stim_label 11 dL_02 \
-stim_file 12 motion_deriv.1D'[5]' -stim_base 12 -stim_label 12 dP_02