[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