Lecture 6
Lecture 6
probability distributions:
μ X
The Normal Distribution:
as mathematical function
(pdf)
1 x 2
1 (
)
f ( x) e 2
2
This is a bell shaped
Note constants: curve with different
=3.14159 centers and spreads
e=2.71828 depending on and
The Normal PDF
It’s a probability function, so no matter what the
values of and , must integrate to 1!
1 x 2
1 ( )
2
e 2 dx 1
Normal distribution is
defined by its mean and
standard dev.
1 x 2
E(X)= = 1 ( )
x 2
e 2 dx
1 x 2
1 ( )
Var(X)=2 = (
x2
2
e 2 dx) 2
Standard Deviation(X)=
**The beauty of the normal curve:
68% of
the data
2 1 x 2
1 ( )
2 2
e 2 dx .95
3 1 x 2
1 ( )
3 2
e 2 dx .997
How good is rule for real data?
25
20
P
e 15
r
c
e
n 10
t
0
80 90 100 110 120 130 140 150 160
POUNDS
95% of 120 = .95 x 120 = ~ 114 runners
In fact, 115 runners fall within 2-SD’s of the mean.
25
20
P
e 15
r
c
e
n 10
t
0
80 90 100 110 120 130 140 150 160
POUNDS
99.7% of 120 = .997 x 120 = 119.6 runners
In fact, all 120 runners fall within 3-SD’s of the
mean.
25
20
P
e 15
r
c
e
n 10
t
0
80 90 100 110 120 130 140 150 160
POUNDS
Example
Suppose SAT scores roughly follows a normal
distribution in the U.S. population of college-
bound students (with range restricted to 200-
800), and the average math SAT is 500 with a
standard deviation of 50, then:
68% of students will have scores between 450 and
550
95% will be between 400 and 600
99.7% will be between 350 and 650
Example
BUT…
What if you wanted to know the math SAT
1 Z 0 2 1
1 ( ) 1 ( Z )2
p( Z ) e 2 1
e 2
(1) 2 2
The Standard Normal Distribution (Z)
All normal distributions can be converted into
the standard normal curve by subtracting the
mean and dividing by the standard deviation:
X
Z
0 2.0 Z ( = 0, = 1)
Example
For example: What’s the probability of getting a math SAT score of 575 or less, =500 and =50?
575 500
Z 1.5
50
i.e., A score of 575 is 1.5 standard deviations above the mean
Yikes!
But to look up Z= 1.5 in standard normal chart (or enter
into SAS) no problem! = .9332
Practice problem
If birth weights in a population are
normally distributed with a mean of 109
oz and a standard deviation of 13 oz,
a. What is the chance of obtaining a birth
weight of 141 oz or heavier when
sampling birth records at random?
b. What is the chance of obtaining a birth
weight of 120 or lighter?
Answer
a. What is the chance of obtaining a birth
weight of 141 oz or heavier when
sampling birth records at random?
141 109
Z 2.46
13
From the chart or SAS Z of .85 corresponds to a left tail area of:
P(Z≤.85) = .8023= 80.23%
Looking up probabilities in
the standard normal table
What is the area to
the left of Z=1.51 in
a standard normal
curve?
Area is 93.45%
Z=1.51
Z=1.51
Normal probabilities in
SAS
data _null_;
The “probnorm(Z)” function gives
theArea=probnorm(1.5); you the probability from negative
put theArea; infinity to Z (here 1.5) in a
run; standard normal curve.
0.9331927987
And if you wanted to go the other direction (i.e., from the area to the Z score
(called the so-called “Probit” function
data _null_;
The “probit(p)” function gives you
theZValue=probit(.93); the Z-value that corresponds to a
put theZValue; left-tail area of p (here .93) from a
run; standard normal curve. The probit
1.4757910282 function is also known as the
inverse standard normal function.
Probit function: the
inverse
(area)= Z: gives the Z-value that goes with the probability you want
For example, recall SAT math scores example. What’s the score that corresponds to the
90th percentile?
In Table, find the Z-value that corresponds to area of .90 Z= 1.28
Or use SAS
data _null_;
theZValue=probit(.90);
put theZValue;
run;
1.2815515655
If Z=1.28, convert back to raw SAT score
1.28 =
X – 500 =1.28 (50)
X=1.28(50) X+ 500
500 = 564 (1.28 standard deviations above the mean!)
50
`
Are my data “normal”?
Not all continuous random
variables are normally distributed!!
It is important to evaluate how well
the data are approximated by a
normal distribution
Are my data normally
distributed?
1. Look at the histogram! Does it appear bell
shaped?
2. Compute descriptive summary measures—
are mean, median, and mode similar?
3. Do 2/3 of observations lie within 1 std dev of
the mean? Do 95% of observations lie within
2 std dev of the mean?
4. Look at a normal probability plot—is it
approximately linear?
5. Run tests of normality (such as Kolmogorov-
Smirnov). But, be cautious, highly influenced
by sample size!
Data from our class…
Median = 6
Mean = 7.1
Mode = 0
SD = 6.8
Range = 0 to 24
(= 3.5 σ)
Data from our class…
Median = 5
Mean = 5.4
Mode = none
SD = 1.8
Range = 2 to 9
(~ 4 σ
Data from our class…
Median = 3
Mean = 3.4
Mode = 3
SD = 2.5
Range = 0 to 12
(~ 5 σ
Data from our class…
Median = 7:00
Mean = 7:04
Mode = 7:00
SD = :55
Range = 5:30 to
9:00
(~4 σ
Data from our class…
3.6 7.2
Data from our class…
5:14
8:54 7:04+/-
2*0:55 =
5:14 – 8:54
Data from our class…
4:19
9:49 7:04+/-
2*0:55 =
4:19 – 9:49
The Normal Probability
Plot
Normal probability plot
Order the data.
Find corresponding standardized normal
quantileth values: i
i quantile ( )
n 1
where is the probit function, which gives the Z value
that corresponds to a particular left - tail area
Plot the observed data values against normal
quantile values.
Evaluate the plot for evidence of linearity.
Normal probability plot
coffee…
Right-Skewed!
(concave up)
Normal probability plot
love of writing…
Neither right-
skewed or left-
skewed, but big
gap at 6.
Norm prob. plot Exercise…
Right-Skewed!
(concave up)
Norm prob. plot Wake up
time
Closest to a
straight line…
Formal tests for normality
Results:
Coffee: Strong evidence of non-
normality (p<.01)
Writing love: Moderate evidence of non-
normality (p=.01)
Exercise: Weak to no evidence of non-
normality (p>.10)
Wakeup time: No evidence of non-
normality (p>.25)
Normal approximation to
the binomial
When you have a binomial distribution where n
is large and p is middle-of-the road (not too
small, not too big, closer to .5), then the
binomial starts to look like a normal
distribution in fact, this doesn’t even take a
particularly large n
Recall: smoking
.27
example… Starting to have a normal
shape even with fairly small
n. You can imagine that if n
got larger, the bars would get
thinner and thinner and this
would look more and more
0 1 2 3 4 5 6 7 8 like a continuous function,
with a bell curve shape. Here
np=4.8.
Normal approximation to
binomial
.27
0 1 2 3 4 5 6 7 8
OR Use SAS:
data _null_;
Cohort=cdf('binomial', 120, .25, 500);
put Cohort;
run;
0.323504227
x np(1 p)
Differs
by a
factor
pˆ p of n.
For proportion:
np(1 p ) p (1 p )
pˆ 2 2
n n
P-hat stands for “sample p (1 p )
proportion.” pˆ
n
It all comes back to Z…
Statistics for proportions are based
on a normal distribution, because
the binomial can be approximated
as normal if np>5