def pred(self, x):
"""Prior predictive. Pr(x)"""
return multivariate_t_density(self.nu_0-self.d+1, self.mu_0,
self.Lam_0*(self.kappa_0+1)/(self.kappa_0 - self.d + 1), x)
def pred(self, x):
"""Prior predictive. Pr(x)"""
assert isinstance(x, np.ndarray)
assert x.shape[-1] == 2
# Generalized from InvGamma. Tested numerically.
return multivariate_t_density(2*self.alpha, self.mu, 1./self.beta/self.alpha*np.eye(2), x)
