[R-es] Help MLG(M) con datos longitudinales

[Olivier Nuñez]

Eduardo,

en tu caso, no creo necesario que especifiques un modelo de  
autocorrelación.
Si utilizas un modelo lineal generalizado mixto (ver por ejemplo la  
función glmer() del paquete "lme4"),
por defecto tus medidas repetidas tendrán correlación en el tiempo.

La sintaxis de la función glmer() no presenta en principio dificultad.
Supongamos que cuentas los huevos en distintos nidos  durante 10 años  
y que quieres explicar la variación de este conteo en función de la  
especie
y de la zona geográfica, se podría considerar el siguiente modelo:

glmer(Conteo ~ especie + zona + (1|nido) , datos, family = poisson)

que asume que los nº de huevos en un mismo nido están correlados de  
un año a otro.
Además se contempla la posibilidad de que en promedio este número  
varíe con la especie y la zona geográfica (¡sin interacción entre  
ambos factores!).

Un saludo. Olivier

--  
____________________________________

Olivier G. Nuñez
Email: onunez en iberstat.es
Tel : +34 663 03 69 09
Web: http://www.iberstat.es

____________________________________




El 13/09/2010, a las 15:26, Eduardo Trumper escribió:

> Estimados,
> Necesito estudiar los efectos de dos factores  sobre una variable  
> de conteo
> con un rango estrecho de variación (es decir, probablemente sin  
> distribución
> Normal). El experimento tiene más de 10 (fechas de muestreo)  
> observaciones
> repetidas en el tiempo.
>
> En primer lugar, considero que debería aplicar un modelo lineal  
> generalizado
> (¿Mixto?) usando Poisson como distribución de referencia.
>
> Segundo, creo que se trata de un problema de datos longitudinales.  
> Tengo
> entendido que se debería aplicar un modelo de correlación serial  
> como el
> autoregresivo de orden 1. Ya sea para éste modelo o para otro que  
> me puedan
> aconsejar, agradecería mucho me ayuden con la sintaxis.
>
>
>
> Eduardo Trumper
>
> 	[[alternative HTML version deleted]]
>
> _______________________________________________
> R-help-es mailing list
> R-help-es en r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es