Premium
Lower Bounds on Integrated Risk, Subject to Inequality Constraints
Author(s) -
Kabaila Paul,
Kong Yining
Publication year - 2016
Publication title -
australian and new zealand journal of statistics
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.434
H-Index - 41
eISSN - 1467-842X
pISSN - 1369-1473
DOI - 10.1111/anzs.12159
Subject(s) - mathematics , minimax , upper and lower bounds , bounding overwatch , mathematical optimization , constraint (computer aided design) , univariate , prior probability , wasserstein metric , estimator , statistics , bayesian probability , multivariate statistics , computer science , mathematical analysis , geometry , artificial intelligence
Summary Statistical decision theory can sometimes be used to find, via a least favourable prior distribution, a statistical procedure that attains the minimax risk. This theory also provides, using an ‘unfavourable prior distribution’, a very useful lower bound on the minimax risk. In the late 1980s, Kempthorne showed how, using a least favourable prior distribution, a specified integrated risk can sometimes be minimised, subject to an inequality constraint on a different risk. Specifically, he was concerned with the solution of a minimax‐Bayes compromise problem (‘compromise decision theory’). Using an unfavourable prior distribution, Kabaila & Tuck ([Kabaila, P., 2008]), provided a very useful lower bound on an integrated risk, subject to an inequality constraint on a different risk. We extend this result to the case of multiple inequality constraints on specified risk functions and integrated risks. We also describe a new and very effective method for the computation of an unfavourable prior distribution that leads to a very useful lower bound. This method is simply to maximize the lower bound directly with respect to the unfavourable prior distribution. Not only does this method result in a relatively tight lower bound, it is also fast because it avoids the repeated computation of the global maximum of a function with multiple local maxima. The advantages of this computational method are illustrated using the problems of bounding the performance of a point estimator of (i) the multivariate normal mean and (ii) the univariate normal mean.