CONTENTS

Exercise: Fisher Information and the Cramér-Rao Bound for Exponential Rate

Problem

Let X1,,XnX_1, \ldots, X_n be i.i.d. Exp(λ)\text{Exp}(\lambda): p(x;λ)=λeλxp(x; \lambda) = \lambda e^{-\lambda x} for x0x \ge 0.

  1. Compute the log-likelihood (λ)\ell(\lambda) and derive the MLE λ^=1/xˉ\hat\lambda = 1/\bar x.

  2. Compute the Fisher information per observation: I(λ)=E ⁣[2logp(X1;λ)λ2].I(\lambda) = -\mathbb{E}\!\left[\frac{\partial^2 \log p(X_1; \lambda)}{\partial \lambda^2}\right]. (Hint: 2/λ2(logλλx)=1/λ2\partial^2/\partial\lambda^2\,(\log\lambda - \lambda x) = -1/\lambda^2.)

  3. The Cramér-Rao lower bound for unbiased estimators of λ\lambda is Var(λ^)1/(nI(λ))=λ2/n\text{Var}(\hat\lambda) \ge 1/(n I(\lambda)) = \lambda^2/n. Does λ^=1/Xˉ\hat\lambda = 1/\bar X saturate this bound? (Note: λ^\hat\lambda is biased.)

  4. Asymptotic normality simulation. For λ=2\lambda = 2 and n=100n = 100, generate m=10,000m = 10{,}000 sample sets and compute λ^\hat\lambda for each. Compare the empirical distribution of n(λ^λ)\sqrt n(\hat\lambda - \lambda) to N(0,λ2)\mathcal{N}(0, \lambda^2). Plot a histogram overlay or report mean and variance.

Hint

For part 2: logp(x;λ)=logλλx\log p(x; \lambda) = \log\lambda - \lambda x. First derivative: 1/λx1/\lambda - x. Second derivative: 1/λ2-1/\lambda^2. Fisher information: E[1/λ2]=1/λ2-\mathbb{E}[-1/\lambda^2] = 1/\lambda^2.

Jump to the solution when you're ready.