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)