12345678910111213141516171819202122232425262728293031323334353637 |
- #!/bin/tcsh
- # ======================================================================
- # Create an ideal time series for run 1 only, combining Vis and Aud, as
- # if there there were only a single stimulus condition.
- #
- # This would create a ideal time series like epi_r1_ideal.1D.
- #
- # Note that the epi_r1 files include 4s (2 TRs) of pre-steady time.
- # ======================================================================
- # ----------------------------------------------------------------------
- # - combine timing for the 2 conditions (vis and aud)
- # - add 4 seconds (2 TRs)
- # - sort the times, per run
- # - extract just the first run, line 1
- # - save as regressor, AV_r1.stim.txt
- timing_tool.py -timing stim_AV1_vis.txt -extend stim_AV2_aud.txt \
- -add_offset 4 -sort -write_timing - \
- | column_cat -line 1 - > AV_r1.stim.txt
- # ----------------------------------------------------------------------
- # convert timing to a single regressor, AV_r1.ideal.1D
- # (152 time points matches epi_r1, at TR=2s)
- 3dDeconvolve -nodata 152 2 \
- -polort -1 \
- -num_stimts 1 \
- -stim_times 1 AV_r1.stim.txt 'BLOCK(20,1)' -stim_label 1 VA \
- -x1D AV_r1.ideal.1D -x1D_stop
- echo ""
- echo "++ made AV_r1.ideal.1D; consider running:"
- echo ""
- echo " 1dplot -one epi_r1_ideal.1D AV_r1.ideal.1D"
- echo ""
|