script r

Esta é uma pré-visualização de arquivo. Entre para ver o arquivo original

#Limpando a memória
rm(list = ls())
#Carregandos os dados do pacote labestData
require(labestData)
dados = PimentelEg6.2
#Veriticando a estrutura do arquivo de dados importado
str(dados)
#Estatística Descritiva
summary(dados)
#
attach(dados)
#Nome das colunas do arquivo importado
names(dados)
#Boxplot
boxplot(prod~varied)
boxplot(prod~linha)
boxplot(prod~coluna)
########################################################
#1ª Forma de analisar os dados
#Ajustando o modelo
modelo = lm(prod~varied+linha+coluna)
modelo
#Summary do modelo
summary(modelo)
#Gráfico para analisar as pressuposicoes da ANAVA
par(mfrow=c(2,2))
plot(modelo)
par(mfrow=c(1,1))
#Teste de normalidade
shapiro.test(rstandard(modelo)) 
#ANAVA
anova(modelo)
#Teste de Tukey
require(TukeyC)
modelo.aov = aov(prod~varied+linha+coluna)
TukeyC(modelo.aov)
plot(TukeyC(modelo.aov))
#Método do SK
require(ScottKnott)
teste = SK(modelo.aov)
teste
plot(teste)
#
detach(dados)
########################################################
#2ª Forma de analisar os dados
attach(dados)
require(ExpDes.pt)
#Usando o teste de Tukey
dql(varied,linha,coluna,prod)
#Usando o teste de SK
dql(varied,linha,coluna,prod,mcomp = 'sk')