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)

