bursts {sjemea} | R Documentation |
Burst detection of MEA spike trains.
Description
For a set of spike trains in an MEA recording, find the
bursts independently within each spike train.
Usage
spikes.to.bursts(s, method="si")
Arguments
s |
MEA data structure |
method |
A string, either "si" (surprise method), "mi"
(maxinterval), "logisi" (Log ISI histogram). |
Value
Return the "all bursts" data structure. This is a list of matrices,
giving the burst information for each electrode.
Each matrix stores basic information about each burst. There is one
row for every burst, with the following columns:
beg | index of the first spike in the burst |
len | number of spikes in this burst |
SI | surprise index (calculated only for the surprise method) |
durn | duration (in s) of the burst |
mean.isis | mean of all interspike intervals. |
If no bursts could be found within a spike train, the value NA is used
rather than an empty matrix.
Examples
data.file <- system.file("examples", "TC89_DIV15_A.nexTimestamps",
package = "sjemea")
s <- sanger.read.spikes(data.file)
s$allb <- spikes.to.bursts.surprise(s)
if (interactive()) spikeview(s)
[Package
sjemea version 0.24
Index]