5

If we have some multi-dimensional normal probability distribution, $X \sim \mathcal{N}(0,\Sigma)$, with zero mean and a known covariance matrix then what is the probability that every component of $X$ is greater than or equal to zero?

$$\mathbb{P}(\forall i.X_i \ge 0) = \: ?$$

ezeidman
  • 147

1 Answers1

9

Stuart and Ord (1994, Kendall's Advanced Theory of Statistics, volume 1, Section 15.10) report symbolic solutions for the standardised bivariate Normal orthant probability $P(X>0,Y>0)$ as:

$$\text{P2}=\frac{\sin ^{-1}(\rho )}{2 \pi }+\frac{1}{4}$$

... while the standardised trivariate Normal orthant probability $P(X>0,Y>0,Z>0)$ is:

$$\text{P3}=\frac{\sin ^{-1}\left(\rho _{\text{xy}}\right)+\sin ^{-1}\left(\rho _{\text{xz}}\right)+\sin ^{-1}\left(\rho _{\text{yz}}\right)}{4 \pi }+\frac{1}{8}$$

I have found these exact symbolic solutions to be useful when testing the accuracy of multivariate Normal cdf functions (which use numerical integration to solve ... unless they know about the special cases :) )

Apparently, no such simple result is available for higher than 3 dimensions ... but there are tables, and some special results reduced to smaller numbers of integrals ...

wolfies
  • 5,174
  • 2
    The results for $n\le 5$ are given in here https://math.stackexchange.com/questions/869502/multivariate-gaussian-integral-over-positive-reals/3148280#3148280 . – Przemo Apr 30 '19 at 15:00
  • 2
    Proofs can be found here: https://math.stackexchange.com/q/255368/321264, https://stats.stackexchange.com/q/379255/119261. – StubbornAtom May 18 '20 at 19:15
  • Does this result still work if the covariance matrix is singular? The derivation relies on the fact that a density exists. –  May 04 '23 at 04:48