fix errors in simulation, p.beta = 0, 0.2, 0.5, 0.8