GetDist: a Python package for analysing Monte Carlo samples
暂无分享,去创建一个
Monte Carlo techniques, including MCMC and other methods, are widely used and generate sets of samples from a parameter space of interest that can be used to infer or plot quantities of interest. This note outlines methods used the Python GetDist package to calculate marginalized one and two dimensional densities using Kernel Density Estimation (KDE). Many Monte Carlo methods produce correlated and/or weighted samples, for example produced by MCMC, nested, or importance sampling, and there can be hard boundary priors. GetDist's baseline method consists of applying a linear boundary kernel, and then using multiplicative bias correction. The smoothing bandwidth is selected automatically following Botev et al., based on a mixture of heuristics and optimization results using the expected scaling with an effective number of samples (defined to account for MCMC correlations and weights). Two-dimensional KDE use an automatically-determined elliptical Gaussian kernel for correlated distributions. The package includes tools for producing a variety of publication-quality figures using a simple named-parameter interface, as well as a graphical user interface that can be used for interactive exploration. It can also calculate convergence diagnostics, produce tables of limits, and output in latex.
[1] Ericka Stricklin-Parker,et al. Ann , 2005 .
[2] Andrew G. Glen,et al. APPL , 2001 .
[3] Computing In Science & Engineering: Web Computing - Java and Grande Applications , 2003, IEEE Distributed Syst. Online.
[4] Ralf S. Klessen,et al. American Institute of Physics Conference Series , 2010 .
[5] Thomas L. Griffiths,et al. Advances in Neural Information Processing Systems 21 , 1993, NIPS 2009.