Nested sampling for frequentist computation: fast estimation of small -values
arXiv:2105.13923 · doi:10.1103/PhysRevLett.128.021801
Abstract
We propose a novel method for computing -values based on nested sampling (NS) applied to the sampling space rather than the parameter space of the problem, in contrast to its usage in Bayesian computation. The computational cost of NS scales as , which compares favorably to the scaling for Monte Carlo (MC) simulations. For significances greater than about in both a toy problem and a simplified resonance search, we show that NS requires orders of magnitude fewer simulations than ordinary MC estimates. This is particularly relevant for high-energy physics, which adopts a gold standard for discovery. We conclude with remarks on new connections between Bayesian and frequentist computation and possibilities for tuning NS implementations for still better performance in this setting.
5 pages + refs, 3 figures. added refs + schematic of algorithm. closely matches published version
References in corpus (9)
- Observation of a new particle in the search for the Standard Model Higgs boson with the ATLAS detector at the LHC
- Observation of a new boson at a mass of 125 GeV with the CMS experiment at the LHC
- The global electroweak fit at NNLO and prospects for the LHC and ILC
- PolyChord: nested sampling for cosmology
- Searching for Dark Matter Annihilation in Recently Discovered Milky Way Satellites with Fermi-LAT
- Study of the mass and spin-parity of the Higgs boson candidate via its decays to Z boson pairs
- anesthetic: nested sampling visualisation
- On the significance of new physics in decays
- Search for supersymmetry with razor variables in pp collisions at sqrt(s) = 7 TeV