Hello,

I am figure out how to re-preprocess my data in order to extract out stats and nifti files for each block, instead of each run or each condition within a run. For example, in my initial preprocess, my generating script (below) resulting in stats files that I could 3dBucket into to extract nifti files that were an average beta weight for entire condition (unlabeled in the below):

afni_proc.py \
-subj_id $subject.$label.noblur \
-dsets $subject/unlabeled1+orig \
$subject/unlabeled2+orig \
-blocks tshift align tlrc volreg mask scale regress \
-copy_anat $subject/mprage+orig \
-align_opts_aea -giant_move \
-tlrc_base MNI152_T1_2009c+tlrc \
-volreg_align_e2a \
-volreg_tlrc_warp \
-regress_censor_motion 1.0 \
-regress_censor_outliers 0.1 \
-regress_stim_times ONSETS/unlabeled.txt \
-regress_basis_multi 'BLOCK(18,1)' \
-regress_est_blur_epits \
-regress_est_blur_errts

I am trying to be able to extract a nifti file each block, instead. There are 4 blocks in each run and 8 blocks total. Here is the generating script I am working with:

afni_proc.py \
-subj_id $subject.$label.noblur \
-dsets $subject/unlabeled1+orig \
$subject/unlabeled2+orig \
-blocks tshift align tlrc volreg mask scale regress \
-copy_anat $subject/mprage+orig \
-align_opts_aea -giant_move \
-tlrc_base MNI152_T1_2009c+tlrc \
-volreg_align_e2a \
-volreg_tlrc_warp \
-regress_censor_motion 1.0 \
-regress_censor_outliers 0.1 \
-regress_stim_times ONSETS/unlabeled1.txt \
-regress_stim_times ONSETS/unlabeled2.txt \
-regress_stim_times ONSETS/unlabeled3.txt \
-regress_stim_times ONSETS/unlabeled4.txt \
-regress_stim_times ONSETS/unlabeled5.txt \
-regress_stim_times ONSETS/unlabeled6.txt \
-regress_stim_times ONSETS/unlabeled7.txt \
-regress_stim_times ONSETS/unlabeled8.txt \
-regress_run_clustsim
-regress_basis_multi 'BLOCK(18,1)' \
-regress_est_blur_epits \
-regress_est_blur_errts
The main difference is that I am feeding it 8 onset files, each with 7 * and 1 number, corresponding to the onset of that specific block.

I am not getting an error, but I am not confident I didn't miss something.

Any advice on how to meet this goal would be very much appreciated!



Edited 2 time(s). Last edit at 12/06/2017 02:12PM by daisyburr.