added code to simulate two-locus model