laser (version 2.4-1)

gamStat: gamStat

Description

Computes the Gamma Statistic of Pybus and Harvey (2000), using branching times

Usage

gamStat(x, return.list=TRUE)

Arguments

x

a numeric vector of branching times

return.list

logical: return a list with pvalue and the statistic?

Value

If return.list = TRUE, a list containing the following elements:

gamstat

the calculated gamma statistic

pval

One-tailed p-value

If return.list = FALSE, it simply returns a single number, the gamma statistic

References

Pybus, O. G., Harvey, P. H. 2000. New uses for incomplete phylogenies.

Examples

  data("plethodon")
  pleth <- plethodon
  gamStat(pleth)