Monte Carlo profile confidence intervals for dynamic systems

J R Soc Interface. 2017 Jul;14(132):20170126. doi: 10.1098/rsif.2017.0126.

Abstract

Monte Carlo methods to evaluate and maximize the likelihood function enable the construction of confidence intervals and hypothesis tests, facilitating scientific investigation using models for which the likelihood function is intractable. When Monte Carlo error can be made small, by sufficiently exhaustive computation, then the standard theory and practice of likelihood-based inference applies. As datasets become larger, and models more complex, situations arise where no reasonable amount of computation can render Monte Carlo error negligible. We develop profile likelihood methodology to provide frequentist inferences that take into account Monte Carlo uncertainty. We investigate the role of this methodology in facilitating inference for computationally challenging dynamic latent variable models. We present examples arising in the study of infectious disease transmission, demonstrating our methodology for inference on nonlinear dynamic models using genetic sequence data and panel time-series data. We also discuss applicability to nonlinear time-series and spatio-temporal data.

Keywords: likelihood-based inference; panel data; phylodynamic inference; sequential Monte Carlo; spatio-temporal data; time series.

Publication types

  • Research Support, N.I.H., Extramural
  • Research Support, U.S. Gov't, Non-P.H.S.

MeSH terms

  • Animals
  • Computer Simulation
  • Likelihood Functions
  • Models, Biological*
  • Models, Statistical
  • Monte Carlo Method*
  • Population Dynamics
  • Time Factors