sesamme.mcmc

Functions

log_likelihood(y, yerr, y_model, mask)

Evaluate the likelihood function by comparing the data and chosen model at every wavelength bin.

log_posterior(theta, x, y, yerr, model_cube, ...)

Calculates the (log of the) posterior probability as log(Ppos) = log(Pprior) + log(likelihood).

log_prior(theta)

Calculate a flat prior probability within the bounds prescribed below.

run_sesamme(filename, runname, x, y, yerr, ...)

Initiate an MCMC procedure and write the results to a file/extension name.

set_chain_size(m)

Updates the number of steps (chain length) of the MCMC run.

set_initial_positions(centers)

Updates the initial positions of the walker ensemble.

set_prior_bounds(prior_dict, ...)

Set boundaries on the priors for an MCMC run.

set_walker_size(m)

Updates the size of the walker ensemble.