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.