Block stimulus model for FMRI Data

Create the block BOLD stimulus for a given task indicator function. Current units in terms of volumes. Sort of a bug, should be in TR.

blockStimulus(scans = 1, onsets = c(1), durations = c(1), rt = 3)

Arguments

scans	number of scans
onsets	vector of onset times (in scans)
durations	vector of duration of ON stimulus in scans or seconds (if !is.null(times))
rt	time between scans in seconds (TR)

Value

Vector with dimension c(scans, 1).

Details

None

References

Worsley, K.J., Liao, C., Aston, J., Petre, V., Duncan, G.H., Morales, F., Evans, A.C. (2002). A general statistical analysis for fMRI data. NeuroImage, 15:1-15.

Examples

  # Example 1
  hrf <- blockStimulus(107, c(18, 48, 78), 15, 2)