De passing van logistische regressie kwantificeren

Introductie tot regressie in R

Richie Cotton

Data Evangelist at DataCamp

Vier mogelijke uitkomsten

	feitelijk onwaar	feitelijk waar
voorspeld onwaar	correct	false negative
voorspeld waar	false positive	correct

Confusion matrix: aantallen uitkomsten

mdl_recency <- glm(has_churned ~ time_since_last_purchase, data = churn, family = "binomial")

actual_response <- churn$has_churned

predicted_response <- round(fitted(mdl_recency))

outcomes <- table(predicted_response, actual_response)

                  actual_response
predicted_response   0   1
                 0 141 111
                 1  59  89

Confusionmatrix visualiseren: mozaïekdiagram

library(ggplot2)
library(yardstick)

confusion <- conf_mat(outcomes)

                  actual_response
predicted_response   0   1
                 0 141 111
                 1  59  89

autoplot(confusion)

Een mozaïekdiagram van de uitkomsten van het churn-versus-recency-model. Er zijn 200 observaties voor zowel echte churns als echte niet-churns, dus elke kolom heeft dezelfde breedte.

Prestatiematen

summary(confusion, event_level = "second")

# A tibble: 13 x 3
   .metric              .estimator .estimate
   <chr>                <chr>          <dbl>
 1 accuracy             binary         0.575
 2 kap                  binary         0.150
 3 sens                 binary         0.445
 4 spec                 binary         0.705
 5 ppv                  binary         0.601
 6 npv                  binary         0.560
 7 mcc                  binary         0.155
 8 j_index              binary         0.150
 9 bal_accuracy         binary         0.575
10 detection_prevalence binary         0.37 
11 precision            binary         0.601
12 recall               binary         0.445
13 f_meas               binary         0.511

Accuracy

summary(confusion) %>% 
  slice(1)

# A tibble: 3 x 3
  .metric  .estimator .estimate
  <chr>    <chr>          <dbl>
1 accuracy binary         0.575

Accuracy is het aandeel correcte voorspellingen.

$$ accuracy = \frac{TN + TP}{TN + FN + FP + TP} $$

confusion

                  actual_response
predicted_response   0   1
                 0 141 111
                 1  59  89

(141 + 89) / (141 + 111 + 59 + 89)

0.575

Sensitiviteit

summary(confusion) %>% 
  slice(3)

# A tibble: 1 x 3
  .metric .estimator .estimate
  <chr>   <chr>          <dbl>
1 sens    binary         0.445

Sensitiviteit is het aandeel true positives.

$$ sensitivity = \frac{TP}{FN + TP} $$

confusion

                  actual_response
predicted_response   0   1
                 0 141 111
                 1  59  89

89 / (111 + 89)

0.445

Specificiteit

summary(confusion) %>% 
  slice(4)

# A tibble: 1 x 3
  .metric .estimator .estimate
  <chr>   <chr>          <dbl>
1 spec    binary         0.705

Specificiteit is het aandeel true negatives.

$$ specificity = \frac{TN}{TN + FP} $$

confusion

                  actual_response
predicted_response   0   1
                 0 141 111
                 1  59  89

141 / (141 + 59)

0.705

Laten we oefenen!

Introductie tot regressie in R