tbea (version 1.6.1)

topoFreq: Frequency of topologies in a tree sample

Description

Frequency of topologies in a tree sample

Usage

topoFreq(mphy, output = "index", maxtrees = 10000)

Value

A list with an element containing the the different tree clusters (as multiPhylo) and the absolute, cumulative, and relative frequencies of each topology in the tree sample.

Arguments

mphy

An object of class multiPhylo

output

A character indicating whether the tree indices or the actual trees should be returned. Defaults to "index"

maxtrees

A numeric indicating whether to warn about having more trees than the arbitrary threshold

Details

This function can be used e.w. with a posterior sample of trees from a Bayesian analysis where we want to explore the distribution of topologies in the posterior of trees. This way we can assess topological uncertainty in a more meaningful way than using a majority-rule consensus.

The use of `maxtrees` is actually a convenience for keeping in mind that large amounts of trees can cause memory issues. This can end up in situations which are difficult to debug but that from personal experience have come from exactly that: More trees than memory can fit or which can be processed for calculating similarity. This number will _not_ break the function call but will return a warning. Try to avoid modifying its default value unless you are sure it will not cause any issues under your computing conditions (e.g. when lots of trees are being processed but also large RAM is available).

Examples

# tests
set.seed(1)
library(ape)
trl <- ape::rmtree(10, 4)
tpf <- topoFreq(ape::unroot(trl), output="trees")