bkn <- read.table("~/Desktop/stats//bkn_stats_100.txt", header = T)
til <- read.table("~/Desktop/stats/til_stats_100.txt", header = T)
plot(density(bkn$tajd[which(bkn$exon_overlap == 100)]))
lines(density(til$tajd[which(bkn$exon_overlap == 100)]), col = "red")
plot(density(bkn$tajd[which(bkn$intron_overlap == 100)]))
lines(density(til$tajd[which(bkn$intron_overlap == 100)]), col = "red")
