1 - Efeito da Contaminação por metais pesados em invertebrados
1.1 - Importando os dados para o Experimento
sesseis <- read.csv(file = "sesseis.csv", header = TRUE)
str(sesseis)
## 'data.frame': 60 obs. of 3 variables:
## $ riqueza : int 68 64 64 63 69 63 70 68 68 62 ...
## $ cobre : Factor w/ 3 levels "alta","baixa",..: 3 3 3 3 3 3 3 3 3 3 ...
## $ orientacao: Factor w/ 2 levels "horizontal","vertical": 2 2 2 2 2 2 2 2 2 2 ...
1.2 - Estimação do Modelo
sesseis.aov <- aov(riqueza ~ cobre*orientacao, data = sesseis)
summary(sesseis.aov)
## Df Sum Sq Mean Sq F value Pr(>F)
## cobre 2 3330 1665.0 192.53 < 0.0000000000000002 ***
## orientacao 1 240 240.0 27.75 0.000002463599 ***
## cobre:orientacao 2 571 285.4 33.00 0.000000000434 ***
## Residuals 54 467 8.6
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
1.3 - Teste de Tukey
TukeyHSD(sesseis.aov)
## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = riqueza ~ cobre * orientacao, data = sesseis)
##
## $cobre
## diff lwr upr p adj
## baixa-alta 7.95 5.708825 10.19117 0
## nenhuma-alta 18.20 15.958825 20.44117 0
## nenhuma-baixa 10.25 8.008825 12.49117 0
##
## $orientacao
## diff lwr upr p adj
## vertical-horizontal -4 -5.522314 -2.477686 0.0000025
##
## $`cobre:orientacao`
## diff lwr upr
## baixa:horizontal-alta:horizontal 0.4 -3.48559529 4.2855953
## nenhuma:horizontal-alta:horizontal 14.2 10.31440471 18.0855953
## alta:vertical-alta:horizontal -11.7 -15.58559529 -7.8144047
## baixa:vertical-alta:horizontal 3.8 -0.08559529 7.6855953
## nenhuma:vertical-alta:horizontal 10.5 6.61440471 14.3855953
## nenhuma:horizontal-baixa:horizontal 13.8 9.91440471 17.6855953
## alta:vertical-baixa:horizontal -12.1 -15.98559529 -8.2144047
## baixa:vertical-baixa:horizontal 3.4 -0.48559529 7.2855953
## nenhuma:vertical-baixa:horizontal 10.1 6.21440471 13.9855953
## alta:vertical-nenhuma:horizontal -25.9 -29.78559529 -22.0144047
## baixa:vertical-nenhuma:horizontal -10.4 -14.28559529 -6.5144047
## nenhuma:vertical-nenhuma:horizontal -3.7 -7.58559529 0.1855953
## baixa:vertical-alta:vertical 15.5 11.61440471 19.3855953
## nenhuma:vertical-alta:vertical 22.2 18.31440471 26.0855953
## nenhuma:vertical-baixa:vertical 6.7 2.81440471 10.5855953
## p adj
## baixa:horizontal-alta:horizontal 0.9996312
## nenhuma:horizontal-alta:horizontal 0.0000000
## alta:vertical-alta:horizontal 0.0000000
## baixa:vertical-alta:horizontal 0.0587141
## nenhuma:vertical-alta:horizontal 0.0000000
## nenhuma:horizontal-baixa:horizontal 0.0000000
## alta:vertical-baixa:horizontal 0.0000000
## baixa:vertical-baixa:horizontal 0.1186003
## nenhuma:vertical-baixa:horizontal 0.0000000
## alta:vertical-nenhuma:horizontal 0.0000000
## baixa:vertical-nenhuma:horizontal 0.0000000
## nenhuma:vertical-nenhuma:horizontal 0.0705255
## baixa:vertical-alta:vertical 0.0000000
## nenhuma:vertical-alta:vertical 0.0000000
## nenhuma:vertical-baixa:vertical 0.0000651