wassermanfunction.R

wasserman_rejections = function(by_mail=.8,rejection=.04,
         total_votes=2000000,
         dem_pct=.65,pre_d_2way=.52){
  
  rejections = (((total_votes*by_mail)*rejection))
  
  pre_rej_vote = c(total_votes*pre_d_2way,total_votes*(1-pre_d_2way))
    
  post_rej_vote = c((total_votes*pre_d_2way) - ((rejections * dem_pct)),
    (total_votes*(1-pre_d_2way)) - ((rejections * (1-dem_pct))))
  
  list(pre_rej_vote,
       pre_rej_vote/sum(pre_rej_vote),
       post_rej_vote,
       post_rej_vote/sum(post_rej_vote),
       c(((post_rej_vote/sum(post_rej_vote)) - (pre_rej_vote/sum(pre_rej_vote)))*100)[1]
       )
}