Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Half-normal distribution

  1. Sep 4, 2011 #1
    Hi: I am reading an article that deals with the distribution function associated with the half-normal distribution. The author presents a formula for the c.d.f. as:

    [tex] \left [ G\left ( x \right ) \right ]^{2r-2}=\left ( \frac{2}{\pi } \right )^{r-1}\left \{ \int_{0}^{x}\exp \left [ -\frac{1}{2} \left (w ^{2}\right )\right ]dw \right \}^{2r-2} [/tex]

    Note that r=2,3,...

    The expression above is also equal to:

    [tex] \left [ G\left ( x \right ) \right ]^{2r-2}=\left ( \frac{2}{\pi } \right )^{r-1}\left \{ \int_{0}^{x}\int_{0}^{x} \exp \left [- \frac{1}{2}\left ( x_{1}^{2}+x_{2}^{2} \right ) \right ]dx_{1}dx_{2}\right \}^{r-1} [/tex]

    which I have no problem with. The author then states that the second equation above is also equal to:

    [tex] \left [ G\left ( x \right ) \right ]^{2r-2}=\left ( \frac{4}{\pi } \right )^{r-1}\left \{ \int_{0}^{\frac{\pi }{4}}\left [ 1-\exp \left ( -\frac{1}{2}x^{2}\sec ^{2}\theta \right ) \right ]d\theta \right \}^{r-1} [/tex]

    Can someone explain to me how the author gets from the second equation to the third (or last) equation above.

    Thanks in advance.
  2. jcsd
  3. Sep 5, 2011 #2
    Since an angle has been introduced in the final result, I would try to express the integral using polar coordinates in the x1,x2 plane.
  4. Sep 5, 2011 #3
    You basically want to show this:

    [tex] \int_{0}^{x}\int_{0}^{x} \exp \left [- \frac{1}{2}\left ( x_{1}^{2}+x_{2}^{2} \right ) \right ]dx_{1}dx_{2} = 2 \int_{0}^{\frac{\pi }{4}}\left [ 1-\exp \left ( -\frac{1}{2}x^{2}\sec ^{2}\theta \right ) \right ]d\theta [/tex]

    Since you are integrating over a rectangle from (0,0) to (x,x) here, you can't really use spherical coordinates.

    You can, however, use an angle theta=0..pi/4 and one cartesian coordinate, i.e. x_1 to describe half of your rectangle. The transformation should look kinda like

    [tex] x_1=x_1 [/tex]
    [tex] x_2 = x_1 \tan \theta [/tex]

    For the argument inside the exponential function, you will get
    [tex] x_1^2 + x_2^2 = x_1 (\tan^2\theta + 1) = x_1^2 \sec^2 \theta [/tex]
    For your Jacobian, you will get:
    [tex] x_1 (\tan\theta)' = x_1 \sec^2\theta [/tex]
    Once you do the x_1-Integration, the theta-Part of the Jacobian will go away.
    Last edited: Sep 5, 2011
  5. Sep 5, 2011 #4

    Okay, thanks.

    Note that I think you meant to write:

    [tex] x_1^2 + x_2^2 = x_1^2 (\tan^2\theta + 1) = x_1^2 \sec^2 \theta [/tex]
    Last edited: Sep 5, 2011
  6. Sep 5, 2011 #5
    Yup, of course, sry ;)
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook