Prevedere il punteggio didattico con l'età

Modellazione con i dati nel Tidyverse

Albert Y. Kim

Assistant Professor of Statistical and Data Sciences

Ripasso: retta di regressione

Previsione per un nuovo docente

Ripasso: tabella di regressione

library(ggplot2)
library(dplyr)
library(moderndive)

# Fit regression model using formula of form: y ~ x
model_score_1 <- lm(score ~ age, data = evals)

# Output regression table using wrapper function
get_regression_table(model_score_1)

# A tibble: 2 x 7
  term      estimate std_error statistic p_value lower_ci...
  <chr>        <dbl>     <dbl>     <dbl>   <dbl>    <dbl>...
1 intercept    4.46      0.127     35.2    0        4.21...
2 age         -0.006     0.003     -2.31   0.021   -0.011...

Valore previsto

Modelli di regressione predittivi in generale: $\hat{y} = \hat{f}(x) = \hat{\beta}_0 + \hat{\beta}_1 \cdot x$
Il nostro modello predittivo: $\hat{\text{score}} = 4.46 - 0.006 \cdot \text{age}$
La nostra previsione: $4.46 - 0.006 \cdot 40 = 4.22$

Errore di previsione

Errore di previsione

Residui come errori del modello

Residuo = $y - \hat{y}$
Corrisponde a $\epsilon$ in $y = f(\vec{x}) + \epsilon$
Per il nostro docente di esempio: $y - \hat{y} = 3.5 - 4.22 = -0.72$
Nella regressione lineare, in media sono 0.

Calcolare tutti i valori previsti

# Fit regression model using formula of form: y ~ x
model_score_1 <- lm(score ~ age, data = evals)

# Get information on each point
get_regression_points(model_score_1)

# A tibble: 463 x 5
      ID score   age score_hat residual
   <int> <dbl> <dbl>     <dbl>    <dbl>
 1     1   4.7    36      4.25    0.452
 2     2   4.1    36      4.25   -0.148
 3     3   3.9    36      4.25   -0.348
 4     4   4.8    36      4.25    0.552
 5     5   4.6    59      4.11    0.488

Retta di regressione «miglior adattamento»

Esercitiamoci!

Modellazione con i dati nel Tidyverse

Preparing Video For Download...