Processing math: 100%

1 votos

Cómo modelar medidas repetidas

Me interesa saber cómo afectan las diferentes variables climáticas a la abundancia en abejas en dos hábitats diferentes (prados y edificios). Tenía 10 transectos en zonas de prados y 10 en zonas de edificios. Recorrí cada transecto 15 veces.

Variable de respuesta: Abundancia
Factores fijos: Temperatura del aire, Sol
Factores aleatorios: Transecto anidado en el hábitat, veces (1 a 15)

Por ejemplo, usando R, ajusté este modelo:

library (lme4)
model <- glmer(Abundance ~ Airtemp * Sunshine + (1|Habitat/transect) + (1|Times),
               family = poisson, data = bees)

¿Es ésta una buena manera de abordar mi pregunta de investigación?

1voto

Bruce ONeel Puntos 391

El concepto clave de las medidas repetidas es que hay agrupación. Por ello, los modelos de efectos mixtos se utilizan a menudo para las medidas repetidas, ya que manejan específicamente los datos agrupados. Estos modelos tienen en cuenta la no independencia de los datos dentro de cada conglomerado, es decir, las observaciones de un conglomerado serán más similares entre sí que las de otro conglomerado.

Para ajustar un modelo de efectos mixtos para medidas repetidas, es necesario especificar la variable de agrupación (es decir, lo que define el conglomerado) y ajustar interceptos aleatorios para ella. En su caso, tiene observaciones anidadas en transectos que, a su vez, están anidados dentro del hábitat. Sin embargo, sólo tiene dos hábitats, por lo que no es una buena idea especificar esto como un efecto aleatorio. En su lugar, puede especificarlo como un efecto fijo.

Si también se especifica Times como un intercepto aleatorio, está diciendo que las observaciones también están anidadas dentro de cada ocasión de medición. Esto puede ser apropiado o no. Si hay alguna razón por la que las observaciones en diferentes transectos en la misma ocasión de medición deberían ser más parecidas entre sí que las observaciones en los mismos transectos pero en un momento diferente, y no está interesado en el efecto fijo de la ocasión de medición, entonces puede incluir este término como un efecto aleatorio en el modelo.

Por lo tanto, el punto clave es que en su caso necesita ajustar un modelo con transect como un intercepto aleatorio. En R, un mejor modelo de partida es:

glmer(Abundance ~ Airtemp * Sunshine + Habitat + (1|transect) + (1|Times),
           family = poisson, data = bees)

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X