[R-sig-ME] Hosmer-Lemeshow test in GLMM (binomial family) using lme4 package

Ben Bolker bbo|ker @end|ng |rom gm@||@com
Mon Jul 10 16:53:58 CEST 2023


    I think it should be appropriate.

   However, Frank Harrell (a leading biostatistician)

https://twitter.com/f2harrell/status/1228423023834718208

"The Hosmer-Lemeshow test has been obsolete for more than a decade.  Not 
recommend.  Low power, hard to interpret, very arbitrary to how deciles 
are computed."

  You can use 
https://search.r-project.org/CRAN/refmans/DescTools/html/HosmerLemeshowTest.html 
(which returns both HL and the recommned le Cessie-van 
Houwelingen-Copas-Hosmer statistics ...)

   cheers
    Ben Bolker


On 2023-07-10 10:27 a.m., Camila Deutsch wrote:
> Hi!
> I ran a GLMM (binomial family) using the package lme4:
> library(lme4)
> #GLMM
> m1 <- glmer(Canto ~ T + H + Patm + V + Pp + CU + (1|Evento), data = Datos,
> family = binomial)
> 
> I was wondering if it would be appropriate to evaluate the goodness of fit
> of the model by applying the Hosmer-Lemeshow test.
> 
> Thank you very much.
> Lic. Camila Deutsch
> Becaria Doctoral-CONICET
> Grupo de Estudios sobre Biodiversidad en Agroecosistemas
> Facultad de Ciencias Exactas y Naturales. IEGEBA
> Universidad de Buenos Aires - CONICET
> 
> 	[[alternative HTML version deleted]]
> 
> _______________________________________________
> R-sig-mixed-models using r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models



More information about the R-sig-mixed-models mailing list