#UNIVERSIDAD NACIONAL DEL ALTIPLANO
#DISEæ¼ã¸±OS EXPERIMENTALES
#FINESI
#TEMA: Diseæ¼ã¸±o de Cuadrados Latinos
library(knitr)
## Warning: package 'knitr' was built under R version 3.6.3
# BUSCAR LA RUTA DE EXCEL
library(readxl)
## Warning: package 'readxl' was built under R version 3.6.3
DCL <- read_excel("G:/VII SEMESTRE/DISENOS EXPERIMENTALES/TAREA 6/DCL.xlsx")
DCL
## # A tibble: 36 x 4
## Observacion Lote Concentraciones Tiempo
## <dbl> <chr> <dbl> <chr>
## 1 12 Lote1 1 A
## 2 24 Lote1 2 B
## 3 10 Lote1 3 C
## 4 18 Lote1 4 D
## 5 21 Lote1 5 E
## 6 18 Lote1 6 F
## 7 21 Lote2 1 A
## 8 26 Lote2 2 B
## 9 24 Lote2 3 C
## 10 16 Lote2 4 D
## # ... with 26 more rows
View(DCL)
head(DCL)
## # A tibble: 6 x 4
## Observacion Lote Concentraciones Tiempo
## <dbl> <chr> <dbl> <chr>
## 1 12 Lote1 1 A
## 2 24 Lote1 2 B
## 3 10 Lote1 3 C
## 4 18 Lote1 4 D
## 5 21 Lote1 5 E
## 6 18 Lote1 6 F
#FACTORES
DCL$Lote <- factor(DCL$Lote)
DCL$Concentraciones <- factor(DCL$Concentraciones)
DCL$Tiempo <- factor(DCL$Tiempo)
#ANOVA
mod1 <- aov(Observacion ~ Lote + Concentraciones + Tiempo, data = DCL)
mod1
## Call:
## aov(formula = Observacion ~ Lote + Concentraciones + Tiempo,
## data = DCL)
##
## Terms:
## Lote Concentraciones Residuals
## Sum of Squares 99.5556 70.5556 464.4444
## Deg. of Freedom 5 5 25
##
## Residual standard error: 4.310195
## 5 out of 16 effects not estimable
## Estimated effects may be unbalanced
summary(mod1)
## Df Sum Sq Mean Sq F value Pr(>F)
## Lote 5 99.6 19.91 1.072 0.400
## Concentraciones 5 70.6 14.11 0.760 0.587
## Residuals 25 464.4 18.58