loglkhs = binom.log_prob(exp_values * shots) - binom.log_prob(sim_values * shots)