def logp(self, F, Y):
return logdensities.gamma(Y, self.shape, self.invlink(F))
def test_gamma(x, shape, scale):
gpf = logdensities.gamma(x, shape, scale).numpy()
sps = scipy.stats.gamma(a=shape, loc=0.0, scale=scale).logpdf(x)
np.testing.assert_allclose(gpf, sps)