FYI, scipy also returns nan for:
st.multinomial.pmf(np.array([0, 1]), 1, np.array([0, 1]))
And for dirichlet scipy return error for array containing 0:
ValueError: Each entry in 'x' must be greater than zero.
FYI, scipy also returns nan for:
st.multinomial.pmf(np.array([0, 1]), 1, np.array([0, 1]))
And for dirichlet scipy return error for array containing 0:
ValueError: Each entry in 'x' must be greater than zero.