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)
