StatsBase.sample
— MethodAbstractMCMC.sample(post::Comrade.Posterior, smplr::Dynesty.NestedSampler, args...; kwargs...)
AbstractMCMC.sample(post::Comrade.Posterior, smplr::Dynesty.DynamicNestedSampler, args...; kwargs...)
Sample the posterior post
using Dynesty.jl
NestedSampler/DynamicNestedSampler
sampler. The args/kwargs
are forwarded to Dynesty
for more information see its docs
This returns a PosteriorSamples object. The samplerstats
includes additional information about the samples, like the log-likelihood, evidence, evidence error, and the sample weights. The final element of the tuple is the original dynesty output file.
To create equally weighted samples the user can use
using StatsBase
chain = sample(post, NestedSampler(dimension(post), 1000))
equal_weighted_chain = sample(chain, Weights(stats.weights), 10_000)