St · Statistics & Probability

What Is the Fisher Information of a Normal Distribution?

Published By the 9to5Science editorial team

At a glance

  1. 1. Choose parameters

    Decide whether the second parameter is the variance or the standard deviation; the information matrix depends on that choice.

  2. 2. Write the log-likelihood

    For one observation, express the normal log-density using the chosen parameters.

  3. 3. Compute the scores

    Differentiate with respect to each parameter and take the expected products of the scores.

  4. 4. Scale for the sample

    For n independent observations, multiply the one-observation information matrix by n.

  5. 5. Transform if needed

    Use the appropriate parameterization when reporting information for standard deviation instead of variance.

The Fisher information for a normal model depends on how you parameterize its spread. Using the mean μ\mu and variance v=σ2v=\sigma^2, the information matrix for nn independent observations is diagonal: its entries are n/vn/v and n/(2v2)n/(2v^2). The off-diagonal entries are zero.

Below is the derivation, along with the equivalent result when the parameter is standard deviation σ\sigma.

Fisher information for one observation

Let XX have a normal distribution with unknown mean μ\mu and unknown variance v>0v>0. Its log-density is

ℓ(μ,v;x)=−12log⁡(2πv)−(x−μ)22v.\ell(\mu,v;x)=-\frac{1}{2}\log(2\pi v)-\frac{(x-\mu)^2}{2v}.

The score is the vector of derivatives of the log-density with respect to the parameters:

∂ℓ∂μ=x−μv,∂ℓ∂v=−12v+(x−μ)22v2.\begin{aligned} \frac{\partial\ell}{\partial\mu}&=\frac{x-\mu}{v},\\ \frac{\partial\ell}{\partial v}&=-\frac{1}{2v}+\frac{(x-\mu)^2}{2v^2}. \end{aligned}

Fisher information is the expected outer product of the score with itself. Equivalently, under the usual regularity conditions, it is the negative expected Hessian of the log-likelihood. For this normal model, use the centered normal moments

E[X−μ]=0,E[(X−μ)2]=v,E[(X−μ)3]=0,E[(X−μ)4]=3v2.\mathbb{E}[X-\mu]=0,\quad \mathbb{E}[(X-\mu)^2]=v,\quad \mathbb{E}[(X-\mu)^3]=0,\quad \mathbb{E}[(X-\mu)^4]=3v^2.

The mean-mean entry is

Iμμ=E[(X−μv)2]=1v.I_{\mu\mu}=\mathbb{E}\left[\left(\frac{X-\mu}{v}\right)^2\right]=\frac{1}{v}.

The cross entry is zero. In the score-product calculation, it is proportional to the third centered moment, which is zero for a normal distribution. For the variance-variance entry, write the variance score as ((X−μ)2−v)/(2v2)((X-\mu)^2-v)/(2v^2). Since the variance of (X−μ)2(X-\mu)^2 is 3v2−v2=2v23v^2-v^2=2v^2, this gives

Ivv=2v24v4=12v2.I_{vv}=\frac{2v^2}{4v^4}=\frac{1}{2v^2}.

Therefore, for one observation and parameter order (μ,v)(\mu,v),

I1(μ,v)=(1v0012v2).I_1(\mu,v)=\begin{pmatrix}\frac{1}{v}&0\\[4pt]0&\frac{1}{2v^2}\end{pmatrix}.

Information for n independent observations

For independent observations, the log-likelihood is a sum of the individual log-likelihoods. Expected Fisher information therefore adds across observations. For X1,…,XnX_1,\ldots,X_n drawn independently from the same normal distribution,

In(μ,v)=nI1(μ,v)=(nv00n2v2).I_n(\mu,v)=nI_1(\mu,v)=\begin{pmatrix}\frac{n}{v}&0\\[4pt]0&\frac{n}{2v^2}\end{pmatrix}.

The result does not depend on the observed sample values. That is a feature of expected Fisher information; the observed information, calculated from a particular sample’s Hessian, can depend on the data.

If the parameter is standard deviation

Some sources use σ\sigma rather than v=σ2v=\sigma^2. Reparameterizing gives, for parameter order (μ,σ)(\mu,\sigma),

In(μ,σ)=(nσ2002nσ2).I_n(\mu,\sigma)=\begin{pmatrix}\frac{n}{\sigma^2}&0\\[4pt]0&\frac{2n}{\sigma^2}\end{pmatrix}.

The variance entry changes because v=σ2v=\sigma^2: a change in σ\sigma produces a change of 2σ2\sigma in vv. Thus the variance-parameter entry n/(2v2)n/(2v^2) is not the same as the standard-deviation-parameter entry 2n/σ22n/\sigma^2. Check which parameter a problem calls the variance or scale before using a formula.

Worked example

Suppose there are n=100n=100 independent observations and the model has variance v=4v=4 (so σ=2\sigma=2). With parameter order (μ,v)(\mu,v), the information matrix is

I100(μ,v)=(100/400100/(2⋅42))=(25003.125).I_{100}(\mu,v)=\begin{pmatrix}100/4&0\\0&100/(2\cdot 4^2)\end{pmatrix} =\begin{pmatrix}25&0\\0&3.125\end{pmatrix}.

Its inverse is

I100−1=(0.04000.32).I_{100}^{-1}=\begin{pmatrix}0.04&0\\0&0.32\end{pmatrix}.

For a regular model, the inverse information gives the Cramér–Rao lower-bound matrix for unbiased estimators, and it also gives the asymptotic covariance matrix of the maximum-likelihood estimator under standard conditions. Here the corresponding asymptotic standard errors are 0.04=0.2\sqrt{0.04}=0.2 for the mean and 0.32≈0.566\sqrt{0.32}\approx0.566 for the variance. These are asymptotic approximations, not a claim that every finite-sample estimator has exactly this covariance.

Common mistakes

  • Mixing up variance and standard deviation. If the parameter is vv, use n/(2v2)n/(2v^2); if it is σ\sigma, use 2n/σ22n/\sigma^2.
  • Forgetting the sample size. The displayed one-observation matrix must be multiplied by nn for nn independent observations.
  • Assuming zero cross-information means the parameters are known independently in every sense. It means the Fisher information matrix has no cross term in this parameterization. It does not make the unknown mean or variance known, and it does not guarantee zero covariance for every finite-sample estimator.
  • Treating inverse information as an exact finite-sample covariance. It is a bound for unbiased estimators and, under regularity conditions, an asymptotic covariance approximation for maximum-likelihood estimators.

References

  1. Statistical Inference, George Casella and Roger L. Berger
  2. Theory of Point Estimation, Erich L. Lehmann and George Casella

Spotted an error? Send a correction. Read how we write and check articles in our editorial policy.