Sou o autor do pacote ez para R e estou trabalhando em uma atualização para incluir o cálculo automático das taxas de verossimilhança (LRs) na saída de ANOVAs. A idéia é fornecer uma LR para cada efeito que seja análogo ao teste desse efeito que a ANOVA alcança. Por exemplo, a LR para efeito principal representa a comparação de um modelo nulo com um modelo que inclui o efeito principal, a LR para uma interação representa a comparação de um modelo que inclui os efeitos principais dos componentes versus um modelo que inclui os efeitos principais e a interação deles, etc.
Agora, meu entendimento sobre computação LR vem de Glover & Dixon ( PDF ), que abrange cálculos básicos e correções de complexidade, e o apêndice de Bortolussi & Dixon ( apêndice PDF ), que abrange cálculos envolvendo variáveis de medidas repetidas. Para testar meu entendimento, desenvolvi esta planilha , que pega os dfs & SSs de um exemplo de ANOVA (gerado a partir de um design 2 * 2 * 3 * 4 usando dados falsos) e percorre o cálculo da LR para cada efeito.
Eu realmente apreciaria se alguém com um pouco mais de confiança com esse cálculo pudesse dar uma olhada e garantir que eu fiz tudo corretamente. Para aqueles que preferem código abstrato, aqui está o código R que implementa a atualização para ezANOVA () (consulte as linhas 15-95).
fonte