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  

|
August 25, 2021 03:23AM
Hi,

To create training/testing dataset for MVPA, I'm doing 3dDeconvolve to extract beta values for each condition for each run.
My '-input' in 3dDeconvolve is a preprocessed file that concatenated 16 runs.
I have onset timing txt files for each condition that consists of 16 rows (i.e. for each run) and N columns (i.e. for each trial).

I put '-concat' to let 3dDeconvolve know the concatenated runs, and I used '-stim_times' with row index'{...}' of a timing txt file (e.g., onset_C_test.txt'{1}') to indicate a specific run's onset timing. For example, onset_C_text.txt looks like this:
*
*
*
*
46.18 74.24 102.33 130.36 158.38 186.43 214.45 242.50
*
*
46.12 74.19 102.19 130.19 158.25 186.26 214.32 242.37
46.16 74.17 102.17 130.23 158.23 186.24 214.28 242.29
*
*
46.14 74.14 102.21 130.23 158.24 186.29 214.36 242.36
*
*
*
*




But afni says errors like this:



** ERROR: mri_read_ascii_ragged_fvect: couldn't open file /MAINTASK/beh_data/onset_C_test.txt{1}
** FATAL ERROR: '-stim_times 17' can't read file '/MAINTASK/beh_data/onset_C_test.txt{1}' [nopt=129]





I'd appreciate if anybody let me know what makes error or any recommendations.
Below is a script that I used.


3dDeconvolve \
-input ${inputEPI} \
-mask ${Bmask} \
-polort A \
-float \
-jobs 2 \
-local_times \
-concat '1D: 0 124 248 372 496 620 744 868 992 1116 1240 1364 1488 1612 1736 1860' \
-num_stimts 86 \
-stim_times 1 ${INITvLaL} 'BLOCK5(36,1)' -stim_label 1 init_vLaL \
-stim_times 2 ${INITvLaR} 'BLOCK5(36,1)' -stim_label 2 init_vLaR \
-stim_times 3 ${INITvLaS} 'BLOCK5(36,1)' -stim_label 3 init_vLaS \
-stim_times 4 ${INITvLaN} 'BLOCK5(36,1)' -stim_label 4 init_vLaN \
-stim_times 5 ${INITvRaL} 'BLOCK5(36,1)' -stim_label 5 init_vRaL \
-stim_times 6 ${INITvRaR} 'BLOCK5(36,1)' -stim_label 6 init_vRaR \
-stim_times 7 ${INITvRaS} 'BLOCK5(36,1)' -stim_label 7 init_vRaS \
-stim_times 8 ${INITvRaN} 'BLOCK5(36,1)' -stim_label 8 init_vRaN \
-stim_times 9 ${TOPUPvLaL} 'BLOCK5(12,1)' -stim_label 9 topup_vLaL \
-stim_times 10 ${TOPUPvLaR} 'BLOCK5(12,1)' -stim_label 10 topup_vLaR \
-stim_times 11 ${TOPUPvLaS} 'BLOCK5(12,1)' -stim_label 11 topup_vLaS \
-stim_times 12 ${TOPUPvLaN} 'BLOCK5(12,1)' -stim_label 12 topup_vLaN \
-stim_times 13 ${TOPUPvRaL} 'BLOCK5(12,1)' -stim_label 13 topup_vRaL \
-stim_times 14 ${TOPUPvRaR} 'BLOCK5(12,1)' -stim_label 14 topup_vRaR \
-stim_times 15 ${TOPUPvRaS} 'BLOCK5(12,1)' -stim_label 15 topup_vRaS \
-stim_times 16 ${TOPUPvRaN} 'BLOCK5(12,1)' -stim_label 16 topup_vRaN \
-stim_times 17 ${maeC}'{1}' 'GAM' -stim_label 17 maeC_R1 \
-stim_times 18 ${maeC}'{2}' 'GAM' -stim_label 18 maeC_R2 \
-stim_times 19 ${maeC}'{3}' 'GAM' -stim_label 19 maeC_R3 \
-stim_times 20 ${maeC}'{4}' 'GAM' -stim_label 20 maeC_R4 \
-stim_times 21 ${maeC}'{5}' 'GAM' -stim_label 21 maeC_R5 \
-stim_times 22 ${maeC}'{6}' 'GAM' -stim_label 22 maeC_R6 \
-stim_times 23 ${maeC}'{7}' 'GAM' -stim_label 23 maeC_R7 \
-stim_times 24 ${maeC}'{8}' 'GAM' -stim_label 24 maeC_R8 \
-stim_times 25 ${maeC}'{9}' 'GAM' -stim_label 25 maeC_R9 \
-stim_times 26 ${maeC}'{10}' 'GAM' -stim_label 26 maeC_R10 \
-stim_times 27 ${maeC}'{11}' 'GAM' -stim_label 27 maeC_R11 \
-stim_times 28 ${maeC}'{12}' 'GAM' -stim_label 28 maeC_R12 \
-stim_times 29 ${maeC}'{13}' 'GAM' -stim_label 29 maeC_R13 \
-stim_times 30 ${maeC}'{14}' 'GAM' -stim_label 30 maeC_R14 \
-stim_times 31 ${maeC}'{15}' 'GAM' -stim_label 31 maeC_R15 \
-stim_times 32 ${maeC}'{16}' 'GAM' -stim_label 32 maeC_R16 \
... [and so on]



* ${maeC} is 'onset_C_test.txt'.



Edited 1 time(s). Last edit at 08/25/2021 03:48AM by Mpark.
Subject Author Posted

How to extract run-wise beta value for each condition

Mpark August 25, 2021 03:23AM

Re: How to extract run-wise beta value for each condition

RWCox August 25, 2021 10:27AM

[SOLVED] Re: How to extract run-wise beta value for each condition

Mpark August 26, 2021 02:09AM

Re: [SOLVED] Re: How to extract run-wise beta value for each condition

RWCox August 26, 2021 11:08AM