| Title: | Detect Population Structure Within Phylogenetic Trees |
|---|---|
| Description: | Algorithms for detecting population structure from the history of coalescent events recorded in phylogenetic trees. This method classifies each tip and internal node of a tree into disjoint sets characterized by similar coalescent patterns. The methods are described in Volz, E., Wiuf, C., Grad, Y., Frost, S., Dennis, A., & Didelot, X. (2020) <doi:10.1093/sysbio/syaa009>. |
| Authors: | Erik Volz |
| Maintainer: | Erik Volz <[email protected]> |
| License: | GPL (>= 2) |
| Version: | 0.1.0 |
| Built: | 2024-11-20 04:05:26 UTC |
| Source: | https://github.com/cran/treestructure |
Plot TreeStructure tree with cluster and partition variables
## S3 method for class 'TreeStructure' plot(x, use_ggtree = TRUE, ...)## S3 method for class 'TreeStructure' plot(x, use_ggtree = TRUE, ...)
x |
A TreeStructure object |
use_ggtree |
Toggle ggtree or ape plotting behaviour |
... |
Additional arguments passed to ggtree or ape::plot.phylo |
Detect cryptic population structure in time trees
trestruct(tre, minCladeSize = 25, minOverlap = -Inf, nsim = 1000, level = 0.01, ncpu = 1, verbosity = 1, debugLevel = 0)trestruct(tre, minCladeSize = 25, minOverlap = -Inf, nsim = 1000, level = 0.01, ncpu = 1, verbosity = 1, debugLevel = 0)
tre |
A tree of type ape::phylo. Must be rooted and binary. |
minCladeSize |
All clusters within parititon must have at least this many tips. |
minOverlap |
Threshold time overlap required to find splits in a clade |
nsim |
Number of simulations for computing null distribution of test statistics |
level |
Significance level for finding new split within a set of tips |
ncpu |
If >1 will compute statistics in parallel using multiple CPUs |
verbosity |
If > 0 will print information about progress of the algorithm |
debugLevel |
If > 0 will produce additional data in return value |
Estimates a partition of a time-scaled tree by contrasting coalescent patterns. The algorithm is premised on a Kingman coalescent null hypothesis and a test statistic is formulated based on the rank sum of node times in the tree.
A TreeStructure object which includes cluster and partitition assignment for each tip of the tree.
E.M. Volz, Wiuf, C., Grad, Y., Frost, S., Dennis, A., Didelot, X.D. (2020) Identification of hidden population structure in time-scaled phylogenies.
Erik M Volz <[email protected]>
tree <- ape::rcoal(50) struct <- trestruct( tree )tree <- ape::rcoal(50) struct <- trestruct( tree )