Math 300 , Spring 2004, Day 29, W, April 14 Hit reload to get most current version

 COW Calculus tutor, see day 28
Continuing with E(X), Var(X)
Gamma, wait till n'th hit, parameter lambda
The waiting time X can be thought of as the sum of the waiting time till first hit + interval till 2nd +...+ interval till n'th.
X = X1+X2+...+Xn
 Since it's a Poisson process, all these Xi are independent exponential.  (Cf. Neg. binomial).
So E(X) = n times expected for exponential, Var(X) = n times Var for exponential.
Normal N(0, 1):   E(X) = 0, done.  Guts of E(X2)  are x2 e-x2/2 , and we'll try integrating by parts.
Let v = x, du = x e-x2/2dx,   and  u = -e-x2/2  , dv = 1.  When we tidy up (the uv term evaluates to 0),
we get integral of Ce-x2/2 dx  where C = 1/sqrt (2pi).
But that's just the integral over the pdf for the standard normal, so it = 1!  So Var(X) = 1..

Got to here Wednesday.
Return to Ch. 4, Sec.4.4, (4.5 optional) 4.6.
Normal distributionN(mu, sigma)

 f(x) = normal formulaF(x) = integral of f. There is no "closed form" formula in elementary functions for F(x).

The parameters are the mean and standard deviation for the normal distribution. You've been told that but we have not proved it. We just proved that the mean and the standard deviation for the Standard Normal distribution Z were 0 and 1.  The proofs from the general formula "become" the same proofs if you make the change of variable (substitution) and change the limits of integration too. (Pp. 224-26)
- - - - - - - - - - - - - - - - - - -

If we have a Normal variable X, and we "standardize", you've been told to take on faith that you get a "standard normal" variable Z--the same distribution but with 0, 1 as parameters.  Now we'll prove it.

How do we find the distribution of a function of a continuous random variable? (Sec. 4.6)

Y = g(X). Two methods:

1) CDF method: Find FY(y0) = P(Y < y0) = P(g(X) < y0)
Do what you need to, to rewrite this as P(X <? g-1?(y0) ) (If g doesn't have an inverse over the whole range this could be more complicated. And the inequality might reverse.)

This can (usually) be written as an integral with upper limit x= g-1 (y0), f(x)dx. Change variables to y, inside the integral and in the limits, until you get an integral with y0 as the upper limit. That's your FY(y0).
Example: Normal distribution, standardize it. .= g(X).
 FZ(z0) = P(Z < z0) = P(g(X) < z0) = P( ) = P() =

   Change variables to  and you'll get  We recognize this is the CDF of the standard normal.  You can take the derivative to get the pdf.

2) Density method: Suppose X has density fX(x). We want the density of Y, fY(y). A thin rectangle at xo of width dx and height fX(xo) gets mapped by y = g(x)into a rectangle at yo of width dy and height fY(yo). (If g is one-to-one, that's all that gets mapped into the new rectangle.)

Set the rectangles equal: fY(yo)dy = fX(xo) dx

Solve for fY(yo):
              fY(yo) = fX(xo) dx/dy. You have to find dx/dy, and change variables to get everything in terms of y.
Since y = g(x), dy/dx = g'(x). You can put it under 1 to get dx/dy. Or find x = g-1 (y) and take the derivative of x with respect to y.

Example: Linear transformations first. Normal.
Let X be .  Let .  A rectangle at xo of width dx and height fX(xo) gets mapped by into a rectangle at zo of width dz and height fZ(zo). fZ(zo)dz = fX(xo) dx      fZ(zo) = fX(xo) dx/dz.

, so .

fX(xo) = , rewrite in terms of z, not hard.
Then fX(xo) dx/dz =, the pdf for Z !

Additional techniques:
(a)  CDF Method:  If the CDF has a nice formula FX(x), you can just substitute into FX(x) from P(X < something ), getting FX(something). If FX does not have the same formula everywhere, you'll have to follow each piecewise formula in the transformation.

(b) CDF method:  If you end up with an integral, instead of changing variables inside the integral,  you can find the density by taking the derivative.  This requires a more sophisticated use of the Fundamental Theorem of the Calculus, with the chain rule, since the upper limit of the integral will not be simple.

Practice here: http://www.math.temple.edu/~cow/
Calculus Book II > 1.Integration > 4.Fundamental Theorem > 1.Differentiation and the Fundamental Theorem
Read the Help, then try problems. Try these: 1, 2, 8, then any you like.  Type x3 as x^3. sqrt(x) or x^.5.
(c) Density method--If the transformation function "flips" things (x1< x2 , but y1 > y2), dy/dx will be negative.  Just throw away the negative (take the absolute value.)

(d) If the transformation function is not 1-1, you'll have two (or more?) x-chunks mapping onto the same y.  Be sure to account for both.

HW: Normal review, Ash p. 136 Note, Ash uses X*, instead of the more conventional Z, for std. normal.
Do the bold ones, read the others to make sure you can do them.
# 1, 2, 3, 4, 5, 6, 7, 9, 17, 14(X = # who show: It's binomial--use Normal approx. to binomial.)

A) Practice the Fundamental Theorem using the Temple U. COW system above, doing 1, 2, 8. "Help" will explain the topic well.  You should check your "fundamental theorem" answers to these 3 by doing the integration shown, then taking the derivative.  (For instance, in 8, integrating 2t gives the indefinite integral t2.  Evaluating between x2 and 2 gives (x2 )2 - (2)2= x4  - 4.  Taking the derivative gives 4x3.  Your "fundamental theorem" answer must match this one.
Then do #10.  Then read the Help to know how to attack problems where both limits of integration move, and do #9.  Practice a few more, from anywhere in the list of 20.
Do the above problems for Day29 homework.
  I suggest printing out this page for reference for the next lecture(s)
4.6: Read as much of 4.6 as you can stand. I have misprints: p.142, bottom formulas are fY and FY,  not X as I have.  The axes in the graphs on p. 148 are rubbery.
B)  X is uniform on [0, 1].  Let Y = cX + d. (assume c > 0)   a) On what interval does Y have positive probability?
b)  Write FY(y0) = P(Y < y0), substitute Y = cX + d, and solve to get P(X.....)  Find a formula for this, using the CDF for X. (no integrals.  Use technique (a) above.
c) Take the derivative with respect to y, thus finding fY(yo).  For F and f, be sure to give the intervals where your functions are in force.
d) Now use the Density method to go straight from fX(x)  to fY(y).  Sketch it for c = 2, d = 3, and show where x= 0, .5, and 1 map to.

C)  Repeat B,  only assume c < 0  (so values flip direction.)  For  c = -2, d = 3,  show where x= 0, .5, and 1 map to, before you begin.  Use technique (c) above.

D) Applying technique (c) above: After the fundamental theorem --chain rule version--is understood, look at the derivation above for standardizing X normal.  We got FZ(zo) = .  Take the derivative of this with respect to z (drop the 0-subscript to do so).
You should get, after everything is done, the formula for fZ(z).


Sievers home
 Math300-Sp04/Dayp29.htm 
10pm
4/13/04

This page belongs to Sally Sievers who is solely responsible for its content. Please see our statement of responsibility.