Hola, ¿qué tal? Yo uso lmer habitualmente (y no lme) y me pierdo un poco con la sintaxis. Pero vamos, casi seguro que lme no espera que "visita" figure en el modelo.
Tampoco está claro si quieres que los sitios sean efectos fijos o aleatorios. Casi me inclino más por lo segundo, pero no sé cuál es el contexto exacto. De todos modos, tienes un ejemplo similar al tuyo tratado con lmer aqui <http://www.alexanderdemos.org/Mixed5.html> (sección *3 Level Random intercepts model*). Un saludo, Carlos J. Gil Bellosta http://www.datanalytics.com . El mar, 11 ene 2022 a las 14:24, Jaume Tormo (<jautor...@gmail.com>) 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. > > -- > Jaume Tormo. > https://flipboard.com/@jaumetormo/hallazgos-interesantes-bj8opmboy > https://acercad.wordpress.com/ > > [[alternative HTML version deleted]] > > _______________________________________________ > R-help-es mailing list > R-help-es@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-help-es > [[alternative HTML version deleted]] _______________________________________________ R-help-es mailing list R-help-es@r-project.org https://stat.ethz.ch/mailman/listinfo/r-help-es