# make a "hlaAlleleClass" object
hla.id <- "A"
hla <- hlaAllele(HLA_Type_Table$sample.id,
HLA_Type_Table[, paste(hla.id, ".1", sep="")],
HLA_Type_Table[, paste(hla.id, ".2", sep="")],
locus=hla.id, assembly="hg19")
summary(hla)
hlaSampleAllele(hla)
hlaSampleAllele(hla, allele.limit=c(
"01:01","02:01","02:06", "03:01", "11:01", "23:01"))
Run the code above in your browser using DataLab