-mlphylo(x, phy, model = DNAmodel(), search.tree = FALSE,
- quiet = FALSE, value = NULL, fixed = FALSE)
-DNAmodel(model = "K80", partition = 1,
- ncat.isv = 1, invar = FALSE,
- equal.isv = TRUE, equal.invar = 1)
-sh.test(..., x, model = DNAmodel(), B = 100)