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)