data(db2)
#Checks that frequencies add to 1
lapply(split(db2$Frequency,db2$Marker),sum)
#Finds number of alleles for all markers
unlist(lapply(split(db2$Frequency,db2$Marker),length))
#A closer look at the marker TH01
TH01=db2[db2$Marker=="TH01",]
barplot(TH01$Frequency)
Run the code above in your browser using DataLab