Numerical evaluation of the transition probability of the simple birth-and-death process

09/24/2019
by   Alberto Pessia, et al.
0

The simple (linear) birth-and-death process is a widely used stochastic model for describing the dynamics of a population. When the process is observed discretely over time, despite the large amount of literature on the subject, little is known about formal estimator properties. Here we will show that its application to observed data is further complicated by the fact that numerical evaluation of the well-known transition probability is an ill-conditioned problem. To overcome this difficulty we will rewrite the transition probability in terms of a Gaussian hypergeometric function and subsequently obtain a three-term recurrence relation for its accurate evaluation. We will also study the properties of the hypergeometric function as a solution to the three-term recurrence relation. We will then provide formulas for the gradient and Hessian of the log-likelihood function and conclude the article by applying our methods for numerically computing maximum likelihood estimates in both simulated and real dataset.

READ FULL TEXT
research
02/14/2018

Parameter estimation for discretely-observed linear birth-and-death processes

Birth-and-death processes are widely used to model the development of bi...
research
12/30/2014

Accurate and Conservative Estimates of MRF Log-likelihood using Reverse Annealing

Markov random fields (MRFs) are difficult to evaluate as generative mode...
research
12/26/2022

The fate of the American dream: A first passage under resetting approach to income dynamics

Detailed knowledge of individual income dynamics is crucial for investig...
research
03/29/2023

On local likelihood asymptotics for Gaussian mixed-effects model with system noise

The Gaussian mixed-effects model driven by a stationary integrated Ornst...
research
12/18/2018

On a flexible construction of a negative binomial model

This work presents a construction of stationary Markov models with negat...
research
10/29/2011

Efficient Marginal Likelihood Computation for Gaussian Process Regression

In a Bayesian learning setting, the posterior distribution of a predicti...

Please sign up or login with your details

Forgot password? Click here to reset