REPO = 'http://cran.xl-mirror.nl/'
#install.packages('vioplot', repos=REPO)
library(vioplot)

dat <- read.table("kings_similarities.tsv", sep = '\t' ,header = TRUE)
dat2 <- dat[dat$book_1 == '2_Kings',]
dat3 <- dat2[dat2$book_2 != '2_Kings',]

boxplot(dat3$sim, ylab = 'Similarity')
summary(dat3$sim)

vioplot(dat3$sim, names = '', col = 'white')

jer <- dat3[dat3$book_2 =='Jeremiah',]
jes <- dat3[dat3$book_2 =='Isaiah',]
chr <- dat3[dat3$book_2 =='2_Chronicles',]
boxplot(jes$sim, jer$sim, chr$sim, names = c('Isaiah', 'Jeremiah', '2 Chronicles'),  ylab = 'Similarity')

jer60 <- jer[jer$sim > 59,]
jes60 <- jes[jes$sim > 59,]
chr60 <- chr[chr$sim > 59,]
boxplot(jes60$sim, jer60$sim, chr60$sim, names = c('Isaiah', 'Jeremiah', '2 Chronicles'), ylab = 'Similarity')

jer80 <- jer[jer$sim > 79,]
jes80 <- jes[jes$sim > 79,]
chr80 <- chr[chr$sim > 79,]
boxplot(jes80$sim, jer80$sim, chr80$sim, names = c('Isaiah', 'Jeremiah', '2 Chronicles'), ylab = 'Similarity')

t.test(jes80$sim, jer80$sim)
t.test(jes80$sim, chr80$sim)
t.test(jer80$sim, chr80$sim)