On Bayes estimators with uniform priors on spheres and their comparative performance with maximum likelihood estimators for estimating bounded multivariate normal means

被引:7
作者
Fourdrinier, Dominique [2 ]
Marchand, Eric [1 ]
机构
[1] Univ Sherbrooke, Dept Math, Sherbrooke, PQ J1K 2R1, Canada
[2] Univ Rouen, LITIS, EA 4108, St Etienne, France
基金
加拿大自然科学与工程研究理事会;
关键词
Restricted parameters; Point estimation; Squared error loss; Dominance; Maximum likelihood; Bayes estimators; Multivariate normal; Unbiased estimate of risk; Sign changes; Modified Bessel functions; MODIFIED BESSEL-FUNCTIONS; DISTRIBUTIONS;
D O I
10.1016/j.jmva.2010.01.011
中图分类号
O21 [概率论与数理统计]; C8 [统计学];
学科分类号
020208 ; 070103 ; 0714 ;
摘要
For independently distributed observables: X-i similar to N(theta(i), sigma(2)), i= 1.....p, we consider estimating the vector theta = (theta(1),...,theta(p))' with loss parallel to d - theta parallel to(2) under the constraint Sigma(p)(i=1) (o(i)-t(i))(2)/sigma(2) <= m(2), with known t(1),...,t(p), sigma(2), m. In comparing the risk performance of Bayesian estimators delta(alpha) associated with uniform priors on spheres of radius alpha centered at (t(1),...,t(p)) with that of the maximum likelihood estimator delta(mle), we make use of Stein's unbiased estimate of risk technique, Karlin's sign change arguments, and a conditional risk analysis to obtain for a fixed (m, p) necessary and sufficient conditions on alpha for delta(alpha) to dominate delta(mle). Large sample determinations of these conditions are provided. Both cases where all such delta(alpha)'s and cases where no such delta(alpha)'s dominate delta(mle) are elicited. We establish, as a particular case, that the boundary uniform Bayes estimator delta(m) dominates delta(mle) if and only if m <= k(p) with lim(p ->infinity) k(p)/root p = root 2, improving on the previously known sufficient condition of Marchand and Perron (2001) [3] for which k(p) >= root p. Finally, we improve upon a universal dominance condition due to Marchand and Perron, by establishing that all Bayesian estimators delta(pi) with pi spherically symmetric and supported on the parameter space dominate delta(mle) whenever m <= c(1)(p) with lim(p ->infinity) c1(p)/root p = root 1/3. (C) 2010 Elsevier Inc. All rights reserved.
引用
收藏
页码:1390 / 1399
页数:10
相关论文
共 13 条
[1]  
AMOS DE, 1974, MATH COMPUT, V28, P239, DOI 10.1090/S0025-5718-1974-0333287-7
[2]  
[Anonymous], 1983, STAT SPHERES
[3]   MINIMAX ESTIMATION OF A BOUNDED NORMAL-MEAN VECTOR [J].
BERRY, JC .
JOURNAL OF MULTIVARIATE ANALYSIS, 1990, 35 (01) :130-139
[4]  
BROWN L, 1981, J ACOUST SOC AM, V376, P824
[5]  
GUEYE NR, 2003, THESIS U MONTREAL
[6]   POLYA TYPE DISTRIBUTIONS .1. [J].
KARLIN, S .
ANNALS OF MATHEMATICAL STATISTICS, 1957, 28 (02) :281-308
[7]   Improving on the mle of a bounded location parameter for spherical distributions [J].
Marchand, É ;
Perron, F .
JOURNAL OF MULTIVARIATE ANALYSIS, 2005, 92 (02) :227-238
[8]   On the minimax estimator of a bounded normal mean [J].
Marchand, É ;
Perron, F .
STATISTICS & PROBABILITY LETTERS, 2002, 58 (04) :327-333
[9]  
Marchand É, 2001, ANN STAT, V29, P1078
[10]  
Marchand E., 2004, ISM Lecture notes-Monograph Series, V45, P21