Walter's confidence interval for the ratio of two binomial variables (and the relative risk)
Walter's confidence interval for the ratio of two binomial variables (and the relative risk).
walter.ci(x1, x2, n1, n2, a = 0.05)
x1 |
An integer number, greater than or equal to zero. |
x2 |
A secondinteger number, greater than or equal to zero. |
n1 |
An integer number, greater than or x1. |
n2 |
A secondinteger number, greater than or equal to x2. |
a |
The significance level. The produced confidence interval has a confidence level equal to 1-a. |
This calculates a (1-a)% confidence interval for the ratio of two binomial variables (and hence for the relative risk) using Walter's suggestion (Walter, 1975). That is, to add 0.5 in each number. This not only overcomes the problem of zero values, but produces intervals that are more accurate than the classical asymptotic confidence interval (Alharbi and Tsagris, 2018).
A list including:
rat |
The ratio of the two binomial distributions. |
ci |
Walter's confidence interval. |
Michail Tsagris
R implementation and documentation: Michail Tsagris <mtsagris@uoc.gr>
Walter S. (1975). The distribution of Levin's measure of attributable risk. Biometrika, 62(2): 371-372.
Alharbi N. and Tsagris M. (2018). Confidence Intervals for the Relative Risk. Biostatistics and Biometrics, 4(5). doi:10.19080/BBOAJ.2018.04.555647
https://juniperpublishers.com/bboaj/pdf/BBOAJ.MS.ID.555647.pdf
x1 <- rbinom(1, 20, 0.7) x2 <- rbinom(1, 30, 0.6) n1 <- 20 n2 <- 30 walter.ci(x1,x2,n1,n2)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.