I am trying to calculate Bayesian Credible Intervals for a proportion (disease prevalence values to be more specific) and am having trouble using R to do this. I am working with ncredint() function but have not had success with it. Please help!
Example: Positive samples = 3 Total sampled = 10 Prevalence = 0.3 pvec <- seq(1,10,by=1) npost = dbinom(pvec,10,prob=0.3, log=FALSE) ncredint(pvec, npost, tol=0.01, verbose=FALSE) But I don't get the right Bayesian CI. What am I doing wrong? Thanks! Tina -- View this message in context: http://r.789695.n4.nabble.com/Bayesian-Credible-Intervals-for-a-Proportion-tp3600976p3600976.html Sent from the R help mailing list archive at Nabble.com. ______________________________________________ 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.