Abstract
When models are defined implicitly by systems of differential equations with no closed-form solution, small local errors in finite-dimensional solution approximations can propagate into deviations from the true underlying model trajectory. Some recent perspectives in quantifying this uncertainty are based on Bayesian probability modeling: a prior is defined over the unknown solution and updated by conditioning on interrogations of the forward model. Improvement in accuracy via grid refinement must be considered in order for such Bayesian numerical methods to compete with state-of-the-art numerical techniques. We review the principles of Bayesian statistical design and apply these to develop an adaptive probabilistic method to sequentially select time-steps for state-space probabilistic ODE solvers. We investigate the behavior of local error under the adaptive scheme which underlies numerical variable step-size methods. Numerical experiments are used to illustrate the performance of such an adaptive scheme, showing improved accuracy over uniform designs in terms of local error.
Similar content being viewed by others
References
Alexanderian, A., Gloor, P.J., Ghattas, O.: On Bayesian A- and D-optimal experimental designs in infinite dimensions. Bayesian Anal. 11(3), 671 (2016). https://doi.org/10.1214/15-BA969
Bernardo, J.M., Smith, A.F.M.: Bayesian Theory. Wiley, New York (1994)
Butcher, J.: Numerical Methods for Ordinary Differential Equations. Wiley, Hoboken (2008). https://doi.org/10.1002/9780470753767
Chaloner, K., Verdinelli, I.: Bayesian experimental design: a review. Stat. Sci. 10, 273 (1995). https://doi.org/10.1214/ss/1177009939
Chkrebtii, O., Campbell, D., Calderhead, B., Girolami, M.: Bayesian solution uncertainty quantification for differential equations. Bayesian Anal. 11(4), 1239 (2016a). https://doi.org/10.1214/16-BA1017
Chkrebtii, O., Campbell, D., Calderhead, B., Girolami, M.: Rejoinder. Bayesian Anal. 4(11), 1295 (2016b). https://doi.org/10.1214/16-BA1017REJ
Cockayne, J., Oates, C., Sullivan, T., Girolami, M.: Bayesian probabilistic numerical methods (2017). arXiv:1702.03673
Diaconis, P.: Bayesian Numerical Analysis, pp. 163–175. Springer, New York (1988)
Hennig, P., Hauberg, S.: Probabilistic solutions to differential equations and their application to Riemannian statistics. In: Proceedings of the 17th International Conference on Artificial Intelligence and Statistics (AISTATS), vol. 33 (2014)
Higdon, D., Gattiker, J., Lawrence, E., Jackson, C., Tobis, M., Pratola, M., Habib, S., Heitmann, K., Price, S.: Computer model calibration using the ensemble Kalman filter. Technometrics 55(4), 488 (2013). https://doi.org/10.1080/00401706.2013.842936
Jazwinski, A.: Stochastic Processes and Filtering Theory. Mathematics in Science and Engineering. Elsevier Science, Amsterdam (1970). https://doi.org/10.1016/S0076-5392(09)60368-4
Kassam, A., Trefethen, L.: Fourth-order time-stepping for stiff PDEs. SIAM J. Sci. Comput. 26(4), 1214 (2005). https://doi.org/10.1137/S1064827502410633
Kullback, S., Leibler, R.A.: On information and sufficiency. Ann. Math. Stat. 22, 49 (1951). https://doi.org/10.1214/aoms/1177729694
Kuramoto, Y., Tsuzuki, T.: Persistent propagation of concentration waves in dissipative media far from thermal equilibrium. Prog. Theor. Phys. 55, 356 (1976). https://doi.org/10.1143/PTP.55.356
Lorenz, E.N.: Deterministic nonperiodic flow. J. Atmos. Sci. 20, 130 (1963). https://doi.org/10.1175/1520-0469
O’Hagan, A.O.: Some Bayesian numerical analysis. Bayesian Stat. 4, 345 (1992)
Paciorek, C.J., Schervish, M.J.: In: Thrun, S., Saul, L.K., Schölkopf, B. (eds.) Advances in Neural Information Processing Systems, vol. 16, pp. 273–280. MIT Press, Cambridge (2004)
Paciorek, C., Schervish, M.: Spatial modelling using a new class of nonstationary covariance functions. Environmetrics 17(5), 483 (2006). https://doi.org/10.1002/env.785
Schober, M., Duvenaud, D.K., Hennig, P.: In: Ghahramani, Z., Welling, M., Cortes, C., Lawrence, N., Weinberger, K. (eds.) Advances in Neural Information Processing Systems, vol. 27, pp. 739–747. Curran Associates, Inc. Red Hook (2014)
Sebastiani, P., Wynn, H.: Maximum entropy sampling and optimal Bayesian experimental design. J. R. Stat. Soc. Ser. B (Stat. Methodol.) 62(1), 145 (2000). https://doi.org/10.1111/1467-9868.00225
Shampine, L.F.: Error estimation and control for ODEs. J. Sci. Comput. 25(1), 3 (2005). https://doi.org/10.1007/BF02728979
Skilling, J.: Bayesian Solution of Ordinary Differential Equations, pp. 23–37. Kluwer Academic Publishers, Seattle (1991)
Stuart, A.M.: Inverse problems: a Bayesian perspective. Acta Numer. 19, 451 (2010). https://doi.org/10.1017/S0962492910000061
Wang, J., Cockayne, J., Oates, C.: On the Bayesian solution of differential equations (2018). arXiv:1805.07109
Wikle, C.K., Berliner, L.M.: A Bayesian tutorial for data assimilation. Physica D 230(1), 1 (2007). https://doi.org/10.1016/j.physd.2006.09.017
Acknowledgements
This material was developed, in part, at the Prob Num 2018 workshop hosted by the Lloyd’s Register Foundation programme on Data-Centric Engineering at the Alan Turing Institute, UK, and supported by the National Science Foundation, USA. This material was based upon work partially supported by the National Science Foundation under Grant DMS-1127914 to the Statistical and Applied Mathematical Sciences Institute (SAMSI) and the Isaac Newton Institute for Mathematical Sciences’ program on Uncertainty Quantification for Complex Systems: theory and methodologies. Any opinions, findings, conclusions or recommendations expressed in this material are those of the authors and do not necessarily reflect the views of the above-named funding bodies and research institutions.
Author information
Authors and Affiliations
Corresponding author
Additional information
Handling Editor: T. J. Sullivan.
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
About this article
Cite this article
Chkrebtii, O.A., Campbell, D.A. Adaptive step-size selection for state-space probabilistic differential equation solvers. Stat Comput 29, 1285–1295 (2019). https://doi.org/10.1007/s11222-019-09899-5
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11222-019-09899-5