gsi: gsi
In santiagosnchez/rBt: r BEAST tools (rBt)

Description Usage Arguments Value References Examples

View source: R/gsi.R

Computes genealogical sorting index (Cummings et al. 2008)

1	gsi(phy, groups)

`phy`	a `phylo` object
`groups`	a `list` where each element is a character vector with groups of tips

a vector with gsi values

Cummings, M.P., M.C. Neal, K L. Shaw (2008) A genealogical approach to quantifying lineage divergence. Evolution 62:2411-2422 doi:10.1111/j.1558-5646.2008.00442.x

set.seed(1)
tr = rtree(10)
groups = split(tr$tip.label, factor(rep(c(1,2,3), c(3,3,4))))
tr_gsi = gsi(tr, groups)
tr_gsi
#     1     2     3 
# 1.000 1.000 0.625