teste para significância da interação em modelos mistos lineares em nlme em R
eu usolme função nonlme Pacote R para testar se os níveis de fatoritems tem interação significativa com os níveis de fatorcondition. O fatorcondition tem dois níveis:Control eTreatmente o fatoritems tem 3 níveis:E1,...,E3. Eu uso o seguinte código:
f.lme = lme(response ~ 0 + factor(condition) * factor(items), random = ~1|subject)
Ondesubject é o efeito aleatório. Desta forma, quando eu corro:
summary(f.lme)$tTable
Eu receberei a seguinte saída:
factor(condition)Control
factor(condition)Treatment
factor(items)E2
factor(items)E3
factor(condition)Treatment:factor(items)E2
factor(condition)Treatment:factor(items)E3
junto comValue, Std.Error, DF, t-value, p-value colunas. Eu tenho duas perguntas:
Se eu quiser compararControl vs.Treatment, devo apenas usarestimable() função emgmodels e fazer um contraste de(-1,1,0,0,0,0)?
Eu estou interessado em saber se os níveis deitems, isto éE1, E2, E3 são diferentescondition, por isso estou interessado em saber se os termos de interação são significativos (apenas verificando op-value coluna??):
factor(condition)Treatment:factor(items)E2 factor(condition)Treatment:factor(items)E3
No entanto, como posso saber sefactor(condition)Treatment:factor(items)E1 é significativo ou não? Não é mostrado na saída de resumo e acho que tem algo a ver com o contraste usado em R ... Muito obrigado!