f1 <- c(42, 2, 15, 61)
f2 <- c(69, 5, 31, 145)
replicate.prop.ps(.05, f1, f2)
# Should return:
# Estimate SE z p
# Original: 0.106557377 0.03440159 3.09745539 1.951898e-03
# Follow-up: 0.103174603 0.02358274 4.37500562 1.214294e-05
# Original - Follow-up: 0.003852359 0.04097037 0.09402793 9.250870e-01
# Average: 0.105511837 0.02048519 5.15064083 2.595979e-07
# LL UL
# Original: 0.03913151 0.17398325
# Follow-up: 0.05695329 0.14939592
# Original - Follow-up: -0.06353791 0.07124263
# Average: 0.06536161 0.14566206
Run the code above in your browser using DataLab