source('http://www-stat.stanford.edu/~jtaylo/venezuela/load-data.R') # Assuming "Poisson" mean.YES <- mixture.prob.tie(pmf.COUNTS$y, dpois(0:1000, mean(VOTED.YES))) * ncenter mean.NO <- mixture.prob.tie(pmf.COUNTS$y, dpois(0:1000, mean(VOTED.NO))) * ncenter print(data.frame(mean.YES, mean.NO))