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:
[size=x-small]*
*
*
*
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
*
*
*
*[/size]
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’.