mccrTest.Rd {laser} | R Documentation |
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.
mccrTest(CladeSize, NumberMissing, NumberOfReps, ObservedGamma = NULL, fname=NULL)
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 |
If ObservedGamma
is supplied, mccrTest
returns the p-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 |
Dan Rabosky DLR32@cornell.edu
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.
mccrTest(CladeSize=25, NumberMissing=5, NumberOfReps=50);