Cuadrado Latino
Luis Fernando Delgado
2026-10-09
Cuadrado Latino
El cuadrado latino supone la existencia de dos fuentes sistemáticas extras diferentes a los tratamientos que deben controlarse mediante bloqueo bidireccional. Las fuentes extras se denominan genéricamente “Hileras y”Columnas”.
Debido a que cada tratamiento debe quedar representado una vez en cada hilera y columna; el numero de tratamiento (t) es igual al numero de columnas (c) y de hileras (h), por ende el numero de unidades experimentales requeridas es 𝑡^2.
Carga de librerias
if (!requireNamespace("tidyverse", quietly = TRUE)) install.packages("tidyverse")
if (!requireNamespace("FielDHub", quietly = TRUE)) install.packages("FielDHub")
if (!requireNamespace("agricolae", quietly = TRUE)) install.packages("agricolae")
library(FielDHub)
library(tidyverse)## ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
## ✔ dplyr 1.2.1 ✔ readr 2.2.0
## ✔ forcats 1.0.1 ✔ stringr 1.6.0
## ✔ ggplot2 4.0.3 ✔ tibble 3.3.1
## ✔ lubridate 1.9.5 ✔ tidyr 1.3.2
## ✔ purrr 1.2.2
## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ dplyr::filter() masks stats::filter()
## ✖ dplyr::lag() masks stats::lag()
## ℹ Use the conflicted package (<http://conflicted.r-lib.org/>) to force all conflicts to become errorslibrary(agricolae)Ejecucion de FieldHub
run_app()
Generar un diseno experimental - Cuadrado Latino
latin_1 <- latin_square(t = 10,
reps = 1,
locationName = "Buga",
seed = 12345)
latin_1$fieldBook## ID LOCATION PLOT SQUARE ROW COLUMN TREATMENT
## 1 1 Buga 101 1 Row 1 Column 1 T7
## 2 2 Buga 102 1 Row 1 Column 2 T3
## 3 3 Buga 103 1 Row 1 Column 3 T5
## 4 4 Buga 104 1 Row 1 Column 4 T4
## 5 5 Buga 105 1 Row 1 Column 5 T9
## 6 6 Buga 106 1 Row 1 Column 6 T10
## 7 7 Buga 107 1 Row 1 Column 7 T6
## 8 8 Buga 108 1 Row 1 Column 8 T1
## 9 9 Buga 109 1 Row 1 Column 9 T8
## 10 10 Buga 110 1 Row 1 Column 10 T2
## 11 11 Buga 120 1 Row 2 Column 1 T8
## 12 12 Buga 119 1 Row 2 Column 2 T9
## 13 13 Buga 118 1 Row 2 Column 3 T6
## 14 14 Buga 117 1 Row 2 Column 4 T2
## 15 15 Buga 116 1 Row 2 Column 5 T5
## 16 16 Buga 115 1 Row 2 Column 6 T1
## 17 17 Buga 114 1 Row 2 Column 7 T7
## 18 18 Buga 113 1 Row 2 Column 8 T10
## 19 19 Buga 112 1 Row 2 Column 9 T3
## 20 20 Buga 111 1 Row 2 Column 10 T4
## 21 21 Buga 121 1 Row 3 Column 1 T9
## 22 22 Buga 122 1 Row 3 Column 2 T8
## 23 23 Buga 123 1 Row 3 Column 3 T1
## 24 24 Buga 124 1 Row 3 Column 4 T7
## 25 25 Buga 125 1 Row 3 Column 5 T4
## 26 26 Buga 126 1 Row 3 Column 6 T5
## 27 27 Buga 127 1 Row 3 Column 7 T3
## 28 28 Buga 128 1 Row 3 Column 8 T6
## 29 29 Buga 129 1 Row 3 Column 9 T2
## 30 30 Buga 130 1 Row 3 Column 10 T10
## 31 31 Buga 140 1 Row 4 Column 1 T5
## 32 32 Buga 139 1 Row 4 Column 2 T4
## 33 33 Buga 138 1 Row 4 Column 3 T7
## 34 34 Buga 137 1 Row 4 Column 4 T3
## 35 35 Buga 136 1 Row 4 Column 5 T1
## 36 36 Buga 135 1 Row 4 Column 6 T6
## 37 37 Buga 134 1 Row 4 Column 7 T10
## 38 38 Buga 133 1 Row 4 Column 8 T2
## 39 39 Buga 132 1 Row 4 Column 9 T9
## 40 40 Buga 131 1 Row 4 Column 10 T8
## 41 41 Buga 141 1 Row 5 Column 1 T1
## 42 42 Buga 142 1 Row 5 Column 2 T10
## 43 43 Buga 143 1 Row 5 Column 3 T4
## 44 44 Buga 144 1 Row 5 Column 4 T8
## 45 45 Buga 145 1 Row 5 Column 5 T7
## 46 46 Buga 146 1 Row 5 Column 6 T2
## 47 47 Buga 147 1 Row 5 Column 7 T9
## 48 48 Buga 148 1 Row 5 Column 8 T5
## 49 49 Buga 149 1 Row 5 Column 9 T6
## 50 50 Buga 150 1 Row 5 Column 10 T3
## 51 51 Buga 160 1 Row 6 Column 1 T10
## 52 52 Buga 159 1 Row 6 Column 2 T1
## 53 53 Buga 158 1 Row 6 Column 3 T9
## 54 54 Buga 157 1 Row 6 Column 4 T6
## 55 55 Buga 156 1 Row 6 Column 5 T2
## 56 56 Buga 155 1 Row 6 Column 6 T8
## 57 57 Buga 154 1 Row 6 Column 7 T4
## 58 58 Buga 153 1 Row 6 Column 8 T3
## 59 59 Buga 152 1 Row 6 Column 9 T5
## 60 60 Buga 151 1 Row 6 Column 10 T7
## 61 61 Buga 161 1 Row 7 Column 1 T2
## 62 62 Buga 162 1 Row 7 Column 2 T5
## 63 63 Buga 163 1 Row 7 Column 3 T3
## 64 64 Buga 164 1 Row 7 Column 4 T9
## 65 65 Buga 165 1 Row 7 Column 5 T8
## 66 66 Buga 166 1 Row 7 Column 6 T4
## 67 67 Buga 167 1 Row 7 Column 7 T1
## 68 68 Buga 168 1 Row 7 Column 8 T7
## 69 69 Buga 169 1 Row 7 Column 9 T10
## 70 70 Buga 170 1 Row 7 Column 10 T6
## 71 71 Buga 180 1 Row 8 Column 1 T6
## 72 72 Buga 179 1 Row 8 Column 2 T2
## 73 73 Buga 178 1 Row 8 Column 3 T8
## 74 74 Buga 177 1 Row 8 Column 4 T10
## 75 75 Buga 176 1 Row 8 Column 5 T3
## 76 76 Buga 175 1 Row 8 Column 6 T7
## 77 77 Buga 174 1 Row 8 Column 7 T5
## 78 78 Buga 173 1 Row 8 Column 8 T9
## 79 79 Buga 172 1 Row 8 Column 9 T4
## 80 80 Buga 171 1 Row 8 Column 10 T1
## 81 81 Buga 181 1 Row 9 Column 1 T3
## 82 82 Buga 182 1 Row 9 Column 2 T6
## 83 83 Buga 183 1 Row 9 Column 3 T2
## 84 84 Buga 184 1 Row 9 Column 4 T1
## 85 85 Buga 185 1 Row 9 Column 5 T10
## 86 86 Buga 186 1 Row 9 Column 6 T9
## 87 87 Buga 187 1 Row 9 Column 7 T8
## 88 88 Buga 188 1 Row 9 Column 8 T4
## 89 89 Buga 189 1 Row 9 Column 9 T7
## 90 90 Buga 190 1 Row 9 Column 10 T5
## 91 91 Buga 200 1 Row 10 Column 1 T4
## 92 92 Buga 199 1 Row 10 Column 2 T7
## 93 93 Buga 198 1 Row 10 Column 3 T10
## 94 94 Buga 197 1 Row 10 Column 4 T5
## 95 95 Buga 196 1 Row 10 Column 5 T6
## 96 96 Buga 195 1 Row 10 Column 6 T3
## 97 97 Buga 194 1 Row 10 Column 7 T2
## 98 98 Buga 193 1 Row 10 Column 8 T8
## 99 99 Buga 192 1 Row 10 Column 9 T1
## 100 100 Buga 191 1 Row 10 Column 10 T9latin_1$infoDesign## $treatments
## [1] 10
##
## $squares
## [1] 1
##
## $locationName
## [1] "Buga"
##
## $seed
## [1] 12345
##
## $id_design
## [1] 3latin_1$squares## $rep1
## Column 1 Column 2 Column 3 Column 4 Column 5 Column 6 Column 7 Column 8
## Row 1 "T7" "T3" "T5" "T4" "T9" "T10" "T6" "T1"
## Row 2 "T8" "T9" "T6" "T2" "T5" "T1" "T7" "T10"
## Row 3 "T9" "T8" "T1" "T7" "T4" "T5" "T3" "T6"
## Row 4 "T5" "T4" "T7" "T3" "T1" "T6" "T10" "T2"
## Row 5 "T1" "T10" "T4" "T8" "T7" "T2" "T9" "T5"
## Row 6 "T10" "T1" "T9" "T6" "T2" "T8" "T4" "T3"
## Row 7 "T2" "T5" "T3" "T9" "T8" "T4" "T1" "T7"
## Row 8 "T6" "T2" "T8" "T10" "T3" "T7" "T5" "T9"
## Row 9 "T3" "T6" "T2" "T1" "T10" "T9" "T8" "T4"
## Row 10 "T4" "T7" "T10" "T5" "T6" "T3" "T2" "T8"
## Column 9 Column 10
## Row 1 "T8" "T2"
## Row 2 "T3" "T4"
## Row 3 "T2" "T10"
## Row 4 "T9" "T8"
## Row 5 "T6" "T3"
## Row 6 "T5" "T7"
## Row 7 "T10" "T6"
## Row 8 "T4" "T1"
## Row 9 "T7" "T5"
## Row 10 "T1" "T9"Impresion del layout
plot(latin_1)
Carga de datos descargados por la aplicacion
cuadrado <- read_csv("Latin_Square_2024-04-25.csv") %>%
select(ROW_SQ, COLUMN_SQ, TREATMENT, YIELD)## Rows: 25 Columns: 10
## ── Column specification ────────────────────────────────────────────────────────
## Delimiter: ","
## chr (4): LOCATION, ROW_SQ, COLUMN_SQ, TREATMENT
## dbl (6): ID, PLOT, ROW, COLUMN, SQUARE, YIELD
##
## ℹ Use `spec()` to retrieve the full column specification for this data.
## ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.str(cuadrado)## tibble [25 × 4] (S3: tbl_df/tbl/data.frame)
## $ ROW_SQ : chr [1:25] "Row 1" "Row 1" "Row 1" "Row 1" ...
## $ COLUMN_SQ: chr [1:25] "Column 1" "Column 2" "Column 3" "Column 4" ...
## $ TREATMENT: chr [1:25] "T2" "T4" "T3" "T1" ...
## $ YIELD : num [1:25] 49.6 33.4 66.4 78 40.4 ...Convertir a factor las filas, columnas y tratamientos
cuadrado$ROW_SQ <- as.factor(cuadrado$ROW_SQ)
cuadrado$COLUMN_SQ <- as.factor(cuadrado$COLUMN_SQ)
cuadrado$TREATMENT <- as.factor(cuadrado$TREATMENT)Grafico de tratamientos
cuadrado %>% ggplot(aes(x = TREATMENT, y = YIELD)) +
geom_boxplot()
Grafico de filas
cuadrado %>% ggplot(aes(x = ROW_SQ, y = YIELD)) +
geom_boxplot()
Grafico de columnas
cuadrado %>% ggplot(aes(x = COLUMN_SQ, y = YIELD)) +
geom_boxplot()
Modelacion
model1 <- aov(YIELD ~ TREATMENT + COLUMN_SQ + ROW_SQ, data = cuadrado)
summary(model1)## Df Sum Sq Mean Sq F value Pr(>F)
## TREATMENT 4 6604 1651.1 153.491 3.42e-10 ***
## COLUMN_SQ 4 10 2.6 0.244 0.908
## ROW_SQ 4 59 14.8 1.372 0.301
## Residuals 12 129 10.8
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1Prueba de medias con TUKEY
medias <- HSD.test(model1, "TREATMENT", alpha = 0.05)
medias$groups## YIELD groups
## T1 79.626 a
## T3 66.736 b
## T2 51.694 c
## T5 44.240 d
## T4 33.782 e