Confidence Intervals

A confidence interval is a random interval produced by a procedure with a target long-run coverage rate. For an approximately normal estimator,

is the usual large-sample form. For a normal mean with unknown variance,

The central limit theorem supplies many standard errors; statistical estimation supplies the estimator; hypothesis testing often uses the same sampling distribution.

Worked simulation

This simulation repeatedly builds t confidence intervals for normal samples and measures empirical coverage, average width, and the critical value used.

import numpy as np
from scipy import stats
 
rng = np.random.default_rng(20260711)
reps, n, mu, sigma = 20000, 30, 5, 2
x = rng.normal(mu, sigma, size=(reps, n))
means = x.mean(axis=1)
s = x.std(axis=1, ddof=1)
tcrit = stats.t.ppf(.975, n - 1)
lo = means - tcrit * s / np.sqrt(n)
hi = means + tcrit * s / np.sqrt(n)
print("coverage", round(((lo <= mu) & (mu <= hi)).mean(), 4),
      "avg_width", round((hi - lo).mean(), 4),
      "tcrit", round(tcrit, 4))

Observed output:

coverage 0.9499 avg_width 1.4796 tcrit 2.0452

Across repeated samples, the 95 percent t-interval covers the fixed mean 0.9499 of the time, close to the nominal 0.95 target. The average width, 1.4796, reflects the sample size, noise scale, and critical value 2.0452.

QuantityMeaning in the simulation
coverageFraction of repeated intervals that contain the fixed true mean.
avg_widthTypical uncertainty width produced by the procedure.
tcritCritical value that widens the interval for finite samples and unknown variance.
One realized intervalA random output of the procedure; it either contains or it does not.

Caveats

One realized interval either contains the parameter or it does not. Coverage can fail under biased sampling, dependence, optional stopping, nonresponse, or variance formulas that ignore clustering. In experiments, pair intervals with effect size and statistical significance, not only a binary decision.

References