require(dplyr)
require(ggplot2)
require(tidyr)
require(reshape2)
require(lemon)
require(plotly)
require(data.table)
require(geobr)
require(terra)
require(openxlsx)
require(scales)
require(viridis)

estacoes = read.xlsx('Data/stations.xlsx', detectDates = T)

a = estacoes %>% 
  select(c('state', "chla_ugL",'pheophytin_ugL',  'phycocianin_ugL', 
           'tss_mgL','tso_mgL', 'tsi_mgL', 'secchi_m','dtc_mgL', 'doc_mgL','dic_mgL', 'turbidity_NTU', 'acdom440'))

a$chla_ugL =  (as.numeric(a$chla_ugL))
a$pheophytin_ugL =  (as.numeric(a$pheophytin_ugL))
a$phycocianin_ugL =  (as.numeric(a$phycocianin_ugL))
a$tss_mgL =   (as.numeric(a$tss_mgL))
a$tsi_mgL =   (as.numeric(a$tsi_mgL))
a$tso_mgL =   (as.numeric(a$tso_mgL))

a$acdom440 =  (as.numeric(a$acdom440))
a$secchi_m = (as.numeric(a$secchi_m))
a$dtc_mgL = (as.numeric(a$dtc_mgL))
a$dic_mgL = (as.numeric(a$dic_mgL))

a$doc_mgL = (as.numeric(a$doc_mgL))
a$turbidity_NTU = (as.numeric(a$turbidity_NTU))



library(tidyverse)
library(GGally)


names(a) = c("estado_sigla", "Chl-a", 'Phaeophytin', 'Phycocyanin', "TSS", "TSO", "TSI", "Zsd", "DTC", "DOC", 'DIC', "Turb" ,"aCDOM(440)")

melted = melt(a) %>% na.omit()

melted = filter(melted, value < 10000)

plt = ggplot(melted, aes(x = estado_sigla, y = value, fill = estado_sigla)) +
  geom_violin(trim=FALSE) +
  geom_boxplot(width=0.1, color="grey", alpha=0.2) +
  scale_y_log10(labels = label_number()) +
  
  labs(x = "State", y = 'Value') +
  facet_wrap(~variable, scale = 'free_y', ncol = 2) +
  scale_fill_viridis(discrete = TRUE) +
  
  theme_bw() + 
  theme(panel.grid.major = element_line(colour = "#d3d3d3"),
        panel.grid.minor = element_blank(),
        panel.border = element_blank(),
        panel.background = element_blank(),
        text=element_text(family = "Tahoma", size = 30),
        axis.title = element_text( size = 30),
        axis.text.x = element_text(colour="black", size = 25),
        axis.text.y = element_text(colour="black", size = 30, angle = 0),
        axis.line = element_line(size=2, colour = "black"),
        strip.text = element_text(size=30)) +
  #legend.title = element_text(size=20), #change legend title font size
  #legend.key.size = unit(2, 'cm')
  #change legend text font size) +
  theme(plot.margin = unit(c(6,6,6,6), "lines"))  + 
  guides(color = guide_legend(override.aes = list(size=5))) +
  theme(legend.position="none")


ggsave(plot = plt, filename = 'Outputs/Figures/Figure_07.jpeg', 
       width = 30, height = 35, dpi = 200, units = 'in')