phangorn (version 1.5-0)

SH.test: Shimodaira-Hasegawa Test

Description

This function computes the Shimodaira--Hasegawa test for a set of trees.

Usage

SH.test(..., B = 10000, data=NULL)

Arguments

...

either a series of objects of class "pml" separated by commas, a list containing such objects or an object of class "pmlPart".

B

the number of bootstrap replicates.

data

an object of class "phyDat".

Value

• a numeric vector with the P-value associated with each tree given in ....

References

Shimodaira, H. and Hasegawa, M. (1999) Multiple comparisons of log-likelihoods with applications to phylogenetic inference. Molecular Biology and Evolution, 16, 1114--1116.

See Also

pml, pmlPart, pmlCluster

Examples

data(Laurasiatherian)
dm <- dist.logDet(Laurasiatherian)
tree1 <- NJ(dm)
tree2 <- unroot(upgma(dm))
fit1 <- pml(tree1, Laurasiatherian)
fit2 <- pml(tree2, Laurasiatherian)
fit1 <- optim.pml(fit1) # optimize edge weights
fit2 <- optim.pml(fit2)
SH.test(fit1, fit2)
example(pmlPart)
SH.test(sp, B=1000)