haplo.stats

library(haplo.stats)

data(hla.demo); attach(hla.demo)

geno<-hla.demo[,c(17,18,21:24)]
label<-c("DQB","DRB","B")

geno.desc<-summaryGeno(geno, miss.val=c(0,NA))
print(geno.desc[c(1:10,80:85,135:140),])

table(geno.desc[,3])
miss.all<-which(geno.desc[,3]==3)
hla.demo.updated<-hla.demo[-miss.all,]

save.em<-haplo.em(geno=geno,locus.label=label,miss.val=c(0,NA))

summary(save.em, nlines=7)