Cours regression + code

02_Regression/lm.R

@@ -0,0 +1,29 @@
+library(tidyverse)
+library(broom)
+
+data(iris)
+
+fit <- lm(Sepal.Length ~ Petal.Length, data = iris)
+
+fit
+
+fit %>% summary
+
+fit %>% glance
+
+fit %>% tidy
+
+fit %>% augment
+
+fit %>% plot
+
+tibble(Petal.Length = 1:7,
+       Sepal.Length = predict(fit, tibble(Petal.Length = Petal.Length))) -> pred
+
+iris %>%
+  ggplot() +
+  aes(x = Petal.Length, y = Sepal.Length) +
+  geom_point() +
+  # geom_abline(intercept = fit %>% tidy %>% filter(term == "(Intercept)") %>% pull(estimate),
+              # slope = fit %>% tidy %>% filter(term == "Petal.Length") %>% pull(estimate)) +
+  geom_point(data = pred, aes(x = Petal.Length, y = Sepal.Length), color = "red")