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.

Author

b avants stnava@gmail.com

Examples

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