[R] Error in optim

djbanana karl79264219 at gmail.com
Fri Oct 28 21:25:27 CEST 2011


I am trying to run this code and obtain the MLEs for my parameters. However I
am getting this error at the end.

Error in optim(c(1.4, 1.1, 0.8, 0.92, 0.4), poisson.lik, v = v) : 
  objective function in *optim evaluates to length 100 not 1*

Code:
poisson.lik <- function(theta,v){
v=matrix(c(1,3,2,0),nrow=2)
n<-nrow(v)
alpha1 <- theta[1]
alpha2 <- theta[2]
beta1 <- theta[3]
beta2 <- theta[4]
gamma <- theta[5]
log.lik <-
v[1]*log(alpha1*beta2*gamma)-alpha1*beta2*gamma-log(factorial(v[1]))+y*log(alpha2*beta1)-alpha2*beta1-log(factorial(v[2]))
return(-log.lik)
}
optim(c(1.4,1.1,0.8,0.92,0.4),poisson.lik,v=v)

Can anyone let me know what I am doing wrong please?

Thanks,
DJ

--
View this message in context: http://r.789695.n4.nabble.com/Error-in-optim-tp3948980p3948980.html
Sent from the R help mailing list archive at Nabble.com.



More information about the R-help mailing list