def test_log_pdf(self): mu = np.log(np.array([[0.5,3.0,5.0]])) sigma = 0.5 x = np.log(np.array([[0.5,3.0,5.0]])) p = model.pdf(x,mu,sigma) p_test = model.log_pdf(x,mu,sigma) assert_almost_equal(np.log(p),p_test)
def test_pdf(self): mu = np.log(np.array([[0.5,3.0,5.0]])) sigma = 0.5 x = np.log(np.array([[0.5,3.0,5.0]])) p0 = scipy.stats.norm.pdf(x[0][0],loc=mu[0][0],scale=sigma) p1 = scipy.stats.norm.pdf(x[0][1],loc=mu[0][1],scale=sigma) p2 = scipy.stats.norm.pdf(x[0][2],loc=mu[0][2],scale=sigma) p_test = model.pdf(x,mu,sigma) assert_equal(p0*p1*p2,p_test)