[R] dmultinomial

[Alexander Shenkin]

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