Abstract:In the context of uncertainty quantification, computational models are required to be repeatedly evaluated. This task is intractable for costly numerical models. Such a problem turns out to be even more severe for stochastic simulators, the output of which is a random variable for a given set of input parameters. To alleviate the computational burden, surrogate models are usually constructed and evaluated instead. However, due to the random nature of the model response, classical surrogate models cannot be applied directly to the emulation of stochastic simulators. To efficiently represent the probability distribution of the model output for any given input values, we develop a new stochastic surrogate model called stochastic polynomial chaos expansions. To this aim, we introduce a latent variable and an additional noise variable, on top of the well-defined input variables, to reproduce the stochasticity. As a result, for a given set of input parameters, the model output is given by a function of the latent variable with an additive noise, thus a random variable. In this paper, we propose an adaptive algorithm which does not require repeated runs of the simulator for the same input parameters. The performance of the proposed method is compared with the generalized lambda model and a state-of-the-art kernel estimator on two case studies in mathematical finance and epidemiology and on an analytical example whose response distribution is bimodal. The results show that the proposed method is able to accurately represent general response distributions, i.e., not only normal or unimodal ones. In terms of accuracy, it generally outperforms both the generalized lambda model and the kernel density estimator.

Output error behavior for discretizations of ergodic, chaotic ODE systems

Small-sample Bayesian error estimation for ergodic, chaotic systems of ordinary differential equations

Chaotic Monte Carlo Computation: A Dynamical Effect of Random-Number Generations

Modelling pathwise uncertainty of Stochastic Differential Equations samplers via Probabilistic Numerics

Random time step probabilistic methods for uncertainty quantification in chaotic and geometric numerical integration

Numerical approximation of ergodic BSDEs using non linear Feynman-Kac formulas

Universal quantification for deterministic chaos in dynamical systems

Estimation of ordinary differential equation models with discretization error quantification

Stochastic Differential Equations for Modeling First Order Optimization Methods

Probabilistic evolution of the error of numerical method for linear stochastic differential equation

Probabilistic ODE Solvers for Integration Error-Aware Numerical Optimal Control

The random timestep Euler method and its continuous dynamics

Real-Time Monitoring of Chaotic Systems with Known Dynamical Equations

Orthogonal polynomial approximation and Extended Dynamic Mode Decomposition in chaos

Modeling still matters: a surprising instance of catastrophic floating point errors in mathematical biology and numerical methods for ODEs

Optimization of Algorithmic Errors in Analog Quantum Simulations

Numerical Studies of Statistical Management Decisions in Conditions of Stochastic Chaos

Quantum Simulation of Nonlinear Dynamical Systems Using Repeated Measurement

Hamiltonian Simulation with Random Inputs

Stochastic polynomial chaos expansions to emulate stochastic simulators

Exploiting Ergodicity for the Analysis of Short Time Series