if(require("BSgenome.Hsapiens.UCSC.hg18")){
bsgenome <- BSgenome.Hsapiens.UCSC.hg18
snp_exp <- getExampleSnpExperiment(bsgenome)
is_het <- isHeterozygous(snp_exp[, 1], c(0.4, 0.6))
table(is_het)
}
Run the code above in your browser using DataLab