Hola, gurús estadísticos y asistentes de programación R,
Estoy interesado en modelar capturas de animales en función de las condiciones ambientales y el día del año. Como parte de otro estudio, tengo recuentos de capturas en ~ 160 días durante tres años. En cada uno de estos días tengo temperatura, lluvia, velocidad del viento, humedad relativa, etc. Debido a que los datos se recopilaron repetidamente de las mismas 5 parcelas, utilizo la parcela como un efecto aleatorio.
Tengo entendido que nlme puede explicar fácilmente la autocorrelación temporal en los residuos, pero no maneja funciones de enlace no gaussianas como lme4 (¿cuál no puede manejar la autocorrelación?). Actualmente, creo que podría funcionar usar el paquete nlme en R en log (count). Entonces, mi solución en este momento sería ejecutar algo como:
m1 <- lme(lcount ~ AirT + I(AirT^2) + RainAmt24 + I(RainAmt24^2) + RHpct + windspeed +
sin(2*pi/360*DOY) + cos(2*pi/360*DOY), random = ~1|plot, correlation =
corARMA(p = 1, q = 1, form = ~DOY|plot), data = Data)
donde DOY = Día del año. Puede haber más interacciones en el modelo final, pero esta es mi idea general. También podría intentar modelar aún más la estructura de varianza con algo como
weights = v1Pow
No estoy seguro de si hay una mejor manera de hacer con una regresión de modelo mixto de Poisson o algo así. Acabo de encontrar una discusión matemática en el Capítulo 4 de "Modelos de regresión para el análisis de series temporales" de Kedem y Fokianos. Estaba un poco más allá de mí en este momento, especialmente en la aplicación (codificándolo en R). También vi una solución MCMC en Zuur et al. Libro de modelos de efectos mixtos (Capítulo 23) en el lenguaje BUGS (usando winBUGS o JAG). ¿Esa es mi mejor opción? ¿Hay un paquete MCMC fácil en R que pueda manejar esto? No estoy realmente familiarizado con las técnicas GAMM o GEE, pero estaría dispuesto a explorar estas posibilidades si las personas pensaran que proporcionarían una mejor visión.Mi objetivo principal es crear un modelo para predecir las capturas de animales dadas las condiciones ambientales. En segundo lugar, me gustaría explicar a qué responden los animales en términos de su actividad.
Cualquier idea sobre la mejor manera de proceder (filosóficamente), cómo codificar esto en R o en ERRORES sería apreciada. Soy bastante nuevo en R y BUGS (winBUGS) pero estoy aprendiendo. Esta es también la primera vez que trato de abordar la autocorrelación temporal.
Gracias dan