areaPerPairI: Calculation of the area per pair index for rooted trees

Description Usage Arguments Value Author(s) References Examples

View source: R/areaPerPairI.R

Description

This function calculates the area per pair index APP(T) for a given rooted tree T. The tree must not necessarily be binary. APP(T) is defined as

APP(T)=2/(n(n-1))*∑_{1<=i<j<=n} d_T(i,j)

in which n denotes the number of leaves in T, and d_T(i,j) denotes the number of edges on the path between the two leaves i and j. Note that APP(T) can also be computed from the Sackin index S(T) and the total cophenetic index TCI(T) of T as APP(T)=2/n*S(T)-4/(n(n-1))*TCI(T) enabling efficient computation.

The area per pair index does not fulfill the definition of an (im)balance index given in "Tree balance indices: a comprehensive survey" (Fischer et al., 2021).

Usage

1

Arguments

tree

A rooted tree in phylo format.

Value

areaPerPairI returns the area per pair index of the given tree.

Author(s)

Luise Kuehn

References

T. Araújo Lima, F. M. D. Marquitti, and M. A. M. de Aguiar. Measuring Tree Balance with Normalized Tree Area. arXiv e-prints, art. arXiv:2008.12867, 2020.

Examples

1
2
tree <- ape::read.tree(text="((((,),),(,)),(((,),),(,)));")
areaPerPairI(tree)

treebalance documentation built on Oct. 17, 2021, 5:06 p.m.