laser (version 2.4-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)

Arguments

CladeSize

The TRUE clade diversity

NumberMissing

The number of missing species, e.g., CladeSize minus the number of taxa in your tree

NumberOfReps

The number of Monte Carlo simulations to conduct. Recommend >= 5000

ObservedGamma

The observed gamma statistic value for the empirical tree. Optional.

fname

An optional filename where simulated trees are stored, if you generated the trees in another program, such as phylogen

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);