R help es - Jul 2011 - Comparación de medias por LSD usando glht

Hola R-listeros.

     Tengo la siguiente duda, estoy haciendo una comparación de medias a 
un modelo obtenido de la siguiente forma:

MODELO <- lme(RESPUESTA ~ FACTOR, random = ~1|REPETICION/FACTOR, data = 
DATOS)

Pares <- glht(MODELO, linfct = mcp(FACTOR = "Tukey"))

cld(Pares)

     Esto funciona muy bien, pero quisiera comparar el resultado 
aplicando una comparación de medias por diferencia mínima significativa. 
¿Cómo puedo hacer esto usando glht?

     De antemano muchas gracias por la ayuda.

     Saludos.

--
     Argel.

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


El 27 de julio de 2011 04:49, Argel Gastélum Arellánez <
argel.gastelum@gmail.com> escribió:
>    Hola R-listeros.
>
>    Tengo la siguiente duda, estoy haciendo una comparación de medias a un
> modelo obtenido de la siguiente forma:
>
> MODELO <- lme(RESPUESTA ~ FACTOR, random = ~1|REPETICION/FACTOR, data
> DATOS)
>
> Pares <- glht(MODELO, linfct = mcp(FACTOR = "Tukey"))
>
> cld(Pares)
>
>    Esto funciona muy bien, pero quisiera comparar el resultado aplicando
> una comparación de medias por diferencia mínima significativa. ¿Cómo puedo
> hacer esto usando glht?
>
>    De antemano muchas gracias por la ayuda.
>
>    Saludos.
>
> --
>    Argel.
>
> ______________________________**_________________
> R-help-es mailing list
> R-help-es@r-project.org
>
https://stat.ethz.ch/mailman/**listinfo/r-help-es<https://stat.ethz.ch/mailman/listinfo/r-help-es>
>
	[[alternative HTML version deleted]]

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.