head(HLA_Type_Table)
dim(HLA_Type_Table) # 60 13
# make a "hlaAlleleClass" object
hla.id <- "C"
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)
subhla1 <- hlaAlleleSubset(hla, 1:100)
summary(subhla1)
subhla2 <- hlaAlleleSubset(hla, 201:300)
summary(subhla2)
H <- hlaCombineAllele(subhla1, subhla2)
summary(H)
Run the code above in your browser using DataLab