Walmes, segue abaixo o script que alterei para o látice parc. balanc.. Consta nele a ANOVA com tratamentos não-ajustados e bloco/repetição não-ajustados.
Acontece que eu tenho 5 experimentos em látice 6x6, cada um em 5 municípios na safra agrícola 2008-09. Como automatizar a análise para que eu utilize um mesmo .txt contendo os locais, tratamentos, blocos, repetições e a variável?
rm(list=ls())
ls()
library(lattice)
library(agricolae)
library(contrast)
library(multcomp)
library(doBy)
library(latticeExtra)
library(nlme)
#-----------------------------------------------------------------------
# lendo arquivos de dados
da <- read.table("C:/ADSR/Tese/teste.txt",
header = TRUE,
colClasses = rep(c("factor", "numeric"), c(3, 1)))
str(da)
#-----------------------------------------------------------------------
# ver os dados
xtabs(~repeticao + bloco, da)
xtabs(~repeticao + cultivar, da)
xyplot(producao ~ cultivar, groups = bloco, data = da)
xyplot(producao ~ cultivar | repeticao, groups = bloco, data = da)
#-----------------------------------------------------------------------
# ajuste do modelo de efeitos fixos para bloco dentro de repeticao
m0 <- lm(producao ~ repeticao/bloco + cultivar, data = da)
par(mfrow = c(2, 2))
plot(m0)
layout(1)
# #Teste de normalidade
# shapiro.test(residuals(m0))
# #Teste de Bartlett
# bartlett.test(residuals(m0)~cultivar, data=da)
anova(m0) # cultivar nao ajustado para blocos
#-----------------------------------------------------------------------
#Só arrumar ordem dos termos
m0 <- lm(terms(producao ~ repeticao/bloco + cultivar,
keep.order = TRUE),
data = da)
anova(m0)