Hosmer-Lemeshow goodness-of-fit test for logistic regression

Hosmer-Lemeshow goodness-of-fit test for logistic regression

Usage

hosmer_lemeshow_test(y, y_prob, n_groups = 10L)

Arguments

y

Binary response vector.

y_prob

Predicted probabilities.

n_groups

Number of decile groups (default 10).

Value

A morie_specification_test.

Examples

set.seed(10)
n <- 200
X <- cbind(1, rnorm(n))
y_prob <- plogis(drop(X %*% c(0, 1)))
y <- rbinom(n, 1, y_prob)
r <- hosmer_lemeshow_test(y, y_prob, n_groups = 5L)
r$statistic
#> [1] 2.163695