source('http://www-stat.stanford.edu/~jtaylo/venezuela/load-data.R')

# Assuming "Poisson", conditional on VOTING.STATION, using model

poiss.mean.NO <- numeric(ncenter)
poiss.mean.YES <- numeric(ncenter)

for (j in 1:ncenter) {
  mean <- mu.YES[j]
  count <- nmachine[j]
  d <- dpois(0:(3*mean), mean)
  d <- d / sum(d)
  if (count > 1) {
    poiss.mean.YES[j] <- sum(prob.tie(d, n=count))
  }

  mean <- mu.NO[j]
  d <- dpois(0:(3*mean), mean)
  d <- d / sum(d)
  if (count > 1) {
    poiss.mean.NO[j] <- sum(prob.tie(d, n=count))
  }
}

mean.YES <- sum(poiss.mean.YES)
mean.NO <- sum(poiss.mean.NO)
print(data.frame(mean.YES, mean.NO))