Boosting Monte Carlo sampling with a non-Gaussian fit

• Verfasst von: Amendola, Luca [VerfasserIn]
• Verfasst von: Gómez-Valent, Adrià [VerfasserIn]
• Titel: Boosting Monte Carlo sampling with a non-Gaussian fit
• Verf.angabe: Luca Amendola and Adrià Gómez-Valent
• E-Jahr: 2020
• Jahr: 14 August 2020
• Umfang: 13 S.
• Fussnoten: Gesehen am 10.12.2020
• Titel Quelle: Enthalten in: Royal Astronomical SocietyMonthly notices of the Royal Astronomical Society
• Ort Quelle: Oxford : Oxford Univ. Press, 1827
• Jahr Quelle: 2020
• Band/Heft Quelle: 498(2020), 1, Seite 181-193
• ISSN Quelle: 1365-2966
• DOI: doi:10.1093/mnras/staa2362
• Datenträger: Online-Ressource
• Sprache: eng
• K10plus-PPN: 1742419852

Abstract

We propose a new method, called Monte Carlo Posterior Fit, to boost the Monte Carlo sampling of likelihood (posterior) functions. The idea is to approximate the posterior function by an analytical multidimensional non-Gaussian fit. The many free parameters of this fit can be obtained by a smaller sampling than is needed to derive the full numerical posterior. In the examples that we consider, based on supernovae and cosmic microwave background data, we find that one needs an order of magnitude smaller sampling than in the standard algorithms to achieve comparable precision. This method can be applied to a variety of situations and is expected to significantly improve the performance of the Monte Carlo routines in all the cases in which sampling is very time consuming. Finally, it can also be applied to Fisher matrix forecasts and can help solve various limitations of the standard approach.