laser (version 2.1)

mccrTest.Rd: Conduct monte carlo constant rates (gamma statistic) test

Description

Conducts gamma statistic analysis (Pybus and Harvey 2000) for incompletely sampled phylogenies. Phylogenies are simulated to the full clade size under the null hypothesis (constant rate pure birth diversification process) and taxa are randomly pruned from the tree to mimic incomplete sampling. The null distribution of the gamma statistic is then tabulated from these phylogenies.

Usage

mccrTest(CladeSize, NumberMissing, NumberOfReps, ObservedGamma = NULL, fname=NULL)

Value

• A list with the following components:
• null.gamma
  The null distribution of gamma. You can plot a histogram or otherwise inspect these values...
• critical.value
  The 0.05 percentile of the null distribution. This is the value corresponding to alpha = 0.05
• p.value
  The actual p-value, only returned if ObservedGamma is supplied by user

Details

If ObservedGamma is supplied, mccrTest returns the p-value.

References

Pybus, O. G. and Harvey, P. H. 2000. Testing macro-evolutionary models using incomplete molecular phylogenies. Proceedings of the Royal Society of London. Series B. Biological Sciences, 267, 2267-2272.

See Also

gamStat

Examples

mccrTest(CladeSize=25, NumberMissing=5, NumberOfReps=50);