Parameter Estimation in Nonlinear Multivariate Stochastic Differential Equations Based on Splitting Schemes
arXiv:2211.11884 · doi:10.1214/24-AOS2371
Abstract
The likelihood functions for discretely observed nonlinear continuous-time models based on stochastic differential equations are not available except for a few cases. Various parameter estimation techniques have been proposed, each with advantages, disadvantages, and limitations depending on the application. Most applications still use the Euler-Maruyama discretization, despite many proofs of its bias. More sophisticated methods, such as Kessler's Gaussian approximation, Ozaki's Local Linearization, Aït-Sahalia's Hermite expansions, or MCMC methods, might be complex to implement, do not scale well with increasing model dimension, or can be numerically unstable. We propose two efficient and easy-to-implement likelihood-based estimators based on the Lie-Trotter (LT) and the Strang (S) splitting schemes. We prove that S has convergence rate of order 1, a property already known for LT. We show that the estimators are consistent and asymptotically efficient under the less restrictive one-sided Lipschitz assumption. A numerical study on the 3-dimensional stochastic Lorenz system complements our theoretical findings. The simulation shows that the S estimator performs the best when measured on precision and computational speed compared to the state-of-the-art.
References in corpus (10)
- Strong and weak divergence in finite time of Euler's method for stochastic differential equations with non-globally Lipschitz continuous coefficients
- A fundamental mean-square convergence theorem for SDEs with locally Lipschitz coefficients and its applications
- Long-Run Accuracy of Variational Integrators in the Stochastic Context
- Maximum-likelihood estimation for diffusion processes via closed-form density expansions
- Splitting and composition methods in the numerical integration of differential equations
- Practical Estimation of High Dimensional Stochastic Differential Mixed-Effects Models
- Analysis of Some Splitting Schemes for the Stochastic Allen-Cahn Equation
- On the approximate maximum likelihood estimation for diffusion processes
- Network inference via approximate Bayesian computation. Illustration on a stochastic multi-population neural mass model
- Parameter Estimation with Increased Precision for Elliptic and Hypo-elliptic Diffusions