[R-br] ENC: Látice parcialmente balanceado em vários locais

[Euriann Yamamoto]

Walmes, eu consegui fazer pelo comando lapply a anova individual dos 5 locais. Só que agora não estou conseguindo colocasr os termos da fonte de variação em ordem:

# ajuste do modelo de efeitos fixos para bloco dentro de repeticao


#Anova Individual
mi <- lapply(split(da, f=da$local), aov,
             formula=producao ~ repeticao/bloco + cultivar)
lapply(mi, summary)#nao ajustado pra bloco

#Arrumar ordem dos termos e tornar cultivar ajustado pra bloco
mi <- lapply(split(da, f=da$local), aov,
             terms(formula=producao ~ repeticao/bloco + cultivar,
             keep.order = TRUE))

Ocorre um erro:
Error in terms.default(formula, "Error", data = data) :
  no terms component nor attribute


Creio que seja na fórmula. Qual seria a ordenação correta para poder utilizar esse comando terms?

Att.,

Euriann L. M. Yamamoto
Doutoranda em Agronomia
Universidade Federal da Grande Dourados





________________________________
De: Euriann Yamamoto
Enviado: quinta-feira, 19 de outubro de 2017 17:05
Para: Walmes Zeviani
Assunto: RE: [R-br] Látice parcialmente balanceado em vários locais

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?

Script:
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)


Euriann L. M. Yamamoto
Doutoranda em Agronomia
Universidade Federal da Grande Dourados




________________________________
De: Walmes Zeviani <walmeszeviani em gmail.com>
Enviado: terça-feira, 17 de outubro de 2017 19:53
Para: Euriann Yamamoto
Assunto: Re: [R-br] Látice parcialmente balanceado em vários locais

Não sei te responder com segurança. Testa e vê no que dá, rs. É bem provável que precise de adaptações, ainda que mínimas.

Walmes.
​
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20171019/018f7919/attachment.html>