Hello All, I'm trying to run a maximum likelihood analysis using dmultinomial (i'm avoiding dmultinom as I'd like to run it with vectors for the ML stuff). However, I'm having a hard time getting even the simplest example running. Any help would be greatly appreciated.
library(mc2d) dmultinomial(x=c(0,0,1),prob=c(1,1,1),size=1,log=TRUE)
Error in if (ncol(x) != K) stop("x[] and prob[] must be equal length
vectors or equal col matrix.") :
argument is of length zero
# Once I get the simple stuff above running, I'd like to be able to do
the following:
testmat=rmultinomial(4, 1, prob) testmat
[,1] [,2] [,3] [1,] 1 0 0 [2,] 0 0 1 [3,] 0 0 1 [4,] 0 0 1
dmultinomial(x=testmat,prob=prob,size=1,log=TRUE)
Error in if (!is.null(size) && nrow(size) != n) size <- matrix(t(size), : missing value where TRUE/FALSE needed thanks, allie