Transient dynamics of Pearson diffusions facilitates estimation of rate parameters

Susanne Ditlevsen; Alejandro Cencerrado Rubio; Petr Lansky

doi:10.1016/j.cnsns.2019.105034

Transient dynamics of Pearson diffusions facilitates estimation of rate parameters

Susanne Ditlevsen^*, Alejandro Cencerrado Rubio, Petr Lansky

^*Corresponding author for this work

Department of Mathematical Sciences

1 Citation (Scopus)

Abstract

Estimation of parameters in stochastic processes has been thoroughly investigated for decades and the asymptotic properties of the estimators are known. However, reaching the regime where the asymptotic properties are valid might require such a long time that the well based theoretical results have no practical value. One example of this situation is at the center of our interest. It concerns determination of the time constant in stochastic Langevin equations with additive or multiplicative white-noise terms. Often the number of observations to achieve the asymptotic conditions is beyond the physical limits. Here we show how to overcome the problem by external perturbation of the system. Furthermore, we show that the perturbation is not at the price of deterioration of the estimates of other parameters unless the observation interval is very short compared to the typical time constant of the system. Three processes from the class of Pearson diffusions are studied. They are frequently used in many applications, in particular, they are examples of leaky integrate-and-fire models, which describe the electrical properties of a neuronal membrane. These neuronal models are often used as examples of systems with excitable dynamics. The most commonly investigated process is the Ornstein-Uhlenbeck process, which has additive noise. Furthermore, the square-root process and the Jacobi process are examples of processes with multiplicative noise. The results are illustrated on computer experiments, which show a striking improvement of the estimates of the rate parameter. It has implications for experimental design, where the information about the parameters can be increased for the same amount and cost of data, which is particularly important when samples are expensive or difficult to obtain.

Original language	English
Article number	105034
Journal	Communications in Nonlinear Science and Numerical Simulation
Volume	82
ISSN	1007-5704
DOIs	https://doi.org/10.1016/j.cnsns.2019.105034
Publication status	Published - Mar 2020

Keywords

Estimation of rate parameters
Neuronal models
Pearson diffusions
Perturbation

Access to Document

10.1016/j.cnsns.2019.105034

Cite this

@article{c0987eb52e88400da6092ff5fe62f20f,

title = "Transient dynamics of Pearson diffusions facilitates estimation of rate parameters",

abstract = "Estimation of parameters in stochastic processes has been thoroughly investigated for decades and the asymptotic properties of the estimators are known. However, reaching the regime where the asymptotic properties are valid might require such a long time that the well based theoretical results have no practical value. One example of this situation is at the center of our interest. It concerns determination of the time constant in stochastic Langevin equations with additive or multiplicative white-noise terms. Often the number of observations to achieve the asymptotic conditions is beyond the physical limits. Here we show how to overcome the problem by external perturbation of the system. Furthermore, we show that the perturbation is not at the price of deterioration of the estimates of other parameters unless the observation interval is very short compared to the typical time constant of the system. Three processes from the class of Pearson diffusions are studied. They are frequently used in many applications, in particular, they are examples of leaky integrate-and-fire models, which describe the electrical properties of a neuronal membrane. These neuronal models are often used as examples of systems with excitable dynamics. The most commonly investigated process is the Ornstein-Uhlenbeck process, which has additive noise. Furthermore, the square-root process and the Jacobi process are examples of processes with multiplicative noise. The results are illustrated on computer experiments, which show a striking improvement of the estimates of the rate parameter. It has implications for experimental design, where the information about the parameters can be increased for the same amount and cost of data, which is particularly important when samples are expensive or difficult to obtain.",

keywords = "Estimation of rate parameters, Neuronal models, Pearson diffusions, Perturbation",

author = "Susanne Ditlevsen and {Cencerrado Rubio}, Alejandro and Petr Lansky",

year = "2020",

month = mar,

doi = "10.1016/j.cnsns.2019.105034",

language = "English",

volume = "82",

journal = "Communications in Nonlinear Science and Numerical Simulation",

issn = "1007-5704",

publisher = "Elsevier",

}

TY - JOUR

T1 - Transient dynamics of Pearson diffusions facilitates estimation of rate parameters

AU - Ditlevsen, Susanne

AU - Cencerrado Rubio, Alejandro

AU - Lansky, Petr

PY - 2020/3

Y1 - 2020/3

N2 - Estimation of parameters in stochastic processes has been thoroughly investigated for decades and the asymptotic properties of the estimators are known. However, reaching the regime where the asymptotic properties are valid might require such a long time that the well based theoretical results have no practical value. One example of this situation is at the center of our interest. It concerns determination of the time constant in stochastic Langevin equations with additive or multiplicative white-noise terms. Often the number of observations to achieve the asymptotic conditions is beyond the physical limits. Here we show how to overcome the problem by external perturbation of the system. Furthermore, we show that the perturbation is not at the price of deterioration of the estimates of other parameters unless the observation interval is very short compared to the typical time constant of the system. Three processes from the class of Pearson diffusions are studied. They are frequently used in many applications, in particular, they are examples of leaky integrate-and-fire models, which describe the electrical properties of a neuronal membrane. These neuronal models are often used as examples of systems with excitable dynamics. The most commonly investigated process is the Ornstein-Uhlenbeck process, which has additive noise. Furthermore, the square-root process and the Jacobi process are examples of processes with multiplicative noise. The results are illustrated on computer experiments, which show a striking improvement of the estimates of the rate parameter. It has implications for experimental design, where the information about the parameters can be increased for the same amount and cost of data, which is particularly important when samples are expensive or difficult to obtain.

AB - Estimation of parameters in stochastic processes has been thoroughly investigated for decades and the asymptotic properties of the estimators are known. However, reaching the regime where the asymptotic properties are valid might require such a long time that the well based theoretical results have no practical value. One example of this situation is at the center of our interest. It concerns determination of the time constant in stochastic Langevin equations with additive or multiplicative white-noise terms. Often the number of observations to achieve the asymptotic conditions is beyond the physical limits. Here we show how to overcome the problem by external perturbation of the system. Furthermore, we show that the perturbation is not at the price of deterioration of the estimates of other parameters unless the observation interval is very short compared to the typical time constant of the system. Three processes from the class of Pearson diffusions are studied. They are frequently used in many applications, in particular, they are examples of leaky integrate-and-fire models, which describe the electrical properties of a neuronal membrane. These neuronal models are often used as examples of systems with excitable dynamics. The most commonly investigated process is the Ornstein-Uhlenbeck process, which has additive noise. Furthermore, the square-root process and the Jacobi process are examples of processes with multiplicative noise. The results are illustrated on computer experiments, which show a striking improvement of the estimates of the rate parameter. It has implications for experimental design, where the information about the parameters can be increased for the same amount and cost of data, which is particularly important when samples are expensive or difficult to obtain.

KW - Estimation of rate parameters

KW - Neuronal models

KW - Pearson diffusions

KW - Perturbation

UR - http://www.scopus.com/inward/record.url?scp=85073538492&partnerID=8YFLogxK

U2 - 10.1016/j.cnsns.2019.105034

DO - 10.1016/j.cnsns.2019.105034

M3 - Journal article

AN - SCOPUS:85073538492

SN - 1007-5704

VL - 82

JO - Communications in Nonlinear Science and Numerical Simulation

JF - Communications in Nonlinear Science and Numerical Simulation

M1 - 105034

ER -

Transient dynamics of Pearson diffusions facilitates estimation of rate parameters

Abstract

Keywords

Access to Document

Other files and links

Fingerprint

Cite this