f <- c(125, 82, 92)
ci.pairs.mult(.05, f)
# Should return:
# Estimate SE LL UL
# 1 2 0.14285714 0.04731825 0.05011508 0.23559920
# 1 3 0.10963455 0.04875715 0.01407230 0.20519680
# 2 3 -0.03322259 0.04403313 -0.11952594 0.05308076
Run the code above in your browser using DataLab