Dear all: I am trying to use the response.zigp and est.zigp from the ZIGP package. Firstly, I generated a simple data via response.zigp and then I try to fit them by est.zigp function, this is the code: ---------------------------------> library(ZIGP) > #try-ZIGP-1 > n<-100 > X<-matrix(c(rep(1,n),runif(n*2,0,1)),n,3) > W<-matrix(c(rep(1,n),rbinom(n,1,0.5)),n,2) > Z<-matrix(c(rep(1,n),rbinom(n,1,0.5)),n,2) > beta<-c(1,-2,1) > alpha<-c(1,1.5) > gamma<-c(-1,2) > y<-response.zigp(X,W,Z,beta,alpha,gamma) > x1<-as.vector(X[,1]) > x2<-as.vector(X[,2]) > x3<-as.vector(X[,3]) > w1<-as.vector(W[,1]) > w2<-as.vector(W[,2]) > z1<-as.vector(Z[,1]) > z2<-as.vector(Z[,2]) > fm.X<-~x1+x2+x3 > fm.W<-~1 > fm.Z<-~z1+z2 > est.zigp(y,fm.X,fm.W,fm.Z,Offset=rep(1,length(y)),init=FALSE)----------------------------- Then I get the *error:* ---------------------------- Error in optim(par = start.delta, fn = loglikelihood.zigp, gr = gradient, : initial value in 'vmmin' is not finite -------------------------------- Any help you can give me would be greatly appreciated. -- View this message in context: http://r.789695.n4.nabble.com/error-in-ZIGP-package-and-optim-tp4127053p4127053.html Sent from the R help mailing list archive at Nabble.com.