[R-es] ANOVA de un factor con datos no balanceados y prueba de tukey
Argel Gast�lum Arell�nez
argel.gastelum en gmail.com
Mar Jun 7 19:16:08 CEST 2011
Hola compa�eros de la lista.
Tengo un experimento con un factor y una variable de respuesta,
pero con datos no balanceados. He aplicado un anova con la funci�n
"aov", y tambi�n con la funci�n "Anova" del paquete "car" sobre un
modelo lineal usando "lm". Los resultados me dan los mismos valores para
la suma de cuadrados y el valor de p en las tres anovas tipo I, II y
III, por lo que asumo que la forma de calcular la suma de cuadrados s�lo
ser� importante cuando se trata de dos factores o m�s. Corr�janme si me
equivoco por favor.
Como la diferencia es significativa, para hacer la comparaci�n de
medias, tengo entendido que la prueba de Tukey necesita tener datos
balanceados, que no es el caso de mi experimento. Entonces, he visto que
se puede hacer la comparaci�n por Tukey-Kramer o la prueba de Dunnett,
sin embargo, viendo la ayuda de R encuentro esto:
- TukeyHSD: "Technically the intervals constructed in this way
would only apply to balanced designs where there are the same number of
observations made at each level of the factor. This function
incorporates an adjustment for sample size that produces sensible
intervals for mildly unbalanced designs."
�Esto significa que puedo aplicar esta prueba en R tambi�n para mis
datos no balanceados?
- HSD.test (paquete agricolae): Aqu� s�lo dice "It is necessary
first makes a analysis of variance".
�Entonces este no hace la correcci�n que menciona TukeyHSD?
- LSD.test (paquete agricolae): Aqu� dice "For equal or different
repetition". Asumo que este s� lo podr�a aplicar a mis datos no balanceados.
- DTK.test (paquete DTK): Aqu� menciona "Conducts a pairwise
multiple comparison test (using the C procedure) for mean differences
with unequal sample sizes and no assumption of equal population variances"
�Qu� me recomiendan aplicar para comparar las medias de mis datos?
De antemano muchas gracias por la ayuda.
Saludos.
--
Argel.
M�s informaci�n sobre la lista de distribuci�n R-help-es