Análise de variância ANOVA para experimento em DIC - Exemplo 01, R

Análise de variância ANOVA para experimento em DIC - Exemplo 01

#R version 3.4.4 
  
print("Scrip para fazer anova no R online Rextester")
# inseir os tratamentos
t1=c("racao1","racao1","racao1","racao1","racao1")
peso1=c(2,2,3,3,5)
t2=c("racao2","racao2","racao2","racao2","racao2")
t3=c("racao3","racao3","racao3","racao3","racao3")
t4=c("racao4","racao4","racao4","racao4","racao4")
t5=c("racao5","racao5","racao5","racao5","racao5")
# criar um vetor para tratamentos
trat=c(t1,t2,t3,t4,t5)
peso1=c(2,2,3,3,5)
peso2=c(3,3,5,5,5)
peso3=c(5,5,5,7,5)
peso4=c(4,7,7,7,8)
peso5=c(8,8,9,9,9)
#criar um vetor para peso
peso=c(peso1,peso2,peso3,peso4,peso5)
#criar um data.frame para analisar a atividade
# iremos chamar de ex1 (exercío 1)
ex1=data.frame(trat,peso); ex1
# vamos testar a normalidade dos tratamentos
# vamos usar a função tapply para separar o teste de shapiro para cada tratamento
shapiro=tapply(peso,trat,shapiro.test); shapiro
# vamos testar a homocedasticidade das variâncias
ex1bartlett=bartlett.test(ex1$peso~ex1$trat); ex1bartlett
# vamos realizar a análise de variância
ex1aov=aov(ex1$peso~ex1$trat)
# vamos solicitar o resultado da análise de variância
summary(ex1aov)
# vamos fazer o teste de Tukey
ex1tukey=TukeyHSD(ex1aov); ex1tukey
# vamos fazer um gráfico boxplot
ex1boxplot=boxplot(ex1$peso~ex1$trat); ex1boxplot
# vamos fazer um gráfico de barras
# para isso vamos criar um vetor com as média usando novamente a função tapply
medias=tapply(ex1$peso,ex1$trat,mean); medias
# agora vamos criar o gráfico de barras usando as médias 
ex1barplot=barplot(medias); ex1barplot


run \| edit \| history \| help	0

λ

                                      .NET NoSQL database for rapid development