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]]
