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: mac...@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@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.
______________________________________________
R-help@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.