Je suis l'auteur du package ez pour R, et je travaille sur une mise à jour pour inclure le calcul automatique des rapports de vraisemblance (LR) dans la sortie des ANOVA. L'idée est de fournir un LR pour chaque effet qui est analogue au test de cet effet que l'ANOVA réalise. Par exemple, le LR pour l'effet principal représente la comparaison d'un modèle nul avec un modèle qui inclut l'effet principal, le LR pour une interaction représente la comparaison d'un modèle qui comprend à la fois les effets principaux des composants et un modèle qui comprend à la fois les effets principaux et leur interaction, etc.
Maintenant, ma compréhension du calcul LR vient de Glover & Dixon ( PDF ), qui couvre les calculs de base ainsi que les corrections de complexité, et l'annexe de Bortolussi & Dixon ( annexe PDF ), qui couvre les calculs impliquant des variables à mesures répétées. Pour tester ma compréhension, j'ai développé cette feuille de calcul , qui prend les dfs et les SS à partir d'un exemple d'ANOVA (généré à partir d'une conception 2 * 2 * 3 * 4 utilisant de fausses données) et parcourt le calcul du LR pour chaque effet.
J'apprécierais vraiment que quelqu'un avec un peu plus de confiance avec un tel calcul puisse jeter un œil et s'assurer que j'ai tout fait correctement. Pour ceux qui préfèrent le code abstrait, voici le code R implémentant la mise à jour de ezANOVA () (voir en particulier les lignes 15-95).
la source