powered by
ldmultinom(x, prob)
p <- rep(0.2, 5) x <- rmultinom(1, 10, p) dmultinom(x, prob=p, log=TRUE) ldmultinom(x, p)
Run the code above in your browser using DataLab