Estimating the Mean of a Multivariate Normal Population with General Quadratic Loss Function

Let X be a p-variate random vector, p 3, having a normal distribution with mean vector and covariance matrix ²I. For estimating with loss function L₁ (, ) = (-) ' (-) = \| - \|², and known ², Stein 2 has shown the inadmissibility of the usual estimator X by considering an alternative estimator with uniformly smaller risk than that of X, the improvement being substantial for close to the origin. The problem of estimating with the same loss function when ² is unknown has been treated similarly by James and Stein 1, when an observation is available on another random variable which is distributed as ²²ₙ independently of X. James and Stein have also demonstrated the inadmissibility of the usual estimator for under the loss function equation*1L₂ (, ) = (-) 'D (-) equation* for the case of known ², where D is a diagonal matrix with unequal diagonal elements. This result has been proved even without the normality assumption, but no explicit formulas for alternative estimators have been given which improve on the usual estimator in some parts of the parameter space under normality. In this note an estimator for is obtained for the case when ² is unknown and the loss function is L₂. An upper bound for the risk function of this estimator is given, which always remains below the risk function of the usual estimator and is substantially smaller for close to the origin. This estimator coincides with the estimator given by James and Stein when the diagonal elements of D become equal. An application of this estimator gives an improvement on the least squares estimator of the parameter vector in a usual linear observational model with normal errors.