Fitting Effective Diffusion Models to Data Associated with a "Glassy Potential": Estimation, Classical Inference Procedures and Some Heuristics

Abstract

A variety of researchers have successfully obtained the parameters of low dimensional diffusion models using the data that comes out of atomistic simulations. This naturally raises a variety of questions about efficient estimation, goodness-of-fit tests, and confidence interval estimation. The first part of this article uses maximum likelihood estimation to obtain the parameters of a diffusion model from a scalar time series. I address numerical issues associated with attempting to realize asymptotic statistics results with moderate sample sizes in the presence of exact and approximated transition densities. Approximate transition densities are used because the analytic solution of a transition density associated with a parametric diffusion model is often unknown.I am primarily interested in how well the deterministic transition density expansions of Ait-Sahalia capture the curvature of the transition density in (idealized) situations that occur when one carries out simulations in the presence of a "glassy" interaction potential. Accurate approximation of the curvature of the transition density is desirable because it can be used to quantify the goodness-of-fit of the model and to calculate asymptotic confidence intervals of the estimated parameters. The second part of this paper contributes a heuristic estimation technique for approximating a nonlinear diffusion model. A "global" nonlinear model is obtained by taking a batch of time series and applying simple local models to portions of the data. I demonstrate the technique on a diffusion model with a known transition density and on data generated by the Stochastic Simulation Algorithm.

30 pages 10 figures Submitted to SIAM MMS (typos removed and slightly shortened)

References in corpus (8)

• Coarse Molecular Dynamics of a Peptide Fragment: Free Energy, Kinetics, and Long-Time Dynamics Computations
• Coarse Stability and Bifurcation Analysis Using Stochastic Simulators: Kinetic Monte Carlo Examples
• Calculating potentials of mean force and diffusion coefficients from nonequilibirum processes without Jarzynski's equality
• Fluctuation-Driven Molecular Transport in an Asymmetric Membrane Channel
• Coarse-Grained Kinetic Computations for Rare Events: Application to Micelle Formation
• Closed-form likelihood expansions for multivariate diffusions
• A test of non-equilibrium thermodynamics in glassy systems: the soft-sphere case
• On the Use of Local Diffusion Models for Path Ensemble Averaging in Potential of Mean Force Computations

Cited by in corpus (8)

• iMapD: intrinsic Map Dynamics exploration for uncharted effective free energy landscapes
• Motion Blur Filtering: A Statistical Approach for Extracting Confinement Forces and Diffusivity from a Single Blurred Trajectory
• Semi-Parametric Drift and Diffusion Estimation for Multiscale Diffusions
• Inferring Latent States and Refining Force Estimates via Hierarchical Dirichlet Process Modeling in Single Particle Tracking Experiments
• Parameter estimation for rough differential equations
• Approximating nonequilibrium processes using a collection of surrogate diffusion models
• Correcting for Bias of Molecular Confinement Parameters Induced by Small Time Series Sample Sizes in Single-Molecule Trajectories Containing Measurement Noise
• A Random Force is a Force, of Course, of Coarse: Decomposing Complex Enzyme Kinetics with Surrogate Models