def test_step_elliptical_slice(self):
start, model, (K, L, mu, std, noise) = mv_prior_simple()
unc = noise ** 0.5
check = (("x", np.mean, mu, unc / 10.0), ("x", np.std, std, unc / 10.0))
with model:
steps = (EllipticalSlice(prior_cov=K), EllipticalSlice(prior_chol=L))
for step in steps:
trace = sample(
5000, tune=0, step=step, start=start, model=model, random_seed=1, chains=1
)
self.check_stat(check, trace, step.__class__.__name__)
def test_step_elliptical_slice(self):
start, model, (K, mu, noise) = mv_prior_simple()
unc = noise ** 0.5
check = (('x', np.mean, mu, unc / 10.),
('x', np.std, unc, unc / 10.))
with model:
steps = (
EllipticalSlice(prior_cov=K),
)
for step in steps:
trace = sample(8000, step=step, start=start, model=model, random_seed=1)
yield self.check_stat, check, trace, step.__class__.__name__