Buenos días:
estoy trabajando con modelos mixtos lineales generalizados mediante la función glmmadmb (paquete glmmADMB) y cuando quiero comparar los modelos mediante el QAIC, resultado entregado por qAIC (paquete bbmle) no se si el valor de dispersión con que tengo que completar es el “Negative binomial dispersion parameter ”, valor que devuelve el summary de glmmadmb o si tengo que poner el de c-hat.
Ejemplo:

install.packages("glmmADMB", repos="http://R-Forge.R-project.org")
library(glmmADMB)
install.packages(“bbmle”)
library(bbmle)

epil2$subject <- factor(epil2$subject)
  (fm <- glmmadmb(y~Base*trt+Age+Visit+(Visit|subject), data=epil2, family="nbinom"))
summary(fm)

#esta es parte de la salida.

Number of observations: total=236, subject=59 
Random effect variance(s):
Group=subject
             Variance    StdDev
(Intercept) 2.172e-01 4.660e-01
Visit       8.504e-09 9.222e-05
Negative binomial dispersion parameter: 7.4702 (std. err.: 1.7605)

Log-likelihood: -624.551 
________________________________________

qAIC(fm,dispersion=XXXXX?)

En caso de ser c-hat ¿cómo calculo ese valor?




Darío Ezequiel Manzoli
Laboratorio de Ecología de Enfermedades
Facultad de Ciencias Veterinarias
Universidad Nacional del Litoral
Kreder 2805 - Esperanza (S3080HOF) Santa Fe
03496-426304(int 111-353)

web: www.fcv.unl.edu.ar/archivos/servicios/laboratorios/lecen/lecen1.htm






	[[alternative HTML version deleted]]

