```
x1 <- 56; n1 <- 70; x2 <- 48; n2 <- 80
xci <- BinomDiffCI(x1, n1, x2, n2, method=c("wald", "waldcc", "ac", "score",
"scorecc", "mn", "mee", "blj", "ha"))
Format(xci[,-1], digits=4)
x1 <- 9; n1 <- 10; x2 <- 3; n2 <- 10
yci <- BinomDiffCI(x1, n1, x2, n2, method=c("wald", "waldcc", "ac", "score",
"scorecc", "mn", "mee", "blj", "ha"))
Format(yci[, -1], digits=4)
# https://www.lexjansen.com/wuss/2016/127_Final_Paper_PDF.pdf, page 9
SetNames(round(
BinomDiffCI(56, 70, 48, 80,
method=c("wald", "waldcc", "hal",
"jp", "mee",
"mn", "score", "scorecc",
"ha", "ac", "blj"))[,-1], 4),
rownames=c("1. Wald, no CC", "2. Wald, CC", "3. Haldane", "4. Jeffreys-Perks",
"5. Mee", "6. Miettinen-Nurminen", "10. Score, no CC", "11. Score, CC",
"12. Hauck-Andersen", "13. Agresti-Caffo", "16. Brown-Li"))
```

Run the code above in your browser using DataLab