Displaying 1 result from an estimated 1 matches for "marp1".
Did you mean:
map1
2012 May 04
0
Converting code from gee() to geeglm()
...was to use geeglm()
I hope someone can help me to solve a problem that arises when converting a
code from gee to geeglm.
*Here is a code that I wrote with the original data, not weighted: *
> m1 <- gee( Bin ~ educ+agemean+ residencysize + yearx , id = rad09 ,
data = Males, subset = marp1 == 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."]
> ef...