Proof of Stein's example

From Wikipedia, the free encyclopedia
Jump to: navigation, search

Stein's example is an important result in decision theory which can be stated as

The ordinary decision rule for estimating the mean of a multivariate Gaussian distribution is inadmissible under mean squared error risk in dimension at least 3.

The following is an outline of its proof. The reader is referred to the main article for more information.

Sketched proof[edit]

The risk function of the decision rule d(\mathbf{x}) = \mathbf{x} is

R(\theta,d) = \mathbb{E}_\theta[ |\mathbf{\theta - X}|^2]
=\int (\mathbf{\theta - x})^T(\mathbf{\theta - x}) \left( \frac{1}{2\pi} \right)^{n/2} e^{(-1/2) (\mathbf{\theta - x})^T (\mathbf{\theta - x}) } m(dx)
 = n.\,

Now consider the decision rule

d'(\mathbf{x}) = \mathbf{x} - \frac{\alpha}{|\mathbf{x}|^2}\mathbf{x}

where \alpha = n-2. We will show that d' is a better decision rule than d. The risk function is

R(\theta,d') = \mathbb{E}_\theta\left[ \left|\mathbf{\theta - X} + \frac{\alpha}{|\mathbf{X}|^2}\mathbf{X}\right|^2\right]
 = \mathbb{E}_\theta\left[ |\mathbf{\theta - X}|^2 + 2(\mathbf{\theta - X})^T\frac{\alpha}{|\mathbf{X}|^2}\mathbf{X} + \frac{\alpha^2}{|\mathbf{X}|^4}|\mathbf{X}|^2 \right]
 = \mathbb{E}_\theta\left[ |\mathbf{\theta - X}|^2 \right] + 2\alpha\mathbb{E}_\theta\left[\frac{\mathbf{(\theta-X)^T X}}{|\mathbf{X}|^2}\right] + \alpha^2\mathbb{E}_\theta\left[\frac{1}{|\mathbf{X}|^2} \right]

— a quadratic in \alpha. We may simplify the middle term by considering a general "well-behaved" function h:\mathbf{x} \mapsto h(\mathbf{x}) \in \mathbb{R} and using integration by parts. For 1\leq i \leq n, for any continuously differentiable h growing sufficiently slowly for large x_i we have:

\mathbb{E}_\theta [ (\theta_i - X_i) h(\mathbf{X}) | X_j=x_j (j\neq i) ]= \int (\theta_i - x_i) h(\mathbf{x}) \left( \frac{1}{2\pi} \right)^{n/2} e^{ -(1/2)\mathbf{(x-\theta)}^T \mathbf{(x-\theta)} } m(dx_i)
= \left[ h(\mathbf{x}) \left( \frac{1}{2\pi} \right)^{n/2} e^{-(1/2) \mathbf{(x-\theta)}^T \mathbf{(x-\theta)} } \right]^\infty_{x_i=-\infty}  
- \int  \frac{\partial h}{\partial x_i}(\mathbf{x}) \left( \frac{1}{2\pi} \right)^{n/2} e^{-(1/2)\mathbf{(x-\theta)}^T \mathbf{(x-\theta)} } m(dx_i)
 = - \mathbb{E}_\theta \left[ \frac{\partial h}{\partial x_i}(\mathbf{X})  | X_j=x_j (j\neq i) \right].

Therefore,

\mathbb{E}_\theta [ (\theta_i - X_i) h(\mathbf{X})]=  - \mathbb{E}_\theta \left[ \frac{\partial h}{\partial x_i}(\mathbf{X}) \right].

(This result is known as Stein's lemma.)

Now, we choose


h(\mathbf{x})  =  \frac{x_i}{|\mathbf{x}|^2}.

If h met the "well-behaved" condition (it doesn't, but this can be remedied -- see below), we would have

\frac{\partial h}{\partial x_i} = \frac{1}{|\mathbf{x}|^2} - \frac{2 x_i^2}{|\mathbf{x}|^4}

and so


\mathbb{E}_\theta\left[\frac{\mathbf{(\theta-X)^T X}}{|\mathbf{X}|^2}\right] = \sum_{i=1}^n \mathbb{E}_\theta \left[ (\theta_i - X_i) \frac{X_i}{|\mathbf{X}|^2} \right]
 = - \sum_{i=1}^n \mathbb{E}_\theta \left[ \frac{1}{|\mathbf{X}|^2} - \frac{2 X_i^2}{|\mathbf{X}|^4} \right]
 = -(n-2)\mathbb{E}_\theta \left[\frac{1}{|\mathbf{X}|^2}\right].

Then returning to the risk function of d' :


R(\theta,d') =  n - 2\alpha(n-2)\mathbb{E}_\theta\left[\frac{1}{|\mathbf{X}|^2}\right] + \alpha^2\mathbb{E}_\theta\left[\frac{1}{|\mathbf{X}|^2} \right].

This quadratic in \alpha is minimized at

\alpha = n-2,\,

giving

R(\theta,d') = R(\theta,d) - (n-2)^2\mathbb{E}_\theta\left[\frac{1}{|\mathbf{X}|^2} \right]

which of course satisfies:


R(\theta,d') < R(\theta,d).

making d an inadmissible decision rule.

It remains to justify the use of


h(\mathbf{X})= \frac{\mathbf{X}}{|\mathbf{X}|^2}.

This function is not continuously differentiable since it is singular at \mathbf{x}=0. However the function


h(\mathbf{X}) = \frac{\mathbf{X}}{\epsilon + |\mathbf{X}|^2}

is continuously differentiable, and after following the algebra through and letting \epsilon \to 0 one obtains the same result.