Approximating the first passage time density from data using generalized Laguerre polynomials

Abstract

This paper analyzes a method to approximate the first passage time probability density function which turns to be particularly useful if only sample data are available. The method relies on a Laguerre-Gamma polynomial approximation and iteratively looks for the best degree of the polynomial such that the fitting function is a probability density function. The proposed iterative algorithm relies on simple and new recursion formulae involving first passage time moments. These moments can be computed recursively from cumulants, if they are known. In such a case, the approximated density can be used also for the maximum likelihood estimates of the parameters of the underlying stochastic process. If cumulants are not known, suitable unbiased estimators relying on k-statistics are employed. To check the feasibility of the method both in fitting the density and in estimating the parameters, the first passage time problem of a geometric Brownian motion is considered.

arXiv comment22 pages, 6 figures

References in corpus (4)

• Long- and short-time asymptotics of the first-passage time of the Ornstein-Uhlenbeck and other mean-reverting processes
• Approximation to Distribution of Product of Random Variables Using Orthogonal Polynomials for Lognormal Density
• Characterization of the cubic exponential families by orthogonality of polynomials
• Symbolic Calculus in Mathematical Statistics: A Review