rgsl/randist/
multinomial.rs

1//
2// A rust binding for the GSL library by Guillaume Gomez (guillaume1.gomez@gmail.com)
3//
4
5/// This function computes the probability P(n_1, n_2, ..., n_K) of sampling `n[K]` from a
6/// multinomial distribution with parameters `p[K]`, using the formula given above.
7#[doc(alias = "gsl_ran_multinomial_pdf")]
8pub fn multinomial_pdf(p: &[f64], n: &[u32]) -> f64 {
9    unsafe { sys::gsl_ran_multinomial_pdf(p.len() as _, p.as_ptr(), n.as_ptr()) }
10}
11
12/// This function returns the logarithm of the probability for the multinomial
13/// distribution P(n_1, n_2, ..., n_K) with parameters `p[K]`.
14#[doc(alias = "gsl_ran_multinomial_lnpdf")]
15pub fn multinomial_lnpdf(p: &[f64], n: &[u32]) -> f64 {
16    unsafe { sys::gsl_ran_multinomial_lnpdf(p.len() as _, p.as_ptr(), n.as_ptr()) }
17}