Documentation |
Multivariate t cumulative distribution function
y = mvtcdf(X,C,DF)
y = mvtcdf(xl,xu,C,DF)
[y,err] = mvtcdf(...)
[...] = mvntdf(...,options)
y = mvtcdf(X,C,DF) returns the cumulative probability of the multivariate t distribution with correlation parameters C and degrees of freedom DF, evaluated at each row of X. Rows of the n-by-d matrix X correspond to observations or points, and columns correspond to variables or coordinates. y is an n-by-1 vector.
C is a symmetric, positive definite, d-by-d matrix, typically a correlation matrix. If its diagonal elements are not 1, mvtcdf scales C to correlation form. mvtcdf does not rescale X. DF is a scalar, or a vector with n elements.
The multivariate t cumulative probability at X is defined as the probability that a random vector T, distributed as multivariate t, will fall within the semi-infinite rectangle with upper limits defined by X, i.e., Pr{T(1)≤X(1),T(2)≤X(2),...T(d)≤X(d)}.
y = mvtcdf(xl,xu,C,DF) returns the multivariate t cumulative probability evaluated over the rectangle with lower and upper limits defined by xl and xu, respectively.
[y,err] = mvtcdf(...) returns an estimate of the error in y. For bivariate and trivariate distributions, mvtcdf uses adaptive quadrature on a transformation of the t density, based on methods developed by Genz, as described in the references. The default absolute error tolerance for these cases is 1e-8. For four or more dimensions, mvtcdf uses a quasi-Monte Carlo integration algorithm based on methods developed by Genz and Bretz, as described in the references. The default absolute error tolerance for these cases is 1e-4.
[...] = mvntdf(...,options) specifies control parameters for the numerical integration used to compute y. This argument can be created by a call to statset. Choices of statset parameters are:
'TolFun' — Maximum absolute error tolerance. Default is 1e-8 when d < 4, or 1e-4 when d ≥ 4.
'MaxFunEvals' — Maximum number of integrand evaluations allowed when d ≥ 4. Default is 1e7. 'MaxFunEvals' is ignored when d < 4.
'Display' — Level of display output. Choices are 'off' (the default), 'iter', and 'final'. 'Display' is ignored when d < 4.
[1] Genz, A. "Numerical Computation of Rectangular Bivariate and Trivariate Normal and t Probabilities." Statistics and Computing. Vol. 14, No. 3, 2004, pp. 251–260.
[2] Genz, A., and F. Bretz. "Numerical Computation of Multivariate t Probabilities with Application to Power Calculation of Multiple Contrasts." Journal of Statistical Computation and Simulation. Vol. 63, 1999, pp. 361–378.
[3] Genz, A., and F. Bretz. "Comparison of Methods for the Computation of Multivariate t Probabilities." Journal of Computational and Graphical Statistics. Vol. 11, No. 4, 2002, pp. 950–971.