library(ggplot2)
install.packages("tidyquant")
library(tidyquant)
pira<-read.csv("piracena.csv", header = T)
pira$data<-as.Date(pira$data,"%m/%d/%Y")
pira$mes<-as.factor(pira$mes)
pira$ano<-as.factor(pira$ano)
pira$cond<-as.numeric(pira$cond)
str(pira)

pira<-group_by(pira,local, mes, ano)
piracena<-summarise(pira, OD = mean(OD_mg, na.rm=T),
                  Q = mean(Qdia, na.rm=T),
                  cond = mean(cond, na.rm=T),
                  resp = mean(resp, na.rm=T),
                  DIC = mean(DIC, na.rm=T),
                  NO3 = mean(NO3, na.rm=T))
piracena                  
write.csv(piracena, "piracena_mean_ano.csv")

pira<-read.csv("piracena_mean.csv", header =T)                  

gr1<-ggplot() + 
  geom_line(pira, mapping = aes (mes_number, OD), color = "red") +
  geom_point(pira, mapping = aes(mes_number,OD), color = "red", size = 0.8) +
  scale_x_continuous(name="Mês", breaks = seq(1,12 ,1), limits = c(1, 12)) +
  facet_wrap(~local) +
  theme_bw()
gr1

artemis<-filter(pira, local == "Copersucar")
gr2<-ggplot() + 
  geom_line(artemis, mapping = aes (mes_number, OD), color = "red") +
  geom_point(artemis, mapping = aes(mes_number,OD), color = "red", size = 0.8) +
  scale_x_continuous(name="Mês", breaks = seq(1,12 ,1), limits = c(1, 12), ) +
  facet_wrap(~local) +
  theme_bw()
gr2

atibaia<-filter(pira, local == "Desembargador")
gr3<-ggplot() + 
  geom_line(atibaia, mapping = aes (mes_number, OD), color = "red") +
  geom_point(atibaia, mapping = aes(mes_number,OD), color = "red", size = 0.8) +
  scale_x_continuous(name="Mês", breaks = seq(1,12 ,1), limits = c(1, 12)) +
  facet_wrap(~local) +
  theme_bw()
gr3

gr4<-ggplot() + 
  geom_point(pira, mapping = aes(Q,OD), color = "red", size = 2) +
  geom_smooth(pira, mapping = aes(Q,OD),method='lm', size = 0.5) + 
  facet_wrap(~local) +
  theme_bw()
gr4

gr5<-ggplot() + 
  geom_point(pira, mapping = aes(Q, cond), color = "red", size = 2) +
  geom_smooth(pira, mapping = aes(Q, cond),method='lm', size = 0.5) + 
  facet_wrap(~local) +
  theme_bw()
gr5

gr6<-ggplot() + 
  geom_point(pira, mapping = aes(Q, DIC), color = "red", size = 2) +
  geom_smooth(pira, mapping = aes(Q, DIC),method='lm', size = 0.5) + 
  facet_wrap(~local) +
  theme_bw()
gr6

gr7<-ggplot() + 
  geom_point(pira, mapping = aes(OD, DIC), color = "red", size = 2) +
  geom_smooth(pira, mapping = aes(OD, DIC),method='lm', size = 0.5) + 
  facet_wrap(~local) +
  theme_bw()
gr7

gr8<-ggplot() + 
  geom_point(pira, mapping = aes(Q, NO3), color = "red", size = 2) +
  geom_smooth(pira, mapping = aes(Q, NO3),method='lm', size = 0.5) + 
  facet_wrap(~local) +
  theme_bw()
gr8

gr9<-ggplot() + 
  geom_point(pira, mapping = aes(OD, NO3), color = "red", size = 2) +
  geom_smooth(pira, mapping = aes(OD, NO3),method='lm', size = 0.5) + 
  facet_wrap(~local) +
  theme_bw()
gr9