DAP_Final/model.R at main · maggie-lyman/DAP_Final · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
## Regression Analysis

#Run a regression of ADU count on median gross rent of each census tract within an ADU zone. Control for census tract population is included.
print(summary(lm(count ~ median_gross_rent + total_pop,
           data = df_tract_counts_census)))

#Plot regression of rents on number of ADUs per census tract
rent_regression <- ggplot(data = df_tract_counts_census,
                          aes(x = median_gross_rent,
                              y = count)) +
  geom_point(fill = "#33756D",
             color = "#33756D",
             shape = 21,
             size = 2) +
  geom_smooth(method = "lm",
              se = FALSE,
              color = "red") +
  scale_y_continuous(breaks = c(0, 5, 10, 15),
                     expand = expansion(mult = c(0.1))) +
  scale_x_continuous(breaks = c(0, 1000, 2000, 3000),
                     labels = c("$0", "$1,000", "$2,000", "$3,000"),
                     expand = expansion(mult = c(0.1))) +
  labs(title = "Median gross rent is correlated with number of ADU permits",
       subtitle = "Number of ADU permits compared to median gross rent with linear regression model, by ADU zone census tract",
       y = "Number of ADU Permits",
       x = "Median Gross Rent",
       caption = "Source: City of Chicago Data Portal | 2022 American Community Survey") +
  theme_minimal() +
  theme(
    text = element_text(color = "#22211d"),
    plot.background = element_rect(fill = "#f5f5f2",
                                   color = NA),
    panel.background = element_rect(fill = "#f5f5f2",
                                    color = NA),
    plot.title = element_text(size= 16,
                              face = "bold",
                              color = "#4e4d47"),
    axis.text = element_text(size= 10),
    axis.title = element_text(size= 12,
                              face = "bold"),
    plot.subtitle = element_text(size= 12,
                                 color = "#4e4d47"),
    plot.caption = element_text(size=10,
                                color = "#4e4d47"))

print(rent_regression)