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