Hi,
I am fitting a GEE model using gee R package, but I am not sure how to get
OR and its CI?
Could anyone give me some hints?
Here are some output:
> gee.obj <- gee(Affection~Sibsex+Probandsex,id = FAMID,family =
> binomial,corstr = "independence",data =seldata)
Beginning Cgee S-function, @(#) geeformula.q 4.13 98/01/27
running glm to get initial regression estimate
(Intercept) Sibsex2 Probandsex2
-1.6067213 0.3495830 0.5645009
> summary(gee.obj)
GEE: GENERALIZED LINEAR MODELS FOR DEPENDENT DATA
gee S-function, version 4.13 modified 98/01/27 (1998)
Model:
Link: Logit
Variance to Mean Relation: Binomial
Correlation Structure: Independent
Call:
gee(formula = Affection ~ Sibsex + Probandsex, id = FAMID, data = seldata,
family = binomial, corstr = "independence")
Summary of Residuals:
Min 1Q Median 3Q Max
-0.3334466 -0.3334466 -0.2607218 0.6665534 0.8329557
Coefficients:
Estimate Naive S.E. Naive z Robust S.E. Robust z
(Intercept) -1.6067213 0.5265207 -3.0515824 0.5475068 -2.9346145
Sibsex2 0.3495830 0.4043623 0.8645293 0.4045141 0.8642049
Probandsex2 0.5645009 0.5543902 1.0182374 0.5580999 1.0114691
Estimated Scale Parameter: 1.027439
Number of Iterations: 1
Working Correlation
[,1] [,2] [,3] [,4]
[1,] 1 0 0 0
[2,] 0 1 0 0
[3,] 0 0 1 0
[4,] 0 0 0 1
> exp(cbind(coef(gee.obj),confint(gee.obj,level = 0.9)))
Waiting for profiling to be done...
Error in profile.glm(object, which = parm, alpha = (1 - level)/4, trace =
trace) :
subscript out of bounds
Error in exp(cbind(coef(gee.obj), confint(gee.obj, level = 0.9))) :
Non-numeric argument to mathematical function
> exp(coef(gee.obj)
+ )
(Intercept) Sibsex2 Probandsex2
0.2005441 1.4184760 1.7585698
Thank you for your help,
pingzhao
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
