R/msaBP.R

In heibl/polenta: Multiple Sequence Alignment with PASTA and GUIDANCE

#' @title Bootstrap a Multiple Sequence Alignnment
#' @description Create a bootstrap replicate of multiple sequence alignment.
#' @param msa An object of class \code{\link{DNAbin}} or \code{\link{AAbin}}.
#' @export

msaBP <- function(msa){

  if (!inherits(msa, c("DNAbin", "AAbin")))
    stop("'msa' not of class DNAbin or AAbin (ape)")

  msa[, sample(ncol(msa), replace = TRUE)]
}