47
0

Estimation of a multivariate normal mean with a bounded signal to noise ratio

Abstract

For normal canonical models with XNp(θ,σ2Ip),    S2σ2χk2,  independentX \sim N_p(\theta, \sigma^{2} I_{p}), \;\; S^{2} \sim \sigma^{2}\chi^{2}_{k}, \;{independent}, we consider the problem of estimating θ\theta under scale invariant squared error loss dθ2σ2\frac{\|d-\theta \|^{2}}{\sigma^{2}}, when it is known that the signal-to-noise ratio θσ\frac{\|\theta\|}{\sigma} is bounded above by mm. Risk analysis is achieved by making use of a conditional risk decomposition and we obtain in particular sufficient conditions for an estimator to dominate either the unbiased estimator δUB(X)=X\delta_{UB}(X)=X, or the maximum likelihood estimator δmle(X,S2)\delta_{\hbox{mle}}(X,S^2), or both of these benchmark procedures. The given developments bring into play the pivotal role of the boundary Bayes estimator δBU\delta_{BU} associated with a prior on (θ,σ)(\theta,\sigma) such that θσ\theta|\sigma is uniformly distributed on the (boundary) sphere of radius mm and a non-informative 1σ\frac{1}{\sigma} prior measure is placed marginally on σ\sigma. With a series of technical results related to δBU\delta_{BU}; which relate to particular ratios of confluent hypergeometric functions; we show that, whenever mpm \leq \sqrt{p} and p2p \geq 2, δBU\delta_{BU} dominates both δUB\delta_{UB} and δmle\delta_{\hbox{mle}}. The finding can be viewed as both a multivariate extension of p=1p=1 result due to Kubokawa (2005) and a unknown variance extension of a similar dominance finding due to Marchand and Perron (2001). Various other dominance results are obtained, illustrations are provided and commented upon. In particular, for mp2m \leq \sqrt{\frac{p}{2}}, a wide class of Bayes estimators, which include priors where θσ\theta|\sigma is uniformly distributed on the ball of radius mm, are shown to dominate δUB\delta_{UB}.

View on arXiv
Comments on this paper