Displaying 1 result from an estimated 1 matches for "fedstatus1".
2008 Oct 10
1
Coefficients in a polynomial glm with family poisson/binomial
...and have resorted to re-estimating them by running a second polynomial (lm() this time) on the predicted values from predict() of the glm. This is clearly not a nice way of doing things.
Could anyone please inform me of why this is happening and of a better way around this?
Code:
glm2<-glm(FEDSTATUS1~AGE+I(AGE^2), family=binomial(link="probit"))
summary(glm2) ### first set of "wrong coefficients"
nd1<-expand.grid(AGE=c(1:70))
Pred.Fed1<-predict(glm2,nd1,type="response")
points(predict(glm2,nd1,type="response")~nd1$AGE, col=2)
AGE11<-c(11:70)
P...