This is now merged into allele_model.R