Smooth functional tempering for nonlinear differential equation models

Authors:
David Campbell;Russell J. Steele
Affiliations:
Department of Statistics and Actuarial Science, Simon Fraser University, Surrey, Canada V3T 0A3;Department of Mathematics and Statistics, McGill University, Montreal, Canada H3A 2K6
Venue:
Statistics and Computing
Year:
2012

Citing 11
Cited 0

Scientific Computing with Ordinary Differential Equations

Scientific Computing with Ordinary Differential Equations
Deterministic Global Optimization in Nonlinear Optimal Control Problems

Journal of Global Optimization
Decoupling dynamical systems for pathway identification from metabolic profiles

Bioinformatics
Parameter estimation using Simulated Annealing for S-system models of biochemical networks

Bioinformatics
On population-based simulation for static inference

Statistics and Computing
Gaussian process modelling of latent chemical species

Bioinformatics
GNU MCSim

Bioinformatics
Monte Carlo Strategies in Scientific Computing

Monte Carlo Strategies in Scientific Computing
Structural and practical identifiability analysis of partially observed dynamical models by exploiting the profile likelihood

Bioinformatics
Estimating Bayes factors via thermodynamic integration and population MCMC

Computational Statistics & Data Analysis
Importance tempering

Statistics and Computing

Quantified Score

Hi-index	0.00

Visualization

Abstract

Differential equations are used in modeling diverse system behaviors in a wide variety of sciences. Methods for estimating the differential equation parameters traditionally depend on the inclusion of initial system states and numerically solving the equations. This paper presents Smooth Functional Tempering, a new population Markov Chain Monte Carlo approach for posterior estimation of parameters. The proposed method borrows insights from parallel tempering and model based smoothing to define a sequence of approximations to the posterior. The tempered approximations depend on relaxations of the solution to the differential equation model, reducing the need for estimating the initial system states and obtaining a numerical differential equation solution. Rather than tempering via approximations to the posterior that are more heavily rooted in the prior, this new method tempers towards data features. Using our proposed approach, we observed faster convergence and robustness to both initial values and prior distributions that do not reflect the features of the data. Two variations of the method are proposed and their performance is examined through simulation studies and a real application to the chemical reaction dynamics of producing nylon.