check in updated R code from Christine for single optional prevalence value