def PY(self, t, xp, x): return dists.IndepProd(*[dists.Poisson(np.exp(self.a[k] + np.sum(x*self.b[k], axis=1))) for k in range(self.dy)])
def PY(self, t, xp, x): return dists.Poisson(rate=np.exp(x))
def PY(self, t, xp, x): return dists.IndepProd(dists.Poisson(rate=np.exp(x)), dists.Poisson(rate=np.exp(x)))