[R-es] Comparación de medias por LSD usando glht
Argel Gastélum Arellánez
argel.gastelum en gmail.com
Mie Jul 27 18:19:11 CEST 2011
El 27/07/11 08:11, Luciano Selzer escribió:
> Hola Argerl,
> Buscando un poco en internet encontré que solo tenes que cambiar el
> método para dar los p:
> cld(summary(Pares, method = univariate())
>
> Espero que sirva
>
> Luciano
Hola Luciano, muchas gracias por tu respuesta. Creí que había que
modificar la opción "Tukey" en la función glht para obtener la
comparación de medias por LSD. Hice la prueba con lo que me recomendaste
"method = univariate()" en la función cld, y obtuve valores de p
ligeramente diferentes, por ejemplo, con univariate resulta algo como esto:
Linear Hypotheses:
Estimate Std. Error z value Pr(>|z|)
4 - 2 == 0 0.001700 0.031188 0.055 1.0000
6 - 2 == 0 0.077867 0.031188 2.497 0.3762
8 - 2 == 0 0.099933 0.031188 3.204 0.0686 .
10 - 2 == 0 0.114618 0.035123 3.263 0.0579 .
12 - 2 == 0 0.120933 0.031188 3.878 <0.01 **
24 - 2 == 0 0.125933 0.031188 4.038 <0.01 **
36 - 2 == 0 0.178933 0.031188 5.737 <0.01 ***
etc...
y como lo hacía antes, con tukey da lo siguiente:
Linear Hypotheses:
Estimate Std. Error z value Pr(>|z|)
4 - 2 == 0 0.001700 0.031188 0.055 1.0000
6 - 2 == 0 0.077867 0.031188 2.497 0.3758
8 - 2 == 0 0.099933 0.031188 3.204 0.0685 .
10 - 2 == 0 0.114618 0.035123 3.263 0.0581 .
12 - 2 == 0 0.120933 0.031188 3.878 <0.01 **
24 - 2 == 0 0.125933 0.031188 4.038 <0.01 **
36 - 2 == 0 0.178933 0.031188 5.737 <0.01 ***
etc...
Por lo que al obtener el resumen compacto con letras da el mismo
resultado. ¿Es correcto? ¿Existe una opción para que en este caso sea
menos estricta para hacer la comparación de medias?
Saludos y de nuevo muchas gracias por tu ayuda.
--
Argel.
Más información sobre la lista de distribución R-help-es