Abstract:This work considers the computation of risk measures for quantities of interest governed by PDEs with Gaussian random field parameters using Taylor approximations. While efficient, Taylor approximations are local to the point of expansion, and hence may degrade in accuracy when the variances of the input parameters are large. To address this challenge, we approximate the underlying Gaussian measure by a mixture of Gaussians with reduced variance in a dominant direction of parameter space. Taylor approximations are constructed at the means of each Gaussian mixture component, which are then combined to approximate the risk measures. The formulation is presented in the setting of infinite-dimensional Gaussian random parameters for risk measures including the mean, variance, and conditional value-at-risk. We also provide detailed analysis of the approximations errors arising from two sources: the Gaussian mixture approximation and the Taylor approximations. Numerical experiments are conducted for a semilinear advection-diffusion-reaction equation with a random diffusion coefficient field and for the Helmholtz equation with a random wave speed field. For these examples, the proposed approximation strategy can achieve less than $1\%$ relative error in estimating CVaR with only $\mathcal{O}(10)$ state PDE solves, which is comparable to a standard Monte Carlo estimate with $\mathcal{O}(10^4)$ samples, thus achieving significant reduction in computational cost. The proposed method can therefore serve as a way to rapidly and accurately estimate risk measures under limited computational budgets.

An Application of the EM-algorithm to Approximate Empirical Distributions of Financial Indices with the Gaussian Mixtures

An Effective EM Algorithm for Mixtures of Gaussian Processes Via the MCMC Sampling and Approximation.

A Robust and Flexible EM Algorithm for Mixtures of Elliptical Distributions with Missing Data

The generalized scale mixtures of asymmetric generalized normal distributions with application to stock data

An Mcmc Based Em Algorithm For Mixtures Of Gaussian Processes

On the Behavior of the Expectation-Maximization Algorithm for Mixture Models

Statistical Convergence of the EM Algorithm on Gaussian Mixture Models

Comparative Analysis on Convergence Rates of the Em Algorithm and Its Two Modifications for Gaussian Mixtures

The Competitive EM Algorithm for Gaussian Mixtures with BYY Harmony Criterion

Mixtures of Gaussians and Minimum Relative Entropy Techniques for Modeling Continuous Uncertainties

Mixture of Gaussian Processes and its Applications

Gaussian mixture Taylor approximations of risk measures constrained by PDEs with Gaussian random field inputs

Approximating the Differential Entropy of Gaussian Mixtures.

A BYY Scale-Incremental EM Algorithm for Gaussian Mixture Learning

Distributed EM Algorithm for Gaussian Mixtures in Sensor Networks

Differentiability and Approximation of Probability Functions under Gaussian Mixture Models: A Bayesian Approach

A Novel Split and Merge EM Algorithm for Gaussian Mixture Model

Gaussian Mixture Models with Rare Events.

Value at Risk and Expected Shortfall for Normal Variance Mean Mixtures of Finite Weighted Inverse Gaussian Distributions

Portfolio Risk Measurement Using a Mixture Simulation Approach

Research on Correct Convergence of the EM Algorithm for Gaussian Mixtures