Displaying 1 result from an estimated 1 matches for "weiht".
Did you mean:
weight
2012 May 04
0
Converting code from gee() to geeglm()
..., "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 == 1 , family = binomial, weights =
Vigtpan , corstr ="unstructured" )
*I get the following error message and not sure how to work on that. Any
suggestions appreciated*...