R/update.raw.geno.R

In zhangh12/ARTP3: Pathway and Gene-Level Association Test

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