Estimating a State-Space Model from Point Process Observations

2003 ◽  
Vol 15 (5) ◽  
pp. 965-991 ◽  
Author(s):  
Anne C. Smith ◽  
Emery N. Brown

A widely used signal processing paradigm is the state-space model. The state-space model is defined by two equations: an observation equation that describes how the hidden state or latent process is observed and a state equation that defines the evolution of the process through time. Inspired by neurophysiology experiments in which neural spiking activity is induced by an implicit (latent) stimulus, we develop an algorithm to estimate a state-space model observed through point process measurements. We represent the latent process modulating the neural spiking activity as a gaussian autoregressive model driven by an external stimulus. Given the latent process, neural spiking activity is characterized as a general point process defined by its conditional intensity function. We develop an approximate expectation-maximization (EM) algorithm to estimate the unobservable state-space process, its parameters, and the parameters of the point process. The EM algorithm combines a point process recursive nonlinear filter algorithm, the fixed interval smoothing algorithm, and the state-space covariance algorithm to compute the complete data log likelihood efficiently. We use a Kolmogorov-Smirnov test based on the time-rescaling theorem to evaluate agreement between the model and point process data. We illustrate the model with two simulated data examples: an ensemble of Poisson neurons driven by a common stimulus and a single neuron whose conditional intensity function is approximated as a local Bernoulli process.

2015 ◽  
Vol 39 (3-4) ◽  
pp. 1400-1408 ◽  
Author(s):  
Qing Qin ◽  
Jiang Wang ◽  
Ming Xue ◽  
Bin Deng ◽  
Xile Wei

2005 ◽  
Vol 93 (2) ◽  
pp. 1074-1089 ◽  
Author(s):  
Wilson Truccolo ◽  
Uri T. Eden ◽  
Matthew R. Fellows ◽  
John P. Donoghue ◽  
Emery N. Brown

Multiple factors simultaneously affect the spiking activity of individual neurons. Determining the effects and relative importance of these factors is a challenging problem in neurophysiology. We propose a statistical framework based on the point process likelihood function to relate a neuron's spiking probability to three typical covariates: the neuron's own spiking history, concurrent ensemble activity, and extrinsic covariates such as stimuli or behavior. The framework uses parametric models of the conditional intensity function to define a neuron's spiking probability in terms of the covariates. The discrete time likelihood function for point processes is used to carry out model fitting and model analysis. We show that, by modeling the logarithm of the conditional intensity function as a linear combination of functions of the covariates, the discrete time point process likelihood function is readily analyzed in the generalized linear model (GLM) framework. We illustrate our approach for both GLM and non-GLM likelihood functions using simulated data and multivariate single-unit activity data simultaneously recorded from the motor cortex of a monkey performing a visuomotor pursuit-tracking task. The point process framework provides a flexible, computationally efficient approach for maximum likelihood estimation, goodness-of-fit assessment, residual analysis, model selection, and neural decoding. The framework thus allows for the formulation and analysis of point process models of neural spiking activity that readily capture the simultaneous effects of multiple covariates and enables the assessment of their relative importance.


2020 ◽  
Vol 0 (0) ◽  
Author(s):  
Ji Chol ◽  
Ri Jun Il

Abstract The modeling of counter-current leaching plant (CCLP) in Koryo Extract Production is presented in this paper. Koryo medicine is a natural physic to be used for a diet and the medical care. The counter-current leaching method is mainly used for producing Koryo medicine. The purpose of the modeling in the previous works is to indicate the concentration distributions, and not to describe the model for the process control. In literature, there are no nearly the papers for modeling CCLP and especially not the presence of papers that have described the issue for extracting the effective components from the Koryo medicinal materials. First, this paper presents that CCLP can be shown like the equivalent process consisting of two tanks, where there is a shaking apparatus, respectively. It allows leachate to flow between two tanks. Then, this paper presents the principle model for CCLP and the state space model on based it. The accuracy of the model has been verified from experiments made at CCLP in the Koryo Extract Production at the Gang Gyi Koryo Manufacture Factory.


1994 ◽  
Vol 20 (2) ◽  
pp. 143-148 ◽  
Author(s):  
Siddhartha Chib ◽  
Ram C. Tiwari

2010 ◽  
Vol 40-41 ◽  
pp. 27-33 ◽  
Author(s):  
Yi Hui Lin ◽  
Hai Bo Zhang

The method of state space model fitting is carried out by using the linear relation of the variable of the differential equations and separating the steady process and instant process to eliminate the steady errors course by instant errors. The improved fitting method is without solving the linear differential equations or using any iterative methods. The coefficient of the state space model can be solve simply using matrix operation under the premise of high accuracy, so it has a higher computational efficiency than former least square method. And this method can also be used with other fitting method. Finally, to illustrate the validity and accuracy of the improved method, a small perturbation state space model of a certain turboshaft engine model has been established by this method, and the simulation result between state space model and nonlinear model was also compared. Also, the state space model could be applied to fault diagnosis and control system design for aeroengines.


2009 ◽  
Vol 10 (2) ◽  
pp. 117-138 ◽  
Author(s):  
Wai-Yuan Tan ◽  
Weiming Ke ◽  
G. Webb

We develop a state space model documenting Gompertz behaviour of tumour growth. The state space model consists of two sub-models: a stochastic system model that is an extension of the deterministic model proposed by Gyllenberg and Webb (1991), and an observation model that is a statistical model based on data for the total number of tumour cells over time. In the stochastic system model we derive through stochastic equations the probability distributions of the numbers of different types of tumour cells. Combining with the statistic model, we use these distribution results to develop a generalized Bayesian method and a Gibbs sampling procedure to estimate the unknown parameters and to predict the state variables (number of tumour cells). We apply these models and methods to real data and to computer simulated data to illustrate the usefulness of the models, the methods, and the procedures.


Mathematics ◽  
2021 ◽  
Vol 9 (12) ◽  
pp. 1445
Author(s):  
Rodi Lykou ◽  
George Tsaklidis

Observational errors of Particle Filtering are studied over the case of a state-space model with a linear observation equation. In this study, the observational errors are estimated prior to the upcoming observations. This action is added to the basic algorithm of the filter as a new step for the acquisition of the state estimations. This intervention is useful in the presence of missing data problems mainly, as well as sample tracking for impoverishment issues. It applies theory of Homogeneous and Non-Homogeneous closed Markov Systems to the study of particle distribution over the state domain and, thus, lays the foundations for the employment of stochastic control against impoverishment. A simulating example is quoted to demonstrate the effectiveness of the proposed method in comparison with existing ones, showing that the proposed method is able to combine satisfactory precision of results with a low computational cost and provide an example to achieve impoverishment prediction and tracking.


Sign in / Sign up

Export Citation Format

Share Document