spikes.to.bursts: Burst detection of MEA spike trains.

In sje30/sjemea: Multi electrode analysis

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)

sje30/sjemea documentation built on May 21, 2024, 5:44 a.m.