Análise de resíduos

Com as linhas de comando abaixo, pode-se efetuar a análise de resíduos do experimento em faixas.

## esse modelo armazena os resíduos no objeto ex10.av1

ex10.av1<- aov(resp ~ bloco + bloco*espA+bloco*denB+espA*denB,data=ex10)

par(mfrow=c(2,2))
plot(ex10.av1)

shapiro.test(ex10.av1$res) # normalidade
plot(ex10$espA,ex10.av1$res) # homocedasticidade
plot(ex10$denB,ex10.av1$res) # homocedasticidade

# verificação de outliers

residuos <- resid(ex10.av1)
s2 <- sum(resid(ex10.av1)^2)/ex10.av1$df.res
 respad <- residuos/sqrt(s2)
 boxplot(respad)

outlier<-c(max(respad),min(respad))
outlier

Como nesse experimento não houve nenhum efeito significativo, não há necessidade de um teste de comparações múltiplas.