Re: [R] residuals of glm

2014-09-29 Thread Ben Bolker
Yuan Jian yahoo.com> writes: > > Hello, > I found the residuals of glm is not the same as calculated manually. > >y=c(12,14,33,50,67,74,123,141,165,204,253,246,240) > > t=1:13 > > m1=glm(y~t+I(t^2),family=poisson(link="log")) > > coefficients(m1)[1]+coefficients(m1)[2]*log(t)+ > coefficients(m1

[R] residuals of glm

2014-09-29 Thread Yuan Jian
Hello, I found the residuals of glm is not the same as calculated manually. >y=c(12,14,33,50,67,74,123,141,165,204,253,246,240) > t=1:13 > m1=glm(y~t+I(t^2),family=poisson(link="log")) > coefficients(m1)[1]+coefficients(m1)[2]*log(t)+coefficients(m1)[3]*log(t^2) [1] 1.901459 2.257258 2.465388 2.6