5.05a, 5.05b - Sample means and unbiased estimates

5.05a, 5.05b - Sample means and unbiased estimates

Averages are steadier than single observations, but averaging does not turn the original population into a normal distribution. In 5.05a you will distinguish the distribution of a sample mean from the population distribution; in 5.05b you will use sample data to estimate an unknown population mean and variance without systematic bias.

The centre and spread of a sample mean

A sensor has mean reading 100100 and standard deviation 1212. Can the mean of 3636 independent readings be six times as steady as one reading? The sensor itself does not improve. What changes is the random quantity you record: an average instead of one reading.

For example, averaging the readings 9797, 104104 and 102102 gives 101101. A new sample will usually give a different average. Before you take the sample, that average is a random variable in its own right.

You already know

For independent observations, expectations and variances add. Multiplying a random variable by aa multiplies its variance by a2a^2. A normal distribution is written N(μ,σ2)\mathrm{N}(\mu, \sigma^2), with variance in the second position.

Keep the random average and its observed value separate in your notation.

Sample mean and sampling distribution

For a randomly and independently selected sample X1,…,XnX_1,\ldots,X_n from one population, the sample mean is X‾=(X1+⋯+Xn)/n\overline X = (X_1 + \cdots + X_n)/n. Its sampling distribution describes the values of X‾\overline X across repeated samples of that same size. After observing values x1,…,xnx_1,\ldots,x_n, the calculated value is x‾\overline x.

Independent fluctuations can offset one another. Averaging therefore preserves the population centre while reducing spread. The following moment results apply whether the population is normal, skewed or discrete, provided its mean and finite variance exist.

Moments of a sample mean

For independent observations from a population with mean μ\mu and finite variance σ2\sigma^2,

E(X‾)=μ,Var(X‾)=σ2n\mathrm{E}(\overline X)=\mu, \qquad \mathrm{Var}(\overline X)=\frac{\sigma^2}{n}

These results hold for every positive sample size nn; normality is not required.

In the formula booklet.

The standard deviation of the sampling distribution has its own useful name. Take the square root of the variance, rather than dividing the population standard deviation by nn.

Standard error of the mean

Under the same independent-sampling conditions,

sd(X‾)=σn\mathrm{sd}(\overline X)=\frac{\sigma}{\sqrt n}

This is the standard error of the sample mean. It is the square root of the sample-mean variance printed in the booklet, rather than a separately printed formula.

Not in the formula booklet: learn it.

Calculate the variance first if there is any risk of confusing the two denominators.

Worked example 1: averaging thirty-six readings

Question. Independent readings have population mean 100100 and variance 144144. Find the expectation, variance and standard error of the mean of 3636 readings.

Plan. Keep the mean, divide the variance by the sample size, then take its square root.

E(X‾)=100same centreVar(X‾)=14436=4independent samplesd(X‾)=4=2square root\begin{aligned} \mathrm{E}(\overline X)&=100 &&\text{same centre}\\ \mathrm{Var}(\overline X)&=\frac{144}{36}=4 &&\text{independent sample}\\ \mathrm{sd}(\overline X)&=\sqrt4=2 &&\text{square root} \end{aligned}

Why this step? The spread needed is the spread of repeated averages, not the spread of individual readings.

Answer. The expectation is 100100, the variance is 44, and the standard error is 22.

Check. 12/36=212/\sqrt{36}=2. The standard error is one sixth of the individual standard deviation.

Worked example 2 changes only the sample size. Watch the difference between reducing the variance and reducing the standard error.

Worked example 2: quadrupling the sample size

Question. The same population has mean 100100 and variance 144144. Find the expectation, variance and standard error of the mean of 144144 independent readings.

Plan. Use the new sample size in the same sampling formulae.

E(X‾)=100unchanged centreVar(X‾)=144144=1divide by nsd(X‾)=1=1standard error\begin{aligned} \mathrm{E}(\overline X)&=100 &&\text{unchanged centre}\\ \mathrm{Var}(\overline X)&=\frac{144}{144}=1 &&\text{divide by }n\\ \mathrm{sd}(\overline X)&=\sqrt1=1 &&\text{standard error} \end{aligned}

Why this step? Four times as many readings divide the variance by 44, so they divide its square root by 22.

Answer. The expectation is 100100, the variance is 11, and the standard error is 11.

Check. Relative to 3636 readings, 144144 readings halve the standard error from 22 to 11. Twice the precision needs four times the data.

That square-root relationship answers the opening puzzle. It also explains why taking a few extra observations gives a modest improvement rather than a dramatic one.

Watch out

Wrong: sd(X‾)=σ/n\mathrm{sd}(\overline X)=\sigma/n. Right: Var(X‾)=σ2/n\mathrm{Var}(\overline X)=\sigma^2/n, so sd(X‾)=σ/n\mathrm{sd}(\overline X)=\sigma/\sqrt n. Repeating one reading nn times does not create an independent sample and does not reduce its uncertainty.

The next calculation gives you the parameters only. Do not invent a distributional shape from them.

Exact normality and the central limit theorem

The mean and variance do not tell you the shape of a distribution. Normality needs a separate reason: either the population is normal, or a sufficiently large sample makes a normal approximation appropriate.

If individual waiting times are exponentially distributed, large samples still contain short waits and occasional long waits. The distribution of individual waits stays skewed. Repeated averages are more balanced, because a sample can combine several short waits with a long one.

[DIAGRAM: asset_name: A skew population and the distributions of its averages; asset_slug: h245_5_05_clt_means; caption: Averaging changes the distribution of the mean while the original population remains skewed; description: Three exact density curves for an exponential variable with mean 4 and the averages of 4 and 36 independent copies, all on the same horizontal scale. Their centres are all 4, their standard deviations are 4, 2 and two thirds, and the final curve has a close normal approximation drawn as a dashed line.]
Averaging changes the distribution of the mean while the original population remains skewed

The picture shows two effects. The averages concentrate around the same centre, and their distribution becomes more nearly symmetric. Only the second effect needs the central limit theorem.

When a sample mean is normal

For a random independent sample of size nn:

  • If X∼N(μ,σ2)X\sim\mathrm{N}(\mu,\sigma^2), then X‾∼N(μ,σ2/n)\overline X\sim\mathrm{N}(\mu,\sigma^2/n) exactly, for every nn.
  • For any population with mean μ\mu and finite positive variance σ2\sigma^2, X‾\overline X is approximately normal when nn is sufficiently large. Use X‾≈N(μ,σ2/n)\overline X\approx\mathrm{N}(\mu,\sigma^2/n).

The second statement is the central limit theorem. A working guide is approximately n>25n>25; this is not a guarantee for every severely skewed population. The results, rather than their proofs, are required.

In the formula booklet.

State which reason you are using before calculating a probability. With a known normal population, even a small sample gives an exactly normal mean. If the population variance is zero, all observations and their average are constant; a normal approximation is unnecessary.

Worked example 3: an exactly normal average

Question. A random variable has distribution N(24,100)\mathrm{N}(24,100). For 100100 independent observations, find the probability that the sample mean exceeds 2626.

Plan. Write the exact distribution of the mean and standardise the threshold.

X‾∼N(24,100100)normal populationP(X‾>26)=P(Z>26−241)Z∼N(0,1)=P(Z>2)=0.022750…normal probability\begin{aligned} \overline X&\sim\mathrm{N}\left(24,\frac{100}{100}\right) &&\text{normal population}\\ \mathrm{P}(\overline X>26)&=\mathrm{P}\left(Z>\frac{26-24}{1}\right) &&Z\sim\mathrm{N}(0,1)\\ &=\mathrm{P}(Z>2)=0.022750\ldots &&\text{normal probability} \end{aligned}

Why this step? The standard deviation in the denominator is 11, the standard error, rather than 1010, the individual standard deviation.

Answer. The probability is 0.02280.0228, to 33 significant figures.

Check. The threshold is two standard errors above the mean, so a small upper-tail probability is sensible.

Worked example 4 changes only what is known about the population shape. The parameter calculations remain exact; the normal probability becomes an approximation.

Worked example 4: a large sample from a non-normal population

Question. A population has mean 2424 and variance 100100, but is not normally distributed. For a random independent sample of size 100100, estimate the probability that the sample mean exceeds 2626.

Plan. Use the central limit theorem for the sample mean, then use its exact mean and variance in the normal approximation.

E(X‾)=24Var(X‾)=1sampling momentsX‾≈N(24,1)large n=100P(X‾>26)≈P(Z>2)=0.022750…standardise the mean\begin{aligned} \mathrm{E}(\overline X)&=24 \quad \mathrm{Var}(\overline X)=1 &&\text{sampling moments}\\ \overline X&\approx\mathrm{N}(24,1) &&\text{large }n=100\\ \mathrm{P}(\overline X>26)&\approx\mathrm{P}(Z>2)=0.022750\ldots &&\text{standardise the mean} \end{aligned}

Why this step? The central limit theorem justifies the normal shape of X‾\overline X. It does not justify the variance formula, which already holds for independent samples of every size.

Answer. The estimated probability is 0.02280.0228.

Check. The same mean, variance and sample size give the same numerical normal calculation as before. Only its exact or approximate status has changed.

Continue this lesson.

You’re previewing NovaLearn without an account. Create a free account to read the rest.

Create free account