
Last chance! 50% off unlimited learning
Sale ends in
This statistic returns ancestral tress in NEWICK format.
sumstat_trees(name = "trees")
The name of the summary statistic. When simulating a model, the value of the statistics are written to an entry of the returned list with this name. Summary statistic names must be unique in a model.
To create a demographic model: coal_model
To calculate this statistic from data: calc_sumstats_from_data
Other summary statistics:
sumstat_dna()
,
sumstat_file()
,
sumstat_four_gamete()
,
sumstat_ihh()
,
sumstat_jsfs()
,
sumstat_mcmf()
,
sumstat_nucleotide_div()
,
sumstat_omega()
,
sumstat_seg_sites()
,
sumstat_sfs()
,
sumstat_tajimas_d()
# 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