nes_counts <- nes %>% group_by(income) %>% summarise(rvote_count = sum(rvote), n=n()) %>% mutate(frac_r = rvote_count/n, frac_r_err = sqrt(frac_r*(1-frac_r)/n) )
