Parameter estimation and associated uncertainty quantification is an important problem in dynamical systems characterized by ordinary differential equation (ODE) models that are often nonlinear. Typically, such models have analytically intractable trajectories which result in likelihoods and posterior distributions that are similarly intractable. Bayesian inference for ODE systems via simulation methods require numerical approximations to produce inference with high accuracy at a cost of heavy computational power and slow convergence. At the same time, Artificial Neural Networks (ANN) offer tractability that can be utilized to construct an approximate but tractable likelihood and posterior distribution. In this paper we propose a hybrid approach, where Laplace-based Bayesian inference is combined with an ANN architecture for obtaining approximations to the ODE trajectories as a function of the unknown initial values and system parameters. Suitable choices of a collocation grid and customized loss functions are proposed to fine tune the ODE trajectories and Laplace approximation. The effectiveness of our proposed methods is demonstrated using an epidemiological system with non-analytical solutions, the Susceptible-Infectious-Removed (SIR) model for infectious diseases, based on simulated and real-life influenza datasets. The novelty and attractiveness of our proposed approach include (i) a new development of Bayesian inference using ANN architectures for ODE based dynamical systems, and (ii) a computationally fast posterior inference by avoiding convergence issues of benchmark Markov Chain Monte Carlo methods. These two features establish the developed approach as an accurate alternative to traditional Bayesian computational methods, with improved computational cost.
翻译:参数估计和相关的不确定性量化是动态系统中的一个重要问题,其特点是普通差异方程模型(ODE)往往非线性分布。通常,这些模型具有分析上难解的轨迹特征,导致可能性和类似难解的近似分布。通过模拟方法对ODE系统的巴伊斯推断要求数字近似值产生高度准确的推断,以沉重的计算力和缓慢的趋同为代价。与此同时,人工神经网络(ANN)提供了可感应性,可用于构建近似但可感应的可能性和后向分布。在本文中,我们提出了一种混合方法,即基于Laplace的Bayes推论与一个ANNE(ANN)结构相结合,以获得对ODE轨迹的近似近似值和系统参数的函数。提出了合合用电网和定制损失功能的适宜选择,以调整ODE的轨迹和Labeferal的近似性方法。我们提出的方法的实效是通过非分析性解决方案、基于Labet-I-Remoal-dealimal-deal ASlimal ASyal ASyal ASyal 方法,以模拟模型为基础,以模拟模型模型为模型的模型的模型的模型和新模型的模型的模型的模型的模型的模型的模型和模型的模型的模型的模型的模型的模型和模型的模型的模型的模型。