Prevalence from inverse binomial sampling

Haldane's estimate (k - 1) / (k + x - 1) and the limits (qbeta(alpha/2, k, x

• 1), qbeta(1 - alpha/2, k, x)) (Hedderich, Sachs & Reynarowych 2023, eqs 6.30-6.31, as its R code).

Usage

invbpr(k, x, conf_level = 0.95)

Arguments

k

Cases observed (>= 2).

x

Non-cases before the k-th case.

conf_level

Confidence level.

Value

Named list: estimate, lower, upper.

References

Haldane, J. B. S. (1945). Biometrika 33, 222-225.

Examples

invbpr(20, 100)
#> $estimate
#> [1] 0.1596639
#> 
#> $lower
#> [1] 0.1048785
#> 
#> $upper
#> [1] 0.2380531
#>