Lecture 11 Fisher Information
Since the data are random, the score function from lecture 4 is a random variable, and the Fisher information, given by the variance of the score, measures how much information the sample contains about the parameter. In lecture 13 we will use the Fisher information to obtain a lower bound for the variance of any unbiased estimator, and in lecture 17 we will see that the reciprocal of the Fisher information is the asymptotic variance of the maximum likelihood estimator. In this lecture we will define the Fisher information, prove two identities which allow us to compute the Fisher information, and we will compute the Fisher information for the standard families of models.
11.1 Setting and Assumptions
Throughout this lecture, are i.i.d. with density or probability weights , where the parameter space is an open interval. As in lecture 4, the log-likelihood is the sum of the terms and the score of definition 4.4 is its derivative . In this section we study the distribution of the score under , for a fixed value of .
The results of this lecture are obtained by differentiating the identity
with respect to . Throughout the section we make the following assumptions:
-
1.
The support of the distribution does not depend on .
-
2.
For every the function is twice continuously differentiable on .
-
3.
The integral of over the support may be differentiated under the integral sign twice, so that
equation (11.1) (11.1)where the integrals on the right-hand side exist, with sums in place of integrals for a discrete model.
-
4.
The score of a single observation has finite variance, for all , so that the quantities defined below are finite.
These four assumptions are referred to as the conditions of this section. The one-parameter models from tables 1.1 and 1.2, and the two-parameter models with one parameter held fixed (with known for the binomial distribution), satisfy these conditions, with the exception of the uniform distribution (see section 11.6). So do the exponential families of lecture 7, as explained after lemma 7.10 there.
11.2 The Score Has Mean Zero
The first identity is that, if is the correct parameter value, the data do not systematically deviate from this value: the score at the true parameter value is zero on average.
Using the chain rule, we find the derivative of as
Since lies in the support with probability one, the denominator is positive and the ratio is well defined. Cancelling the density against the denominator we find
where the integral is over the support of the distribution, which does not depend on (with a sum in place of the integral for a discrete model), and the third equality comes from the case of (11.1). Summing over gives the statement for the whole sample. This completes the proof. ∎
If the parameter value is wrong, the score is not mean zero: for the exponential model from example 4.5, the score has expectation if the true parameter value is , positive for and negative for . The score on average points towards the correct parameter value from both sides. In lecture 16 we will see the corresponding property of the log-likelihood, which on average is larger at the true parameter value than at any other value, and use it to prove that the maximum likelihood estimator is consistent.
11.3 Fisher Information
The next question we could ask is about the fluctuations of the score around zero. If the score is typically large in magnitude, this indicates that the data allow us to distinguish well between nearby parameter values. In contrast, if the score is typically close to zero, this indicates that the log-likelihood is nearly flat at the truth. The variance of the score is a measure for how informative the sample is.
Under the conditions of section 11.1, the Fisher information of the sample about is given by
where the variance of the score is taken at . The Fisher information of a single observation is denoted by and is given by the variance of .
The subscript in is the sample size. Where we compare two different parametrisations of the same model, we write , or for the information of a single observation about the parameter named in the subscript.
The variance of the score is often awkward to compute directly, and the second expression of the following lemma is nearly always the quickest route to the information.
The first equality follows from the fact that has mean zero by lemma 11.1, and thus the variance equals the second moment. The second equality is obtained by taking the derivative of the ratio (11.2) using the quotient rule:
Similar to the proof of lemma 11.1, the expectation of the first term vanishes, using the case of (11.1):
Since the second moment of is finite by the conditions of section 11.1, we can take expectations term by term and get for all .
To go from to we also need to consider the cross terms: are independent of each other, since the are independent. Thus the variance of the sum is the sum of the variances. Since each term also has mean zero we find
where we used in the last step. This completes the proof. ∎
The second expression motivates the name Fisher Information: the quantity is the expected curvature of the log-likelihood at the true parameter value. Large information means that the log-likelihood has a sharp maximum, well-defined by the data. Small information means the log-likelihood is flat. Note that the identity holds in expectation only: for a given sample the curvature is random and in general differs from .
The proof also shows how the information of the sample relates to the information of one observation: each observation contributes the same amount, and observations carry times the information of one.
The last display of the proof of lemma 11.3 writes as the sum of the variances , which are all equal to since the observations are identically distributed. This completes the proof. ∎
11.4 Examples
We now compute the information for four of the standard families of tables 1.1 and 1.2; the exercises cover the binomial, gamma and geometric distributions. By proposition 11.4 it suffices to work with a single observation .
For the Bernoulli distribution with success probability we have and thus
Using we find
Thus, the information of the sample is . Exercise 11.1 checks these identities for this model directly.
For the Poisson distribution with mean we have and thus
In the case of the normal distribution with unknown mean and known variance , the parameter is and . Thus we have
Since the second derivative does not depend on the data, we do not need to consider the expectation. Note that by lemma 2.2, and in lecture 13 we will see that no unbiased estimator can have a smaller variance.
For the exponential distribution with rate , we have and thus
Again, the second derivative is non-random. Section 11.5 shows how the information changes if we describe the model by the mean instead of the rate.
11.5 Reparametrisation
The information depends on which parameter we use to describe a model. The exponential distribution, for example, can be described by its rate or by its mean, and the following proposition gives the general rule for how the information changes.
Let where is a twice continuously differentiable bijection with for all . Let denote the Fisher information for a single observation when the model is parametrised by and let be the information when the model is parametrised by . Then,
With the new parametrisation, the density is and the log-likelihood for a single observation is . Since is continuous and never vanishes on the interval , it has the same sign throughout, so that is strictly monotonic and is twice continuously differentiable. Using the chain rule, the conditions of section 11.1 for the parametrisation by carry over to the parametrisation by ; we omit the verification. The derivative of the inverse function is for . This gives
Since is fixed for fixed , is a constant and taking variances under gives , as claimed. ∎
For the exponential distribution, reparametrising from the rate to the mean changes the information from example 11.8 to ; see exercise 11.6. Thus, the information is not a property of the model in isolation, but depends on the parametrisation. In lecture 28 we will use the transformation rule of the proposition to construct a prior distribution, the Jeffreys prior, which does not depend on the parametrisation chosen.
11.6 When the Identities Fail
The uniform distribution from example 1.2 shows what can go wrong if the conditions from section 11.1 are violated.
Let be uniformly distributed on , i.e. for . The support of the distribution depends on , but we can still compute the derivatives on the support to see whether the identities hold. For we have and thus and . The score is a negative constant and thus we have , in contradiction to lemma 11.1. The variance is zero, whereas and : the three expressions for lemma 11.3 take three different values, one of which is negative. The proof breaks down at the interchange between differentiation and integration: for this model, the two sides of (11.1) with differ by the term from the moving upper limit, which is missed by the differentiation under the integral sign. Exercise 11.7 can be used to work out the two computations.
For this model the Fisher information is not defined, and the two results based on it, in lectures 13 and 17, do not apply.
-
•
Under the conditions of section 11.1 (a support free of , twice continuously differentiable density, differentiation under the integral sign and a score with finite variance), the score at the true parameter value has mean zero.
-
•
The Fisher information is the variance of the score, and it equals and , the expected curvature of the log-likelihood at the truth.
-
•
Information from independent observations can be combined, . The information depends on the parametrisation, and under we get .
-
•
We have learned the following standard values: for the Bernoulli distribution, for the Poisson distribution, for the mean of the normal distribution and for the rate parameter in the exponential distribution. For the uniform distribution these relations do not hold, since the support of the distribution moves with the parameter.
- •
Let be i.i.d. Bernoulli with success probability , and let be the score of a single observation from example 11.5.
-
1.
Verify lemma 11.1 for this model by computing directly from the two possible values of .
-
2.
Compute directly, and check that it agrees with the value of found in example 11.5.
-
3.
Compare with the variance of the estimator for .
-
4.
Now let be a single observation from the binomial distribution with known number of trials and success probability . Show that the information of about is , and explain the relation to the Bernoulli case using proposition 11.4.
Let be an exponential family in its natural parametrisation from definition 7.9, where ranges over an open interval inside the natural parameter space , and let be a single observation.
-
1.
Show that the Fisher information about is given by .
-
2.
Now let the family be written as as in definition 7.1, where is twice continuously differentiable and takes values in the interior of . Show that and that the information about is given by .
- 3.
Let be i.i.d. with density for , where , the gamma distribution with known shape and scale from exercise 5.2, i.e. in the rate parametrisation from table A.2.
-
1.
Compute the score of a single observation and verify that it has mean zero.
-
2.
Show that , once via and once via the variance of the score, using .
- 3.
-
4.
In exercise 5.2 we found that the MLE is unbiased with . Compare this with .
Let be i.i.d. from the geometric distribution with parameter , i.e. for , where and , as in exercises 4.2 and 5.1.
-
1.
Compute the score of a single observation and verify that it has mean zero.
-
2.
Show that and write down .
-
3.
Verify the value from the previous part by taking the variance of the score.
Let be i.i.d. Poisson with mean .