MH

class sherpa.sim.mh.MH(fcn, sigma, mu, dof, *args, rng=None)

Bases: Sampler

The Metropolis Hastings Sampler

Changed in version 4.16.0: The rng parameter was added.

Random number generation is controlled by the rng argument. If set to None (the default) then the routines from numpy.random are used, and so can be controlled by calling numpy.random.seed, otherwise it takes a numpy.random.Generator object (or a numpy.random.RandomState object which should only be used for testing or checking against old code).

Methods

accept(current, current_stat, proposal, ...)	Should the proposal be accepted (using the Cash statistic and the t distribution)?
draw(current)	Create a new set of parameter values using the t distribution.
mh(current)	MH jumping rule
update(stat, mu[, init])	include prior

accept_mh
calc_fit_stat
calc_stat
dmvt
init
reject
tear_down

Methods Summary

accept(current, current_stat, proposal, ...)	Should the proposal be accepted (using the Cash statistic and the t distribution)?
accept_mh(current, current_stat, proposal, ...)
calc_fit_stat(proposed_params)
calc_stat(proposed_params)
dmvt(x[, log, norm])
draw(current)	Create a new set of parameter values using the t distribution.
init([log, inv, defaultprior, priorshape, ...])
mh(current)	MH jumping rule
reject()
tear_down()
update(stat, mu[, init])	include prior

Methods Documentation

accept(current, current_stat, proposal, proposal_stat, **kwargs)

Should the proposal be accepted (using the Cash statistic and the t distribution)?

accept_mh(current, current_stat, proposal, proposal_stat)

calc_fit_stat(proposed_params)

calc_stat(proposed_params)

dmvt(x, log=True, norm=False)

draw(current)

Create a new set of parameter values using the t distribution.

Given the best-guess (mu) and current (current) set of parameters, along with the covariance matrix (sigma), return a new set of parameters.

init(log=False, inv=False, defaultprior=True, priorshape=False, priors=(), originalscale=True, scale=1, sigma_m=False)

mh(current)

MH jumping rule

reject()

tear_down()

update(stat, mu, init=True)

include prior