Figure 3.1 code does not work

[RaymondBalise]

The code for figure 3.1 does not render properly.

There is a known bug in boom::augment.lm() that causes this:
tidymodels/broom#937

This shows the guts of the problem:

library(broom)
aModel <- augment(lm(log(Sale_Price) ~ Year_Built, data = ames_train))
names(aModel)
aModel <- augment(lm(Sale_Price ~ Year_Built, data = ames_train))
names(aModel)

This code works:

library(tidyverse)
ames <- AmesHousing::make_ames() %>% 
  mutate(log_Sale_Price = log(Sale_Price))

# Load and split the Ames housing data using stratified sampling
set.seed(123)  # for reproducibility
split  <- rsample::initial_split(ames, prop = 0.7, strata = "Sale_Price")
ames_train  <- rsample::training(split)
ames_test   <- rsample::testing(split)

models <- c("Non-log transformed model residuals", 
            "Log transformed model residuals")
list(
  m1 = lm(Sale_Price ~ Year_Built, data = ames_train) ,
  m2 = lm(log_Sale_Price ~ Year_Built, data = ames_train)
) %>%
  map2_dfr(models, ~ broom::augment(.x) %>% mutate(model = .y)) %>%
  ggplot(aes(.resid)) +
  geom_histogram(bins = 75) +
  facet_wrap(~ model, scales = "free") +
  ylab(NULL) +
  xlab("Residuals")