simulate_dirichlet_multinomial function simulate_dirichlet_multinomial(n:Integer, α:Real[_]) -> Integer[_] Simulate a Dirichlet-multinomial distribution. n: Number of trials. α: Concentrations.