Análise de Variâncias DIC e teste de Tukey 1, R

Análise de Variâncias DIC e teste de Tukey 1

#R version 3.4.4

#Análise de Variâncias DIC e teste de Tukey 1

# Análise de Variâncias para experimento em DIC e teste de Tukey  
# Criar os vetores para valores cada resultado (peso de camarões)
peso = c (21,23,26,21,19,
    31,26,27,26,33,
    41,43,44,41,36,
    54,56,57,58,56)

# Criar vetor alfa numérico para cada tratamento-repetição

trat = c ("ração1", "ração1", "ração1", "ração1", "ração1",
    "ração2", "ração2", "ração2", "ração2", "ração2",
    "ração3", "ração3", "ração3", "ração3", "ração3",
    "ração4", "ração4", "ração4", "ração4", "ração4")

# Criar um vetor em formato data.frame
# nome desse vetor eu chamarei de ex1 (Que será o exemplo 1)
ex1 = data.frame (trat,peso)

# Verificar a normalidade dos dados dos tratamentos
normalidade = tapply(ex1$peso, ex1$trat, shapiro.test)

# Verificar a homocedasticidade das variâncias
homocedasticidade =  bartlett.test(ex1$peso, ex1$trat)

# Realizar a ANOVA
ex1_anova = aov (ex1$peso ~ ex1$trat)

# Para verificar a diferença entre os tratamentos
# Vamos utilizar o teste de Tukey
ex1_tukey = TukeyHSD(ex1_anova)

print("data.frame")
ex1
print("normalidade dos dados - teste de Shapiro Wilk")
normalidade
print("homocedasticidade das variâncias - teste de Bartlett")
homocedasticidade
print("análise de variância - ANOVA")
summary(ex1_anova)
print("Teste de comparação de médias - Teste de Tukey")
ex1_tukey
print("**** FINAL DA ANÁLISE ****")


run \| edit \| history \| help	0