# Normal Models¶

The theoretical foundation of inference is to model our data with some distribution. By far the most common such distribution that arises in practice is the normal distribution.

## Normal curves¶

### The standard normal¶

The standard normal distribution is a specific bell shaped curve; its graph is shown at the top of the page above. It has a maximum as it crosses the $y$-axis; its mean is zero. It decreases as we move away from the mean and tapers off to the $x$-axis. It's scaled so that the total area under its graph is one. The area under the curve is distributed in a specific way; it's standard deviation is one.

### More normals¶

The standard normal is one of a specific family of bell shaped curves. We can change the mean by shifting the graph to the left or to the right. We can change the standard deviation by compressing or dilating along the $x$-axis while simultaneously doing the opposite along the $y$-axes. Several such curves are shown in the figure below.

Note that the mean and standard deviation are denoted $m$ and $s$ in the figure above. Often, we denote these with the Greek letters $\mu$ and $\sigma$.

### A formula¶

It's worth mentioning that there is a specific formula that generates the normal curves, namely $$f_{\mu,\sigma}(x) = \frac{1}{\sqrt{2\pi}\sigma} e^{-(x-\mu)^2/(2\sigma^2)}.$$

You can explore the graphs of these functions on Desmos.

### Modeling histograms with normal curves¶

Often, if we scale a histogram so that the total area of its rectangles is one, then the normal curve with the same mean and standard deviation will match the histogram fairly closely. An example is shown below.

We've got a dynamic illustration on our webpage.

## The 68-95-99.7 rule¶

There is a statistical rule of thumb called the 68-95-99.7 rule that states that

• 68% of the population lies within 1 standard deviation of the mean,
• 95% of the population lies within 2 standard deviations of the mean, and
• 99.7% of the population lies within 3 standard deviations of the mean.

For this to work, the property being measured should be normally distributed. If so, then the rule of thumb follows from the fact that this property holds for the standard normal:

### Example¶

Let's suppose that the mean life expectancy of a cat is 14 years with a standard deviation of 2.5 years. Assuming that the cats' life spans are normally distributed, is it reasonable to expect a cat to live to 22 years old?

Solution: Well, the $z$-score for a 22 year old kitty cat would be $$Z = \frac{22-14}{2.5} = 3.2.$$ As we know, only $0.3\%$ of cats live beyond a $z$-score of 3, so a 22 year old cat would be quite rare indeed.

## Standard deviation as a ruler¶

The mean and standard deviation can be used to standardize data. A principle tool to do so is the $z$-score.

### The heptathlon¶

I've got a CSV on my website that contains the results of the 2008 Olympic Heptathlon. Here are the first few rows:

In [4]:
import pandas as pd

Out[4]:
Name Country Score javelin eight_hundred shot_put two_hundred long_jump hundred_hurdles high_jump
0 Nataliia Dobrynska Ukraine 6733 48.60 137.72 17.29 24.39 6.63 13.44 1.80
1 Hyleas Fountain United States 6619 41.93 135.45 13.36 23.21 6.38 12.78 1.89
2 Tatyana Chernova Russia 6591 48.37 126.50 12.88 23.95 6.47 13.65 1.83
3 Kelly Sotherton Great Britain 6517 37.66 127.34 13.87 23.39 6.33 13.18 1.83
4 Jessica Zelinka Canada 6490 43.91 127.95 13.79 23.64 6.12 12.97 1.77

As it turns out, Nataliia Dobrynska won the long jump by about half a meter over the average jump; Hyleas Fountain won the 200m by a second and a half faster than the mean. Which of these two results deserves more points?

One way to compare these results involving different units is to use the standard deviation as a ruler. To do so, we standardize. Given a data point $X$ chosen from data that is normally distributed with mean $\mu$ and standard deviation $\sigma$ we compute $$Z = \frac{x-\mu}{\sigma}.$$ With this choice of $Z$, note that $$X = \mu + Z\sigma.$$ Thus, $Z$ measures exactly how many standard deviations $X$ is from the mean. This computation is sometimes called the $Z$-score.

Here's how we apply this to the 200m vs long jump question:

In [5]:
two_hundred_x = df.two_hundred.min()
two_hundred_mean = df.two_hundred.mean()
two_hundred_std = df.two_hundred.std()
(two_hundred_x-two_hundred_mean)/two_hundred_std

Out[5]:
-2.0357396116514703
In [6]:
long_jump_x = df.long_jump.max()
long_jump_mean = df.long_jump.mean()
long_jump_std = df.long_jump.std()
(long_jump_x-long_jump_mean)/long_jump_std

Out[6]:
2.2511914941134483

It looks like the long jump is the more impressive result.

### SAT Example¶

The SAT is designed to have mean score of 500 with a standard deviation of 100. Let's suppose you score a 650.

a. What is your $z$-score? b. What is your percentile, according to the normal model?

Solution: For part (a), the $z$-score is simply $$Z = \frac{650-500}{100} = 1.5.$$ We can answer part (b) by computing an area under the standard normal curve as shown below:

This type of area is typically computed using either software or a table. We will learn how to do this in Python soon but, for now, let's take a look at this table.

## Assessing normality¶

In the previous problem, we assumed that the lifespan of cats was normally distributed. Is this a valid assumption? I have no idea, but given a data set there are ways to check.

### A histogram¶

Human heights are classically known to be normal. Let's take a look at the heights of the over 9500 men chosen from our CDC data set.

In [8]:
import pandas as pd
df.height.hist(bins = 20, grid=False, edgecolor='black');


Looks kinda normal.

### The normal probability plot¶

There's a more sensitive tool, though: the normal probability plot - more generally called a quantile-quantile plot.

Here’s the basic idea: Suppose we have some sample data consisting of $n$ numerical values. We’ll make a plot with $n$ points - one point for each of the $n$ values. The vertical component for a given data point is just the value of the data point itself. The horizontal component for that given data point is determined by the $z$-score of that percentile.

For example, $z=1$ is greater than about 84% of all values produced by the normal distribution. You can look that up in a normal table or use the following Python command:

In [9]:
from scipy.stats import norm
norm.cdf(1)

Out[9]:
0.8413447460685429

Furthermore, the $84^{\text{th}}$ quantile for the men's CDC height data is 72 inches, as the following computation reveals:

In [10]:
df.height.quantile(0.84)

Out[10]:
72.0

Thus, the point $(1,72)$ should be on the normal probability plot for this data. Here's the normal probability plot of that height data, together with the point $(1,72)$ shown in yellow.

By contrast, income is classically not normally distributed. Here's a normal probability plot for a random sample of incomes taken from the US census: