[R-br] Soma de quadrados tipo III

[Thiago de Paula Protásio]

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
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130409/f245b989/attachment.html>