A Nonparametric Bayesian Framework for Uncertainty Quantification in Stochastic Simulation

When we use simulation to assess the performance of stochastic systems, the input models used to drive simulation experiments are often estimated from finite real-world data. There exist both input model and simulation estimation uncertainties in the system performance estimates. Without strong prior information on the input models and the system mean response surface, in this paper, we propose a Bayesian nonparametric framework to quantify the impact from both sources of uncertainty. Specifically, since the real-world data often represent the variability caused by various latent sources of uncertainty, Dirichlet Processes Mixtures (DPM) based nonparametric input models are introduced to model a mixture of heterogeneous distributions, which can faithfully capture the important features of real-world data, such as multi-modality and skewness. Bayesian posteriors of flexible input models characterize the input model estimation uncertainty, which automatically accounts for both model selection and parameter value uncertainty. Then, input model estimation uncertainty is propagated to outputs by using direct simulation. Thus, under very general conditions, our framework delivers an empirical credible interval accounting for both input and simulation uncertainties. A variance decomposition is further developed to quantify the relative contributions from both sources of uncertainty. Our approach is supported by rigorous theoretical and empirical study.

[1]  Canan G. Corlu,et al.  Simulation of inventory systems with unknown input models: a data-driven approach , 2017, Int. J. Prod. Res..

[2]  T. Ferguson A Bayesian Analysis of Some Nonparametric Problems , 1973 .

[3]  M. Escobar,et al.  Bayesian Density Estimation and Inference Using Mixtures , 1995 .

[4]  Wei Xie,et al.  An Efficient Budget Allocation Approach for Quantifying the Impact of Input Uncertainty in Stochastic Simulation , 2017, ACM Trans. Model. Comput. Simul..

[5]  Alexander J. Smola,et al.  Kernel methods and the exponential family , 2006, ESANN.

[6]  Wei Xie,et al.  A Bayesian Framework for Quantifying Uncertainty in Stochastic Simulation , 2014, Oper. Res..

[7]  E. Parzen On Estimation of a Probability Density Function and Mode , 1962 .

[8]  Wei Xie,et al.  A Factor-Based Bayesian Framework for Risk Analysis in Stochastic Simulations , 2017, ACM Trans. Model. Comput. Simul..

[9]  J. Ghosh,et al.  POSTERIOR CONSISTENCY OF DIRICHLET MIXTURES IN DENSITY ESTIMATION , 1999 .

[10]  S. MacEachern Estimating normal means with a conjugate style dirichlet process prior , 1994 .

[11]  P. Müller,et al.  Bayesian curve fitting using multivariate normal mixtures , 1996 .

[12]  Michael P. Wiper,et al.  Mixtures of Gamma Distributions With Applications , 2001 .

[13]  Heng Lian Cross-validation for comparing multiple density estimation procedures , 2007 .

[14]  L. Wasserman,et al.  Consistency of Bernstein polynomial posteriors , 2002 .

[15]  Wei Xie,et al.  Stochastic simulation under input uncertainty: A Review , 2020 .

[16]  Barry L. Nelson,et al.  Stochastic kriging for simulation metamodeling , 2008, 2008 Winter Simulation Conference.

[17]  Bahar Biller,et al.  Copula-Based Multivariate Input Models for Stochastic Simulation , 2009, Oper. Res..

[18]  Carl E. Rasmussen,et al.  Dirichlet Process Gaussian Mixture Models: Choice of the Base Distribution , 2010, Journal of Computer Science and Technology.

[19]  Wei Xie,et al.  Quantifying Input Uncertainty via Simulation Confidence Intervals , 2014, INFORMS J. Comput..

[20]  J. Ghosh,et al.  On convergence of posterior distributions , 1995 .

[21]  Bradley P. Carlin,et al.  Bayesian measures of model complexity and fit , 2002 .

[22]  Albert Y. Lo,et al.  On a Class of Bayesian Nonparametric Estimates: I. Density Estimates , 1984 .

[23]  Stephen E. Chick,et al.  Input Distribution Selection for Simulation Experiments: Accounting for Input Uncertainty , 2001, Oper. Res..

[24]  Lianming Wang,et al.  Fast Bayesian Inference in Dirichlet Process Mixture Models , 2011, Journal of computational and graphical statistics : a joint publication of American Statistical Association, Institute of Mathematical Statistics, Interface Foundation of North America.