MathGroup Archive 2006

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: Integrate the Multivariate normal distribution


Miguel Lejeune wrote:
> I would like to compute (numerically) the following integrale (PDF is the
> probability density function of the bivariate normal distribution):
>
> Integrale of  [PDF[ndist, {x1, x2}] with respect to x2 (dx2) and the
> integration bounds are -Infinity and 1.
> with
> ndist = MultinormalDistribution[{0, 0}, r];
> r = {{1, 0.2}, {0.2, 1}};
>
> Could you indicate me how I should do?

In a bivariate normal distribution, both marginal distributions are
normal, and the conditional distribution of either variable given the
other is also normal. When the mean vector = {0,0} and the covariance
matrix = {{1,r},{r,1}}, the bivariate pdf can be written as

f[x,y,r] = f[x] * f[(y-r*x)/Sqrt[1-r^2]]/Sqrt[1-r^2],

where f[z] = Exp[-z^2/2]/Sqrt[2Pi] is the standard normal pdf. Then

Integrate[f[x,y,r],{y,-Infinity,u}] = f[z] * F[(u-r*x)/Sqrt[1-r^2]],

where F[z] = Erfc[-z/Sqrt[2]]/2 is the standard normal cdf.


  • Prev by Date: solving an equation with sums
  • Next by Date: Re: matrix substitution
  • Previous by thread: Re: Integrate the Multivariate normal distribution
  • Next by thread: Re: Re: Integrate the Multivariate normal distribution