glm convergence warning
Hi Sophie It helps if you do some detective work Try fit1 <- glm(formula = Y ~ offset(cst) + beta + tau, family = poisson,trace = T, maxit = 200) and compare Regards Duncan Duncan Mackay Department of Agronomy and Soil Science University of New England Armidale NSW 2351 Email: home: mackay at northnet.com.au
At 05:35 18/05/2012, you wrote:
Hi, When I run the following code : Y <- c(rep(0,35),1,2,0,6,8,16,43) cst <- log(choose(42, 42:1)) beta <- 42:1 tau <- (beta^2)/2 fit <- glm(formula = Y ~ offset(cst) + beta + tau, family = poisson) fit fit$converged glm prints a warning saying that the algorithm did not converge. However, fit$converged takes the value TRUE. I don't understand why fit$converged is not always FALSE when the warning "algorithm did not converge" is produced. Could someone help me understand why I get this result? Thanks a lot, Sophie
______________________________________________ R-help at r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.