Tarea - Varianzas desiguales

Modelación de cultivos - 2020-II - Universidad Nacional de Colombia

set.seed(2011)
gradosbrix<-replicate(1000,rnorm(48,22,0.8))
l1=list()
l2=list()
var<-gl(3,16,48,paste0("V",1:3))
for(i in 1:dim(gradosbrix)[2]){
  l1[[i]]=(aov(gradosbrix[,i]~var))$residuals
  l2[[i]]=shapiro.test(l1[[i]])$p.value
}
table(unlist(l2)>0.05) #Resultados de la prueba shapiro a los residuales

## 
## FALSE  TRUE 
##    35   965

#lapply(l1,shapiro.test)

gradosbrix<-c(runif(16,20,24),runif(16,15,35),runif(16,20,21)) #Varianzas desiguales 
var<-gl(3,16,48,paste0("V",1:3))
boxplot(gradosbrix~var,main="Varianzas desiguales",col=rainbow(3))

set.seed(2011)
gradosbrix1<-replicate(1000,rnorm(16,20,0.8))
gradosbrix2<-replicate(1000,rnorm(16,22,1.6))
gradosbrix3<-replicate(1000,rnorm(16,24,2.4))
gradosbrixT<-rbind(gradosbrix1,gradosbrix2,gradosbrix3)
var<-gl(3,16,48,paste0("V",1:3))
l1=list()
l2=list()
for(i in 1:dim(gradosbrixT)[2]){
  l1[[i]]=(aov(gradosbrixT[,i]~var))$residuals
  l2[[i]]=shapiro.test(l1[[i]])$p.value
}
table(unlist(l2)>0.05) #Resultados de la prueba shapiro a los residuales

## 
## FALSE  TRUE 
##   253   747

#lapply(l1,shapiro.test)