¿Hay algún problema con el resultado de la corrección de Huynh-Feldt cuando se ejecuta un ANOVA tanto para factores entre y dentro de los factores (diseño mixto) en R? Mi problema es que cuando ejecuto un ANOVA "dentro del sujeto (sólo)" usando las funciones "Anova" o "ezANOVA", la salida de la corrección de Huynh-Feldt es exactamente la misma que obtengo cuando uso SPSS. Sin embargo, cuando ejecuto un ANOVA en un diseño mixto, y considero tanto los efectos entre y dentro del sujeto usando las funciones 'Anova' o 'ezANOVA', mis correcciones de Greenhouse-Geisser coinciden con SPSS, pero las correcciones de Huynh-Feldt están ligeramente fuera comparadas con la salida de SPSS? ¿Alguna razón por la que esto podría ser? ¡Agradezco cualquier consejo/ayuda que alguien pueda dar!
Respuesta
¿Demasiados anuncios?Si hubiera un error en Anova también lo habría en ezANOVA porque Mike Lawrence construyó ez sobre el paquete Anova de John Fox. Y el paquete Anova, a su vez, utiliza el estándar de R anova
de mando. R está consiguiendo que esto sea correcto. Peter Dalgaard (pp.3-4) señalaron que SAS y SPSS obtienen el cálculo H-F equivocado y lo han hecho durante décadas. Es porque utilizan la fórmula original publicada y no la ecuación corregida (Lecoutre, 1991) que utiliza R.
(Mi especulación sobre por qué no lo han arreglado es similar a la razón por la que Excel sigue teniendo errores... tienen una mayor obligación con los principales clientes de ser compatibles con los errores que de ser correctos. Arreglar esto sería admitir que también se equivocaron. Aunque si se hubiera hecho en 1992 habría sido razonable. Ahora es un gran argumento a favor del código abierto. Otra posibilidad es que no les importe. Me pregunto si uno de ellos se limita a probar su resultado contra el otro).