Soy el autor de la paquete ez para R, y estoy trabajando en una actualización para incluir el cálculo automático de los cocientes de probabilidad (LRs) en la salida de los ANOVAs. La idea es proporcionar un LR para cada efecto que sea análogo a la prueba de ese efecto que el ANOVA logra. Por ejemplo, el LR para el efecto principal representa la comparación de un modelo nulo con un modelo que incluye el efecto principal, el LR para una interacción representa la comparación de un modelo que incluye ambos efectos principales componentes contra un modelo que incluye ambos efectos principales y su interacción, etc.
Ahora, mi comprensión del cálculo de LR viene de Glover y Dixon ( PDF ), que abarca los cálculos básicos así como las correcciones de la complejidad, y el apéndice de Bortolussi y Dixon ( apéndice PDF ), que abarca los cálculos que implican variables de medidas repetidas. Para probar mi comprensión, desarrollé esta hoja de cálculo que toma los dfs & SSs de un ejemplo de ANOVA (generado a partir de un diseño de 2*2*3*4 utilizando datos falsos) y pasa por el cálculo del LR para cada efecto.
Apreciaría mucho que alguien con un poco más de confianza en este tipo de cálculo echara un vistazo y se asegurara de que hice todo correctamente. Para aquellos que prefieren el código abstracto, aquí está el código R implementando la actualización de ezANOVA() (ver esp. líneas 15-95).
