if (FALSE) {
library(phyclust, quiet = TRUE)
# For nucleotide
set.seed(1234)
X <- seq.data.toy$org
X[sample(c(T, F), length(X), replace = TRUE, prob = c(0.05, 0.95))] <-
.nucleotide$nid[.nucleotide$code == "-"]
plotgaps(X)
}
Run the code above in your browser using DataLab