A seguir ajusta-se o modelo e inspecionam-se os resíduos.
##Anova
tr.av <- aov(resp ~ trat, data=tr)
tr.av
anova(tr.av)
##Homocedasticidade
plot.default(tr$trat,tr.av$res)
##Normalidade
qqnorm(tr.av$res,ylab="Residuos", main=NULL)
qqline(tr.av$res)
title("Grafico Normal de Probabilidade dos Resíduos")
O gráfico de resíduos vs valores preditos mostra claramente uma heterogeneidade de variâncias e o
mostra um comportamento dos dados que se afasta muito da distribuição normal. A menssagem
é clara mas testes podem ser feitos para verificar o desvio dos pressupostos.
bartlett.test(tr$resp, tr$trat)
Bartlett test for homogeneity of variances
data: tr$resp and tr$trat
Bartlett's K-squared = 29.586, df = 4, p-value = 5.942e-06
shapiro.test(tr.av$res)
Shapiro-Wilk normality test
data: tr.av$res
W = 0.8961, p-value = 0.006742
Nos resultados acima a homogeneidade de variâncias foi rejeitada e também a normalidade dos resíduos.