qmvnorm  R Documentation 
Computes the equicoordinate quantile function of the multivariate normal
distribution for arbitrary correlation matrices
based on inversion of pmvnorm
, using a stochastic root
finding algorithm described in Bornkamp (2018).
qmvnorm(p, interval = NULL, tail = c("lower.tail", "upper.tail", "both.tails"),
mean = 0, corr = NULL, sigma = NULL, algorithm = GenzBretz(),
ptol = 0.001, maxiter = 500, trace = FALSE, seed = NULL, ...)
p 
probability. 
interval 
optional, a vector containing the endpoints of the interval to be searched. Does not need to contain the true quantile, just used as starting values by the rootfinder. If equal to NULL a guess is used. 
tail 
specifies which quantiles should be computed.

mean 
the mean vector of length n. 
corr 
the correlation matrix of dimension n. 
sigma 
the covariance matrix of dimension n. Either 
algorithm 
an object of class 
ptol , maxiter , trace 
Parameters passed to the stochastic rootfinding
algorithm. Iteration stops when the 95% confidence interval
for the predicted quantile is inside [pptol, p+ptol]. 
seed 
an object specifying if and how the random number generator
should be initialized, see 
... 
additional parameters to be passed to

Only equicoordinate quantiles are computed, i.e., the quantiles in each dimension coincide. The result is seed dependend.
A list with two components: quantile
and f.quantile
give the location of the quantile and the difference between the distribution
function evaluated at the quantile and p
.
Bornkamp, B. (2018). Calculating quantiles of noisy distribution functions using local linear regressions. Computational Statistics, 33, 487–501.
pmvnorm
, qmvt
qmvnorm(0.95, sigma = diag(2), tail = "both")
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.