powered by
dmnom(x, size=sum(x), prob=stop("no prob arg"))
dmnom(c(1,1,4,4),10,c(.2,.2,.3,.3))
Run the code above in your browser using DataLab