12345678910111213141516171819202122232425262728293031 |
- # create random order stim files
- RSFgen -nt 300 -num_stimts 3 \
- -nreps 1 50 -nreps 2 50 -nreps 3 50 \
- -seed 1234568 -prefix RSF.stim.001.
- # from the random order stim files, make ideal reference functions
- waver -GAM -dt 1.0 -input RSF.stim.001.1.1D > wav.hrf.001.1.1D
- waver -GAM -dt 1.0 -input RSF.stim.001.2.1D > wav.hrf.001.2.1D
- waver -GAM -dt 1.0 -input RSF.stim.001.3.1D > wav.hrf.001.3.1D
- # now evaluate the stimulus timings
- 3dDeconvolve \
- -nodata \
- -nfirst 4 \
- -nlast 299 \
- -polort 1 \
- -num_stimts 3 \
- -stim_file 1 "wav.hrf.001.1.1D" \
- -stim_label 1 "stim_A" \
- -stim_file 2 "wav.hrf.001.2.1D" \
- -stim_label 2 "stim_B" \
- -stim_file 3 "wav.hrf.001.3.1D" \
- -stim_label 3 "stim_C" \
- -glt 1 contrasts/contrast_AB \
- -glt 1 contrasts/contrast_AC \
- -glt 1 contrasts/contrast_BC \
- >& 3dD.nodata.001
|