Hello Gang,
thank you very much for your reply. I tried your suggestion but it still does not output the R^2 stats. Maybe -rout does not output the R^2 stats for the error in resting state? Below is the script that I used:

afni_proc.py -subj_id ${subj}\
-script proc.py_scripts_remltest/proc_${subj}.sh\
-out_dir ${top_dir}${subj}/${subj}.3dDeconvolve_remltest\
-dsets ${func_dir}pb04.P1002.scale_masked_Despiked_wds.nii.gz\
-blocks regress -regress_reml_exec -regress_opts_reml -rout -regress_motion_file /tmp/mri2/hierarchical/BOLD_variability/data/P1/G1/P1002/P1002.preprocessedtest/dfile_rall.1D -regress_apply_mot_types demean deriv -regress_extra_ortvec ${stim_times_dir}cwm.1D ${stim_times_dir}ccsf.1D -regress_extra_ortvec_labels cwm ccsf

The only bucket datasets that I get are: TSNR.P1002 and stats.P1002_REMLvar. The last one contains the reml variance parameters. Do you have any other suggestion that I can try?

Best
Carolin