code and data for initial submission and preprint

v1.0.0

improve fig 3 and get upper/lower for counties above p_hat