Hi,
I am trying to create a plot which shows the time course of percent signal change within an ROI (e.g., amygdala) for 3 different events (viewing pleasant, unpleasant, neutral images). Each image was on the screen for 2.5 s with 12.5s fixation cross btw events. I have 3 runs in total, and in each run the 3 different events were randomly presented. I have concatenated the 3 runs into one file for each subject "sbj_allScaled" which is motion corrected, timeshifted, filtered and scaled.
I have run a voxel-wise analysis (script attached below) and so already have a 1D file for each condition (pleasant, unpleasant, neutral) for each subject which details the onset of the stimulus type within each run. So each subject has 3 1D files (pleasant, unpleasant, neutral), each with 3 rows. I have pasted 1 file below as an example:
---------------
42.5 132.5 147.5 177.5 192.5 222.5 267.5 357.5 417.5 432.5
12.5 57.5 72.5 102.5 147.5 207.5 267.5 342.5 402.5 432.5
57.5 72.5 102.5 177.5 192.5 237.5 297.5 357.5 372.5 417.5
---------------
I also have the mask ready, and I have the timing and order of the events established, but am unsure of the best way to achieve my goal.
Thanks,
Steve
---------------------------------
foreach sub (17593 17600 17675 17769 17781 17824 17850 17858 17879 17883 17969 17990 18058)
3dDeconvolve -input ${sub}_allScaled+orig \
-concat '1D: 0 185 370' \
-nfirst 0 \
-num_stimts 3 \
-stim_times 1 ${sub}_pleasant.1D 'BLOCK(2.5,1)' -stim_label 1 pleasant \
-stim_times 2 ${sub}_unpleasant.1D 'BLOCK(2.5,1)' -stim_label 2 unpleasant \
-stim_times 3 ${sub}_neutral.1D 'BLOCK(2.5,1)' -stim_label 3 neutral \
-polort A \
-glt 1 P-N.txt -glt_label 1 P-N \
-glt 1 U-N.txt -glt_label 2 U-N \
-glt 1 U-P.txt -glt_label 3 U-P \
-fitts fitts_${sub}_glm \
-tout -bucket ${sub}_glm
end