MC_Gauss: MC estimate for the remainder

Description Usage Arguments Value References

View source: R/MC0_2.R

Description

Standard Monte Carlo estimate for P(max X^{-q} >threshold | max X^{q}≤ threshold) or P(min X^{-q} <threshold | min X^{q}≥ threshold) where X is a normal vector. It is used for the bias correction in ProbaMax and ProbaMin.

Usage

1
2
3
MC_Gauss(compBdg, problem, delta = 0.1, type = "M",
  trmvrnorm = trmvrnorm_rej_cpp, typeReturn = 0, verb = 0,
  params = NULL)

Arguments

compBdg

total computational budget in seconds.

problem

list defining the problem with mandatory fields:

  • muEq = mean vector of X^{q};

  • sigmaEq = covariance matrix of X^q;

  • threshold = threshold;

  • muEmq = mean vector of X^{-q};

  • wwCondQ = “weights” for X^{-q} | X^q [ the vector Σ^{-q,q}(Σ^q)^{-1}];

  • sigmaCondQChol = Cholesky factorization of the conditional covariance matrix Σ^{-q | q}.

delta

total proportion of budget assigned to initial estimate (default 0.1), the actual proportion used might be smaller.

type

type of excursion: "m", for minimum below threshold or "M", for maximum above threshold.

trmvrnorm

function to generate truncated multivariate normal samples, it must have the following signature trmvrnorm(n,mu,sigma,upper,lower,verb), where

  • n: number of simulations;

  • mu: mean vector of the Normal variable of dimension d;

  • sigma: covariance matrix of dimension d x d;

  • upper: vector of upper limits of length d;

  • lower: vector of lower limits of length d;

  • verb: the level of verbosity 3 basic, 4 extended.

It must return a matrix d x n of realizations. If not specified, the rejection sampler trmvrnorm_rej_cpp is used.

typeReturn

integer: 0 (only the estimate) or 1 (heavy return with variance of the estimate, parameters of the estimator and computational times).

verb

the level of verbosity, also sets the verbosity of trmvrnorm (to verb-1).

params

system dependent parameters (if NULL they are estimated).

Value

A list containing the estimated probability of excursion, see typeReturn for details.

References

Azzimonti, D. and Ginsbourger, D. (2018). Estimating orthant probabilities of high dimensional Gaussian vectors with an application to set estimation. Journal of Computational and Graphical Statistics, 27(2), 255-267. Preprint at hal-01289126

Azzimonti, D. (2016). Contributions to Bayesian set estimation relying on random field priors. PhD thesis, University of Bern.

Genz, A. (1992). Numerical computation of multivariate normal probabilities. Journal of Computational and Graphical Statistics, 1(2):141–149.


anMC documentation built on Oct. 30, 2019, 11:41 a.m.