gamStat: gamStat

Description Usage Arguments Value Author(s) References Examples

Description

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

Usage

1
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

Author(s)

Dan Rabosky [email protected]

References

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

Examples

1
2
3
  data("plethodon")
  pleth <- plethodon
  gamStat(pleth)

Example output

Loading required package: ape
Loading required package: geiger
------------------------------
Calculated gamma: 1.877344
pvalue: 0.9697645 
test: one-tailed; Ho: rates have not decreased over time 
*assumes complete taxon sampling.
$gamstat
[1] 1.877344

$pval
[1] 0.9697645

$test
[1] "one-tailed; Ho: rates have not decreased over time"

laser documentation built on Jan. 15, 2017, 4:49 p.m.