powered by
hoslem.test(x, y, g = 10)
"htest"
sum((observed - expected)^2 / expected)
g - 2
g
set.seed(123) n <- 500 x <- rnorm(n) y <- rbinom(n, 1, plogis(0.1 + 0.5*x)) m <- glm(y ~ x, family=binomial) hoslem.test(m$y, fitted(m))
Run the code above in your browser using DataLab