Constructing Confidence Regions of Optimal Expected Size

• 成果类型: Article
• 署名作者: Schafer, Chad M.; Stark, Philip B.
• 署名单位: Carnegie Mellon University; University of California System; University of California Berkeley
• 刊物名称: JOURNAL OF THE AMERICAN STATISTICAL ASSOCIATION
• ISSN/ISSBN: 0162-1459
• DOI: 10.1198/jasa.2009.tm07420
• 发表日期: 2009
• 页码: 1080-1089
• 关键词: dark energy; interval estimation; scale-parameters; sets; constraints; Admissibility; supernovae; location

摘要

This article presents a Monte Carlo method for approximating the minimax expected size (MES) confidence set for a parameter known to lie in a compact set. The algorithm is motivated by problems in the physical sciences in which parameters are unknown physical constants related to the distribution of observable phenomena through complex numerical models. The method repeatedly draws parameters at random from the parameter space and simulates data as if each of those values were the true value of the parameter. Each set of simulated data is compared to the observed data using a likelihood ratio test. Inverting the likelihood ratio test minimizes the probability of including false values in the confidence region, which in turn minimizes the expected size of the confidence region. We prove that as the size of the Simulations grows, this Monte Carlo confidence set estimator converges to the Gamma-minimax procedure, where Gamma is a polytope of priors. Fortran-90 implementations of the algorithm for both serial and parallel computers are available. We apply the method to an inference problem in cosmology.