Gráfico de líneas para experimento de larga duración

Gráficos de líneas creados en R, usando ggplot2 y plotly. Los datos se encuentran en la página 191 y el ejemplo en la 192 de Experimentación en Agricultura.

# Gráfico (opción 2)
library(plotly)

p <- plot_ly(olivoT, x = ~FECHA, y = ~AS, color= ~TRAT) %>% add_lines() %>%
  layout(title = "Gráfico de líneas",
         xaxis = list(title = "Años del Experimento"),
         yaxis = list (title = "asin(sqrt(Afectadas/100))*(180/pi)")); p

Códigos completos para cargar datos, transformarlos y graficarlos en R

#***********************************************************
# EXPERIMENTACIÓN EN AGRICULTURA                              
# Fernández Escobar, R.; Trapero, A.; Domínguez, J. 2010       
# CAPÍTULO 13: ANÁLISIS DE LA VARIANZA COMBINADO 
# ANÁLISIS DE EXPERIMENTOS DE LARGA DURACIÓN (páginas 190 - 192)
#************************************************************

#### Importar Tabla 13.5
olivo <- read.csv("C:/Users/Administrator/Desktop/tabla 13.6.csv", header=TRUE, sep=",")
attach(olivo); olivo

##    FECHA TRAT Afectadas
## 1      1   TC      0.20
## 2      2   TC     32.60
## 3      3   TC     57.30
## 4      4   TC     40.80
## 5      5   TC     16.90
## 6      6   TC      0.40
## 7      7   TC      5.10
## 8      8   TC     16.70
## 9      9   TC     26.90
## 10    10   TC     10.20
## 11    11   TC      6.30
## 12    12   TC      1.00
## 13     1   TD      0.16
## 14     2   TD     12.40
## 15     3   TD     18.00
## 16     4   TD     15.20
## 17     5   TD      2.50
## 18     6   TD      0.20
## 19     7   TD      2.30
## 20     8   TD      4.80
## 21     9   TD      8.10
## 22    10   TD      3.60
## 23    11   TD      1.80
## 24    12   TD      0.16

#### Creando variables independientes como factores
TRAT   <- factor(TRAT, levels = c('TC','TD'), labels = c("Tratamiento Convencional", "Tratamiento Dirigido"))
AS     <- asin(sqrt(Afectadas/100))*(180/pi) # datos ajustados a la normalidad mediante la transformación angular
olivoT <- data.frame(FECHA, TRAT, Afectadas, AS); olivoT

##    FECHA                     TRAT Afectadas        AS
## 1      1 Tratamiento Convencional      0.20  2.563200
## 2      2 Tratamiento Convencional     32.60 34.817483
## 3      3 Tratamiento Convencional     57.30 49.197596
## 4      4 Tratamiento Convencional     40.80 39.698582
## 5      5 Tratamiento Convencional     16.90 24.273708
## 6      6 Tratamiento Convencional      0.40  3.626123
## 7      7 Tratamiento Convencional      5.10 13.051796
## 8      8 Tratamiento Convencional     16.70 24.120456
## 9      9 Tratamiento Convencional     26.90 31.241880
## 10    10 Tratamiento Convencional     10.20 18.625095
## 11    11 Tratamiento Convencional      6.30 14.536577
## 12    12 Tratamiento Convencional      1.00  5.739170
## 13     1     Tratamiento Dirigido      0.16  2.292443
## 14     2     Tratamiento Dirigido     12.40 20.618039
## 15     3     Tratamiento Dirigido     18.00 25.104090
## 16     4     Tratamiento Dirigido     15.20 22.946521
## 17     5     Tratamiento Dirigido      2.50  9.097436
## 18     6     Tratamiento Dirigido      0.20  2.563200
## 19     7     Tratamiento Dirigido      2.30  8.722994
## 20     8     Tratamiento Dirigido      4.80 12.655533
## 21     9     Tratamiento Dirigido      8.10 16.535238
## 22    10     Tratamiento Dirigido      3.60 10.937416
## 23    11     Tratamiento Dirigido      1.80  7.710285
## 24    12     Tratamiento Dirigido      0.16  2.292443

attach(olivoT); str(olivoT)

# Gráfico (opción 1)
library(ggplot2)

ggplot(data = olivoT, aes(x        = FECHA, 
                          y        = AS, 
                          group    = TRAT, 
                          color    = TRAT, 
                          linetype = TRAT, 
                          shape    = TRAT ))  + 
  ggtitle("Gráfico de líneas")                + 
  xlab("Años del Experimento")                + 
  ylab("asin(sqrt(Afectadas/100))*(180/pi)")  + 
  geom_point(size = 5.0)                      + 
  geom_line(size  = 2.0, linetype = 'solid')  + 
  scale_color_brewer(palette = "Set2")

# Gráfico (opción 2)
library(plotly)

p <- plot_ly(olivoT, x = ~FECHA, y = ~AS, color= ~TRAT) %>% add_lines() %>%
  layout(title = "Gráfico de líneas",
         xaxis = list(title = "Años del Experimento"),
         yaxis = list (title = "asin(sqrt(Afectadas/100))*(180/pi)")); pGráfico de líneas

# No contamos con los datos completos de las 5 repeticiones, por eso no agregamos tabla de ANOVA

Datos y cógios adjuntos: tabla 13.6.csv, capítulo 13f.r