gamStat: gamStat

In laser: Likelihood Analysis of Speciation/Extinction Rates from Phylogenies

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

Author(s)

Dan Rabosky [e-mail]

References

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

Examples

  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.