It would help us to help you if you told us which package ncredint is in (it is not in any of the packages that is installed on the computer I am presently using, but could be in multiple others).
Does this interval match what you are expecting? > library(TeachingDemos) > hpd(qbeta, shape1=4, shape2=8) [1] 0.09337233 0.58795256 -----Original Message----- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of tinazilla Sent: Wednesday, June 15, 2011 6:42 PM To: r-help@r-project.org Subject: [R] Bayesian Credible Intervals for a Proportion 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. ______________________________________________ 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.
