setwd("/Volumes/GoogleDrive/Mi unidad/Agrosavia/Informes/Informe 8")
info7 <- read.table("curv_1.csv", header=T, sep=",")
attach(info7)
library(nlme) #activates the nlme library
library(ggplot2)
## pH loess-time-course figure
pht<- ggplot(info7, aes(x = time.num)) +
  facet_grid(gen~.) +
  geom_point(aes(y = ph.testa, colour = "ph.testa")) +
  geom_smooth(aes(y = ph.testa, colour = "ph.testa"), method = "loess", se=T) +
  scale_y_continuous(name = "pH") +  # Etiqueta de la variable continua
  scale_x_continuous(name = "day", breaks=seq(0,7,1)) + # Etiqueta de los grupos
  ggtitle("pH time course") +   # Título del plot
  theme(axis.line = element_line(colour = "black", # Personalización del tema
                                 size = 0.25))
pht + geom_smooth(aes(y = ph.grano, colour = "ph.grano"), method = "loess", se=T) +
  geom_point(aes(y=ph.grano, colour = "ph.grano"))
## `geom_smooth()` using formula 'y ~ x'
## `geom_smooth()` using formula 'y ~ x'

## Acidity loess-time-course figure
act1<- ggplot(info7, aes(x = time.num)) +
  facet_grid(gen~.) +
  geom_point(aes(y=acidez.testa, colour = "acidez.testa")) +
  geom_smooth(aes(y = acidez.testa, colour = "acidez.testa"), method = "loess", se=T) +
  scale_y_continuous(name = "Total acidity") +  # Etiqueta de la variable continua
  scale_x_continuous(name = "day", breaks=seq(0,7,1)) + # Etiqueta de los grupos
  ggtitle("Acidity time course") +   # Título del plot
  theme(axis.line = element_line(colour = "black", # Personalización del tema
                                 size = 0.25))
act1 + geom_smooth(aes(y = acidez.grano, colour = "acidez.grano"), method = "loess", se=T) +
  geom_point(aes(y=acidez.grano, colour = "acidez.grano"))
## `geom_smooth()` using formula 'y ~ x'
## `geom_smooth()` using formula 'y ~ x'