- Para este ejempo usaremos los siguientes paquetes:
library(agricolae)
library(desplot)
library(tidyverse)- Como ejemplo usaremos el data melon de agricolae
| row | col | variety | yield |
|---|---|---|---|
| 1 | 1 | V1 | 45 |
| 1 | 2 | V4 | 29 |
| 1 | 3 | V2 | 37 |
| 1 | 4 | V3 | 38 |
| 2 | 1 | V2 | 50 |
| 2 | 2 | V3 | 53 |
| 2 | 3 | V4 | 41 |
| 2 | 4 | V1 | 40 |
| 3 | 1 | V3 | 43 |
| 3 | 2 | V2 | 41 |
| 3 | 3 | V1 | 41 |
| 3 | 4 | V4 | 35 |
| 4 | 1 | V4 | 35 |
| 4 | 2 | V1 | 63 |
| 4 | 3 | V3 | 63 |
| 4 | 4 | V2 | 41 |
Este fue Una evaluación en cuadrado latino del sistema de riego por exudación utilizando cuatro variedades de melón, bajo la modalidad de siembra, FILA SIMPLE. El objetivo es analizar el comportamiento de tres variedades híbridas de melón y un estándar.
Graficamos
desplot(melon, variety ~ col+row,
main = "Diseño experimetal de 4 variedades de Melon",
col.regions = c('#8dd3c7','#ffffb3','#bebada','#fb8072'),
out1 = variety, cex = 1.5, show.key = F, text = variety, flip = T, ticks = T
)
- Una vez se termine el experimento es posible visualizar los resultados de la siguiente manera, graficando la variable de respuesta produccion
desplot(melon, yield ~ col+row,
main = "Resultados diseño experimetal de 4 variedades de Melon",
out1 = variety, cex = 1.5, show.key = F, flip = T, ticks = T,
xlab = "kg/parcela", ylab = ""
)
- Hacemos un ANOVA para ver los resultados en datos.
| term | df | sumsq | meansq | statistic | p.value |
|---|---|---|---|---|---|
| row | 1 | 227.8125 | 227.8125 | 3.7847323 | 0.0803536 |
| col | 1 | 46.5125 | 46.5125 | 0.7727292 | 0.4000040 |
| variety | 3 | 483.6875 | 161.2292 | 2.6785591 | 0.1037840 |
| Residuals | 10 | 601.9250 | 60.1925 | NA | NA |
Como se aprecia en la tabla de resultado del anova no hubo diferencias.
Grafiquemos otro ejemplo
desplot(cannabis,
block ~ col+row, main = "Diseño experimetal de 3 especies de Cannabis y 4 concentraciones de K",
col=K, text=Gen, cex=1, out1=block,
out2=Gen, out2.gpar=list(col = "black", lwd = 1, lty = 1),
col.regions = c('#8dd3c7','#ffffb3','#bebada','#fb8072','#80b1d3','#fdb462'))
Adaptacion de agridat::data(yates.oats)