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
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.