Hi Afni experts,
I am using McRetroTS for the first time and I got an error when processing the oba.slibase.1D file in afni_proc.py. I have 160 volumes with a TR=3s and I remove 4 TRs. Also my puls/resp sampling rate is 50 Hz. The scans are interleaved in +z direction and even #slices on a Siemens Trio, therefore alt+z2. Do you know what the error is related to and how to correct it? My thoughts are it is with the generation of oba.slibase.1D based on the message
afni_proc.py -subj_id $sub -script $scrpt_file1 \
-out_dir $out_dir \
-dsets $rest_file \
-copy_anat $anat_file \
-blocks despike ricor tshift volreg regress \
-tcat_remove_first_trs 4 \
-ricor_regs_nfirst 4 \
-ricor_regs $ricor_file \
-tshift_interp -heptic \
-tshift_opts_ts -tpattern alt+z2 \
-volreg_interp -heptic \
-volreg_base_ind 0 159 \
-regress_censor_motion 0.5 \
-regress_censor_outliers 0.1 \
-regress_apply_mot_types demean deriv
My output is the following:
** warning: removing first 4 TRs from beginning of each run
--> the stimulus timing files must reflect the removal of these TRs
-- tcat: reps is now 156
-- ricor: have unscaled short input, will revert back to it
** ERROR: ricor nsliregs x nslices != nvec (49,42,2080)
(# slice 0 labels found = 13)
** script creation failure for block 'ricor'
-- volreg: base/sub indices are 0, 155
++ updating polort to 4, from run len 468.0 s
-- will apply motion types: demean, deriv
++ adding labels: []
-- no 3dClustSim (since no blur estimation)
----------------------------------------------------------------------
** failed command (create_script):
1) I created a truncated .resp and .puls file in matlab to correspond to the timing of the scans in column format
24000 points*(1sec/50pts)*(1vol/3sec) = 160 volumes which matches my acquisition.
Files and output are text/plain; charset=us-ascii
2)I used the following commands to make oba.slicbase.1D after installing McRetroTS:
McRetroTS resp.1D puls.1D 3 160 50 1 Opt.SliceOrder=alt+z2 which created oba.slibase.1D
3)I run afni_proc.py and get an error.
Thanks for the help!
-Ajay
