Dear Mr.Zeileis & all. (1) Thx for your reply. Yes, I am talk about the function zeroinfl() from the package "pscl". I want to use Newton Raphson to get parameter estimation ZINB, so I try this: ---------------------------------------------------------------------------------------------------------------------------------- > zinb <- zeroinfl(y ~x1+x2+x3+x4+x5|x1+x2+x3+x4+x5, data = dat, link="logit",dist = "negbin", method="Newton Raphson") Error in match.arg(method) : 'arg' should be one of “Nelder-Mead”, “BFGS”, “CG”, “L-BFGS-B”, “SANN” ----------------------------------------------------------------------------------------------------------------------------------- So, what is the solution (R syntax) to get estimation by Newton Raphson? Thanks Regards, Krist. [[alternative HTML version deleted]]