Displaying 1 result from an estimated 1 matches for "printco".
Did you mean:
printcp
2012 May 04
0
Converting code from gee() to geeglm()
...cept) 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+agemean+ residencysize + yearx , id = rad09 ,
data = Males, subset = marp1...