








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

