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
_______________________________________________
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es