Hola a todos/as, Quería preguntaros por una cosa bastante sencilla, entiendo, pero que no logro solucionar. Se trata de un análisis planeado. Por ejemplo, si tenemos una variable intrasujetos y una intersujetos, muchas veces interesa saber la significación por niveles puesto que la significación podría restringirse a ciertos nievles de una determinada variable. Algunas hipótesis de los trabajos pueden no referirse en exclusiva a los efectos principales, sino a esto y ahí surge mi problema. Si tuviéramos una tabla como la siguiente por ejemplo Participante factor1 factor2 puntuacion Sujeto1 a1 b1 2 Sujeto1 a1 b2 3 Sujeto1 a1 b3 2 Sujeto1 a1 b4 5 Sujeto2 a1 b1 5 Sujeto2 a1 b2 6 Sujeto2 a1 b3 3 Sujeto2 a1 b4 2 Sujeto2 a2 b1 3 Sujeto3 a2 b2 4 Sujeto3 a2 b3 3 Sujeto3 a2 b4 4 Sujeto4 a2 b1 4 Sujeto4 a2 b2 5 Sujeto4 a2 b3 4 Sujeto4 a2 b4 2 e introdujera>resultado=aov(puntuación~factor1*factor2+Error(subject/factor2), data =a) >summary(resultado)Tendría los resultados de las anovas inter e intrasujetos y la interacción, pero no sabría si existe significación de un factor tan solo en uno de los niveles del otro, por ejemplo si solo existiría significación para el nivel 1 y 2 del factor2 pero no de los niveles 3 y 4. No logro encontrar en donde busco los comandos para esta comparación planeada. ¿Alguien podría ayudarme? Gracias Un saludo cordial Miguel [[alternative HTML version deleted]]
Prueba require(nlme) a$interact=interaction(a$factor1,a$factor2) resultado=lme(puntuación~ interact, random=~1|subject/factor2, data =a) summary(glht(resultado, linfct=mcp(interact="Tukey"))) Un saludo. Olivier -- ____________________________________ Olivier G. Nuñez Email: onunez en iberstat.es Tel : +34 663 03 69 09 Web: http://www.iberstat.es ____________________________________ El 23/05/2011, a las 23:12, Miguel Lázaro escribió:> Hola a todos/as, > Quería preguntaros por una cosa bastante sencilla, entiendo, pero > que no logro solucionar. Se trata de un análisis planeado. Por > ejemplo, si tenemos una variable intrasujetos y una intersujetos, > muchas veces interesa saber la significación por niveles puesto que > la significación podría restringirse a ciertos nievles de una > determinada variable. Algunas hipótesis de los trabajos pueden no > referirse en exclusiva a los efectos principales, sino a esto y ahí > surge mi problema. Si tuviéramos una tabla como la siguiente por > ejemplo > > Participante factor1 factor2 puntuacion > Sujeto1 a1 b1 2 > Sujeto1 a1 b2 3 > Sujeto1 a1 b3 2 > Sujeto1 a1 b4 5 > Sujeto2 a1 b1 5 > Sujeto2 a1 b2 6 > Sujeto2 a1 b3 3 > Sujeto2 a1 b4 2 > Sujeto2 a2 b1 3 > Sujeto3 a2 b2 4 > Sujeto3 a2 b3 3 > Sujeto3 a2 b4 4 > Sujeto4 a2 b1 4 > Sujeto4 a2 b2 5 > Sujeto4 a2 b3 4 > Sujeto4 a2 b4 2 > > > e introdujera >> resultado=aov(puntuación~factor1*factor2+Error(subject/factor2), >> data =a) >> summary(resultado) > > Tendría los resultados de las anovas inter e intrasujetos y la > interacción, pero no sabría si existe significación de un factor > tan solo en uno de los niveles del otro, por ejemplo si solo > existiría significación para el nivel 1 y 2 del factor2 pero no de > los niveles 3 y 4. > No logro encontrar en donde busco los comandos para esta > comparación planeada. ¿Alguien podría ayudarme? > Gracias > Un saludo cordial > Miguel > [[alternative HTML version deleted]] > > _______________________________________________ > R-help-es mailing list > R-help-es en r-project.org > https://stat.ethz.ch/mailman/listinfo/r-help-es
La librería es multcomp:
install.packages("multcomp")
require(multcomp)
Un saludo. Olivier
--
____________________________________
Olivier G. Nuñez
Email: onunez@iberstat.es
Tel : +34 663 03 69 09
Web: http://www.iberstat.es
____________________________________
El 25/05/2011, a las 1:42, Miguel Lázaro escribió:
> Muchas gracias a ambos,
> he llegado hasta casi el final pero me indica que la función glht
> no está disponible. ¿dónde la puedo descargar, de qué paquete es?
> es un pquete glht?
> intall.packages(c("??")
>
> Muchas gracias también por la bibliografía. Yo empiezo en esto muy
> paso a paso como se ve claramente... y utilizo el de Harald Baayen
> Analyzing linguistic data: a practical introduction to statistics.
> Me parece bastante bueno y tiene un draft gratuito por si queréis
> mirarlo en internet.
> http://www.ualberta.ca/~baayen/publications/baayenCUPstats.pdf
>
> un saludo
> Miguel
[[alternative HTML version deleted]]