def check_exponential(inverse_scale): rng = check_random_state(1) p_carl = Exponential(inverse_scale=inverse_scale) p_scipy = st.expon(scale=1. / inverse_scale) X = rng.rand(50, 1) assert_array_almost_equal(p_carl.pdf(X), p_scipy.pdf(X.ravel())) assert_array_almost_equal(p_carl.cdf(X), p_scipy.cdf(X.ravel())) assert_array_almost_equal(-np.log(p_carl.pdf(X)), p_carl.nll(X))
def check_exponential(inverse_scale): rng = check_random_state(1) p_carl = Exponential(inverse_scale=inverse_scale) p_scipy = st.expon(scale=1. / inverse_scale) X = rng.rand(50, 1) assert_array_almost_equal(p_carl.pdf(X), p_scipy.pdf(X.ravel())) assert_array_almost_equal(p_carl.cdf(X), p_scipy.cdf(X.ravel())) assert_array_almost_equal(-np.log(p_carl.pdf(X)), p_carl.nnlf(X))