Interval.R

setwd("~/Google Drive/Agrosavia/Informes/Informe 2")
datos<-read.table("ph_regalias.csv", header=T, sep=',')
library(ggplot2)
library(Rmisc)

## Loading required package: lattice

## Loading required package: plyr

datos2<-summarySE (datos, measurevar = "testa", groupvars = c("día", "muestra"))
datos3<-summarySE (datos, measurevar = "almendra", groupvars = c("día", "muestra"))
datos2$ID <- seq.int(nrow(datos2))
datos3$ID <- seq.int(nrow(datos3))
total <- merge (datos2, datos3, by = "ID")

#Graficación por muestra
pht<- ggplot(total, aes(x = día.x)) +
  facet_grid(~muestra.x) +
  geom_errorbar(aes(ymin=testa-ci.x, ymax=testa+ci.x, colour="testa"), width=.1) +
  geom_line(aes(y=testa, colour="testa")) +
  geom_point(aes(y=testa, colour = "testa")) +
  scale_y_continuous(name = "pH") +  # Etiqueta de la variable continua
  scale_x_continuous(name = "día", breaks=seq(0,7,1)) + # Etiqueta de los grupos
  ggtitle("pH en el tiempo") +   # Título del plot
  theme(axis.line = element_line(colour = "black", # Personalización del tema
                                 size = 0.25))
pht+ geom_errorbar(aes(ymin=almendra-ci.y, ymax=almendra+ci.y, colour="almendra"), width=.1) +
  geom_line(aes(y=almendra, colour="almendra")) +
  geom_point(aes(y=almendra, colour = "almendra"))

# Graficación de toda la población
datos4<-summarySE (total, measurevar = "testa", groupvars = c("día.x"))
datos5<-summarySE (total, measurevar = "almendra", groupvars = c("día.x"))
datos4$ID <- seq.int(nrow(datos4))
datos5$ID <- seq.int(nrow(datos5))
total2 <- merge (datos4, datos5, by = "ID")

pht<- ggplot(total2, aes(x = día.x.x)) +
  geom_errorbar(aes(ymin=testa-ci.x, ymax=testa+ci.x, colour="testa"), width=.1) +
  geom_line(aes(y=testa, colour="testa")) +
  geom_point(aes(y=testa, colour = "testa")) +
  scale_y_continuous(name = "pH") +  # Etiqueta de la variable continua
  scale_x_continuous(name = "día", breaks=seq(0,7,1)) + # Etiqueta de los grupos
  ggtitle("pH en el tiempo") +   # Título del plot
  theme(axis.line = element_line(colour = "black", # Personalización del tema
                                 size = 0.25))
pht+ geom_errorbar(aes(ymin=almendra-ci.y, ymax=almendra+ci.y, colour="almendra"), width=.1) +
  geom_line(aes(y=almendra, colour="almendra")) +
  geom_point(aes(y=almendra, colour = "almendra"))