R/removeSequences.R

In strataG: Summaries and Population Structure Analyses of Genetic Data

#' @title Remove Sequences
#' @description Remove sequences not used by samples.
#'   
#' @param g a \linkS4class{gtypes} object.
#' 
#' @return a new \linkS4class{gtypes} object with unused sequences removed.
#' 
#' @author Eric Archer \email{eric.archer@@noaa.gov}
#' 
#' @export
#' 
removeSequences <- function(g) {
  if(is.null(getSequences(g))) return(g)
  haps <- getAlleleNames(g)
  g@sequences <- g %>% 
    getLociNames() %>% 
    purrr::map(function(x) getSequences(g)[[x]][haps[[x]]]) %>% 
    stats::setNames(names(haps)) %>% 
    as.multidna()
  g
}