mcmc_runner

zdm.MCMC.mcmc_runner(logpf, outfile, state, params, surveys, nwalkers=10, nsteps=100, nthreads=1, Pn=False, pNreps=True, ptauw=False, log_halo=False, lin_host=False, ind_surveys=False, g0info=None, reset=False)[source]

Handles the MCMC running.

Inputs:

logpf (function) = Log posterior function handle outfile (string) = Name of the output file (excluding .h5 extension) state (params.state) = State object to modify params (dictionary) = Parameter names, min and max values surveys (list) = surveys_sep[0] : list of non-repeater surveys

surveys_sep[1] : list of repeater surveys

grid_params (dictionary) = nz, ndm, dmmax nwalkers (int) = Number of walkers nsteps (int) = Number of steps nthreads (int) = Number of threads (currently not implemented - uses default) Pn (bool) = Include Pn or not pNreps (bool) = Include pNreps or not ptauw (bool) = Include ptauw or not log_halo (bool) = Use a log uniform prior on DMhalo ind_surveys (bool) = Return individual survey data g0info (list) = List of [zDMgrid, zvals, DMvals] Passed to use as speedup if needed

Outputs:

posterior_sample (emcee.EnsembleSampler) = Final sample outfile.h5 (HDF5 file) = HDF5 file containing the sampler