if (FALSE) {
model <- coal_model(20, 1, 50000) +
feat_recombination(50) +
feat_mutation(1000) +
feat_selection(strength_A = 1000, time = 0.03) +
sumstat_omega()
stats <- simulate(model)
plot(stats$omega$omega, type = "l")}
Run the code above in your browser using DataLab