Calculates the Multivariate Normal Distribution with Product Correlation Structure published by Charles Dunnett, Algorithm AS 251.1 Appl.Statist. (1989), Vol.38, No.3, tools:::Rd_expr_doi("10.2307/2347754").
mvnprd(..., A, B, BPD, EPS = 1e-06, INF, IERC = 1, HINC = 0)Ensures that all arguments (starting from the "...") are to be named and that a warning will be displayed if unknown arguments are passed.
Upper limits of integration. Array of N dimensions
Lower limits of integration. Array of N dimensions
Values defining correlation structure. Array of N dimensions
desired accuracy. Defaults to 1e-06
Determines where integration is done to infinity. Array of N dimensions. Valid values for INF(I): 0 = c(B(I), Inf), 1 = c(-Inf, A(I)), 2 = c(B(I), A(I))
error control. If set to 1, strict error control based on fourth derivative is used. If set to zero, error control based on halving intervals is used
Interval width for Simpson's rule. Value of zero caused a default .24 to be used
This is a wrapper function for the original Fortran 77 code. For a multivariate normal vector with correlation structure defined by RHO(I,J) = BPD(I) * BPD(J), computes the probability that the vector falls in a rectangle in n-space with error less than eps.