def logp(self, F, Y):
return logdensities.bernoulli(Y, self.invlink(F))
def predict_density(self, Fmu, Fvar, Y):
p = self.predict_mean_and_var(Fmu, Fvar)[0]
return logdensities.bernoulli(Y, p)
def test_bernoulli(x, p):
gpf = logdensities.bernoulli(x, p).numpy()
sps = scipy.stats.bernoulli.logpmf(k=x, p=p)
np.testing.assert_allclose(gpf, sps)
