Множественная логистическая регрессия

Промежуточная регрессия в R

Richie Cotton

Data Evangelist at DataCamp

Набор данных об оттоке клиентов банка

has_churned	time_since_first_purchase	time_since_last_purchase
0	0.3993247	-0.5158691
1	-0.4297957	0.6780654
0	3.7383122	0.4082544
0	0.6032289	-0.6990435
...	...	...
response	length of relationship	recency of activity

¹ https://www.rdocumentation.org/packages/bayesQR/topics/Churn

glm()

glm(response ~ explanatory, data = dataset, family = binomial)

glm(response ~ explanatory1 + explanatory2, data = dataset, family = binomial)

glm(response ~ explanatory1 * explanatory2, data = dataset, family = binomial)

Процесс прогнозирования

explanatory_data <- expand_grid(
  explanatory1 = some_values,
  explanatory2 = some_values
)
prediction_data <- explanatory_data %>% 
  mutate(
    has_churned = predict(mdl, explanatory_data, type = "response")
  )

Четыре возможных исхода

	фактически ложь	фактически истина
предсказано: ложь	верно	ложноотрицательный
предсказано: истина	ложноположительный	верно

¹ https://campus.datacamp.com/courses/introduction-to-regression-in-r/simple-logistic-regression?ex=10

Матрица ошибок

actual_response <- dataset$response
predicted_response <- round(fitted(mdl))

outcomes <- table(predicted_response, actual_response)

confusion <- conf_mat(outcomes)

autoplot(confusion)

summary(confusion, event_level = "second")

Визуализация

Для категориальных переменных используйте разбивку на панели.
Для двух числовых объясняющих переменных кодируйте отклик цветом.
Значениям отклика ниже 0.5 назначьте один цвет, выше 0.5 — другой.

scale_color_gradient2(midpoint = 0.5)

Давайте потренируемся!

Промежуточная регрессия в R