dados<- read.csv2("eucalipto.csv")
dados
str(dados)
mod.l1<- lm(altura ~ Bloco + Trat, data=dados, subset=c(Local=="L1"))
anova(mod.l1)
mod.l2<- lm(altura ~ Bloco + Trat, data=dados, subset=c(Local=="L2"))
anova(mod.l2)
mod.l3<- lm(altura ~ Bloco + Trat, data=dados, subset=c(Local=="L3"))
anova(mod.l3)
# razão entre os quadrados médios dos resíduos
QMResiduo1<- anova(mod.l1)$"Mean Sq"[3]
QMResiduo2<- anova(mod.l2)$"Mean Sq"[3]
QMResiduo3<- anova(mod.l3)$"Mean Sq"[3]
QMResiduo<- c(QMResiduo1, QMResiduo2, QMResiduo3)
max(QMResiduo)/min(QMResiduo)
#análise conjunta - exploratória
require(lattice)
xyplot(altura ~ Trat, group=Local, type="a", 
       auto.key=list(space = "right", points = FALSE, lines = TRUE),
       data=dados)
# teste F para efeito da interaçao Local:Trat
summary(aov(altura~Local+Local:Bloco+Trat+ Local:Trat, data=dados))
# teste F para efeito de Tratamentos
modelo <- aov(altura~Local+Local:Bloco+Trat+Error(Local:(Bloco+Trat)), data=dados)
summary(modelo)[[2]]
#desdobramento
dados$LT<- as.factor(dados$Local:dados$Trat)
#efeito de tratamento dentro de cada nível de local
mod.conj<- aov(altura ~ Local + Local:Bloco + LT, data=dados)
summary(mod.conj, split=list(LT=list(TdL1=1:5,TdL2=6:10, TdL3=11:15)))
#comparações múltiplas
require(agricolae)
#dentro de local 1
(tukey.l1<- HSD.test(dados$altura[dados$Local=="L1"], 
                    dados$Trat[dados$Local=="L1"], 
                    45, 3.0))
#dentro de local 2
(tukey.l2<- HSD.test(dados$altura[dados$Local=="L2"], 
                     dados$Trat[dados$Local=="L2"],
                     45, 3.0))
#dentro de local 3
(tukey.l3<- HSD.test(dados$altura[dados$Local=="L3"], 
                     dados$Trat[dados$Local=="L3"],
                     45, 3.0))

par(mfrow=c(1,3))
bar.group(tukey.l1$groups, ylim=c(0,30), 
          main="Local 1", xlab="Tratamentos", ylab="Altura (m)")
bar.group(tukey.l2$groups, ylim=c(0,30), 
          main="Local 2", xlab="Tratamentos", ylab="Altura (m)")
bar.group(tukey.l3$groups, ylim=c(0,30), 
          main="Local 3", xlab="Tratamentos", ylab="Altura (m)")