11 Normal Distributions
- Normal distributions are probably the most important distributions in probability and statistics.
- A continuous random variable \(X\) has a Normal (a.k.a., Gaussian) distribution with mean \(\mu\in (-\infty,\infty)\) and standard deviation \(\sigma>0\) if its pdf is \[ f_X(x) = \frac{1}{\sigma\sqrt{2\pi}}\,\exp\left(-\frac{1}{2}\left(\frac{x-\mu}{\sigma}\right)^2\right), \quad -\infty<x<\infty \]
- If \(X\) has a Normal(\(\mu\), \(\sigma\)) distribution then \[\begin{align*} \text{E}(X) & = \mu\\ \text{SD}(X) & = \sigma \end{align*}\]
- In R:
rnorm(N_rep, mean, sd)
to simulate valuesdnorm(x, mean, sd)
to compute the probability density functionpnorm(x, mean, sd)
to compute the cumulative distribution function \(\text{P}(X \le x)\).qnorm(p, mean, sd)
to compute the quantile function which returns \(x\) for which \(\text{P}(X\le x) = p\).
- A Normal density is a particular “bell-shaped” curve which is symmetric about its mean \(\mu\). The mean \(\mu\) is a location parameter: \(\mu\) indicates where the center and peak of the distribution is.
- The standard deviation \(\sigma\) is a scale parameter: \(\sigma\) indicates the distance from the mean to where the concavity of the density changes. That is, there are inflection points at \(\mu\pm \sigma\).
Example 11.1 The pdfs in the plot below represent the distribution of hypothetical test scores in three classes. The test scores in each class follow a Normal distribution. Identify the mean and standard deviation for each class.
- Any Normal distribution follows the “empirical rule” which determines the percentiles that give a Normal distribution its particular bell shape.
- For a Normal distribution:
- 38% of values are within 0.5 standard deviations of the mean
- 68% of values are within 1 standard deviation of the mean
- 87% of values are within 0.5 standard deviations of the mean
- 95% of values are within 2 standard deviations of the mean
- 99% of values are within 2.6 standard deviations of the mean
- 99.7% of values are within 3 standard deviations of the mean
- The table below lists some percentiles of a Normal distribution.
Percentile | SDs away from the mean |
---|---|
0.1% | 3.09 SDs below the mean |
0.5% | 2.58 SDs below the mean |
1% | 2.33 SDs below the mean |
2.5% | 1.96 SDs below the mean |
10% | 1.28 SDs below the mean |
15.9% | 1 SDs below the mean |
25% | 0.67 SDs below the mean |
30.9% | 0.5 SDs below the mean |
50% | 0 SDs above the mean |
69.1% | 0.5 SDs above the mean |
75% | 0.67 SDs above the mean |
84.1% | 1 SDs above the mean |
90% | 1.28 SDs above the mean |
97.5% | 1.96 SDs above the mean |
99% | 2.33 SDs above the mean |
99.5% | 2.58 SDs above the mean |
99.9% | 3.09 SDs above the mean |
- The Normal(0, 1) distribution (mean 0, SD 1) is called the “Standard Normal” distribution
- Standardization:
- If \(X\) has a Normal(\(\mu\), \(\sigma\)) distribution then \(Z = \frac{X-\mu}{\sigma}\) has a Normal(0, 1) distribution
- If \(Z\) has a Normal(0, 1) distribution then \(X = \mu + \sigma Z\) has a Normal(\(\mu\), \(\sigma\)) distribution
Example 11.2 The wrapper of a package of candy lists a weight of 47.9 grams. Naturally, the weights of individual packages vary somewhat. Suppose package weights have an approximate Normal distribution with a mean of 49.8 grams and a standard deviation of 1.3 grams.
- Sketch the distribution of package weights. Carefully label the variable axis. It is helpful to draw two axes: one in the measurement units of the variable, and one in standardized units.
- Why wouldn’t the company print the mean weight of 49.8 grams as the weight on the package?
- Estimate the probability that a package weighs less than the printed weight of 47.9 grams.
- Estimate the probability that a package weighs between 47.9 and 53.0 grams.
- Suppose that the company only wants 1% of packages to be underweight. Find the weight that must be printed on the packages.
- Find the 25th percentile (a.k.a., first (lower) quartile) of package weights.
- Find the 75th percentile (a.k.a., third (upper) quartile) of package weights. How can you use the work you did in the previous part?
= 10000
N_rep
= rnorm(N_rep, mean = 49.8, sd = 1.3)
x
head(x) |>
kbl()
x |
---|
49.14364 |
50.99079 |
49.90919 |
48.84746 |
51.00235 |
50.64904 |
hist(x,
freq = FALSE,
main = "")
sum(x < 47.9) / N_rep
[1] 0.0731
pnorm(47.9, 49.8, 1.3)
[1] 0.07193386
pnorm((47.9 - 49.8) / 1.3)
[1] 0.07193386
sum((x > 47.9) & (x < 53)) / N_rep
[1] 0.9202
pnorm(53, 49.8, 1.3) - pnorm(47.9, 49.8, 1.3)
[1] 0.921149
quantile(x, 0.01)
1%
46.84973
qnorm(0.01, 49.8, 1.3)
[1] 46.77575
qnorm(0.01)
[1] -2.326348
49.8 + 1.3 * qnorm(0.01)
[1] 46.77575
Example 11.3 Daily high temperatures (degrees Fahrenheit) in San Luis Obispo in August follow (approximately) A Normal distribution with a mean of 76.9 degrees F. The temperature exceeds 100 degrees Fahrenheit on about 1.5% of August days.
- What is the standard deviation?
- Suppose the mean increases by 2 degrees Fahrenheit. On what percentage of August days will the daily high temperature exceed 100 degrees Fahrenheit? (Assume the standard deviation does not change.)
- A mean of 78.9 is 1.02 times greater than a mean of 76.9. By what (multiplicative) factor has the percentage of 100-degree days increased? What do you notice?
Example 11.4 In a large class, scores on midterm 1 follow (approximately) a Normal\((\mu_1, \sigma)\) distribution and scores on midterm 2 follow (approximately) a Normal\((\mu_2, \sigma)\) distribution. Note that the SD \(\sigma\) is the same on both exams. The 40th percentile of midterm 1 scores is equal to the 70th percentile of midterm 2 scores. Compute
\[ \frac{\mu_1-\mu_2}{\sigma} \]
(This is one statistical measure of effect size.)
Example 11.5 A bank uses an applicant’s score on some criteria to decide whether or not to approve their loan application. Based on past history, the bank has determined that
- Scores for those who repay the loan follow a Normal distribution with mean 60 and SD 8
- Scores for those who do NOT repay the loan follow a Normal distribution with mean 40 and SD 12
When someone applies for a loan the bank does not know whether the applicant will eventually repay the loan. How can the bank use the applicant’s score to decide whether or not to approve the loan?
- Draw sketches of these two normal curves on the same axis.
- Suggest a decision rule, based on an applicant’s score, for deciding whether or not to approve the applicant’s loan.
- Describe the two kinds of classification errors that could be made in this situation.
- Determine the probability that we incorrectly reject the application of someone who would repay. Shade in the corresponding region under the Normal curve, and interpret this probability.
- Determine the probability that we incorrectly approve the application of someone who would NOT repay. Shade in the corresponding region under the Normal curve, and interpret this probability.
- In which direction — smaller or larger — would you need to change the decision rule’s cutoff value in order to decrease the probability that an applicant who would repay the loan is incorrectly rejected?
- Would the probability of the other kind of error — incorrectly approving a loan for an applicant who would not repay it — increase or decrease with this new cutoff value?
- Determine the cutoff value needed to decrease the probability that an applicant who would repay the loan is incorrectly rejected to 0.05.
- Determine the other error probability with this new cut-off rule.
- Now repeat the two previous parts with the goal of decreasing to 0.05 the probability of incorrectly approving an applicant who would not repay the loan.
- If you consider the two kinds of errors to be equally serious, how might you decide which of the three decision rules considered so far is the best?
- Which error do you think the bank would consider more serious? In which direction would that move the threshold?
- Making decisions/predictions based on data often involves trade-offs.
- Decreasing the probability of one kind of error often comes at the expense of increasing the probability of another kind of error.