Rough R script for drift simulation two trait, four locus population change