Services & Resources / Wolfram Forums
-----
 /
MathGroup Archive
2000
*January
*February
*March
*April
*May
*June
*July
*August
*September
*October
*November
*December
*Archive Index
*Ask about this page
*Print this page
*Give us feedback
*Sign up for the Wolfram Insider

MathGroup Archive 2000

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

Search the Archive

Integrating Probability Surface_HOW?

  • To: mathgroup at smc.vnet.net
  • Subject: [mg22176] Integrating Probability Surface_HOW?
  • From: beintinc at aol.com (BEIntInc)
  • Date: Thu, 17 Feb 2000 01:24:08 -0500 (EST)
  • Sender: owner-wri-mathgroup at wolfram.com

I need to find the cumulative probability density contained within an
elliptical region on the bivariate Gaussian surface.  Everything outside the
elliptical region is zero, the interior of the elliptical region is magnitude
1; i.e., the ellipse "windows" the PDF.

To see this, run the script below (Mathematica 4.0) , and you will see the
rotated 4-sigma contour of the PDF, along with an ellipse (dashed).  I used the
additional PDF (ndist1) simply to form the shape of the windowing ellipse.

My problem is, how to find the volume under the dashed ellipse; in other words,
the volume contained within the ellipse is the desired CDF of the rotated PDF. 
(I am probably not precise here, but hopefully it makes sense).


I was hoping I could do

ff = EllipsoidQuantile[ndist1, .3935]
RegionProbability[ndist, ff]

but that does not work.

Any help greatly appreciated,

Regards,

John Brooks

%***************************************************************************
<< Statistics`MultinormalDistribution`

r = {{0.14606165662744, -0.08743998612085}, {-0.08743998612085, 
          0.08233731794636}}*10^-4;
ndist = MultinormalDistribution[{0, 0}, r]

pdf = PDF[ndist, {x1, x2}]

BW1 = .00912442396313364^2

BW2 = .01446829793425295^2

r1 = {{BW2, 0}, {0, BW1}};
ndist1 = MultinormalDistribution[{0, 0}, r1]

pdf1 = PDF[ndist1, {x1, x2}]

Show[Graphics[{EllipsoidQuantile[ndist, .9996], Dashing[{.04, .03}], 
      EllipsoidQuantile[ndist1, .3935]}, Axes -> True, GridLines -> Automatic,
     AspectRatio -> Automatic, Frame -> True]]

%***********************************************************************


  • Prev by Date: Re: Graphics and animation
  • Next by Date: Re: FindMinimum
  • Previous by thread: Some more questions Part II: Mathematica
  • Next by thread: Re: Linux: how compatible with add-ons ?