Hi experts,
I am using afni_proc.py to process some resting-state data. I encountered an issue when runing the script generated by afni_proc.py. The error message says:
** ERROR: *** Censoring has made regression impossible
***
** FATAL ERROR: 3dDeconvolve dies: Insufficient data (61) for estimating 124 parameters
** Program compile date = Apr 16 2021
The 3dDeconvolve part in the script:
# ------------------------------
# run the regression analysis
3dDeconvolve -input pb05.$subj.r*.scale+tlrc.HEAD \
-censor censor_${subj}_combined_2.1D \
-ortvec bandpass_rall.1D bandpass \
-ortvec motion_demean.1D mot_demean \
-ortvec motion_deriv.1D mot_deriv \
-polort 4 -float \
-num_stimts 0 \
-fout -tout -x1D X.xmat.1D -xjpeg X.jpg \
-x1D_uncensored X.nocensor.xmat.1D \
-fitts fitts.$subj \
-errts errts.${subj} \
-x1D_stop \
-bucket stats.$subj
# -- use 3dTproject to project out regression matrix --
# (make errts like 3dDeconvolve, but more quickly)
3dTproject -polort 0 -input pb05.$subj.r*.scale+tlrc.HEAD \
-censor censor_${subj}_combined_2.1D -cenmode ZERO \
-ort X.nocensor.xmat.1D -prefix errts.${subj}.tproject
# if 3dDeconvolve fails, terminate the script
if ( $status != 0 ) then
echo '---------------------------------------'
echo '** 3dDeconvolve error, failing...'
echo ' (consider the file 3dDeconvolve.err)'
exit
endif
The script is generated by bellow codes.
# script generated by the command:
#
# afni_proc.py -subj_id S002 -dsets S002.rs.02aptopup.nii.gz -copy_anat \
# S002.T1.nii -blocks despike tshift align tlrc volreg blur mask scale \
# regress -tcat_remove_first_trs 0 -tlrc_base MNI152_T1_2009c+tlrc \
# -tlrc_NL_warp -volreg_align_to MIN_OUTLIER -volreg_align_e2a \
# -volreg_tlrc_warp -mask_epi_anat yes -blur_size 6 \
# -regress_censor_motion 0.2 -regress_censor_outliers 0.05 \
# -regress_bandpass 0.01 0.1 -regress_apply_mot_types demean deriv \
# -regress_est_blur_epits -regress_est_blur_errts
Could you please help me to have a look? Thank you!
Xiyue