Pessoal; Novamente recorro a ajuda de vcs... Preciso analisar um experimento em esquema fatorial duplo em blocos casualizados com diferentes números de repetições. Mas estou tendo dificuldades em obter as médias ajustadas para cada fator, pois a interação foi não significativa. Por favor, vejam a rotina abaixo: anava1=aov(QEA~especie*adesivo+Bloco) Anova(anava1,type="III") # ajuste do modelo aditivo com teste F marginal m1 <- lm(QEA~especie*adesivo+Bloco) drop1(m1, test="F") # estimativas dos efeitos (solução) sob a restrição do R summary(m1) # obtenção das médias ajustadas dos níveis de tratamento require(contrast) lapply(levels(especie), function(i){ contrast(m1, type="average", list(especie=i, adesivo=levels(adesivo))) } ) Erro em generateData(fit = list(coefficients = c(9.49304920888927, -3.79433616881114, : not enough factors Meus sinceros agradecimentos. -- Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Departamento de Ciências Florestais Ciência e Tecnologia da Madeira (035) 9183 - 2246