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