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]]