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