raster: Generate a Raster Plot

In STAR: Spike Train Analysis with R

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

## 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 [e-mail]

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

## 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")

Example output

Loading required package: survival
Loading required package: mgcv
Loading required package: nlme
This is mgcv 1.8-17. For overview type 'help("mgcv-package")'.
Loading required package: R2HTML
Loading required package: gss
Loading required package: codetools

STAR documentation built on May 2, 2019, 11:44 a.m.