¿Cómo obtener el parámetro de dispersión de un modelo mixto binomial?

Question

Estoy modelando datos con un modelo mixto generalizado con distribución de error binomial y estoy preocupado por la sobre dispersión. Sé que el parámetro de dispersión se puede medir como devianza/df, pero para modelos mixtos (= con efecto aleatorio), el número de grados de libertad no se puede extraer (estoy utilizando la función lmer de R). ¿Hay alguna manera de averiguar si mis datos están sobre dispersos?

Answer 1

Puede usar esta función atribuida a D. Bates para obtener el parámetro de escala:

dispersion_glmer<- function(modelglmer)
{   
# calculando la escala estimada (modelo binomial)
#siguiendo a D. Bates:
#Esa cantidad es la raíz cuadrada de la suma de los residuos penalizados al cuadrado dividido por n, el número de observaciones, evaluado como:

    n <- length(modelglmer@resid)

    return(  sqrt( sum(c(modelglmer@resid, modelglmer@u) ^2) / n ) )
}

Este es un enlace a más información sobre el parámetro de escala.