[R-br] two-way anova - violação das pressuposições

[Marcelo Laia]

Colegas,

Estamos com um conjunto de dados que viola pressuposições da ANOVA,
principalmente a normalidade. Eu já li, por diversas vezes, que essa violação
pode não ser tão danosa assim, devido a robustez do teste.

Por favor, tire um tempo para ver os resultados e depois peço uma ajuda.

Genotipo e Isolado são fatores
Area está em centímetros quadrados

> leveneTest(Area ~ Genotipo*Isolado, data = cerato.desc)
Levene's Test for Homogeneity of Variance (center = median)
       Df F value    Pr(>F)    
group  41  3.4755 4.718e-08 ***
      126                      
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
> 

cat("Normality p-values by Factor Genotipo: ")
for (i in unique(factor(cerato.desc$Genotipo))){
  cat(shapiro.test(cerato.desc[cerato.desc$Genotipo==i, ]$Area)$p.value," ")
}
7.074459e-10  3.200422e-06

#Shapiro-Wilk normality tests by Isolado
for (i in unique(factor(cerato.desc$Isolado))){
  cat(shapiro.test(cerato.desc[cerato.desc$Isolado==i, ]$Area)$p.value," ")
}
0.09534117  0.4006495  0.6065291  0.2093362  0.6138097  0.5604402  0.1302976  0.3135567  0.905537  0.7294285  0.0966383  0.1512716  0.8469947  0.1226855  0.2713435  0.9695489  0.2747097  0.5476302  0.0008750702  0.03693436  0.4197769 

> bartlett.test(Area~Genotipo,data = cerato.desc )

	Bartlett test of homogeneity of variances

data:  Area by Genotipo
Bartlett's K-squared = 19.769, df = 1, p-value = 8.738e-06

> bartlett.test(Area~Isolado,data = cerato.desc )

	Bartlett test of homogeneity of variances

data:  Area by Isolado
Bartlett's K-squared = 171.26, df = 20, p-value < 2.2e-16

A pergunta é: mesmo com esses resultados eu poderia afirmar que o teste, neste
caso, será robusto o suficiente para essas violações?

Obrigado!

-- 
Marcelo