# source("THC.R")
# install.packages("cowplot")

library("cowplot")
library(ggplot2)
library(ggpubr)

thc<-read.table("THC.csv", dec = ".", header = TRUE, sep = "$", quote = "&", na.strings = "NA")

thc0<-thc[thc$tempo=="T0",]

thc1<-thc[thc$tempo=="T1",]

t1 <- ggboxplot(thc, "trattamento", "THC", ggtheme = theme_light(base_size=10), add = "jitter", xlab = "temperatura", ylab = "cells/µL")

t2 <- ggboxplot(thc, "tempo", "THC", ggtheme = theme_light(base_size=10), add = "jitter", xlab = "tempo (h)", ylab = "cells/µL")

t3 <- ggboxplot(thc0, "trattamento", "THC", ggtheme = theme_light(base_size=10), add = "jitter", xlab = "temperatura", ylab = "cells/µL")

t4 <- ggboxplot(thc1, "trattamento", "THC", ggtheme = theme_light(base_size=10), add = "jitter", xlab = "temperatura", ylab = "cells/µL")

plot_grid(t1, t2, t3, t4, labels = c("Tot", "Tot", "T0", "T1"), nrow = 2, ncol = 2, label_size = 12)

ggsave("THC.pdf", plot = last_plot(), device = "pdf")

# confronto tempo 0, parametrico
t.test(thc0$THC ~ thc0$trattamento)

# confronto tempo 1, parametrico
t.test(thc1$THC ~ thc1$trattamento)

# confronto tempo 0, non parametrico
wilcox.test(thc0$THC ~ thc0$trattamento)

# confronto tempo 1, non parametrico
wilcox.test(thc1$THC ~ thc1$trattamento)

# confronto tempi, parametrico
t.test(thc$THC ~ thc$tempo)

# confronto tempi, non parametrico
wilcox.test(thc$THC ~ thc$tempo)