
Comparto con Olivier que nlme es la herramienta adecuada.
Si tienes tiempo, consulta el excelente libro de Gelman y Hill "Data Analysis Using Regression and Multilevel/Hierarchical Models" - la segunda parte explica distintos modelos anidados/jerárquicos/multinivel utilizando R.

José


-----Original Message-----
From: r-help-es-bounces@r-project.org on behalf of Olivier Nuñez
Sent: Tue 24/05/2011 11:35
To: Miguel Lázaro
Cc: R-es CRAN
Subject: Re: [R-es] comparaciones planeadas
 
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@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@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es

_______________________________________________
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es


	[[alternative HTML version deleted]]

