bursts                package:sjemea                R Documentation

_B_u_r_s_t _d_e_t_e_c_t_i_o_n _o_f _M_E_A _s_p_i_k_e _t_r_a_i_n_s.

_D_e_s_c_r_i_p_t_i_o_n:

     For a set of spike trains in an MEA recording, find the bursts
     independently within each spike train.

_U_s_a_g_e:

     spikes.to.bursts(s, method="si")

_A_r_g_u_m_e_n_t_s:

       s: MEA data structure

  method: A string, either "si" (surprise method), "mi" (maxinterval),
          "logisi" (Log ISI histogram).

_V_a_l_u_e:

     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.

_E_x_a_m_p_l_e_s:

     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)

