Generate a Raster Plot

Description

Given a list of spike trains (or a repeatedTrain object) where each train was acquired during, say, one presentation of a given stimulus, a raster plot is generated. If stimulus time properties are specified, the stimulus application time also appears on the plot.

Usage

1
2
3
4
5
## S3 method for class 'repeatedTrain'
plot(x, stimTimeCourse = NULL,
       colStim = "grey80", xlim, pch, xlab, ylab, main, ...)
raster(x, stimTimeCourse = NULL, colStim = "grey80",
       xlim, pch, xlab, ylab, main, ...)

Arguments

x

a repeatedTrain object or a list which can be coerced to such an object.

stimTimeCourse

NULL (default) or a two elements vector specifying the time boundaries (in s) of a stimulus presentation.

colStim

the background color used for the stimulus.

xlim

a numeric (default value supplied). See plot.

pch

data symbol used for the spikes. See plot.

xlab

a character (default value supplied). See plot.

ylab

a character (default value supplied). See plot.

main

a character (default value supplied). See plot.

...

see plot.

Details

Basic raster plot stuff.

Value

Nothing is returned raster is used for its side effect, a plot is generated on the current graphical device.

Note

Brillinger (1992) calls these plots "rastor" instead of raster...

Author(s)

Christophe Pouzat christophe.pouzat@gmail.com

References

Brillinger, David R. (1992) Nerve Cell Spike Train Data Analysis: A Progression of Technique. JASA 87: 260–271.

See Also

as.repeatedTrain, is.repeatedTrain, print.repeatedTrain, summary.repeatedTrain, psth

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
## Load Vanillin responses data (first cockroach data set)
data(CAL1V)
## convert them into repeatedTrain objects
## The stimulus command is on between 4.49 s and 4.99s
CAL1V <- lapply(CAL1V,as.repeatedTrain)
## look at the individual raster plots
raster(CAL1V[["neuron 1"]],stimTimeCourse=c(4.49,4.99),main="N1")
plot(CAL1V[["neuron 2"]],stimTimeCourse=c(4.49,4.99),main="N2")
plot(CAL1V[["neuron 3"]],stimTimeCourse=c(4.49,4.99),main="N3")
plot(CAL1V[["neuron 4"]],stimTimeCourse=c(4.49,4.99),main="N4")

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.