The Pmc sampler
Adaptive Importance Sampling
Author(s): CosmoSIS Team
Cite: MNRAS 405.4 2381-2390 (2010)
Population Monte-Carlo uses importance sampling with an initial distribution that is gradually adapted as more samples are taken and their likelihood found.
At each iteration some specified number of samples are drawn from a mixed Gaussian distribution. Their posteriors are then evaluated and importance weights calculated. This approximate distribution is then used to update the Gaussian mixture model so that it more closely mirrors the underlying distribution.
Components are dropped if they are found not to be necessary.
This is a python re-implementation of the CosmoPMC alogorithm in the cited paper.
No special installation required; everything is packaged with CosmoSIS
These parameters can be set in the sampler's section in the ini parameter file.
If no default is specified then the parameter is required. A listing of "(empty)" means a blank string is the default.
|components||integer||Number of components in the Gaussian mixture||5|
|iterations||integer||Number of iterations (importance updates) of PMC||30|
|student||boolean||Do not use this. It is a not yet functional attempt to use a Student t mixture.||F|
|final_samples||integer||Samples to take after the updating of the mixture is complete||5000|
|samples_per_iteration||integer||Number of samples per iteration of PMC||1000|
|nu||float||Do not use this. It is the nu parameter for the non-function Student t mode.||2.0|