Please don't post in HTML. Your post is almost unreadable Also,lease have a look at http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example and/or http://adv-r.had.co.nz/Reproducibility.html for some suggestions about what you should include in your question.
Welcome to R. John Kane Kingston ON Canada > -----Original Message----- > From: rosit...@gmail.com > Sent: Thu, 14 Apr 2016 19:07:04 +0100 > To: pd....@cbs.dk, r-help@r-project.org > Subject: [R] Odds Ratio and OR CI > > Howdy everyone > > > > I’m trying to get Odds ratio and OR confidence intervals using a probit > model, but I'm not getting. > > > > Do you think you can help me? > > > > I’m new with R L > > > > naive = > summary(glm(pcr.data[,7]~boldBeta_individual+pcr.data$age,family=binomial(link=probit))) > > naive_answer = c(naive$coefficients[,1:3]) > #naive estimates for > > > #alpha (first 4 collumns: intercept; beta_intercept, beta_slope and age) > and > > > #and SE(last 4 collumns: intercept; beta_intercept, beta_slope and age) > > > > OR.naive = exp(1.6*coef(naive)) > > > > (till here works, the problem is with the confidence interval) > > > > I tried to get the Standard error from the variance, but I’m not sure if > this can be done as I’ve done. > > > > > > Var_coef <- 1.6^2*var(coef(naive)) > > SE_coef <- Var_coef/sqrt(nsample) ########## I thi k > this is correct > > > > OR.naive.inf <- exp(OR.naive - (1.96 * SE_coef)) > > OR.naive.sup <- exp(OR.naive + (1.96 * SE_coef)) > > > > if I used logit link I would get the CI with confint(naïve) command, but > with probit I don't think so. Is there a way? > > > > What should I do? > > > > > > > > > Atenciosamente, > Rosa Oliveira > > -- > ____________________________________________________________________________ > > > > Rosa Celeste dos Santos Oliveira, > > E-mail: rosit...@gmail.com <mailto:rosit...@gmail.com> > Tlm: +351 939355143 > Linkedin: https://pt.linkedin.com/in/rosacsoliveira > <https://pt.linkedin.com/in/rosacsoliveira> > ____________________________________________________________________________ > "Many admire, few know" > Hippocrates > > ______________________________________________ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > 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. ____________________________________________________________ FREE 3D EARTH SCREENSAVER - Watch the Earth right on your desktop! ______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see 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.
