Open Access
Translator Disclaimer
May, 1980 Empirical Bayes Estimation of the Multivariate Normal Covariance Matrix
L. R. Haff
Ann. Statist. 8(3): 586-597 (May, 1980). DOI: 10.1214/aos/1176345010

Abstract

Let $\mathbf{S}_{p \times p}$ have a Wishart distribution with scale matrix $\Sigma$ and $k$ degrees of freedom. Estimators of $\Sigma$ are given for each of the loss functions $L_1(\hat{\Sigma}, \Sigma) = \operatorname{tr} (\hat{\Sigma}\Sigma^{-1}) - \log \det (\hat{\Sigma}\Sigma^{-1}) - p$ and $L_2(\hat{\Sigma}, \Sigma) = \operatorname{tr} (\hat{\Sigma}\Sigma^{-1} - I)^2$. The obvious estimators of $\Sigma$ are the scalar multiples of $\mathbf{S}$, i.e., $a\mathbf{S}$ where $0 < a \leqslant 1/k$. (Recall that $(1/k)\mathbf{S}$ is unbiased.) For each problem $(\Sigma, \hat{\Sigma}, L_i), i = 1, 2$, we provide empirical Bayes estimators which dominate $a\mathbf{S}$ by a substantial amount. It is seen that the uniform reduction in the risk function determined by $L_2$ is at least $100(p + 1)/(k + p + 1){\tt\%}$. Dominance results for $L_1$ and $L_2$ were first given by James and Stein.

Citation

Download Citation

L. R. Haff. "Empirical Bayes Estimation of the Multivariate Normal Covariance Matrix." Ann. Statist. 8 (3) 586 - 597, May, 1980. https://doi.org/10.1214/aos/1176345010

Information

Published: May, 1980
First available in Project Euclid: 12 April 2007

zbMATH: 0441.62045
MathSciNet: MR568722
Digital Object Identifier: 10.1214/aos/1176345010

Subjects:
Primary: 62F10
Secondary: 62C99

Rights: Copyright © 1980 Institute of Mathematical Statistics

JOURNAL ARTICLE
12 PAGES


SHARE
Vol.8 • No. 3 &bull