def test_sample_results_object_initializes(self):
s = SamplingResults(self.sampler)
assert s.samples.shape[0] == self.nwalkers * self.niter
assert s.acceptance > 0.25
assert np.isclose(s.L,
s.acceptance * self.nwalkers * self.niter)
def test_infer_computes_correct_values(self):
s = SamplingResults(self.sampler)
test_mean = 2.00190793
test_std = 0.00195719
test_ci = [[1.99435539], [1.9971502]]
assert np.isclose(test_mean, s.mean[0], atol=0.01, rtol=0.01)
assert np.isclose(test_std, s.std[0], atol=0.01, rtol=0.01)
assert np.all(np.isclose(test_ci, s.ci, atol=0.01, rtol=0.01))
def test_infer_computes_correct_values(self):
s = SamplingResults(self.sampler)
def test_rhat_computes_correct_answer(self):
s = SamplingResults(self.sampler)
rhat_test = 3.81886815e-06
assert np.isclose(rhat_test, s.rhat[0], atol=0.001, rtol=0.001)
def test_sampling_results_acceptance_ratio(self):
s = SamplingResults(self.sampler)
assert s.acceptance > 0.25
assert np.isclose(s.L, s.acceptance * self.nwalkers * self.niter)
def test_sample_results_produces_attributes(self):
s = SamplingResults(self.sampler)
assert s.samples.shape[0] == self.nwalkers * self.niter
def test_sample_results_object_initializes(self):
SamplingResults(self.sampler)