def init(self,
log=False,
inv=False,
defaultprior=True,
priorshape=False,
priors=(),
originalscale=True,
scale=1,
sigma_m=False,
p_M=.5,
simarf=None,
p_M_arf=.5,
sigma_arf=0.1):
# nsubiters is missing from init() to indicate that nubiters=1 for
# full bayes
if isinstance(simarf, PCA1DAdd):
self.simarf = simarf
else:
self.simarf = ARFSIMFactory()(simarf)
if not isinstance(self.simarf, PCA1DAdd):
raise TypeError("Simulation ARF must be PCA for FullBayes" +
" not %s" % type(self.simarf).__name__)
self.accept_arfs = [0]
self.p_M_arf = p_M_arf
self.rrsig = sigma_arf
return MetropolisMH.init(self, log, inv, defaultprior, priorshape,
priors, originalscale, scale, sigma_m, p_M)
