[R] Bayesian Credible Intervals for a Proportion

[tinazilla]

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.