Detection of bursts in extracellular spike trains using hidden semi-Markov point process models.
Neurons in vitro and in vivo have epochs of bursting or "up state" activity during which firing rates are dramatically elevated. Various methods of detecting bursts in extracellular spike trains have appeared in the literature, the most widely used apparently being Poisson Surprise (PS). A natural description of the phenomenon assumes (1) there are two hidden states, which we label "burst" and "non-burst," (2) the neuron evolves stochastically, switching at random between these two states, and (3) within each state the spike train follows a time-homogeneous point process. If in (2) the transitions from non-burst to burst and burst to non-burst states are memoryless, this becomes a hidden Markov model (HMM). For HMMs, the state transitions follow exponential distributions, and are highly irregular. Because observed bursting may in some cases be fairly regular-exhibiting inter-burst intervals with small variation-we relaxed this assumption. When more general probability distributions are used to describe the state transitions the two-state point process model becomes a hidden semi-Markov model (HSMM). We developed an efficient Bayesian computational scheme to fit HSMMs to spike train data. Numerical simulations indicate the method can perform well, sometimes yielding very different results than those based on PS.
Tokdar
Surya
S
Department of Statistical Science, Duke University, Durham, NC, 27708, USA. st118@stat.duke.edu.
Xi
Peiyi
P
Department of Statistics, Carnegie Mellon University, Pittsburgh, PA, 15213, USA.
Kelly
Ryan C
RC
Department of Computer Science, Center for the Neural Basis of Cognition, Carnegie Mellon University, Pittsburgh, PA, 15213, USA.
Kass
Robert E
RE
Department of Statistics, Machine Learning Department, and Center for the Neural Basis of Cognition, Carnegie Mellon University, Pittsburgh, PA, 15213, USA.
D000200:Action Potentials / Q000502:physiology*; D000818:Animals; D003198:Computer Simulation; D008390:Markov Chains; D008959:Models, Neurological; D009474:Neurons / Q000502:physiology*; D012165:Retinal Ganglion Cells / Q000502:physiology; D012680:Sensitivity and Specificity; D012815:Signal Processing, Computer-Assisted; D013997:Time Factors
Neural Comput. 2002 Feb;14(2):325-46
11802915
Exp Brain Res. 1991;83(2):317-28
2022242
Trends Neurosci. 2003 Mar;26(3):161-7
12591219
Biol Cybern. 1991;65(6):459-67
1958731
Biophys J. 1972 Jan;12(1):1-24
4332108
PLoS Biol. 2005 Jun;3(6):e175
15857153
Neurosci Biobehav Rev. 2002 Mar;26(2):127-85
11856557
J Neurosci Methods. 2005 Jun 30;145(1-2):23-35
15922023
Proc Natl Acad Sci U S A. 1995 Sep 12;92(19):8616-20
7567985
J Neurosci Methods. 2007 Mar 15;160(2):349-58
17070926
Nature. 2003 Jan 30;421(6922):539-43
12556894
Biol Cybern. 1997 Oct;77(4):289-95
9394447
Brain Res. 1997 Mar 7;750(1-2):129-40
9098537
J Neurophysiol. 1985 Apr;53(4):926-39
3998798
Neural Comput. 2009 Jul;21(7):1797-862
19323637