allelenames <- paste("A",11:13,sep="")
allele1 <- factor(c("A11","A11","A12","A12","A13","A12"),levels=allelenames)
allele2 <- factor(c("A11","A12","A12","A13","A13","A11"),levels=allelenames)
GC <- table(allele1,allele2)
GC <- as.matrix(unclass(GC))
GCf <- fold(GC)
Run the code above in your browser using DataLab