# NOT RUN {
library("data.table")
n <- 100
pop <- 1000
MoE_Y(P = 0.5, n = n, pop = pop)
DT <- data.table(P = seq(0, 1, 0.01))
DT[, Y := round(pop * P)]
DT[, AMoE := MoE_Y(P, n = 100, pop = 1000)]
DT[Y > 0, RMoE := AMoE / Y]
DT
# }
Run the code above in your browser using DataLab