Computes the distribution function of the truncated multivariate normal
distribution for arbitrary limits and correlation matrices
based on the pmvnorm()
implementation of the algorithms by Genz and Bretz.
ptmvnorm(lowerx, upperx, mean=rep(0, length(lowerx)), sigma,
lower = rep(-Inf, length = length(mean)),
upper = rep( Inf, length = length(mean)),
maxpts = 25000, abseps = 0.001, releps = 0)
the vector of lower limits of length n.
the vector of upper limits of length n.
the mean vector of length n.
the covariance matrix of dimension n. Either corr
or
sigma
can be specified. If sigma
is given, the
problem is standardized. If neither corr
nor
sigma
is given, the identity matrix is used
for sigma
.
Vector of lower truncation points,\
default is rep(-Inf, length = length(mean))
.
Vector of upper truncation points,\
default is rep( Inf, length = length(mean))
.
maximum number of function values as integer.
absolute error tolerance as double.
relative error tolerance as double.
The evaluated distribution function is returned with attributes
estimated absolute error and
status messages.
The computation of truncated multivariate normal probabilities and densities is done using conditional probabilities
from the standard/untruncated multivariate normal distribution.
So we refer to the documentation of the mvtnorm
package and the methodology is described in
Genz (1992, 1993) and Genz/Bretz (2009).
For properties of the truncated multivariate normal distribution see for example Johnson/Kotz (1970) and Horrace (2005).
Genz, A. (1992). Numerical computation of multivariate normal probabilities. Journal of Computational and Graphical Statistics, 1, 141--150
Genz, A. (1993). Comparison of methods for the computation of multivariate normal probabilities. Computing Science and Statistics, 25, 400--405
Genz, A. and Bretz, F. (2009). Computation of Multivariate Normal and t Probabilities. Lecture Notes in Statistics, Vol. 195, Springer-Verlag, Heidelberg.
Johnson, N./Kotz, S. (1970). Distributions in Statistics: Continuous Multivariate Distributions Wiley & Sons, pp. 70--73
Horrace, W. (2005). Some Results on the Multivariate Truncated Normal Distribution. Journal of Multivariate Analysis, 94, 209--221
# NOT RUN {
sigma <- matrix(c(5, 0.8, 0.8, 1), 2, 2)
Fx <- ptmvnorm(lowerx=c(-1,-1), upperx=c(0.5,0), mean=c(0,0),
sigma=sigma, lower=c(-1,-1), upper=c(1,1))
# }
Run the code above in your browser using DataLab