advanced32 min studyLesson 705 in Mathematics

Local asymptotic normality

Approximate a regular experiment through local likelihood ratios.

On this page

What you’ll understand

  • State a differentiability-in-quadratic-mean route to iid LAN.
  • Identify central sequences and information in a local likelihood expansion.
  • Separate local parameter shifts from fixed alternatives and estimator normality.

Asymptotic normality describes one estimator. Local asymptotic normality describes how an entire statistical experiment changes near a fixed parameter through its likelihood ratios. The relevant alternatives approach the reference parameter at inverse-root-sample-size distance. Under regularity, their log likelihood ratios resemble those of a Gaussian shift experiment. This supplies a shared foundation for local testing, efficiency bounds, and changes of measure. The approximation is local and probabilistic; it is not a claim that the original observations become normally distributed or that every estimator is efficient.

A quadratic-mean derivative of the model

Consider an iid family dominated by a common measure, with an interior reference parameter. For simplicity, assume its densities are positive on the same measurable set throughout a parameter neighborhood, ensuring local mutual absolute continuity. Differentiability in quadratic mean means that the square root of the density admits an integrated squared-error linear approximation, with remainder smaller than the squared parameter displacement. Its derivative is one half of the score times the square-root density. Require a mean-zero score with finite second moment and, for the nondegenerate version used here, a positive-definite information matrix. This is a model condition, distinct from smoothness of a numerical objective.

The information matrix is the score covariance. Summing scores and dividing by the square root of sample size gives the central sequence, whose limit is normal by the multivariate central limit theorem. For this iid dominated quadratic-mean differentiable model, the local log likelihood ratio has a linear central-sequence term and a deterministic quadratic correction. The statement holds at each fixed local displacement; stronger uniform assertions over a growing parameter region would require additional work.

log⁡dPθ0+h/n ndPθ0 n=hTΔn−12hTIh+oPθ0n(1),Δn=1n∑iℓ˙θ0(Xi)⇒N(0,I).\log\frac{dP_{\theta_0+h/\sqrt n}^{\,n}}{dP_{\theta_0}^{\,n}}=h^{\mathsf T}\Delta_n-\tfrac12h^{\mathsf T}Ih+o_{P_{\theta_0}^{n}}(1),\qquad \Delta_n=\frac1{\sqrt n}\sum_i\dot\ell_{\theta_0}(X_i)\Rightarrow N(0,I).
LAN expands likelihood ratios under the reference product measure, with h fixed.↔ Scroll wide equations horizontally.

The Gaussian shift is a local experiment

The quadratic correction centers the exponential of the Gaussian limit at mean one. It is not an arbitrary penalty inserted into a fitting procedure. Under the local alternative, a suitable change-of-measure argument shifts the central sequence's limiting mean by information times the local displacement, while preserving its covariance. The next lesson proves the relevant transfer principle through contiguity and exponential tilting. Keeping the probability measure explicit avoids confusing the reference distribution with the local-alternative distribution.

A fixed alternative stays a positive distance from the reference parameter and generally lies outside this expansion. Likewise, a normal limit for a maximum-likelihood estimator does not establish quadratic-mean differentiability or LAN for the underlying model. Parameter-dependent support can introduce different scales: estimating a uniform endpoint uses sample extremes and an inverse-sample-size rate. One must examine the model's likelihood geometry rather than assuming that every familiar estimator supplies the same local Gaussian experiment.

Δn under Pθ0+h/n n ⇒ N(Ih,I),n(θ^n−θ0)≈I−1Δn for an efficient estimator.\Delta_n\ \text{under }P_{\theta_0+h/\sqrt n}^{\,n}\ \Rightarrow\ N(Ih,I),\qquad \sqrt n(\widehat\theta_n-\theta_0)\approx I^{-1}\Delta_n\ \text{for an efficient estimator}.
The shifted central sequence describes local information; the estimator approximation needs its own efficiency conditions.↔ Scroll wide equations horizontally.

What the local scale buys

An inverse-root-sample-size change is small enough to remain statistically comparable to the reference law, yet large enough to produce nontrivial local power in regular models. Smaller changes are often asymptotically invisible at first order; fixed changes can become perfectly distinguishable. These statements concern regular models at the chosen point, not a universal hierarchy for singular or boundary parameters. Information determines the meaningful local distance, so a displacement should be interpreted together with the model's score variance.

For calculations, first identify the per-observation score and information, then write the central sequence with all variance factors intact. Expand the actual likelihood ratio to verify the sign of the quadratic term. In a normal location model this expansion is exact, which provides a useful calibration for the asymptotic notation. For a Bernoulli family at an interior probability the expansion is approximate, while approaching a degenerate probability may invalidate the fixed-interior argument altogether.

An exact local normal likelihood ratio

Observe n=100n=100 independent N(θ,9)N(\theta,9) variables. Compare θ0=2\theta_0=2 with θ0+3/n=2.3\theta_0+3/\sqrt n=2.3. If the observed mean is 2.62.6, calculate the central sequence and log likelihood ratio.

  1. The score at the reference point is (X−2)/9(X-2)/9 and information is I=1/9I=1/9. Thus Δn=n(X‾−2)/9=2/3\Delta_n=\sqrt n(\overline X-2)/9=2/3 for the observed mean.

  2. With local displacement h=3h=3, the normal log likelihood ratio is exactly hΔn−h2I/2=2−1/2=1.5h\Delta_n-h^2I/2=2-1/2=1.5. The likelihood ratio is therefore e1.5e^{1.5}.

  3. Under the reference law, Δn\Delta_n is exactly normal with variance 1/91/9, so the log likelihood ratio is N(−1/2,1)N(-1/2,1). Under the alternative, n(X‾−2)/3\sqrt n(\overline X-2)/3 is exactly N(1,1)N(1,1).

The result

The same likelihood ratio has different laws under the two product measures. The Gaussian-shift interpretation is exact in this model and local in the general theorem.

Explore the model

Interactive lab

Interactive lab · sampling

Inspect an interior Bernoulli local displacement

Pause/Reset and set BOTH controls to p=0.25p=0.25, n=100n=100. Predict mean 0.250.25, variance 0.0018750.001875. For a fresh comparison Pause/Reset and set BOTH to p=0.30p=0.30, n=100n=100: the mean shifts by 0.05=0.5/1000.05=0.5/\sqrt{100} and variance is 0.00210.0021. Changing either control retains the paused clock.

Inspect an interior Bernoulli local displacementExact binomial masses for the mean of n independent Bernoulli trials, compared with normal mass between adjacent half-count boundaries; normal tails outside the finite support remain outside the bars. Playback highlights outcomes, not random samples. The Berry–Esseen bound applies to standardized CDFs, not individual bars or continuity-corrected errors. At p = 0 or 1, normal comparison and that bound are disabled. Normal-CDF evaluation uses a numerical approximation with absolute error below 7.5×10⁻⁸ (at most 1.5×10⁻⁷ per interval, before rounding). Current readings: Mean sample proportion: 0.5; Sample-mean variance: 0.0125; Uniform CDF error bound: ≤ 0.1062. All parameters are adjustable with the labeled controls below.Exact distribution of the mean of 20 independent Bernoulli variables with success probability 0.5; playback highlights outcome 0.00.250.50.75100.0580.1150.1730.23Sample mean · k / nProbability massNormal bin masses are an approximationk = 0 · exact probability = 9.54e-7
Exact binomial massNormal interval massHighlighted outcome / mean
Mean sample proportion
0.5
Sample-mean variance
0.0125
Uniform CDF error bound
≤ 0.1062
t = 0 s

Exact binomial masses for the mean of n independent Bernoulli trials, compared with normal mass between adjacent half-count boundaries; normal tails outside the finite support remain outside the bars. Playback highlights outcomes, not random samples. The Berry–Esseen bound applies to standardized CDFs, not individual bars or continuity-corrected errors. At p = 0 or 1, normal comparison and that bound are disabled. Normal-CDF evaluation uses a numerical approximation with absolute error below 7.5×10⁻⁸ (at most 1.5×10⁻⁷ per interval, before rounding). Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.

Look forThe exact binomial laws illustrate one finite local displacement. They do not display likelihood ratios or prove LAN for a parameter family.

Check your understanding

Which object does LAN approximate under the reference measure?

Try it yourself

For a Bernoulli family at p0=1/4p_0=1/4, use local displacement h=1/2h=1/2. Calculate the limiting mean and variance of the log likelihood ratio under the reference law, and identify the local probability at sample size 100100.

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

Bernoulli information is 1/[p0(1−p0)]1/[p_0(1-p_0)]. The LAN log ratio has mean −h2I/2-h^2I/2 and variance h2Ih^2I.

Compare with the solution

Information is I=16/3I=16/3, so h2I=4/3h^2I=4/3. The limiting log likelihood ratio is normal with mean −2/3-2/3 and variance 4/34/3. Its exponent has expectation exp⁡(−2/3+(4/3)/2)=1\exp(-2/3+(4/3)/2)=1, as required for this Gaussian tilt. At sample size 100100, the local probability is 1/4+(1/2)/10=0.301/4+(1/2)/10=0.30. This is an asymptotic likelihood statement, unlike the exact normal-location calculation.

Take it with you

  • LAN describes local likelihood geometry under a named reference measure.
  • Quadratic-mean differentiability and finite score information provide a regular iid route.
  • A local Gaussian experiment does not establish global normality or automatic estimator efficiency.

Read further

Original sources and deeper explanations to follow at your own pace.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.