Dear Lian, You might be interested to hear that a new package is available on CRAN to perform Bayesian estimation of true prevalence from apparent prevalence: http://cran.r-project.org/package=prevalence <http://cran.r-project.org/package=prevalence>
The function 'truePrev()', that estimates true prevalence from individual samples, is also implemented as an online Shiny application: http://users.ugent.be/~bdvleess/R/prevalence/shiny/ <http://users.ugent.be/~bdvleess/R/prevalence/shiny/> Best wishes, Brecht LianD wrote > Thanks Bert > > I've been through my variables again and have managed to get the code > working - shouldn't have tried to deal with it at the end of the day > yesterday! > > all the best > > Lian -- View this message in context: http://r.789695.n4.nabble.com/Bayesian-estimate-of-prevalence-with-an-imperfect-test-tp4265595p4671014.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.