simulation.r

N < - 8619170             # population size of voters
ss <- c(1000, 1200, 3000) # samples size
p <- .27                 # .27 Serra .25 Haddad .19 Russomano true population proportion
nsim <- 100              # number of simulations

pop.prob <- list()
for (i in 1:length(ss)) {
  n <- ss[i]
  x <- rhyper(nsim, N * p, N * (1 - p), n)
  pop.prob[[i]] <- x / n
}
names(pop.prob) <- paste("n =", ss)