mvorpb: Depends: R(>= 3.0.3)

Description Usage Arguments Details Value Author(s) References Examples

View source: R/mvorpb.R

Description

Orthant Probability of the Multivariate Normal Distribution

Usage

1
mvorpb( dim.p , m.tgt , v.tgt , n.itr , it.rg )

Arguments

dim.p

Dimension to evaluate (Integer, Scalar)

m.tgt

Mean vector(Real vector, length: dim.p).

v.tgt

Covariance Matrix (Real square matrix of size dim.p).

n.itr

Number of intervals for numerical integration. Too small n.itr increases the error of the returned value, but too large n.itr requires large memory size as well as computation time.

it.rg

Maximum point of the numerical integration range. Too small it.rg increases the error of the returned value due tu integration range. Estimated upper bound of caused by finiteness of it.rg will be given by the attribute with returned value. On the other hand, too large it.rg also increases the error because widths between grid points in numerical integration is inverse proportional to it.rg, and it.rg should be set with the consideration of balance it.rg and n.itr.

Details

This program evaluates the orthant probability of multivariate normal distribution based on orthogonal projections to subspaces. The covariance matrix is assumed to be nonsingular. The implemented procedure is described in Nomura (2014). The value is evaluated through applying numerical integration iteratatively. The number of the intervals for numerical integration is given by n.itr and grid points are placed equidistantly between zero and it.rg in the transformed space.

Value

The evaluated orthant probability is returned with attributes

error-itg-rg

Estimated upper bound of error caused by finiteness of the interval of the integration

.

Author(s)

Noboru Nomura

References

Nomura N. (2014). Evaluation of Gaussian orthant probabilities based on orthogonal projections to subspaces. Statistics and Computing, Online First.

Examples

1
2
3
4
5
6
7
dim.p <- 8
mu   <- c( rep(- 0.5   , dim.p -1) , 3 )
cov  <-  diag( dim.p )  * 0.5 + 0.5
n.itr <- 800
integ.range <- 10
res.val <-  mvorpb( dim.p , mu , cov  ,n.itr , integ.range )
print(res.val)

Example output

[1] 0.03887031
attr(,"error-itg-rg")
[1] 9.030753e-49

mvprpb documentation built on May 29, 2017, 7:55 p.m.