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]