R/update.ref.geno.R

update.ref.geno <- function(ref.geno, exc.snps){
  
  if(length(exc.snps) == 0){
    return(ref.geno)
  }
  
  id <- which(!(colnames(ref.geno) %in% exc.snps))
  if(length(id) == 0){
    msg <- "All SNPs excluded in update.ref.geno"
    stop(msg)
  }
  
  ref.geno <- ref.geno[, id, drop = FALSE]
  ref.geno
  
}

Try the ARTP2 package in your browser

Any scripts or data that you put into this service are public.

ARTP2 documentation built on May 2, 2019, 3:38 p.m.