modify a code to simulate data from binomial with overdispersion