[R-es] Espcificar analisis de medidas repetidas en lme()

[Marcelino de la Cruz Rot]

Hola Jaume:

En teoría en cada fila (cada observación)  de la tabla de datos tienes 
los valores adecuados de "id" y "visita", por lo que R "entiende" de 
sobra la estructura individual-temporal de las observaciones.

Yo pondría:

lme(x ~ sitio + visita, random = ~ 1|id, data = datos)


o

lme(x ~ sitio + visita, random = ~ 1 + sitio + visita |id, data = datos)
  

o

lme(x ~ sitio * visita, random = ~ 1|id, data = datos)


etc, dependiendo de los efectos fijos y aleatorios que quieras testar.


Un saludo,

Marcelino


El 11/01/2022 a las 14:23, Jaume Tormo escribió:
> Hola,
>
> Estoy tratando de especificar un modelo de medidas repetidas que es como
> sigue:
> - 4 sitios de estudio. (A-D)
> - Dentro de cada sitio hay 5 puntos donde se mide la variable X
> - Repetimos las medidas en el tiempo. Los lugares y puntos son los mismos
> en cada visita.
>
> Tenemos las siguientes variables:
> id - identificador del punto de medida
> sitio - Los cuatro lugares donde se mide sitio1, sitio2, ... sitio4
> visita - La visitas que se hacen para medir. visita1, visita2,... visitaN
> x - variable medida
>
> Lo voy a analizar con lme pero no estoy seguro de cómo especificar esto en
> la fórmula. Yo creo que es así:
>
> lme(x~sitio, random = ~ 1|id/visita, data = datos)
>
> Porque entiendo que hay que decir id/visita para que R entienda que a lo
> largo de las visitas los id son el mismo punto de medida.
> Pero no estoy seguro de si debe ser al revés ~1|visita/id
>
> ¿Que opinais?¿Voy por buen camino?
>
> Muchas gracias.
>

-- 
Marcelino de la Cruz Rot
Depto. de Biología y Geología
Física y Química Inorgánica
Universidad Rey Juan Carlos
Móstoles España