Punto 1

library(readxl)
## Warning: package 'readxl' was built under R version 4.2.2
Parcial_1 <- read_excel("C:/Users/fcecursos/Downloads/Parcial 1.xlsx")
View(Parcial_1)
datos = Parcial_1

trt = datos$Tratamiento
sul = datos$Suelo
pend = datos$Pendiente
da = datos$Area
library(lattice)
library(ggplot2)
## Warning: package 'ggplot2' was built under R version 4.2.2
library(lattice)
bwplot(trt ~ da | sul * pend)

mod1 = aov(da ~ trt + sul + pend,
           datos)
summary(mod1)
##             Df Sum Sq Mean Sq F value Pr(>F)  
## trt          2  665.4   332.7  29.918 0.0323 *
## sul          2   28.5    14.3   1.282 0.4383  
## pend         2   14.1     7.1   0.634 0.6119  
## Residuals    2   22.2    11.1                 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
res_mod = mod1$residuals
shapiro.test(res_mod)
## 
##  Shapiro-Wilk normality test
## 
## data:  res_mod
## W = 0.81562, p-value = 0.03074
bartlett.test(res_mod,
              trt)
## 
##  Bartlett test of homogeneity of variances
## 
## data:  res_mod and trt
## Bartlett's K-squared = 7.2669e-16, df = 2, p-value = 1
plot(datos$Area | res_mod)

library(TukeyC)
## Warning: package 'TukeyC' was built under R version 4.2.3
tt = TukeyHSD(mod1, da)
## Warning in TukeyHSD.aov(mod1, da): 'which' specified some non-factors which
## will be dropped
plot(tt)