Then the distribution of X is called the standard normal distribution, denoted by N(0, 1).

Transcription

1 Normal distributions Normal distributions are commonly used because distributions of averages of IID random variables can often be approximated well by normal distributions The standard normal distribution ( 標準常態分布 ) N(0, 1) Suppose that X has the following PDF: f X (x) = 1 2π e x2 2, < x < Then the distribution of X is called the standard normal distribution, denoted by N(0, 1) Normal distribution N(µ, σ 2 ) Suppose that µ and σ > 0 are constants, and X is a random variable such that X µ σ N(0, 1) The distribution of X is called the normal distribution with mean µ and variance σ 2 (see Page 4 of this handout for mean and variance of N(µ, σ 2 )), denoted by N(µ, σ 2 ) X has the following PDF: f X (x) = µ)2 1 e (x2σ 2, < x < 2πσ 2 ( a µ P (a < N(µ, σ 2 ) < b) is P < N(0, 1) < b µ ) The notation σ σ P (a < N(µ, σ 2 ) < b) means P (a < X < b), where X N(µ, σ 2 ) P (a < N(0, 1) < b) can be found using software that computes or using normal tables b a 1 2π e x2 /2 dx R codes for computing P (0 < N(0, 1) < 15) Use integrate integrate(f, 0, 15) computes 15 0 f(x)dx f <- function(x){ exp(-x^2/2)/sqrt(2*pi) } integrate(f, 0, 15) 1

2 Use pnorm pnorm(a) gives P (N(0, 1) a) pnorm(15)-pnorm(0) gives P (0 < N(0, 1) 15) = P (0 < N(0, 1) < 15) Example 1 Use the table Normal probabilities at to find P (0 < N(0, 1) < 15) and P (0 < N(0, 1) < 151) z P (0 < N(0, 1) < 15) = P (0 < N(0, 1) < 151) = P (0 < N(0, 1) < ( )) = The N(0, 1) PDF is symmetric about 0 Therefore, for a > 0, P ( a < N(0, 1) < 0) = P (0 < N(0, 1) < a) and P (N(0, 1) < a) = P (N(0, 1) > a) (a) (b) Figure 1: P ( 15 < N(0, 1) < 0) = P (0 < N(0, 1) < 15) =

3 Example 2 Find P ( 15 < N(0, 1) < 13) Sol P ( 15 < N(0, 1) < 13) = P ( 15 < N(0, 1) < 0) + P (0 N(0, 1) < 13) = P (0 < N(0, 1) < 15) + P (0 < N(0, 1) < 13) = = Example 3 Sol Find P (N(0, 1) > 13) and P (N(0, 1) < 15) P (N(0, 1) > 13) = P (N(0, 1) > 0) P (0 < N(0, 1) 13) = P (N(0, 1) > 0) P (0 < N(0, 1) < 13) = = 00968, and P (N(0, 1) < 15) = P (N(0, 1) > 15) = P (N(0, 1) > 0) P (0 < N(0, 1) 15) = P (N(0, 1) > 0) P (0 < N(0, 1) < 15) = = Example 4 Suppose that X N(1000, (100) 2 ) (a) Find P (1000 < X < 1100) (b) Find P (X < 790) Sol (a) (b) P (1000 < X < 1100) = ( P < X 1000 < = P (0 < N(0, 1) < 1) = P (X < 790) = ( ) X P < = P (N(0, 1) < 21) = P (N(01) > 21) = P (N(0, 1) > 0) P (0 < N(0, 1) < 21) = = )

4 Mean and variance for N(µ, σ 2 ) Suppose that X N(µ, σ 2 ) Then E(X) = µ and V ar(x) = σ 2 (1) One may derive (1) using the following facts 1 (X µ)/σ N(0, 1) 2 The mean and variance for N(0, 1) are 0 and 1 respectively 3 For any constants a and b, E(a + bx) = a + be(x) and V ar(a + bx) = b 2 V ar(x) (2) The proof of the second equality in (2) is left as an exercise (Problem 14) The mean and variance for N(0, 1) can be computed using R and the following formulas: and V ar(x) = where f is the N(0, 1) PDF R codes E(X) = xf(x)dx (x E(X)) 2 f(x)dx, Define g(x) = xf(x), where f is the N(0, 1) PDF, and then compute g(x)dx, which is the mean for N(0, 1) g <- function(x){ x*exp(-x^2/2)/sqrt(2*pi) } integrate(g, -Inf, Inf) The result shows that the mean for N(0, 1) is 0 Define g(x) = x 2 f(x), where f is the N(0, 1) PDF, and then compute g(x)dx, which is the variance for N(0, 1) g <- function(x){ (x^2)*exp(-x^2/2)/sqrt(2*pi) } integrate(g, -Inf, Inf) The result shows that the variance for N(0, 1) is 1 Shift foremen income example (Page 231 in the 14th Ed or Page 235 in the 15th Ed) The weekly incomes of shift foremen in the glass industry follow the normal probability distribution with a mean of $1,000 and a standard deviation of $100 4

5 (a) What is the probability of selecting a shift foreman in the glass industry whose income is between $1,000 and $1,100? (b) What is the probability of selecting a shift foreman in the glass industry whose income is less than $790? Sol See the solution to Example 4 The empirical rule Suppose that X N(µ, σ 2 ), σ > 0 Then = 068 if k = 1; P (µ kσ < X < µ + kσ) = 095 if k = 2; = 1 if k = 3 Recall that for a random variable X such that E(X) = µ and V ar(x) = σ, Chebyshev s theorem states that P (µ kσ < X < µ + kσ) if k = 1; k 2 = 075 if k = 2; (3) = 089 if k = 3 If X N(µ, σ 2 ), then from (1), (3) holds However, the guaranteed coverage probabilities from (3) are not very close to the exact coverage probability (see the empirical rule for probability comparison) Example in the text (Page 228 in the 14th Ed or Page 232 in the 15th Ed) As part of its quality assurance program, the Autolite Battery Company conducts tests on battery life For a particular D-cell alkaline battery, the mean life is 19 hours The useful life of the battery follows a normal distribution with a standard deviation of 12 hours Answer the following questions (a) About 95 percent of the batteries failed between what two values? (b) Virtually all of the batteries failed between what two values? Ans About 95 percent of the battery failed between 166 hours and 214 hours (19±2(12)) Virtually all of the batteries failed between 154 hours and 226 hours (19 ± 3(12)) The six-sigma problem Suppose that in the manufacturing process of some product, the measurement of a quality characteristic can be viewed as a random variable with mean µ and standard deviation σ Then we can use µ ± kσ to indicate a reasonable range for the measurement Suppose that k = 6 What can we say about the coverage probability for the 5

6 range µ ± 6σ? What can we say about the coverage probability if the measurement is normally distributed? Ans The coverage probability is at least 1 1/ If the measurement is normally distributed, the coverage probability is P ( 6 < N(0, 1) < 6), which is very close to 1 Suppose that X N(µ, σ 2 ), then a + bx N(E(a + bx), V ar(a + bx)) Example 5 Suppose that X N(1, 4) and Y = 5 2X Find P (3 < Y < 9) Sol E(Y ) = 5 2E(X) = 3 and V ar(y ) = 4V ar(x) = 16, so ( P (3 < Y < 9) = P 0 < Y 3 ) < 15 = P (0 < N(0, 1) < 15) =