AFNI wizards,

I was fiddling around with anaticor on multiecho data, and I ran into an issue at the 3dREML stage. Error here:

** ERROR: -dsort dataset './Local_WMe_rall+tlrc.HEAD' doesn't match input dataset 'pb06.501_con_gam_anaticor.r01.scale+tlrc.HEAD pb06.501_con_gam_anaticor.r02.scale+tlrc.HEAD'
 +       nt[dsort]=672  nt[input] = 224
 +       nx[dsort]=55  nx[input] = 55
 +       ny[dsort]=65  ny[input] = 65
 +       nz[dsort]=55  nz[input] = 55
** FATAL ERROR: Can't continue after -dsort mismatch error

The original data is two runs of 112 timepoints, each with 3 echoes. It looks to me that the Local_WMe_rall was created from 2 runs * 3 echoes (672 timepoints), rather than the optimally combined data (224 timepoints).

afni_proc code here:

# ==========================================================================
# script generated by the command:
#
# afni_proc.py -subj_id 501_act_gam_anaticor -blocks despike tshift align     \
#     tlrc volreg mask combine blur scale regress -copy_anat                  \
#     ./../derivatives/anatomical.nii.gz -tcat_remove_first_trs 3             \
#     -tshift_opts_ts -tpattern                                               \
#     @/mnt/data1/dose_response/data/slice_times.txt -dsets_me_echo           \
#     sub-501_ses-main_task-act_run-1_echo-1_bold.nii.gz                      \
#     sub-501_ses-main_task-act_run-2_echo-1_bold.nii.gz -dsets_me_echo       \
#     sub-501_ses-main_task-act_run-1_echo-2_bold.nii.gz                      \
#     sub-501_ses-main_task-act_run-2_echo-2_bold.nii.gz -dsets_me_echo       \
#     sub-501_ses-main_task-act_run-1_echo-3_bold.nii.gz                      \
#     sub-501_ses-main_task-act_run-2_echo-3_bold.nii.gz -echo_times 11 28 45 \
#     -volreg_base_dset sub-501_ses-main_task-act_run-1_echo-1_sbref.nii.gz   \
#     -mask_epi_anat yes -combine_method OC -align_opts_aea -cost lpc+ZZ      \
#     -giant_move -tlrc_base /home/dowdlelt/abin/MNI152_T1_2009c+tlrc         \
#     -tlrc_NL_warp -volreg_tlrc_warp -regress_motion_per_run                 \
#     -regress_stim_times                                                     \
#     /mnt/data1/dose_response/model_params/afni/tms_onsets.txt               \
#     -regress_anaticor_fast -regress_stim_labels tms -regress_basis GAM      \
#     -regress_apply_mot_types demean deriv -regress_est_blur_errts           \
#     -regress_reml_exec -blur_size 10 -bash -scr_overwrite -execute

I think I could hunt through the generated script and correct the portion that extracts the Anaticor regressor (Thanks Bootcamp!) - specfically, this section here:

# --------------------------------------------------
# fast ANATICOR: generate local WMe time series averages
# create catenated volreg dataset
3dTcat -prefix rm.all_runs.volreg pb03.$subj.r*.volreg+tlrc.HEAD

# mask white matter before blurring
3dcalc -a rm.all_runs.volreg+tlrc -b mask_WMe_resam+tlrc                      \
       -expr "a*bool(b)" -datum float -prefix rm.all_runs.volreg.mask

# generate ANATICOR voxelwise regressors via blur
3dmerge -1blur_fwhm 30 -doall -prefix Local_WMe_rall                          \
    rm.all_runs.volreg.mask+tlrc

but I wanted to make sure I wasn't making some obvious mistake in attempting this. Notice any incorrect input? If I run without anaticor, everything works perfeclty (including the reml portion) Thank you for your help.