powered by
The multinomial function and its logarithm
multinomial(x) lmultinomial(x)
Numeric vector
Robin K. S. Hankin
Function multinomial() returns
multinomial()
n n_1,n_2,...,n_k = n!n_1!n_2! n_k!omitted; see PDF
where _i n_i=nn_1+...+n_k=n, and function lmultinomial() returns the natural logarithm of this.
lmultinomial()
x <- runif(10) exp(lmultinomial(x)) - multinomial(x) #should be small
Run the code above in your browser using DataLab