AFNI Message Board

Dear AFNI users-

We are very pleased to announce that the new AFNI Message Board framework is up! Please join us at:

https://discuss.afni.nimh.nih.gov

Existing user accounts have been migrated, so returning users can login by requesting a password reset. New users can create accounts, as well, through a standard account creation process. Please note that these setup emails might initially go to spam folders (esp. for NIH users!), so please check those locations in the beginning.

The current Message Board discussion threads have been migrated to the new framework. The current Message Board will remain visible, but read-only, for a little while.

Sincerely, AFNI HQ

History of AFNI updates  

Previous Message Next Message
Message List New Topic Print View
|
xiyuezh
February 10, 2022 11:49PM
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:

Quote

** 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
Reply Quote
RSS
Subject Author Posted

afni_proc.py error: FATAL ERROR: 3dDeconvolve Insufficient data

xiyuezh February 10, 2022 11:49PM