if (FALSE) {
library(phyclust, quiet = TRUE)
Tt <- 0.5
Q <- list(pi = c(0.25, 0.25, 0.25, 0.25), kappa = 0.5)
phyclust.Pt(Q, Tt, "HKY85")
Q <- list(pi = c(0.5, 0.5), kappa = 0.5)
phyclust.Pt(Q, Tt, "SNP_JC69")
}
Run the code above in your browser using DataLab