def __init__(self, alpha):
        self._alpha = np.array(alpha)
        if self.alpha.ndim != 1:
            raise ValueError('The input alpha must be a 1D list of concentration parameters.')

        self._dist = st.dirichlet(alpha=self.alpha)