# homemade Venn code
KK<-readLines("KaengKrachan_babblers.txt")
KY<-readLines("KhaoYai_babblers.txt")
DI<-readLines("DoiInthanon_babblers.txt")


plot(NULL, xlim=c(0,2), ylim=c(0,2), axes=F, xlab="", ylab="", main="Babbler overlap between three Thai N.P.s")
cx1<-0.7
cx2<-1.3
cx3<-1
cy1<-0.8
cy2<-0.8
cy3<-0.8+(((cx2-cx1)*3^0.5)/2)

theta<-seq(0,2*pi,0.01) 
radius<-0.5
x<-radius*cos(theta) 
y<-radius*sin(theta)
lines(x+cx1,y+cy1,lwd=3,col= "blue")
lines(x+cx3,y+cy3,lwd=3,col= "darkgreen")
lines(x+cx2,y+cy2,lwd=3,col= "red")

text(0.7,0.22,"Kaeng Krachan N.P.",col="blue")
text(0.45,1.8,"Khao Yai N.P.", col="darkgreen")
text(cy3, 0.22,"Doi Inthanon N.P.", col="red")

all<-intersect(intersect(KK,KY), DI)
text(1, 1.04, length(all)) 

kk_ky_only<-length(intersect(KK,KY))-length(all)
kk_di_only<-length(intersect(KK,DI))-length(all)
di_ky_only<-length(intersect(DI,KY))-length(all)
text(1, 0.65, kk_di_only)
text(0.75, 1.1, kk_ky_only)
text(1.25, 1.1, di_ky_only)

kk<-length(KK)-length(all)-kk_ky_only-kk_di_only
ky<-length(KY)-length(all)-di_ky_only-kk_ky_only
di<-length(DI)-length(all)-di_ky_only-kk_di_only

text(0.6,0.65,kk)
text(1.4,0.65,ky)
text(1,1.4,di)