[R-es] comparaciones planeadas

[Olivier Nuñez]

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