# Without recombination:
model <- coal_model(4, 2) + sumstat_trees()
stats <- simulate(model)
print(stats$trees)
# With recombination:
model <- model + feat_recombination(5)
stats <- simulate(model)
print(stats$trees)
Run the code above in your browser using DataLab