Hello Gang,
I tried to run the deconvolution script I prepared (see below), but I keep having problems calculating the contrasts. If I have spaces in the gltsym lines, it won"t do the contrasts at all, if I take out the spaces, it calculates nonsense and gives the following message:
++ '-stim_times_AM2 1 onsettimes_pleasant_S1.txt' has 1 auxiliary values per time point
++ '-stim_times_AM2 1 onsettimes_pleasant_S1.txt' will have 2 regressors
++ '-stim_times_AM2 2 onsettimes_neutral_S1.txt' has 1 auxiliary values per time point
++ '-stim_times_AM2 2 onsettimes_neutral_S1.txt' will have 2 regressors
++ '-stim_times_AM2 3 onsettimes_unpleasant_S1.txt' has 1 auxiliary values per time point
++ '-stim_times_AM2 3 onsettimes_unpleasant_S1.txt' will have 2 regressors
++ 3dDeconvolve: AFNI version=AFNI_2010_10_19_1028 (Jan 3 2011) [64-bit]
++ Authored by: B. Douglas Ward, et al.
++ current memory malloc-ated = 113,462 bytes (about 113 thousand)
++ loading dataset psc.ssm.tsm.cl.dsp.hpf.reg.ts.1+orig psc.ssm.tsm.cl.dsp.hpf.reg.ts.2+orig
++ current memory malloc-ated = 275,373,774 bytes (about 275 million)
++ Auto-catenated datasets start at: 0 210
*+ WARNING: Input polort=1; Longest run=630.0 s; Recommended minimum polort=5
++ -stim_times using TR=3 s for stimulus timing conversion
++ -stim_times using TR=3 s for any -iresp output datasets
++ [you can alter the -iresp TR via the -TR_times option]
++ -stim_times_AM2 1 using GLOBAL times
++ '-stim_times_AM2 1' average amplitude#1=-0.157895
++ -stim_times_AM2 2 using GLOBAL times
++ '-stim_times_AM2 2' average amplitude#1=-0.588235
++ -stim_times_AM2 3 using GLOBAL times
++ '-stim_times_AM2 3' average amplitude#1=-0.15
++ -stim_times 4 using GLOBAL times
------------------------------------------------------------
GLT matrix from 'SYM:pleasant[0]-neutral[0]':
0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0
------------------------------------------------------------
GLT matrix from 'SYM:pleasant[1]-unpleasant[1]':
0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0
------------------------------------------------------------
GLT matrix from 'SYM:neutral[1]-unpleasant[1]':
0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0
++ Number of time points: 420 (no censoring)
+ Number of parameters: 17 [10 baseline ; 7 signal]
++ Memory required for output bricks = 298,844,160 bytes (about 299 million)
++ Wrote matrix values to file Xmat.x1D
Fatal Signal 11 (SIGSEGV) received
calculate_results
3dDeconvolve main
Bottom of Debug Stack
** AFNI version = AFNI_2010_10_19_1028 Compile date = Jan 3 2011
** [[Precompiled binary macosx_10.5_Intel_64: Jan 3 2011]]
** Program Abort **
** If you report this crash to the AFNI message board,
** please copy the error messages EXACTLY, and give
** the command line you used to run the program, and
** any other information needed to repeat the problem.
** You may later be asked to upload data to help debug.
Can you see what the problem is? The experiment consists of the presentation of pictures of different valence (pleasant, neutral, unpleasant, filler items [others]) for 18 sec. The corresponding valence ratings the participants gave for these pictures (except for others) were used as modulation parameters.
It would be great if you could help me with this problem.
All the best,
Jana
---------------------------------------------------------------------------------
3dDeconvolve\
-input psc.ssm.tsm.cl.dsp.hpf.reg.ts.1+orig psc.ssm.tsm.cl.dsp.hpf.reg.ts.2+orig\
-polort 1\
-num_stimts 10\
-global_times\
-stim_times_AM2 1 onsettimes_pleasant_S1.txt 'BLOCK(18,1)'\
-stim_label 1 pleasant\
-stim_times_AM2 2 onsettimes_neutral_S1.txt 'BLOCK(18,1)'\
-stim_label 2 neutral\
-stim_times_AM2 3 onsettimes_unpleasant_S1.txt 'BLOCK(18,1)'\
-stim_label 3 unpleasant\
-stim_times 4 '1D: 9 27 339 507' 'BLOCK(18,1)'\
-stim_label 4 other\
-stim_file 5 mot[0]\
-stim_base 5\
-stim_label 5 Roll\
-stim_file 6 mot[1]\
-stim_base 6\
-stim_label 6 Pitch\
-stim_file 7 mot[2]\
-stim_base 7\
-stim_label 7 Yaw\
-stim_file 8 mot[3]\
-stim_base 8\
-stim_label 8 dS\
-stim_file 9 mot[4]\
-stim_base 9\
-stim_label 9 dL\
-stim_file 10 mot[5]\
-stim_base 10\
-stim_label 10 dP\
-iresp 1 iresp_modulation/pleasant.irf\
-iresp 2 iresp_modulation/neutral.irf\
-iresp 3 iresp_modulation/unpleasant.irf\
-iresp 4 iresp_modulation/other.irf\
-num_glt 6\
-gltsym 'SYM:pleasant[0]-neutral[0]'\
-glt_label 1 'pleasant_vs_neutral\
-gltsym 'SYM:pleasant[0]-unpleasant[0]'\
-glt_label 2 'pleasant_vs_unpleasant'\
-gltsym 'SYM:neutral[0]-unpleasant[0]'\
-glt_label 3 'neutral_vs_unpleasant'\
-gltsym 'SYM:pleasant[1]-neutral[1]'\
-glt_label 4 'pleasant_vs_neutral_modulation\
-gltsym 'SYM:pleasant[1]-unpleasant[1]'\
-glt_label 5 'pleasant_vs_unpleasant_modulation'\
-gltsym 'SYM:neutral[1]-unpleasant[1]'\
-glt_label 6 'neutral_vs_unpleasant_modulation'\
-fout\
-rout\
-tout\
-fitts fitts_modulation_S1\
-x1D Xmat.x1D\
-bucket decon_modulation_S1\