Displaying 1 result from an estimated 1 matches for "efrinedri".
2012 May 04
0
Converting code from gee() to geeglm()
...p1 == 1 , + family = binomial, corstr
="unstructured" )
(Intercept) educ agemean residencysize yearx
-0.23875 -0.17931 -0.01470 -0.07418 -0.15200
> se <- summary(m1)$coefficients["yearx", "Robust S.E."]
> efrinedri <- coef(m1)["yearx"] + c(-1, 1) * se * qnorm(0.975)
> printco( y1 = summary(m1)$coefficients["yearx", "Estimate"] , uppdown =
efrinedri )
[1] 0.85 ( 0.59 , 1.24 )
*
*
*Trying to convert it to geeglm() with sample weiht: *
m1 <- geeglm( Bin ~ educ...