Exceedance probability for PoolPrev, HierPoolPrev, or getPrevalence

[AngusMcLure]

Can we perhaps have an option to return the exceedance probability for a given threshold to PoolPrev and/or HierPoolPrev or even getPrevalence?

Something like PoolPrev(Data, Result, NumInPool , Threshold = X/100) would return the same as PoolPrev(Data, Result, NumInPool), but with an additional column giving the posterior probability that prevalence exceeds X %.

Similarly something like

mod <- PoolRegBayes(Result ~ Cov1 + Cov2 + (1|Cov3), Data, NumInPool)
getPrevalence(mod, Threshold = X)

would return the same as getPrevalence(mod), but with an additional column in the data.frames at each level giving the posterior probability that prevalence exceeds X %.

Exceedance probability isn't a properly defined concept for the frequentist analyses (but similar-ish to one-side p-values), so would have to be restricted to the Bayesian versions of outputs