Hello, I want to compute OR from a simulation, I simulate outcome-y and binary covariate, and I repeat it 1000 times. I built a logistic model, and want to compute OR for every sample At the end I'm looking for the mean(OR) of the 1000 repetitions. Can I use cc (from epicalc)? example: if yall is the outcome and unit is the exposure (together its one sample of 500 observations) for(i in 1:1000:{ or_all[i]<-cc(yall, unit) } I need graph and confidence interval only for the mean(OR), not for every iteration. Thank You, Sigalit. [[alternative HTML version deleted]]