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@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.