-
Notifications
You must be signed in to change notification settings - Fork 0
Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
- Loading branch information
Daniel Berry
authored and
Daniel Berry
committed
Nov 22, 2016
1 parent
5c6b580
commit be43b30
Showing
11 changed files
with
688 additions
and
26 deletions.
There are no files selected for viewing
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,133 @@ | ||
| --- | ||
| title: "Food Deserts in Chicago" | ||
| author: "Daniel Berry" | ||
| output: revealjs::revealjs_presentation | ||
|
|
||
| --- | ||
|
|
||
| # Introduction | ||
|
|
||
| ## What is a Food Desert? | ||
| - In general: a place where it is more difficult to access healthy food | ||
| - Defiition used for this work: a city block located more than 1 mile from a supermarket | ||
| - Supermarket defined as a grocery store larger than 10000 sq ft | ||
| - Distance is great circle distance between center of grocery store and center of city block | ||
| - Other definitions exist that also cover rural areas and take into account car ownership (harder to travel w/o a car) | ||
|
|
||
| ## Where are food deserts in Chicago | ||
|
|
||
| ```{r, out.width = "600px", echo = FALSE} | ||
| knitr::include_graphics("deserts_plot.png") | ||
| ``` | ||
|
|
||
| # Chicago Demographics | ||
|
|
||
| ## Distribution of Black People | ||
|
|
||
| ```{r, out.width = "600px", echo = FALSE} | ||
| knitr::include_graphics("pct_black_plot.png") | ||
| ``` | ||
|
|
||
| ## Distribution of White People | ||
|
|
||
| ```{r, out.width = "600px", echo = FALSE} | ||
| knitr::include_graphics("pct_white_plot.png") | ||
| ``` | ||
|
|
||
| ## Income | ||
|
|
||
| ```{r, out.width = "600px", echo = FALSE} | ||
| knitr::include_graphics("income_plot.png") | ||
| ``` | ||
|
|
||
| ## Vacancy | ||
| ```{r, out.width = "600px", echo = FALSE} | ||
| knitr::include_graphics("vacant_plot.png") | ||
| ``` | ||
|
|
||
| # Data | ||
|
|
||
| ## Source | ||
| All data from the [Chicago Open Data Portal](https://data.cityofchicago.org). Several files: | ||
|
|
||
| - Crimes 2001 - Present | ||
| - 311 Service Requests - Vacant Buildings | ||
| - CTA Ridership Avg Weekly Boardings Oct 2010 | ||
| - City Block Population | ||
| - Public Health Statistics Selected Indicators | ||
| - Census Data: Selected Socioeconomic Indicators | ||
| - Race by Community Area | ||
|
|
||
| # Models | ||
|
|
||
| ## Complete Pooling | ||
|
|
||
| $$ y_{ij} = \text{logit}^{-1}\left( \alpha + X_{B}\beta_{B} \right) $$ | ||
|
|
||
| ## No Pooling | ||
|
|
||
| $$ y_{i} = \text{logit}^{-1}\left( \alpha + X_{B}\beta_{B} + \gamma_j \right) $$ | ||
|
|
||
| ## Partial Pooling | ||
|
|
||
| $$ y_{i} = \text{logit}^{-1}\left( \alpha_{j[i]} + X_{B}\beta_{B} \right) $$ | ||
|
|
||
| Where $\alpha_j \sim N(\mu_\alpha, \sigma^2_\alpha)$ | ||
|
|
||
| ## Hierarchical | ||
|
|
||
| $$ y_{i} = \text{logit}^{-1}\left( \alpha_{j[i]} + X_{B}\beta_{B} + \epsilon_i \right) $$ | ||
| Where $\alpha_j \sim N(X_N \beta_N, \sigma^2_\alpha)$ | ||
|
|
||
| # Model summaries | ||
|
|
||
| ## Complete Pooling | ||
|
|
||
| ```{r, echo = F} | ||
| library(lme4) | ||
| load('cp') | ||
| summary(cp) | ||
| ``` | ||
|
|
||
| ## No Pooling | ||
|
|
||
| ```{r, echo = F} | ||
| load('np') | ||
| summary(np) | ||
| ``` | ||
|
|
||
| ## Partial Pooling | ||
|
|
||
| ```{r, echo = F} | ||
| load('pp') | ||
| summary(pp) | ||
| ``` | ||
|
|
||
| ## Hierarchical | ||
|
|
||
| ```{r, echo = F} | ||
| load('mlm') | ||
| summary(mlm) | ||
| ``` | ||
|
|
||
| # Results | ||
|
|
||
| ## Was pooling effective? | ||
|
|
||
| - Variance ratio $\approx 3$ indicates much higher variability within a neighborhood than between neighborhoods. | ||
|
|
||
| - AICs for random intercept models were higher than no pooling model. | ||
|
|
||
| - Cross validated MSEs (or [Brier Scores](https://en.wikipedia.org/wiki/Brier_score)) were a way to quantify how accurate we are on previously unseen city blocks within a neighborhood: | ||
| - Complete Pooling: 0.07382216 | ||
| - No Pooling: 0.05328587 | ||
| - Partial Pooling: 0.05329956 | ||
| - Hierarchical: 0.05323632 | ||
|
|
||
| ## Thoughts | ||
|
|
||
| - Models are an improvement over just using neighborhood level variables, but I'm not convinced that the hierarchical model is an improvement over simpler no pooling or partial pooling. | ||
|
|
||
| - Unfortunately this project doesn't really give us any more information about the causes of food deserts that we didn't really know before. More importantly this project doesn't help resolve the issue at all. | ||
|
|
||
| # Questions? |
Large diffs are not rendered by default.
Oops, something went wrong.