Hello, I have set up an analysis using afni_proc.py, but I would like to use only simple linear regression, not random effects. This is partly for time and space considerations, but also because I would like to compare my results to results from another package that does not have a random effects method. I understand that I could do ONLY the REMLfit by using the -x1D_stop option, but I want to do the opposite. Also, since the command and matrix for 3DREMLfit will be produced, I can go back and do the RE analysis later (which I may elect to do, but at this time do not want to). How can I do this? Do I just omit "-regress_reml_exec"?

Thank you,
Priya

Here is my afni_proc.py:


afni_proc.py -subj_id ${subj} \
-copy_anat ${subj}_structural+orig \
-dsets 00005+orig 00006+orig 00007+orig 00008+orig 00009+orig 00010+orig\
-blocks despike tshift align tlrc volreg blur mask \
scale regress \
-tcat_remove_first_trs 5 \
-align_opts_aea -cost lpc+ZZ \
-volreg_align_to MIN_OUTLIER \
-volreg_align_e2a \
-volreg_tlrc_warp \
-mask_epi_anat yes \
-regress_stim_times \
${subj}FarFracFrac.txt \
${subj}NearFracFrac.txt \
${subj}MedFracFrac.txt \
${subj}FarFracLine.txt \
${subj}NearFracLine.txt \
${subj}MedFracLine.txt \
${subj}FarLineLine.txt \
${subj}NearLineLine.txt \
${subj}MedLineLine.txt \
-regress_stim_labels farff nearff meDFf \
farfl \
nearfl \
meDFl \
farll \
nearll \
medll \
-regress_basis 'BLOCK(6,1)' \
-regress_motion_per_run \
-regress_censor_motion .3 \
-regress_reml_exec \
-regress_opts_3dD \
-gltsym 'SYM: +nearff +nearfl +nearll -farff -farfl -farll' \
-glt_label 1 all_near_far \
-gltsym 'SYM: +nearff -farff' \
-glt_label 1 FF_near_far \
-gltsym 'SYM: + +nearll -farll' \
-glt_label 1 LL_near_far \
-gltsym 'SYM: +nearfl -farfl' \
-glt_label 1 FL_near_far \
-regress_est_blur_epits \
-regress_est_blur_errts