R/set.excluded.sites.R

In phraber/lassie: Tools to select sites and sequences that represent variation in a longitudinally sampled protein alignment

#' @keywords internal
set.excluded.sites <- function(S, site_list) {

    if (class(S) != "swarmtools")
        stop("ERROR: Unexpected status for arguments to set.excluded.sites()")

    ### TO DO: Sanitize then parse site_list

#    if (is.null(S$excluded_sites))
        S$excluded_sites <- sort(unique(site_list))
#    else
#        S$excluded_sites <- sort(unique(c(S$excluded_sites, site_list)))

#    if (is.null(S$aas_aln) | is.null(S$refseq_lut))
#        S <- prep.aln(S)

### TO DO: ensure workflow will proceed if this is the last thing set by ui
#    if (!is.null(S$tf_loss_cutoff))
#        S <- select.sites(S)

    return ( S )
}