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

Marcelo Laia marcelolaia em gmail.com
Qui Mar 28 12:18:23 -03 2019


Na minha busca por informação (para o meu conhecimento, mesmo) a respeito do
tema, recebi a inestimável ajuda de vocês e também encontrei algumas coisas na
internet.

Daquilo que encontrei na net, eis um artigo que achei interessante:

https://www.r-bloggers.com/normality-tests-don%E2%80%99t-do-what-you-think-they-do/

ou

https://wp.me/pMm6L-fpD

Abraços a todos.

Laia ML



> > 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

-- 
Marcelo


Mais detalhes sobre a lista de discussão R-br