AFNI Message Board

Dear AFNI users-

We are very pleased to announce that the new AFNI Message Board framework is up! Please join us at:

https://discuss.afni.nimh.nih.gov

Existing user accounts have been migrated, so returning users can login by requesting a password reset. New users can create accounts, as well, through a standard account creation process. Please note that these setup emails might initially go to spam folders (esp. for NIH users!), so please check those locations in the beginning.

The current Message Board discussion threads have been migrated to the new framework. The current Message Board will remain visible, but read-only, for a little while.

Sincerely, AFNI HQ

History of AFNI updates  

|
July 30, 2014 04:18PM
Hi,
I am a little confused here. We collect 2 separate EPIs for 2 separate events (with two visual conditions, vision and no vision). We collect 2 separate EPIs. I was under the impression that if STIM files are correctly defined, both the EPIs can be run in one afni_poc.py. I did a test:
1. Created afni_proc.py with both the runs and stim files with 2 lines (one for each task) and performed BLOCK(30,1) regression and got beta-coefficients for the tasks

2. Created 2 separate afni_proc.py for each task and modified the stim files (one line only). and performed same regression : model: BLOCK(30,1).

Technically, both methods should give same beta coefficients but they don't. Can you please take a look at the proc.py and let me know what is wrong or what is wrong in my understanding.

AFNI proc.py for 2 EPIs:
afni_proc.py -subj_id $sub \
-script $scrpt_file -scr_overwrite -out_dir $out_dir \
-blocks tshift align tlrc volreg blur mask scale regress \
-copy_anat $anat_file \
-tcat_remove_first_trs $TR \
-tshift_opts_ts -tpattern $pattern \
-dsets \
$same_file \
$diff_file \
-align_opts_aea -cost $costopt \
-tlrc_base MNI_avg152T1+tlrc \
-volreg_align_to last \
-volreg_align_e2a \
-volreg_tlrc_warp \
-blur_to_fwhm \
-blur_size $blur_size \
-blur_in_automask \
-regress_stim_times \
$stim_dir/times_same_vis.1D \
$stim_dir/times_same_nv.1D \
$stim_dir/times_diff_vis.1D \
$stim_dir/times_diff_nv.1D \
-test_stim_files yes \
-regress_stim_labels \
same_vis same_nv diff_vis diff_nv \
-regress_basis_multi \
'BLOCK(30,1)' \
-regress_censor_motion ${motion} \
-regress_opts_3dD \
-tout \
-regress_make_ideal_sum sum_ideal.1D \
-regress_est_blur_epits \
-regress_est_blur_errts
-----------------------------------------------------------

For 1 EPI at a time
afni_proc.py -subj_id $sub \
-script $scrpt_file -scr_overwrite -out_dir $out_dir \
-blocks tshift align tlrc volreg blur mask scale regress \
-copy_anat $anat_file \
-tcat_remove_first_trs $TR \
-tshift_opts_ts -tpattern $pattern \
-dsets \
$same_file \
-align_opts_aea -cost $costopt \
-tlrc_base MNI_avg152T1+tlrc \
-volreg_align_to last \
-volreg_align_e2a \
-volreg_tlrc_warp \
-blur_to_fwhm \
-blur_size $blur_size \
-blur_in_automask \
-regress_stim_times \
$stim_dir/only_same_vis.1D \
$stim_dir/only_same_nv.1D \
-test_stim_files yes \
-regress_stim_labels \
same_vis same_nv \
-regress_basis_multi \
'BLOCK(30,1)' \
-regress_censor_motion ${motion} \
-regress_opts_3dD \
-tout \
-regress_make_ideal_sum sum_ideal.1D \
-regress_est_blur_epits \
-regress_est_blur_errts
---------------------------------------------------------------------------------------------------
Subject Author Posted

2 tasks, 2 EPIs, need 2 separate Proc.py?

pshuk July 30, 2014 04:18PM

Re: 2 tasks, 2 EPIs, need 2 separate Proc.py?

rick reynolds July 31, 2014 11:08AM

Re: 2 tasks, 2 EPIs, need 2 separate Proc.py?

pshuk July 31, 2014 04:39PM

Re: 2 tasks, 2 EPIs, need 2 separate Proc.py?

rick reynolds July 31, 2014 09:04PM