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

[Olivier Nuñez]

Eduardo,

si no te conformas con una estructura de correlación sencilla  
("simetría compuesta"),
puedes sofisticar el modelo, de manera que la correlación vaya  
decreciendo conforme se alejan la fechas.
En este caso, es mejor utilizar la función glmmPQL() del paquete MASS.
Utilizando el ejemplo de mi anterior mensaje con una correlación de  
tipo AR(1), tendríamos la sintaxis:

fit=glmmPQL(Conteo ~ especie + zona, random= ~ 1 | nido, family =  
poisson)
update(fit , correlation=corAR1(rho, form = ~ 1 | nido))

donde "rho" es un valor inicial para el parámetro de correlación (por  
ejemplo, rho=0.5).

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