diff --git a/INSTALLATION-es.md b/INSTALLATION-es.md
index 39796e06..46e0aa64 100644
--- a/INSTALLATION-es.md
+++ b/INSTALLATION-es.md
@@ -55,6 +55,6 @@ Si bien se puede usar cualquier IDE, describimos cómo configurar VS Code
 
 Para el desarrollo de front-end, lea el [Client README](client/README.md).
 
-Para el desorrollo de back-end, lea el [Score README](score/README.md).
+Para el desorrollo de back-end, lea el [Data Pipeline README](data/data-pipeline/README.md).
 
 Para la implementación de software, lea el [Infrastructure README](infrastructure/README.md).
\ No newline at end of file
diff --git a/INSTALLATION.md b/INSTALLATION.md
index 850d69db..9b19fea5 100644
--- a/INSTALLATION.md
+++ b/INSTALLATION.md
@@ -78,7 +78,7 @@ There is also a docker-compose file that will eventually include everything you
 If you wish to check out our client-side code (i.e. the GatsbyJS statically generated website with the map), check out the `client` directory and its [README](client/README.md).
 
 **Running the backend**
-If you want to run the data pipeline with ETL and score generation, check out the `score` directory and its [README](score/README.md).
+If you want to run the data pipeline with ETL and score generation, check out the `data/data-pipeline` directory and its [README](data/data-pipeline/README.md).
 
 **Deployment**
 For core team contributors working on deployment, check out the `infrastructure` directory and its [README](infrastructure/README.md) for information on deploying the backend to AWS.
diff --git a/data/data-pipeline/README.md b/data/data-pipeline/README.md
index 6f1e1fb8..e55fa2c8 100644
--- a/data/data-pipeline/README.md
+++ b/data/data-pipeline/README.md
@@ -34,7 +34,7 @@ This application is used to compare experimental versions of the Justice40 score
 _**NOTE:** These scores **do not** represent final versions of the Justice40 scores and are merely used for comparative purposes. As a result, the specific input columns and formulas used to calculate them are likely to change over time._
 
 
-### Score comparison workflow
+### Score generation and comparison workflow
 
 The descriptions below provide a more detailed outline of what happens at each step of ETL and score calculation workflow.
 
@@ -58,8 +58,8 @@ TODO add mermaid diagram
 1. The `etl-run` command will execute the corresponding ETL script for each data source in `etl/sources/`. For example, `etl/sources/ejscreen/etl.py` is the ETL script for EJSCREEN data.
 1. Each ETL script will extract the data from its original source, then format the data into `.csv` files that get stored in the relevant folder in `data/dataset/`. For example, HUD Housing data is stored in `data/dataset/hud_housing/usa.csv`
 
-_**NOTE:** You have the option to pass the name of a specific data source to the `etl-run` command, which will limit the execution of the ETL process to that specific data source._
-_For example: `poetry run python application.py etl-run ejscreen` would only run the ETL process for EJSCREEN data._
+_**NOTE:** You have the option to pass the name of a specific data source to the `etl-run` command using the `-d` flag, which will limit the execution of the ETL process to that specific data source._
+_For example: `poetry run python application.py etl-run -d ejscreen` would only run the ETL process for EJSCREEN data._
 
 #### Step 2: Calculate the Justice40 score experiments
 
@@ -74,7 +74,31 @@ _For example: `poetry run python application.py etl-run ejscreen` would only run
 
 #### Step 3: Compare the Justice40 score experiments to other indices
 
-1. TODO: Describe the steps for this
+We are building a comparison tool to enable easy (or at least straightforward) comparison of the Justice40 score with other existing indices. The goal of having this is so that as we experiment and iterate with a scoring methodology, we can understand how our score overlaps with or differs from other indices that communities, nonprofits, and governmentss use to inform decision making.
+
+Right now, our comparison tool exists simply as a python notebook in `data/data-pipeline/ipython/scoring_comparison.ipynb`.
+
+To run this comparison tool:
+
+1. Make sure you've gone through the above steps to run the data ETL and score generation.
+1. From this directory (`data/data-pipeline`), navigate to the `ipython` directory: `cd ipython`.
+1. Ensure you have `pandoc` installed on your computer. If you're on a Mac, run `brew install pandoc`; for other OSes, see pandoc's [installation guide](https://pandoc.org/installing.html).
+1. Install the extra dependencies:
+```
+  pip install pypandoc
+  pip install requests
+  pip install us
+  pip install tqdm
+  pip install dynaconf
+  pip instal xlsxwriter
+```
+1. Start the notebooks: `jupyter notebook`
+1. In your browser, navigate to one of the URLs returned by the above command.
+1. Select `scoring_comparison.ipynb` from the options in your browser.
+1. Run through the steps in the notebook. You can step through them one at a time by clicking the "Run" button for each cell, or open the "Cell" menu and click "Run all" to run them all at once.
+1. Reports and spreadsheets generated by the comparison tool will be available in `data/data-pipeline/data/comparison_outputs`.
+
+*NOTE:* This may take several minutes or over an hour to fully execute and generate the reports.
 
 ### Data Sources
 
diff --git a/data/data-pipeline/ipython/scoring_comparison.ipynb b/data/data-pipeline/ipython/scoring_comparison.ipynb
index 125cbb83..76c47c45 100644
--- a/data/data-pipeline/ipython/scoring_comparison.ipynb
+++ b/data/data-pipeline/ipython/scoring_comparison.ipynb
@@ -9,11 +9,6 @@
    },
    "outputs": [],
    "source": [
-    "# Before running this script as it currently stands, you'll need to run these notebooks (in any order):\n",
-    "# * score_calc.ipynb\n",
-    "# * calenviroscreen_etl.ipynb\n",
-    "# * hud_recap_etl.ipynb\n",
-    "\n",
     "import collections\n",
     "import functools\n",
     "import IPython\n",
@@ -92,7 +87,243 @@
    "metadata": {
     "scrolled": false
    },
-   "outputs": [],
+   "outputs": [
+    {
+     "name": "stderr",
+     "output_type": "stream",
+     "text": [
+      "/usr/local/lib/python3.9/site-packages/IPython/core/interactiveshell.py:3169: DtypeWarning: Columns (87,88,90) have mixed types.Specify dtype option on import or set low_memory=False.\n",
+      "  has_raised = await self.run_ast_nodes(code_ast.body, cell_name,\n"
+     ]
+    },
+    {
+     "data": {
+      "text/html": [
+       "<div>\n",
+       "<style scoped>\n",
+       "    .dataframe tbody tr th:only-of-type {\n",
+       "        vertical-align: middle;\n",
+       "    }\n",
+       "\n",
+       "    .dataframe tbody tr th {\n",
+       "        vertical-align: top;\n",
+       "    }\n",
+       "\n",
+       "    .dataframe thead th {\n",
+       "        text-align: right;\n",
+       "    }\n",
+       "</style>\n",
+       "<table border=\"1\" class=\"dataframe\">\n",
+       "  <thead>\n",
+       "    <tr style=\"text-align: right;\">\n",
+       "      <th></th>\n",
+       "      <th>GEOID10</th>\n",
+       "      <th>Housing burden (percent)</th>\n",
+       "      <th>Total population</th>\n",
+       "      <th>Air toxics cancer risk</th>\n",
+       "      <th>Respiratory hazard index</th>\n",
+       "      <th>Diesel particulate matter</th>\n",
+       "      <th>Particulate matter (PM2.5)</th>\n",
+       "      <th>Ozone</th>\n",
+       "      <th>Traffic proximity and volume</th>\n",
+       "      <th>Proximity to RMP sites</th>\n",
+       "      <th>...</th>\n",
+       "      <th>Score D (top 25th percentile)</th>\n",
+       "      <th>Score E (percentile)</th>\n",
+       "      <th>Score E (top 25th percentile)</th>\n",
+       "      <th>GEOID</th>\n",
+       "      <th>State Abbreviation</th>\n",
+       "      <th>County Name</th>\n",
+       "      <th>State Code</th>\n",
+       "      <th>State Name</th>\n",
+       "      <th>GEOID10_TRACT</th>\n",
+       "      <th>GEOID10_STATE</th>\n",
+       "    </tr>\n",
+       "  </thead>\n",
+       "  <tbody>\n",
+       "    <tr>\n",
+       "      <th>0</th>\n",
+       "      <td>010010201001</td>\n",
+       "      <td>0.15</td>\n",
+       "      <td>692</td>\n",
+       "      <td>49.38</td>\n",
+       "      <td>0.79</td>\n",
+       "      <td>0.28</td>\n",
+       "      <td>10.00</td>\n",
+       "      <td>40.12</td>\n",
+       "      <td>91.02</td>\n",
+       "      <td>0.09</td>\n",
+       "      <td>...</td>\n",
+       "      <td>False</td>\n",
+       "      <td>0.35</td>\n",
+       "      <td>False</td>\n",
+       "      <td>1001</td>\n",
+       "      <td>AL</td>\n",
+       "      <td>Autauga County</td>\n",
+       "      <td>1.00</td>\n",
+       "      <td>Alabama</td>\n",
+       "      <td>01001020100</td>\n",
+       "      <td>01</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>1</th>\n",
+       "      <td>010010201002</td>\n",
+       "      <td>0.15</td>\n",
+       "      <td>1153</td>\n",
+       "      <td>49.38</td>\n",
+       "      <td>0.79</td>\n",
+       "      <td>0.28</td>\n",
+       "      <td>10.00</td>\n",
+       "      <td>40.12</td>\n",
+       "      <td>2.62</td>\n",
+       "      <td>0.07</td>\n",
+       "      <td>...</td>\n",
+       "      <td>False</td>\n",
+       "      <td>0.11</td>\n",
+       "      <td>False</td>\n",
+       "      <td>1001</td>\n",
+       "      <td>AL</td>\n",
+       "      <td>Baldwin County</td>\n",
+       "      <td>2.00</td>\n",
+       "      <td>Alaska</td>\n",
+       "      <td>01001020100</td>\n",
+       "      <td>01</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>2</th>\n",
+       "      <td>010010202001</td>\n",
+       "      <td>0.25</td>\n",
+       "      <td>1020</td>\n",
+       "      <td>50.32</td>\n",
+       "      <td>0.81</td>\n",
+       "      <td>0.30</td>\n",
+       "      <td>10.07</td>\n",
+       "      <td>40.22</td>\n",
+       "      <td>4.68</td>\n",
+       "      <td>0.08</td>\n",
+       "      <td>...</td>\n",
+       "      <td>False</td>\n",
+       "      <td>0.51</td>\n",
+       "      <td>False</td>\n",
+       "      <td>1001</td>\n",
+       "      <td>AL</td>\n",
+       "      <td>Barbour County</td>\n",
+       "      <td>4.00</td>\n",
+       "      <td>Arizona</td>\n",
+       "      <td>01001020200</td>\n",
+       "      <td>01</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>3</th>\n",
+       "      <td>010010202002</td>\n",
+       "      <td>0.25</td>\n",
+       "      <td>1152</td>\n",
+       "      <td>50.32</td>\n",
+       "      <td>0.81</td>\n",
+       "      <td>0.30</td>\n",
+       "      <td>10.07</td>\n",
+       "      <td>40.22</td>\n",
+       "      <td>218.65</td>\n",
+       "      <td>0.09</td>\n",
+       "      <td>...</td>\n",
+       "      <td>False</td>\n",
+       "      <td>0.59</td>\n",
+       "      <td>False</td>\n",
+       "      <td>1001</td>\n",
+       "      <td>AL</td>\n",
+       "      <td>Bibb County</td>\n",
+       "      <td>5.00</td>\n",
+       "      <td>Arkansas</td>\n",
+       "      <td>01001020200</td>\n",
+       "      <td>01</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>4</th>\n",
+       "      <td>010010203001</td>\n",
+       "      <td>0.21</td>\n",
+       "      <td>2555</td>\n",
+       "      <td>50.77</td>\n",
+       "      <td>0.82</td>\n",
+       "      <td>0.36</td>\n",
+       "      <td>10.12</td>\n",
+       "      <td>40.31</td>\n",
+       "      <td>69.64</td>\n",
+       "      <td>0.08</td>\n",
+       "      <td>...</td>\n",
+       "      <td>False</td>\n",
+       "      <td>0.47</td>\n",
+       "      <td>False</td>\n",
+       "      <td>1001</td>\n",
+       "      <td>AL</td>\n",
+       "      <td>Blount County</td>\n",
+       "      <td>6.00</td>\n",
+       "      <td>California</td>\n",
+       "      <td>01001020300</td>\n",
+       "      <td>01</td>\n",
+       "    </tr>\n",
+       "  </tbody>\n",
+       "</table>\n",
+       "<p>5 rows × 93 columns</p>\n",
+       "</div>"
+      ],
+      "text/plain": [
+       "        GEOID10  Housing burden (percent)  Total population  \\\n",
+       "0  010010201001                      0.15               692   \n",
+       "1  010010201002                      0.15              1153   \n",
+       "2  010010202001                      0.25              1020   \n",
+       "3  010010202002                      0.25              1152   \n",
+       "4  010010203001                      0.21              2555   \n",
+       "\n",
+       "   Air toxics cancer risk  Respiratory hazard index  \\\n",
+       "0                   49.38                      0.79   \n",
+       "1                   49.38                      0.79   \n",
+       "2                   50.32                      0.81   \n",
+       "3                   50.32                      0.81   \n",
+       "4                   50.77                      0.82   \n",
+       "\n",
+       "   Diesel particulate matter  Particulate matter (PM2.5)  Ozone  \\\n",
+       "0                       0.28                       10.00  40.12   \n",
+       "1                       0.28                       10.00  40.12   \n",
+       "2                       0.30                       10.07  40.22   \n",
+       "3                       0.30                       10.07  40.22   \n",
+       "4                       0.36                       10.12  40.31   \n",
+       "\n",
+       "   Traffic proximity and volume  Proximity to RMP sites  ...  \\\n",
+       "0                         91.02                    0.09  ...   \n",
+       "1                          2.62                    0.07  ...   \n",
+       "2                          4.68                    0.08  ...   \n",
+       "3                        218.65                    0.09  ...   \n",
+       "4                         69.64                    0.08  ...   \n",
+       "\n",
+       "   Score D (top 25th percentile)  Score E (percentile)  \\\n",
+       "0                          False                  0.35   \n",
+       "1                          False                  0.11   \n",
+       "2                          False                  0.51   \n",
+       "3                          False                  0.59   \n",
+       "4                          False                  0.47   \n",
+       "\n",
+       "   Score E (top 25th percentile)  GEOID  State Abbreviation     County Name  \\\n",
+       "0                          False   1001                  AL  Autauga County   \n",
+       "1                          False   1001                  AL  Baldwin County   \n",
+       "2                          False   1001                  AL  Barbour County   \n",
+       "3                          False   1001                  AL     Bibb County   \n",
+       "4                          False   1001                  AL   Blount County   \n",
+       "\n",
+       "   State Code  State Name  GEOID10_TRACT  GEOID10_STATE  \n",
+       "0        1.00     Alabama    01001020100             01  \n",
+       "1        2.00      Alaska    01001020100             01  \n",
+       "2        4.00     Arizona    01001020200             01  \n",
+       "3        5.00    Arkansas    01001020200             01  \n",
+       "4        6.00  California    01001020300             01  \n",
+       "\n",
+       "[5 rows x 93 columns]"
+      ]
+     },
+     "execution_count": 3,
+     "metadata": {},
+     "output_type": "execute_result"
+    }
+   ],
    "source": [
     "# Load CEJST score data\n",
     "cejst_data_path = DATA_DIR / \"score\" / \"csv\" / \"full\" / \"usa.csv\"\n",
@@ -119,7 +350,228 @@
    "metadata": {
     "scrolled": false
    },
-   "outputs": [],
+   "outputs": [
+    {
+     "data": {
+      "text/html": [
+       "<div>\n",
+       "<style scoped>\n",
+       "    .dataframe tbody tr th:only-of-type {\n",
+       "        vertical-align: middle;\n",
+       "    }\n",
+       "\n",
+       "    .dataframe tbody tr th {\n",
+       "        vertical-align: top;\n",
+       "    }\n",
+       "\n",
+       "    .dataframe thead th {\n",
+       "        text-align: right;\n",
+       "    }\n",
+       "</style>\n",
+       "<table border=\"1\" class=\"dataframe\">\n",
+       "  <thead>\n",
+       "    <tr style=\"text-align: right;\">\n",
+       "      <th></th>\n",
+       "      <th>GEOID10_TRACT</th>\n",
+       "      <th>Total Population</th>\n",
+       "      <th>California County</th>\n",
+       "      <th>ZIP</th>\n",
+       "      <th>Nearby City \\r\\n(to help approximate location only)</th>\n",
+       "      <th>Longitude</th>\n",
+       "      <th>Latitude</th>\n",
+       "      <th>calenviroscreen_score</th>\n",
+       "      <th>calenviroscreen_percentile</th>\n",
+       "      <th>DRAFT CES 4.0\\r\\nPercentile Range</th>\n",
+       "      <th>...</th>\n",
+       "      <th>Poverty</th>\n",
+       "      <th>Poverty Pctl</th>\n",
+       "      <th>Unemployment</th>\n",
+       "      <th>Unemployment Pctl</th>\n",
+       "      <th>Housing Burden</th>\n",
+       "      <th>Housing Burden Pctl</th>\n",
+       "      <th>Pop. Char.</th>\n",
+       "      <th>Pop. Char. Score</th>\n",
+       "      <th>Pop. Char. Pctl</th>\n",
+       "      <th>calenviroscreen_priority_community</th>\n",
+       "    </tr>\n",
+       "  </thead>\n",
+       "  <tbody>\n",
+       "    <tr>\n",
+       "      <th>0</th>\n",
+       "      <td>06019001100</td>\n",
+       "      <td>2760</td>\n",
+       "      <td>Fresno</td>\n",
+       "      <td>93706</td>\n",
+       "      <td>Fresno</td>\n",
+       "      <td>-119.78</td>\n",
+       "      <td>36.71</td>\n",
+       "      <td>94.61</td>\n",
+       "      <td>100.00</td>\n",
+       "      <td>95-100% (highest scores)</td>\n",
+       "      <td>...</td>\n",
+       "      <td>76.60</td>\n",
+       "      <td>98.43</td>\n",
+       "      <td>16.20</td>\n",
+       "      <td>97.15</td>\n",
+       "      <td>30.70</td>\n",
+       "      <td>90.61</td>\n",
+       "      <td>93.73</td>\n",
+       "      <td>9.72</td>\n",
+       "      <td>99.87</td>\n",
+       "      <td>True</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>1</th>\n",
+       "      <td>06077000700</td>\n",
+       "      <td>4177</td>\n",
+       "      <td>San Joaquin</td>\n",
+       "      <td>95206</td>\n",
+       "      <td>Stockton</td>\n",
+       "      <td>-121.29</td>\n",
+       "      <td>37.94</td>\n",
+       "      <td>90.83</td>\n",
+       "      <td>99.99</td>\n",
+       "      <td>95-100% (highest scores)</td>\n",
+       "      <td>...</td>\n",
+       "      <td>70.60</td>\n",
+       "      <td>96.43</td>\n",
+       "      <td>18.50</td>\n",
+       "      <td>98.45</td>\n",
+       "      <td>35.20</td>\n",
+       "      <td>95.61</td>\n",
+       "      <td>93.40</td>\n",
+       "      <td>9.68</td>\n",
+       "      <td>99.84</td>\n",
+       "      <td>True</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>2</th>\n",
+       "      <td>06077000100</td>\n",
+       "      <td>4055</td>\n",
+       "      <td>San Joaquin</td>\n",
+       "      <td>95202</td>\n",
+       "      <td>Stockton</td>\n",
+       "      <td>-121.29</td>\n",
+       "      <td>37.95</td>\n",
+       "      <td>85.75</td>\n",
+       "      <td>99.97</td>\n",
+       "      <td>95-100% (highest scores)</td>\n",
+       "      <td>...</td>\n",
+       "      <td>81.80</td>\n",
+       "      <td>99.50</td>\n",
+       "      <td>17.90</td>\n",
+       "      <td>98.17</td>\n",
+       "      <td>36.40</td>\n",
+       "      <td>96.51</td>\n",
+       "      <td>95.71</td>\n",
+       "      <td>9.92</td>\n",
+       "      <td>99.97</td>\n",
+       "      <td>True</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>3</th>\n",
+       "      <td>06071001600</td>\n",
+       "      <td>5527</td>\n",
+       "      <td>San Bernardino</td>\n",
+       "      <td>91761</td>\n",
+       "      <td>Ontario</td>\n",
+       "      <td>-117.62</td>\n",
+       "      <td>34.06</td>\n",
+       "      <td>83.56</td>\n",
+       "      <td>99.96</td>\n",
+       "      <td>95-100% (highest scores)</td>\n",
+       "      <td>...</td>\n",
+       "      <td>67.10</td>\n",
+       "      <td>94.82</td>\n",
+       "      <td>6.70</td>\n",
+       "      <td>57.20</td>\n",
+       "      <td>32.10</td>\n",
+       "      <td>92.65</td>\n",
+       "      <td>80.59</td>\n",
+       "      <td>8.36</td>\n",
+       "      <td>93.06</td>\n",
+       "      <td>True</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>4</th>\n",
+       "      <td>06037204920</td>\n",
+       "      <td>2639</td>\n",
+       "      <td>Los Angeles</td>\n",
+       "      <td>90023</td>\n",
+       "      <td>Los Angeles</td>\n",
+       "      <td>-118.20</td>\n",
+       "      <td>34.02</td>\n",
+       "      <td>82.90</td>\n",
+       "      <td>99.95</td>\n",
+       "      <td>95-100% (highest scores)</td>\n",
+       "      <td>...</td>\n",
+       "      <td>64.90</td>\n",
+       "      <td>93.51</td>\n",
+       "      <td>5.60</td>\n",
+       "      <td>43.81</td>\n",
+       "      <td>25.00</td>\n",
+       "      <td>77.95</td>\n",
+       "      <td>83.95</td>\n",
+       "      <td>8.70</td>\n",
+       "      <td>95.78</td>\n",
+       "      <td>True</td>\n",
+       "    </tr>\n",
+       "  </tbody>\n",
+       "</table>\n",
+       "<p>5 rows × 59 columns</p>\n",
+       "</div>"
+      ],
+      "text/plain": [
+       "  GEOID10_TRACT   Total Population California County    ZIP  \\\n",
+       "0   06019001100               2760           Fresno   93706   \n",
+       "1   06077000700               4177       San Joaquin  95206   \n",
+       "2   06077000100               4055       San Joaquin  95202   \n",
+       "3   06071001600               5527    San Bernardino  91761   \n",
+       "4   06037204920               2639       Los Angeles  90023   \n",
+       "\n",
+       "  Nearby City \\r\\n(to help approximate location only)  Longitude  Latitude  \\\n",
+       "0                                             Fresno     -119.78     36.71   \n",
+       "1                                           Stockton     -121.29     37.94   \n",
+       "2                                           Stockton     -121.29     37.95   \n",
+       "3                                            Ontario     -117.62     34.06   \n",
+       "4                                        Los Angeles     -118.20     34.02   \n",
+       "\n",
+       "   calenviroscreen_score  calenviroscreen_percentile  \\\n",
+       "0                  94.61                      100.00   \n",
+       "1                  90.83                       99.99   \n",
+       "2                  85.75                       99.97   \n",
+       "3                  83.56                       99.96   \n",
+       "4                  82.90                       99.95   \n",
+       "\n",
+       "  DRAFT CES 4.0\\r\\nPercentile Range  ...  Poverty  Poverty Pctl  Unemployment  \\\n",
+       "0          95-100% (highest scores)  ...    76.60         98.43         16.20   \n",
+       "1          95-100% (highest scores)  ...    70.60         96.43         18.50   \n",
+       "2          95-100% (highest scores)  ...    81.80         99.50         17.90   \n",
+       "3          95-100% (highest scores)  ...    67.10         94.82          6.70   \n",
+       "4          95-100% (highest scores)  ...    64.90         93.51          5.60   \n",
+       "\n",
+       "   Unemployment Pctl  Housing Burden  Housing Burden Pctl  Pop. Char.   \\\n",
+       "0              97.15           30.70                90.61        93.73   \n",
+       "1              98.45           35.20                95.61        93.40   \n",
+       "2              98.17           36.40                96.51        95.71   \n",
+       "3              57.20           32.10                92.65        80.59   \n",
+       "4              43.81           25.00                77.95        83.95   \n",
+       "\n",
+       "   Pop. Char. Score  Pop. Char. Pctl  calenviroscreen_priority_community  \n",
+       "0              9.72            99.87                                True  \n",
+       "1              9.68            99.84                                True  \n",
+       "2              9.92            99.97                                True  \n",
+       "3              8.36            93.06                                True  \n",
+       "4              8.70            95.78                                True  \n",
+       "\n",
+       "[5 rows x 59 columns]"
+      ]
+     },
+     "execution_count": 4,
+     "metadata": {},
+     "output_type": "execute_result"
+    }
+   ],
    "source": [
     "# Load CalEnviroScreen 4.0\n",
     "CALENVIROSCREEN_SCORE_FIELD = \"calenviroscreen_score\"\n",
@@ -146,7 +598,168 @@
    "metadata": {
     "scrolled": true
    },
-   "outputs": [],
+   "outputs": [
+    {
+     "data": {
+      "text/html": [
+       "<div>\n",
+       "<style scoped>\n",
+       "    .dataframe tbody tr th:only-of-type {\n",
+       "        vertical-align: middle;\n",
+       "    }\n",
+       "\n",
+       "    .dataframe tbody tr th {\n",
+       "        vertical-align: top;\n",
+       "    }\n",
+       "\n",
+       "    .dataframe thead th {\n",
+       "        text-align: right;\n",
+       "    }\n",
+       "</style>\n",
+       "<table border=\"1\" class=\"dataframe\">\n",
+       "  <thead>\n",
+       "    <tr style=\"text-align: right;\">\n",
+       "      <th></th>\n",
+       "      <th>FID</th>\n",
+       "      <th>GEOID10_TRACT</th>\n",
+       "      <th>STATE</th>\n",
+       "      <th>STUSAB</th>\n",
+       "      <th>STATE_NAME</th>\n",
+       "      <th>COUNTY</th>\n",
+       "      <th>COUNTY_NAME</th>\n",
+       "      <th>CNTY_FIPS</th>\n",
+       "      <th>TRACT</th>\n",
+       "      <th>RCAP_90</th>\n",
+       "      <th>RCAP_00</th>\n",
+       "      <th>RCAP_10</th>\n",
+       "      <th>hud_recap_priority_community</th>\n",
+       "      <th>SHAPE_Length</th>\n",
+       "      <th>SHAPE_Area</th>\n",
+       "    </tr>\n",
+       "  </thead>\n",
+       "  <tbody>\n",
+       "    <tr>\n",
+       "      <th>0</th>\n",
+       "      <td>29993</td>\n",
+       "      <td>01001020100</td>\n",
+       "      <td>1</td>\n",
+       "      <td>AL</td>\n",
+       "      <td>Alabama</td>\n",
+       "      <td>1</td>\n",
+       "      <td>Autauga</td>\n",
+       "      <td>1001</td>\n",
+       "      <td>20100</td>\n",
+       "      <td>0.00</td>\n",
+       "      <td>0.00</td>\n",
+       "      <td>0.00</td>\n",
+       "      <td>False</td>\n",
+       "      <td>0.15</td>\n",
+       "      <td>0.00</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>1</th>\n",
+       "      <td>30627</td>\n",
+       "      <td>01001020200</td>\n",
+       "      <td>1</td>\n",
+       "      <td>AL</td>\n",
+       "      <td>Alabama</td>\n",
+       "      <td>1</td>\n",
+       "      <td>Autauga</td>\n",
+       "      <td>1001</td>\n",
+       "      <td>20200</td>\n",
+       "      <td>0.00</td>\n",
+       "      <td>0.00</td>\n",
+       "      <td>0.00</td>\n",
+       "      <td>False</td>\n",
+       "      <td>0.09</td>\n",
+       "      <td>0.00</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>2</th>\n",
+       "      <td>29992</td>\n",
+       "      <td>01001020300</td>\n",
+       "      <td>1</td>\n",
+       "      <td>AL</td>\n",
+       "      <td>Alabama</td>\n",
+       "      <td>1</td>\n",
+       "      <td>Autauga</td>\n",
+       "      <td>1001</td>\n",
+       "      <td>20300</td>\n",
+       "      <td>0.00</td>\n",
+       "      <td>0.00</td>\n",
+       "      <td>0.00</td>\n",
+       "      <td>False</td>\n",
+       "      <td>0.10</td>\n",
+       "      <td>0.00</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>3</th>\n",
+       "      <td>30079</td>\n",
+       "      <td>01001020400</td>\n",
+       "      <td>1</td>\n",
+       "      <td>AL</td>\n",
+       "      <td>Alabama</td>\n",
+       "      <td>1</td>\n",
+       "      <td>Autauga</td>\n",
+       "      <td>1001</td>\n",
+       "      <td>20400</td>\n",
+       "      <td>0.00</td>\n",
+       "      <td>0.00</td>\n",
+       "      <td>0.00</td>\n",
+       "      <td>False</td>\n",
+       "      <td>0.12</td>\n",
+       "      <td>0.00</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>4</th>\n",
+       "      <td>30078</td>\n",
+       "      <td>01001020500</td>\n",
+       "      <td>1</td>\n",
+       "      <td>AL</td>\n",
+       "      <td>Alabama</td>\n",
+       "      <td>1</td>\n",
+       "      <td>Autauga</td>\n",
+       "      <td>1001</td>\n",
+       "      <td>20500</td>\n",
+       "      <td>0.00</td>\n",
+       "      <td>0.00</td>\n",
+       "      <td>0.00</td>\n",
+       "      <td>False</td>\n",
+       "      <td>0.16</td>\n",
+       "      <td>0.00</td>\n",
+       "    </tr>\n",
+       "  </tbody>\n",
+       "</table>\n",
+       "</div>"
+      ],
+      "text/plain": [
+       "     FID GEOID10_TRACT  STATE STUSAB STATE_NAME  COUNTY COUNTY_NAME  \\\n",
+       "0  29993   01001020100      1     AL    Alabama       1     Autauga   \n",
+       "1  30627   01001020200      1     AL    Alabama       1     Autauga   \n",
+       "2  29992   01001020300      1     AL    Alabama       1     Autauga   \n",
+       "3  30079   01001020400      1     AL    Alabama       1     Autauga   \n",
+       "4  30078   01001020500      1     AL    Alabama       1     Autauga   \n",
+       "\n",
+       "   CNTY_FIPS  TRACT  RCAP_90  RCAP_00  RCAP_10  hud_recap_priority_community  \\\n",
+       "0       1001  20100     0.00     0.00     0.00                         False   \n",
+       "1       1001  20200     0.00     0.00     0.00                         False   \n",
+       "2       1001  20300     0.00     0.00     0.00                         False   \n",
+       "3       1001  20400     0.00     0.00     0.00                         False   \n",
+       "4       1001  20500     0.00     0.00     0.00                         False   \n",
+       "\n",
+       "   SHAPE_Length  SHAPE_Area  \n",
+       "0          0.15        0.00  \n",
+       "1          0.09        0.00  \n",
+       "2          0.10        0.00  \n",
+       "3          0.12        0.00  \n",
+       "4          0.16        0.00  "
+      ]
+     },
+     "execution_count": 5,
+     "metadata": {},
+     "output_type": "execute_result"
+    }
+   ],
    "source": [
     "# Load HUD data\n",
     "hud_recap_data_path = DATA_DIR / \"dataset\" / \"hud_recap\" / \"usa.csv\"\n",
@@ -164,7 +777,228 @@
    "metadata": {
     "scrolled": true
    },
-   "outputs": [],
+   "outputs": [
+    {
+     "data": {
+      "text/html": [
+       "<div>\n",
+       "<style scoped>\n",
+       "    .dataframe tbody tr th:only-of-type {\n",
+       "        vertical-align: middle;\n",
+       "    }\n",
+       "\n",
+       "    .dataframe tbody tr th {\n",
+       "        vertical-align: top;\n",
+       "    }\n",
+       "\n",
+       "    .dataframe thead th {\n",
+       "        text-align: right;\n",
+       "    }\n",
+       "</style>\n",
+       "<table border=\"1\" class=\"dataframe\">\n",
+       "  <thead>\n",
+       "    <tr style=\"text-align: right;\">\n",
+       "      <th></th>\n",
+       "      <th>GEOID10_TRACT</th>\n",
+       "      <th>Total Population</th>\n",
+       "      <th>California County</th>\n",
+       "      <th>ZIP</th>\n",
+       "      <th>Nearby City \\r\\n(to help approximate location only)</th>\n",
+       "      <th>Longitude</th>\n",
+       "      <th>Latitude</th>\n",
+       "      <th>calenviroscreen_score</th>\n",
+       "      <th>calenviroscreen_percentile</th>\n",
+       "      <th>DRAFT CES 4.0\\r\\nPercentile Range</th>\n",
+       "      <th>...</th>\n",
+       "      <th>COUNTY</th>\n",
+       "      <th>COUNTY_NAME</th>\n",
+       "      <th>CNTY_FIPS</th>\n",
+       "      <th>TRACT</th>\n",
+       "      <th>RCAP_90</th>\n",
+       "      <th>RCAP_00</th>\n",
+       "      <th>RCAP_10</th>\n",
+       "      <th>hud_recap_priority_community</th>\n",
+       "      <th>SHAPE_Length</th>\n",
+       "      <th>SHAPE_Area</th>\n",
+       "    </tr>\n",
+       "  </thead>\n",
+       "  <tbody>\n",
+       "    <tr>\n",
+       "      <th>0</th>\n",
+       "      <td>06019001100</td>\n",
+       "      <td>2760.00</td>\n",
+       "      <td>Fresno</td>\n",
+       "      <td>93706.00</td>\n",
+       "      <td>Fresno</td>\n",
+       "      <td>-119.78</td>\n",
+       "      <td>36.71</td>\n",
+       "      <td>94.61</td>\n",
+       "      <td>100.00</td>\n",
+       "      <td>95-100% (highest scores)</td>\n",
+       "      <td>...</td>\n",
+       "      <td>19</td>\n",
+       "      <td>Fresno</td>\n",
+       "      <td>6019</td>\n",
+       "      <td>1100</td>\n",
+       "      <td>0.00</td>\n",
+       "      <td>1.00</td>\n",
+       "      <td>1.00</td>\n",
+       "      <td>True</td>\n",
+       "      <td>0.09</td>\n",
+       "      <td>0.00</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>1</th>\n",
+       "      <td>06077000700</td>\n",
+       "      <td>4177.00</td>\n",
+       "      <td>San Joaquin</td>\n",
+       "      <td>95206.00</td>\n",
+       "      <td>Stockton</td>\n",
+       "      <td>-121.29</td>\n",
+       "      <td>37.94</td>\n",
+       "      <td>90.83</td>\n",
+       "      <td>99.99</td>\n",
+       "      <td>95-100% (highest scores)</td>\n",
+       "      <td>...</td>\n",
+       "      <td>77</td>\n",
+       "      <td>San Joaquin</td>\n",
+       "      <td>6077</td>\n",
+       "      <td>700</td>\n",
+       "      <td>0.00</td>\n",
+       "      <td>0.00</td>\n",
+       "      <td>0.00</td>\n",
+       "      <td>True</td>\n",
+       "      <td>0.07</td>\n",
+       "      <td>0.00</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>2</th>\n",
+       "      <td>06077000100</td>\n",
+       "      <td>4055.00</td>\n",
+       "      <td>San Joaquin</td>\n",
+       "      <td>95202.00</td>\n",
+       "      <td>Stockton</td>\n",
+       "      <td>-121.29</td>\n",
+       "      <td>37.95</td>\n",
+       "      <td>85.75</td>\n",
+       "      <td>99.97</td>\n",
+       "      <td>95-100% (highest scores)</td>\n",
+       "      <td>...</td>\n",
+       "      <td>77</td>\n",
+       "      <td>San Joaquin</td>\n",
+       "      <td>6077</td>\n",
+       "      <td>100</td>\n",
+       "      <td>1.00</td>\n",
+       "      <td>1.00</td>\n",
+       "      <td>1.00</td>\n",
+       "      <td>True</td>\n",
+       "      <td>0.06</td>\n",
+       "      <td>0.00</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>3</th>\n",
+       "      <td>06071001600</td>\n",
+       "      <td>5527.00</td>\n",
+       "      <td>San Bernardino</td>\n",
+       "      <td>91761.00</td>\n",
+       "      <td>Ontario</td>\n",
+       "      <td>-117.62</td>\n",
+       "      <td>34.06</td>\n",
+       "      <td>83.56</td>\n",
+       "      <td>99.96</td>\n",
+       "      <td>95-100% (highest scores)</td>\n",
+       "      <td>...</td>\n",
+       "      <td>71</td>\n",
+       "      <td>San Bernardino</td>\n",
+       "      <td>6071</td>\n",
+       "      <td>1600</td>\n",
+       "      <td>0.00</td>\n",
+       "      <td>0.00</td>\n",
+       "      <td>0.00</td>\n",
+       "      <td>True</td>\n",
+       "      <td>0.25</td>\n",
+       "      <td>0.00</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>4</th>\n",
+       "      <td>06037204920</td>\n",
+       "      <td>2639.00</td>\n",
+       "      <td>Los Angeles</td>\n",
+       "      <td>90023.00</td>\n",
+       "      <td>Los Angeles</td>\n",
+       "      <td>-118.20</td>\n",
+       "      <td>34.02</td>\n",
+       "      <td>82.90</td>\n",
+       "      <td>99.95</td>\n",
+       "      <td>95-100% (highest scores)</td>\n",
+       "      <td>...</td>\n",
+       "      <td>37</td>\n",
+       "      <td>Los Angeles</td>\n",
+       "      <td>6037</td>\n",
+       "      <td>204920</td>\n",
+       "      <td>0.00</td>\n",
+       "      <td>0.00</td>\n",
+       "      <td>0.00</td>\n",
+       "      <td>False</td>\n",
+       "      <td>0.04</td>\n",
+       "      <td>0.00</td>\n",
+       "    </tr>\n",
+       "  </tbody>\n",
+       "</table>\n",
+       "<p>5 rows × 73 columns</p>\n",
+       "</div>"
+      ],
+      "text/plain": [
+       "  GEOID10_TRACT   Total Population California County      ZIP  \\\n",
+       "0   06019001100            2760.00           Fresno  93706.00   \n",
+       "1   06077000700            4177.00       San Joaquin 95206.00   \n",
+       "2   06077000100            4055.00       San Joaquin 95202.00   \n",
+       "3   06071001600            5527.00    San Bernardino 91761.00   \n",
+       "4   06037204920            2639.00       Los Angeles 90023.00   \n",
+       "\n",
+       "  Nearby City \\r\\n(to help approximate location only)  Longitude  Latitude  \\\n",
+       "0                                             Fresno     -119.78     36.71   \n",
+       "1                                           Stockton     -121.29     37.94   \n",
+       "2                                           Stockton     -121.29     37.95   \n",
+       "3                                            Ontario     -117.62     34.06   \n",
+       "4                                        Los Angeles     -118.20     34.02   \n",
+       "\n",
+       "   calenviroscreen_score  calenviroscreen_percentile  \\\n",
+       "0                  94.61                      100.00   \n",
+       "1                  90.83                       99.99   \n",
+       "2                  85.75                       99.97   \n",
+       "3                  83.56                       99.96   \n",
+       "4                  82.90                       99.95   \n",
+       "\n",
+       "  DRAFT CES 4.0\\r\\nPercentile Range  ...  COUNTY     COUNTY_NAME  CNTY_FIPS  \\\n",
+       "0          95-100% (highest scores)  ...      19          Fresno       6019   \n",
+       "1          95-100% (highest scores)  ...      77     San Joaquin       6077   \n",
+       "2          95-100% (highest scores)  ...      77     San Joaquin       6077   \n",
+       "3          95-100% (highest scores)  ...      71  San Bernardino       6071   \n",
+       "4          95-100% (highest scores)  ...      37     Los Angeles       6037   \n",
+       "\n",
+       "    TRACT  RCAP_90  RCAP_00  RCAP_10  hud_recap_priority_community  \\\n",
+       "0    1100     0.00     1.00     1.00                          True   \n",
+       "1     700     0.00     0.00     0.00                          True   \n",
+       "2     100     1.00     1.00     1.00                          True   \n",
+       "3    1600     0.00     0.00     0.00                          True   \n",
+       "4  204920     0.00     0.00     0.00                         False   \n",
+       "\n",
+       "   SHAPE_Length  SHAPE_Area  \n",
+       "0          0.09        0.00  \n",
+       "1          0.07        0.00  \n",
+       "2          0.06        0.00  \n",
+       "3          0.25        0.00  \n",
+       "4          0.04        0.00  \n",
+       "\n",
+       "[5 rows x 73 columns]"
+      ]
+     },
+     "execution_count": 6,
+     "metadata": {},
+     "output_type": "execute_result"
+    }
+   ],
    "source": [
     "# Join all dataframes that use tracts\n",
     "census_tract_dfs = [calenviroscreen_df, hud_recap_df]\n",
@@ -176,7 +1010,9 @@
     "    census_tract_dfs,\n",
     ")\n",
     "\n",
-    "if census_tract_df[GEOID_TRACT_FIELD_NAME].str.len().unique() != [11]:\n",
+    "tract_values = census_tract_df[GEOID_TRACT_FIELD_NAME].str.len().unique()\n",
+    "if any(tract_values != [11]):\n",
+    "    print(tract_values)\n",
     "    raise ValueError(\"Some of the census tract data has the wrong length.\")\n",
     "\n",
     "if len(census_tract_df) > 74134:\n",
@@ -192,7 +1028,228 @@
    "metadata": {
     "scrolled": false
    },
-   "outputs": [],
+   "outputs": [
+    {
+     "data": {
+      "text/html": [
+       "<div>\n",
+       "<style scoped>\n",
+       "    .dataframe tbody tr th:only-of-type {\n",
+       "        vertical-align: middle;\n",
+       "    }\n",
+       "\n",
+       "    .dataframe tbody tr th {\n",
+       "        vertical-align: top;\n",
+       "    }\n",
+       "\n",
+       "    .dataframe thead th {\n",
+       "        text-align: right;\n",
+       "    }\n",
+       "</style>\n",
+       "<table border=\"1\" class=\"dataframe\">\n",
+       "  <thead>\n",
+       "    <tr style=\"text-align: right;\">\n",
+       "      <th></th>\n",
+       "      <th>GEOID10</th>\n",
+       "      <th>Housing burden (percent)</th>\n",
+       "      <th>Total population</th>\n",
+       "      <th>Air toxics cancer risk</th>\n",
+       "      <th>Respiratory hazard index</th>\n",
+       "      <th>Diesel particulate matter</th>\n",
+       "      <th>Particulate matter (PM2.5)</th>\n",
+       "      <th>Ozone_x</th>\n",
+       "      <th>Traffic proximity and volume</th>\n",
+       "      <th>Proximity to RMP sites</th>\n",
+       "      <th>...</th>\n",
+       "      <th>COUNTY</th>\n",
+       "      <th>COUNTY_NAME</th>\n",
+       "      <th>CNTY_FIPS</th>\n",
+       "      <th>TRACT</th>\n",
+       "      <th>RCAP_90</th>\n",
+       "      <th>RCAP_00</th>\n",
+       "      <th>RCAP_10</th>\n",
+       "      <th>hud_recap_priority_community</th>\n",
+       "      <th>SHAPE_Length</th>\n",
+       "      <th>SHAPE_Area</th>\n",
+       "    </tr>\n",
+       "  </thead>\n",
+       "  <tbody>\n",
+       "    <tr>\n",
+       "      <th>0</th>\n",
+       "      <td>010010201001</td>\n",
+       "      <td>0.15</td>\n",
+       "      <td>692</td>\n",
+       "      <td>49.38</td>\n",
+       "      <td>0.79</td>\n",
+       "      <td>0.28</td>\n",
+       "      <td>10.00</td>\n",
+       "      <td>40.12</td>\n",
+       "      <td>91.02</td>\n",
+       "      <td>0.09</td>\n",
+       "      <td>...</td>\n",
+       "      <td>1.00</td>\n",
+       "      <td>Autauga</td>\n",
+       "      <td>1001.00</td>\n",
+       "      <td>20100.00</td>\n",
+       "      <td>0.00</td>\n",
+       "      <td>0.00</td>\n",
+       "      <td>0.00</td>\n",
+       "      <td>False</td>\n",
+       "      <td>0.15</td>\n",
+       "      <td>0.00</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>1</th>\n",
+       "      <td>010010201002</td>\n",
+       "      <td>0.15</td>\n",
+       "      <td>1153</td>\n",
+       "      <td>49.38</td>\n",
+       "      <td>0.79</td>\n",
+       "      <td>0.28</td>\n",
+       "      <td>10.00</td>\n",
+       "      <td>40.12</td>\n",
+       "      <td>2.62</td>\n",
+       "      <td>0.07</td>\n",
+       "      <td>...</td>\n",
+       "      <td>1.00</td>\n",
+       "      <td>Autauga</td>\n",
+       "      <td>1001.00</td>\n",
+       "      <td>20100.00</td>\n",
+       "      <td>0.00</td>\n",
+       "      <td>0.00</td>\n",
+       "      <td>0.00</td>\n",
+       "      <td>False</td>\n",
+       "      <td>0.15</td>\n",
+       "      <td>0.00</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>2</th>\n",
+       "      <td>010010202001</td>\n",
+       "      <td>0.25</td>\n",
+       "      <td>1020</td>\n",
+       "      <td>50.32</td>\n",
+       "      <td>0.81</td>\n",
+       "      <td>0.30</td>\n",
+       "      <td>10.07</td>\n",
+       "      <td>40.22</td>\n",
+       "      <td>4.68</td>\n",
+       "      <td>0.08</td>\n",
+       "      <td>...</td>\n",
+       "      <td>1.00</td>\n",
+       "      <td>Autauga</td>\n",
+       "      <td>1001.00</td>\n",
+       "      <td>20200.00</td>\n",
+       "      <td>0.00</td>\n",
+       "      <td>0.00</td>\n",
+       "      <td>0.00</td>\n",
+       "      <td>False</td>\n",
+       "      <td>0.09</td>\n",
+       "      <td>0.00</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>3</th>\n",
+       "      <td>010010202002</td>\n",
+       "      <td>0.25</td>\n",
+       "      <td>1152</td>\n",
+       "      <td>50.32</td>\n",
+       "      <td>0.81</td>\n",
+       "      <td>0.30</td>\n",
+       "      <td>10.07</td>\n",
+       "      <td>40.22</td>\n",
+       "      <td>218.65</td>\n",
+       "      <td>0.09</td>\n",
+       "      <td>...</td>\n",
+       "      <td>1.00</td>\n",
+       "      <td>Autauga</td>\n",
+       "      <td>1001.00</td>\n",
+       "      <td>20200.00</td>\n",
+       "      <td>0.00</td>\n",
+       "      <td>0.00</td>\n",
+       "      <td>0.00</td>\n",
+       "      <td>False</td>\n",
+       "      <td>0.09</td>\n",
+       "      <td>0.00</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>4</th>\n",
+       "      <td>010010203001</td>\n",
+       "      <td>0.21</td>\n",
+       "      <td>2555</td>\n",
+       "      <td>50.77</td>\n",
+       "      <td>0.82</td>\n",
+       "      <td>0.36</td>\n",
+       "      <td>10.12</td>\n",
+       "      <td>40.31</td>\n",
+       "      <td>69.64</td>\n",
+       "      <td>0.08</td>\n",
+       "      <td>...</td>\n",
+       "      <td>1.00</td>\n",
+       "      <td>Autauga</td>\n",
+       "      <td>1001.00</td>\n",
+       "      <td>20300.00</td>\n",
+       "      <td>0.00</td>\n",
+       "      <td>0.00</td>\n",
+       "      <td>0.00</td>\n",
+       "      <td>False</td>\n",
+       "      <td>0.10</td>\n",
+       "      <td>0.00</td>\n",
+       "    </tr>\n",
+       "  </tbody>\n",
+       "</table>\n",
+       "<p>5 rows × 165 columns</p>\n",
+       "</div>"
+      ],
+      "text/plain": [
+       "        GEOID10  Housing burden (percent)  Total population  \\\n",
+       "0  010010201001                      0.15               692   \n",
+       "1  010010201002                      0.15              1153   \n",
+       "2  010010202001                      0.25              1020   \n",
+       "3  010010202002                      0.25              1152   \n",
+       "4  010010203001                      0.21              2555   \n",
+       "\n",
+       "   Air toxics cancer risk  Respiratory hazard index  \\\n",
+       "0                   49.38                      0.79   \n",
+       "1                   49.38                      0.79   \n",
+       "2                   50.32                      0.81   \n",
+       "3                   50.32                      0.81   \n",
+       "4                   50.77                      0.82   \n",
+       "\n",
+       "   Diesel particulate matter  Particulate matter (PM2.5)  Ozone_x  \\\n",
+       "0                       0.28                       10.00    40.12   \n",
+       "1                       0.28                       10.00    40.12   \n",
+       "2                       0.30                       10.07    40.22   \n",
+       "3                       0.30                       10.07    40.22   \n",
+       "4                       0.36                       10.12    40.31   \n",
+       "\n",
+       "   Traffic proximity and volume  Proximity to RMP sites  ...  COUNTY  \\\n",
+       "0                         91.02                    0.09  ...    1.00   \n",
+       "1                          2.62                    0.07  ...    1.00   \n",
+       "2                          4.68                    0.08  ...    1.00   \n",
+       "3                        218.65                    0.09  ...    1.00   \n",
+       "4                         69.64                    0.08  ...    1.00   \n",
+       "\n",
+       "   COUNTY_NAME  CNTY_FIPS    TRACT  RCAP_90  RCAP_00  RCAP_10  \\\n",
+       "0      Autauga    1001.00 20100.00     0.00     0.00     0.00   \n",
+       "1      Autauga    1001.00 20100.00     0.00     0.00     0.00   \n",
+       "2      Autauga    1001.00 20200.00     0.00     0.00     0.00   \n",
+       "3      Autauga    1001.00 20200.00     0.00     0.00     0.00   \n",
+       "4      Autauga    1001.00 20300.00     0.00     0.00     0.00   \n",
+       "\n",
+       "   hud_recap_priority_community  SHAPE_Length  SHAPE_Area  \n",
+       "0                         False          0.15        0.00  \n",
+       "1                         False          0.15        0.00  \n",
+       "2                         False          0.09        0.00  \n",
+       "3                         False          0.09        0.00  \n",
+       "4                         False          0.10        0.00  \n",
+       "\n",
+       "[5 rows x 165 columns]"
+      ]
+     },
+     "execution_count": 7,
+     "metadata": {},
+     "output_type": "execute_result"
+    }
+   ],
    "source": [
     "# Join tract indices and CEJST data.\n",
     "# Note: we're joining on the census *tract*, so there will be multiple CBG entries joined to the same census tract row from CES,\n",
@@ -302,7 +1359,383 @@
    "metadata": {
     "scrolled": true
    },
-   "outputs": [],
+   "outputs": [
+    {
+     "name": "stdout",
+     "output_type": "stream",
+     "text": [
+      "Converting calenviroscreen_priority_community to boolean.\n",
+      "Converting hud_recap_priority_community to boolean.\n"
+     ]
+    },
+    {
+     "data": {
+      "application/vnd.jupyter.widget-view+json": {
+       "model_id": "223dcb75c0384fd5b93bc2ac3bc07656",
+       "version_major": 2,
+       "version_minor": 0
+      },
+      "text/plain": [
+       "  0%|          | 0/52 [00:00<?, ?it/s]"
+      ]
+     },
+     "metadata": {},
+     "output_type": "display_data"
+    },
+    {
+     "data": {
+      "text/html": [
+       "<div>\n",
+       "<style scoped>\n",
+       "    .dataframe tbody tr th:only-of-type {\n",
+       "        vertical-align: middle;\n",
+       "    }\n",
+       "\n",
+       "    .dataframe tbody tr th {\n",
+       "        vertical-align: top;\n",
+       "    }\n",
+       "\n",
+       "    .dataframe thead th {\n",
+       "        text-align: right;\n",
+       "    }\n",
+       "</style>\n",
+       "<table border=\"1\" class=\"dataframe\">\n",
+       "  <thead>\n",
+       "    <tr style=\"text-align: right;\">\n",
+       "      <th></th>\n",
+       "      <th></th>\n",
+       "      <th>GEOID10_STATE</th>\n",
+       "      <th>State name</th>\n",
+       "      <th>Total CBGs in state</th>\n",
+       "      <th>Total population in state</th>\n",
+       "      <th>Score A (top 25th percentile) (priority population)</th>\n",
+       "      <th>Score A (top 25th percentile) (total CBGs)</th>\n",
+       "      <th>Score A (top 25th percentile) (percent CBGs)</th>\n",
+       "      <th>Score A (top 25th percentile) (percent population)</th>\n",
+       "      <th>Score B (top 25th percentile) (priority population)</th>\n",
+       "      <th>Score B (top 25th percentile) (total CBGs)</th>\n",
+       "      <th>...</th>\n",
+       "      <th>Score E (top 25th percentile) (percent CBGs)</th>\n",
+       "      <th>Score E (top 25th percentile) (percent population)</th>\n",
+       "      <th>calenviroscreen_priority_community (priority population)</th>\n",
+       "      <th>calenviroscreen_priority_community (total CBGs)</th>\n",
+       "      <th>calenviroscreen_priority_community (percent CBGs)</th>\n",
+       "      <th>calenviroscreen_priority_community (percent population)</th>\n",
+       "      <th>hud_recap_priority_community (priority population)</th>\n",
+       "      <th>hud_recap_priority_community (total CBGs)</th>\n",
+       "      <th>hud_recap_priority_community (percent CBGs)</th>\n",
+       "      <th>hud_recap_priority_community (percent population)</th>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>GEOID10_STATE</th>\n",
+       "      <th></th>\n",
+       "      <th></th>\n",
+       "      <th></th>\n",
+       "      <th></th>\n",
+       "      <th></th>\n",
+       "      <th></th>\n",
+       "      <th></th>\n",
+       "      <th></th>\n",
+       "      <th></th>\n",
+       "      <th></th>\n",
+       "      <th></th>\n",
+       "      <th></th>\n",
+       "      <th></th>\n",
+       "      <th></th>\n",
+       "      <th></th>\n",
+       "      <th></th>\n",
+       "      <th></th>\n",
+       "      <th></th>\n",
+       "      <th></th>\n",
+       "      <th></th>\n",
+       "      <th></th>\n",
+       "      <th></th>\n",
+       "    </tr>\n",
+       "  </thead>\n",
+       "  <tbody>\n",
+       "    <tr>\n",
+       "      <th>01</th>\n",
+       "      <th>0</th>\n",
+       "      <td>01</td>\n",
+       "      <td>Alabama</td>\n",
+       "      <td>3438</td>\n",
+       "      <td>4850771</td>\n",
+       "      <td>1547345</td>\n",
+       "      <td>1326</td>\n",
+       "      <td>0.39</td>\n",
+       "      <td>0.32</td>\n",
+       "      <td>1556417</td>\n",
+       "      <td>1323</td>\n",
+       "      <td>...</td>\n",
+       "      <td>0.23</td>\n",
+       "      <td>0.19</td>\n",
+       "      <td>0</td>\n",
+       "      <td>0</td>\n",
+       "      <td>0.00</td>\n",
+       "      <td>0.00</td>\n",
+       "      <td>235117</td>\n",
+       "      <td>258</td>\n",
+       "      <td>0.08</td>\n",
+       "      <td>0.05</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>02</th>\n",
+       "      <th>0</th>\n",
+       "      <td>02</td>\n",
+       "      <td>Alaska</td>\n",
+       "      <td>534</td>\n",
+       "      <td>738565</td>\n",
+       "      <td>63868</td>\n",
+       "      <td>57</td>\n",
+       "      <td>0.11</td>\n",
+       "      <td>0.09</td>\n",
+       "      <td>63868</td>\n",
+       "      <td>57</td>\n",
+       "      <td>...</td>\n",
+       "      <td>0.14</td>\n",
+       "      <td>0.12</td>\n",
+       "      <td>0</td>\n",
+       "      <td>0</td>\n",
+       "      <td>0.00</td>\n",
+       "      <td>0.00</td>\n",
+       "      <td>6536</td>\n",
+       "      <td>8</td>\n",
+       "      <td>0.01</td>\n",
+       "      <td>0.01</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>04</th>\n",
+       "      <th>0</th>\n",
+       "      <td>04</td>\n",
+       "      <td>Arizona</td>\n",
+       "      <td>4178</td>\n",
+       "      <td>6809946</td>\n",
+       "      <td>1956052</td>\n",
+       "      <td>1230</td>\n",
+       "      <td>0.29</td>\n",
+       "      <td>0.29</td>\n",
+       "      <td>1960856</td>\n",
+       "      <td>1231</td>\n",
+       "      <td>...</td>\n",
+       "      <td>0.30</td>\n",
+       "      <td>0.30</td>\n",
+       "      <td>0</td>\n",
+       "      <td>0</td>\n",
+       "      <td>0.00</td>\n",
+       "      <td>0.00</td>\n",
+       "      <td>560353</td>\n",
+       "      <td>378</td>\n",
+       "      <td>0.09</td>\n",
+       "      <td>0.08</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>05</th>\n",
+       "      <th>0</th>\n",
+       "      <td>05</td>\n",
+       "      <td>Arkansas</td>\n",
+       "      <td>2147</td>\n",
+       "      <td>2977944</td>\n",
+       "      <td>960799</td>\n",
+       "      <td>817</td>\n",
+       "      <td>0.38</td>\n",
+       "      <td>0.32</td>\n",
+       "      <td>975780</td>\n",
+       "      <td>826</td>\n",
+       "      <td>...</td>\n",
+       "      <td>0.20</td>\n",
+       "      <td>0.18</td>\n",
+       "      <td>0</td>\n",
+       "      <td>0</td>\n",
+       "      <td>0.00</td>\n",
+       "      <td>0.00</td>\n",
+       "      <td>101200</td>\n",
+       "      <td>106</td>\n",
+       "      <td>0.05</td>\n",
+       "      <td>0.03</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>06</th>\n",
+       "      <th>0</th>\n",
+       "      <td>06</td>\n",
+       "      <td>California</td>\n",
+       "      <td>23212</td>\n",
+       "      <td>38982847</td>\n",
+       "      <td>12610810</td>\n",
+       "      <td>7102</td>\n",
+       "      <td>0.31</td>\n",
+       "      <td>0.32</td>\n",
+       "      <td>12556846</td>\n",
+       "      <td>7065</td>\n",
+       "      <td>...</td>\n",
+       "      <td>0.40</td>\n",
+       "      <td>0.42</td>\n",
+       "      <td>9610287</td>\n",
+       "      <td>5690</td>\n",
+       "      <td>0.25</td>\n",
+       "      <td>0.25</td>\n",
+       "      <td>1748765</td>\n",
+       "      <td>1013</td>\n",
+       "      <td>0.04</td>\n",
+       "      <td>0.04</td>\n",
+       "    </tr>\n",
+       "  </tbody>\n",
+       "</table>\n",
+       "<p>5 rows × 32 columns</p>\n",
+       "</div>"
+      ],
+      "text/plain": [
+       "                GEOID10_STATE  State name  Total CBGs in state  \\\n",
+       "GEOID10_STATE                                                    \n",
+       "01            0            01     Alabama                 3438   \n",
+       "02            0            02      Alaska                  534   \n",
+       "04            0            04     Arizona                 4178   \n",
+       "05            0            05    Arkansas                 2147   \n",
+       "06            0            06  California                23212   \n",
+       "\n",
+       "                 Total population in state  \\\n",
+       "GEOID10_STATE                                \n",
+       "01            0                    4850771   \n",
+       "02            0                     738565   \n",
+       "04            0                    6809946   \n",
+       "05            0                    2977944   \n",
+       "06            0                   38982847   \n",
+       "\n",
+       "                 Score A (top 25th percentile) (priority population)  \\\n",
+       "GEOID10_STATE                                                          \n",
+       "01            0                                            1547345     \n",
+       "02            0                                              63868     \n",
+       "04            0                                            1956052     \n",
+       "05            0                                             960799     \n",
+       "06            0                                           12610810     \n",
+       "\n",
+       "                 Score A (top 25th percentile) (total CBGs)  \\\n",
+       "GEOID10_STATE                                                 \n",
+       "01            0                                        1326   \n",
+       "02            0                                          57   \n",
+       "04            0                                        1230   \n",
+       "05            0                                         817   \n",
+       "06            0                                        7102   \n",
+       "\n",
+       "                 Score A (top 25th percentile) (percent CBGs)  \\\n",
+       "GEOID10_STATE                                                   \n",
+       "01            0                                          0.39   \n",
+       "02            0                                          0.11   \n",
+       "04            0                                          0.29   \n",
+       "05            0                                          0.38   \n",
+       "06            0                                          0.31   \n",
+       "\n",
+       "                 Score A (top 25th percentile) (percent population)  \\\n",
+       "GEOID10_STATE                                                         \n",
+       "01            0                                               0.32    \n",
+       "02            0                                               0.09    \n",
+       "04            0                                               0.29    \n",
+       "05            0                                               0.32    \n",
+       "06            0                                               0.32    \n",
+       "\n",
+       "                 Score B (top 25th percentile) (priority population)  \\\n",
+       "GEOID10_STATE                                                          \n",
+       "01            0                                            1556417     \n",
+       "02            0                                              63868     \n",
+       "04            0                                            1960856     \n",
+       "05            0                                             975780     \n",
+       "06            0                                           12556846     \n",
+       "\n",
+       "                 Score B (top 25th percentile) (total CBGs)  ...  \\\n",
+       "GEOID10_STATE                                                ...   \n",
+       "01            0                                        1323  ...   \n",
+       "02            0                                          57  ...   \n",
+       "04            0                                        1231  ...   \n",
+       "05            0                                         826  ...   \n",
+       "06            0                                        7065  ...   \n",
+       "\n",
+       "                 Score E (top 25th percentile) (percent CBGs)  \\\n",
+       "GEOID10_STATE                                                   \n",
+       "01            0                                          0.23   \n",
+       "02            0                                          0.14   \n",
+       "04            0                                          0.30   \n",
+       "05            0                                          0.20   \n",
+       "06            0                                          0.40   \n",
+       "\n",
+       "                 Score E (top 25th percentile) (percent population)  \\\n",
+       "GEOID10_STATE                                                         \n",
+       "01            0                                               0.19    \n",
+       "02            0                                               0.12    \n",
+       "04            0                                               0.30    \n",
+       "05            0                                               0.18    \n",
+       "06            0                                               0.42    \n",
+       "\n",
+       "                 calenviroscreen_priority_community (priority population)  \\\n",
+       "GEOID10_STATE                                                               \n",
+       "01            0                                                  0          \n",
+       "02            0                                                  0          \n",
+       "04            0                                                  0          \n",
+       "05            0                                                  0          \n",
+       "06            0                                            9610287          \n",
+       "\n",
+       "                 calenviroscreen_priority_community (total CBGs)  \\\n",
+       "GEOID10_STATE                                                      \n",
+       "01            0                                                0   \n",
+       "02            0                                                0   \n",
+       "04            0                                                0   \n",
+       "05            0                                                0   \n",
+       "06            0                                             5690   \n",
+       "\n",
+       "                 calenviroscreen_priority_community (percent CBGs)  \\\n",
+       "GEOID10_STATE                                                        \n",
+       "01            0                                               0.00   \n",
+       "02            0                                               0.00   \n",
+       "04            0                                               0.00   \n",
+       "05            0                                               0.00   \n",
+       "06            0                                               0.25   \n",
+       "\n",
+       "                 calenviroscreen_priority_community (percent population)  \\\n",
+       "GEOID10_STATE                                                              \n",
+       "01            0                                               0.00         \n",
+       "02            0                                               0.00         \n",
+       "04            0                                               0.00         \n",
+       "05            0                                               0.00         \n",
+       "06            0                                               0.25         \n",
+       "\n",
+       "                 hud_recap_priority_community (priority population)  \\\n",
+       "GEOID10_STATE                                                         \n",
+       "01            0                                             235117    \n",
+       "02            0                                               6536    \n",
+       "04            0                                             560353    \n",
+       "05            0                                             101200    \n",
+       "06            0                                            1748765    \n",
+       "\n",
+       "                 hud_recap_priority_community (total CBGs)  \\\n",
+       "GEOID10_STATE                                                \n",
+       "01            0                                        258   \n",
+       "02            0                                          8   \n",
+       "04            0                                        378   \n",
+       "05            0                                        106   \n",
+       "06            0                                       1013   \n",
+       "\n",
+       "                 hud_recap_priority_community (percent CBGs)  \\\n",
+       "GEOID10_STATE                                                  \n",
+       "01            0                                         0.08   \n",
+       "02            0                                         0.01   \n",
+       "04            0                                         0.09   \n",
+       "05            0                                         0.05   \n",
+       "06            0                                         0.04   \n",
+       "\n",
+       "                 hud_recap_priority_community (percent population)  \n",
+       "GEOID10_STATE                                                       \n",
+       "01            0                                               0.05  \n",
+       "02            0                                               0.01  \n",
+       "04            0                                               0.08  \n",
+       "05            0                                               0.03  \n",
+       "06            0                                               0.04  \n",
+       "\n",
+       "[5 rows x 32 columns]"
+      ]
+     },
+     "execution_count": 10,
+     "metadata": {},
+     "output_type": "execute_result"
+    }
+   ],
    "source": [
     "def get_state_distributions(\n",
     "    df: pd.DataFrame, priority_communities_fields: typing.List[str]\n",
@@ -1132,7 +2565,50 @@
    "metadata": {
     "scrolled": true
    },
-   "outputs": [],
+   "outputs": [
+    {
+     "name": "stdout",
+     "output_type": "stream",
+     "text": [
+      "Running comparisons for Score A against CalEnviroScreen 4.0...\n"
+     ]
+    },
+    {
+     "data": {
+      "application/vnd.jupyter.widget-view+json": {
+       "model_id": "531ec4deb2f54c26ad0f5311fdea0e60",
+       "version_major": 2,
+       "version_minor": 0
+      },
+      "text/plain": [
+       "  0%|          | 0/8057 [00:00<?, ?it/s]"
+      ]
+     },
+     "metadata": {},
+     "output_type": "display_data"
+    },
+    {
+     "name": "stdout",
+     "output_type": "stream",
+     "text": [
+      "Running comparisons for Score A against HUD RECAP...\n"
+     ]
+    },
+    {
+     "data": {
+      "application/vnd.jupyter.widget-view+json": {
+       "model_id": "170da68ae0734892bef4a452b5de45f7",
+       "version_major": 2,
+       "version_minor": 0
+      },
+      "text/plain": [
+       "  0%|          | 0/74001 [00:00<?, ?it/s]"
+      ]
+     },
+     "metadata": {},
+     "output_type": "display_data"
+    }
+   ],
    "source": [
     "# Actually execute the functions\n",
     "file_paths = execute_comparisons(\n",
@@ -1143,6 +2619,14 @@
     "\n",
     "print(file_paths)"
    ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "887ee948",
+   "metadata": {},
+   "outputs": [],
+   "source": []
   }
  ],
  "metadata": {
@@ -1161,7 +2645,7 @@
    "name": "python",
    "nbconvert_exporter": "python",
    "pygments_lexer": "ipython3",
-   "version": "3.7.1"
+   "version": "3.9.5"
   }
  },
  "nbformat": 4,