[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