pride-data-analysis/analysis/analysis1.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Nat Scott"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Research question/interests\n",
    "\n",
    "**Is there a correlation between political alignment & living in neighbourhoods with large quantities of LGBT people?** The obvious answer to this question is \"yes, they are going to mostly be democrats\" but anyone who's ever been around queer people will know that this question is quite a bit more nuanced than that, and this nuance is what we hope to capture in investigating this question.\n",
    "\n",
    "- The gaybourhoods data set does not include data on residents political alignments, however, there is a wealth of electoral data available freely online that we intend on incorporating into this project. The primary difficulty then will be developing a geographic \"compatibility layer\" between the data sets so that the data can be understood in the same context. To build this, we intend on working with the OpenStreetMap API to create an additional column representing observations position space in a more neutral way, such as their coordinates.\n",
    "- Alternatively, we've also considered working with an additional data set that links US zip codes to their longitude and lattitude positions. As such, incorporating this data would be as easy as merging the two tables.\n",
    "\n",
    "\n",
    "**Is there a correlation between geographical stratums & being LGBT?** This question is more abstract, and will serve as a preliminary exploration of the data in hopes of establishing two key details along the way that will shape the rest of the project: how do we quantify queerness, and how do we best represent it visually?\n",
    "\n",
    "- Once again, representing this data visually will require determining the coordinates associated with each observation.\n",
    "- The gaybourhoods data set defines a \"gaybourhood index\" which effectively measures how friendly a given neighbourhood is to queer people. Since this index is entirely subjective, we will need to closely evaluate it's usefulness for our project and investigate different ways to quantify \"queer-friendliness\"\n",
    "- In addition to the last point, since, of course, no matter what choice of observations we make, the measurement will still be subjective, answering this research question will come more so in the form of comparing and contrasting different measurements to see what they tell us.\n",
    "- Obviously, visualizing this among many aspects of the other research questions would involve projecting the data onto a map of the United States, so visualizing this research question would motivate many of the visualizations for other components of this project"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Analysis Pipeline"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import pandas as pd\n",
    "import seaborn as sns\n",
    "import numpy as np\n",
    "from PIL import Image\n",
    "\n",
    "import matplotlib.pyplot as plt\n",
    "import matplotlib.patches as mpatches\n",
    "from matplotlib.colors import LinearSegmentedColormap\n",
    "\n",
    "sns.set_theme(style=\"white\")\n",
    "sns.set(font_scale=1.2)\n",
    "sns.set_palette(\"copper_r\")\n",
    "plt.rcParams[\"axes.labelsize\"] = 12\n",
    "\n",
    "# Absolutely diabolical method of doing relative imports with a package who shares its name with\n",
    "# something in the stdlib in Jupyter Lab because it seems impossible otherwise\n",
    "__import__(\"sys\").path.append(\"./code\")\n",
    "from project_functions1 import *\n",
    "\n",
    "CITIES = [\n",
    "    { \"name\": \"Atlanta\", \"pos\": (-84.791944, -84.001944, 33.478333, 34.0375) },\n",
    "    { \"name\": \"Austin\", \"pos\": (-98.071667, -97.371111, 30.033889, 30.569722) },\n",
    "    { \"name\": \"Boston\", \"pos\": (-71.284444, -70.880278, 42.206389, 42.484444) },\n",
    "    { \"name\": \"Chicago\", \"pos\": (-88.044167, -87.508333, 41.624444, 42.026389) },\n",
    "    { \"name\": \"Denver\", \"pos\": (-105.104444, -104.625556, 39.578611, 39.920278) },\n",
    "    { \"name\": \"Houston\", \"pos\": (-95.858333, -95.013056, 29.472778, 30.134167) },\n",
    "    { \"name\": \"Los Angeles\", \"pos\": (-118.615556, -117.618333, 33.701111, 34.411667) },\n",
    "    { \"name\": \"Miami\", \"pos\": (-80.2375, -80.145833, 25.734722, 25.812222) },\n",
    "    { \"name\": \"New York\", \"pos\": (-74.459722, -73.393333, 40.355556, 41.102222) },\n",
    "    { \"name\": \"New Orleans\", \"pos\": (-90.220833, -89.915833, 29.813056, 30.045556) },\n",
    "    { \"name\": \"Philadelphia\", \"pos\": (-75.344167, -74.94, 39.846667, 40.146389) },\n",
    "    { \"name\": \"Portland\", \"pos\": (-122.795278, -122.493333, 45.465556, 45.654444) },\n",
    "    { \"name\": \"San Francisco\", \"pos\": (-122.521667, -122.352778, 37.690278, 37.812222) },\n",
    "    { \"name\": \"Seattle\", \"pos\": (-122.473056, -122.151667, 47.493333, 47.718611) },\n",
    "    { \"name\": \"Washington DC\", \"pos\": (-77.1125, -76.919722, 38.833333, 38.963889) },\n",
    "]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Data Analysis Pipeline"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Now in one, new-and-improved, non-descript method imported from another file\n",
    "gb, pol, counties, cords = load_and_process()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "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>SS_Index</th>\n",
       "      <th>TOTINDEX</th>\n",
       "      <th>lat</th>\n",
       "      <th>long</th>\n",
       "      <th>kinsey</th>\n",
       "      <th>percent_democrat</th>\n",
       "      <th>neighbourhood_kinsey</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>55.012669</td>\n",
       "      <td>67.077054</td>\n",
       "      <td>34.093828</td>\n",
       "      <td>-118.381697</td>\n",
       "      <td>6</td>\n",
       "      <td>0.456450</td>\n",
       "      <td>1.132075</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>58.412517</td>\n",
       "      <td>61.866815</td>\n",
       "      <td>37.758057</td>\n",
       "      <td>-122.435410</td>\n",
       "      <td>6</td>\n",
       "      <td>0.742633</td>\n",
       "      <td>2.533333</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>30.731192</td>\n",
       "      <td>37.908747</td>\n",
       "      <td>40.742039</td>\n",
       "      <td>-74.000620</td>\n",
       "      <td>6</td>\n",
       "      <td>0.775215</td>\n",
       "      <td>1.091667</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>21.995808</td>\n",
       "      <td>37.530067</td>\n",
       "      <td>40.734012</td>\n",
       "      <td>-74.006746</td>\n",
       "      <td>6</td>\n",
       "      <td>0.794248</td>\n",
       "      <td>1.101695</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>21.284224</td>\n",
       "      <td>35.843573</td>\n",
       "      <td>37.773134</td>\n",
       "      <td>-122.411167</td>\n",
       "      <td>5</td>\n",
       "      <td>0.742633</td>\n",
       "      <td>2.533333</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2323</th>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>45.528666</td>\n",
       "      <td>-122.678981</td>\n",
       "      <td>0</td>\n",
       "      <td>0.753689</td>\n",
       "      <td>1.708333</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2324</th>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>47.606211</td>\n",
       "      <td>-122.333792</td>\n",
       "      <td>0</td>\n",
       "      <td>0.545340</td>\n",
       "      <td>2.148148</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2325</th>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>47.449678</td>\n",
       "      <td>-122.307657</td>\n",
       "      <td>0</td>\n",
       "      <td>0.543670</td>\n",
       "      <td>1.357143</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2326</th>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>47.604569</td>\n",
       "      <td>-122.335359</td>\n",
       "      <td>0</td>\n",
       "      <td>0.545340</td>\n",
       "      <td>2.148148</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2327</th>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>47.649339</td>\n",
       "      <td>-122.310294</td>\n",
       "      <td>0</td>\n",
       "      <td>0.545340</td>\n",
       "      <td>1.714286</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>2328 rows × 7 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "       SS_Index   TOTINDEX        lat        long  kinsey  percent_democrat  \\\n",
       "0     55.012669  67.077054  34.093828 -118.381697       6          0.456450   \n",
       "1     58.412517  61.866815  37.758057 -122.435410       6          0.742633   \n",
       "2     30.731192  37.908747  40.742039  -74.000620       6          0.775215   \n",
       "3     21.995808  37.530067  40.734012  -74.006746       6          0.794248   \n",
       "4     21.284224  35.843573  37.773134 -122.411167       5          0.742633   \n",
       "...         ...        ...        ...         ...     ...               ...   \n",
       "2323   0.000000   0.000000  45.528666 -122.678981       0          0.753689   \n",
       "2324   0.000000   0.000000  47.606211 -122.333792       0          0.545340   \n",
       "2325   0.000000   0.000000  47.449678 -122.307657       0          0.543670   \n",
       "2326   0.000000   0.000000  47.604569 -122.335359       0          0.545340   \n",
       "2327   0.000000   0.000000  47.649339 -122.310294       0          0.545340   \n",
       "\n",
       "      neighbourhood_kinsey  \n",
       "0                 1.132075  \n",
       "1                 2.533333  \n",
       "2                 1.091667  \n",
       "3                 1.101695  \n",
       "4                 2.533333  \n",
       "...                    ...  \n",
       "2323              1.708333  \n",
       "2324              2.148148  \n",
       "2325              1.357143  \n",
       "2326              2.148148  \n",
       "2327              1.714286  \n",
       "\n",
       "[2328 rows x 7 columns]"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "gb.to_csv(\"../data/processed/gaybourhoods-nat.csv\")\n",
    "\n",
    "SS_INDEX_MAX = gb.SS_Index.max()\n",
    "\n",
    "gb"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "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>county</th>\n",
       "      <th>lat</th>\n",
       "      <th>long</th>\n",
       "      <th>percent</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>Autauga AL</td>\n",
       "      <td>32.532237</td>\n",
       "      <td>-86.646439</td>\n",
       "      <td>0.265878</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>Baldwin AL</td>\n",
       "      <td>30.659218</td>\n",
       "      <td>-87.746067</td>\n",
       "      <td>0.215894</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>Barbour AL</td>\n",
       "      <td>31.870253</td>\n",
       "      <td>-85.405103</td>\n",
       "      <td>0.513685</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>Bibb AL</td>\n",
       "      <td>33.015893</td>\n",
       "      <td>-87.127148</td>\n",
       "      <td>0.261520</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>Blount AL</td>\n",
       "      <td>33.977357</td>\n",
       "      <td>-86.566440</td>\n",
       "      <td>0.123719</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2879</th>\n",
       "      <td>Sweetwater WY</td>\n",
       "      <td>41.660328</td>\n",
       "      <td>-108.875677</td>\n",
       "      <td>0.282569</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2880</th>\n",
       "      <td>Teton WY</td>\n",
       "      <td>44.048662</td>\n",
       "      <td>-110.426087</td>\n",
       "      <td>0.541957</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2881</th>\n",
       "      <td>Uinta WY</td>\n",
       "      <td>41.284726</td>\n",
       "      <td>-110.558947</td>\n",
       "      <td>0.190655</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2882</th>\n",
       "      <td>Washakie WY</td>\n",
       "      <td>43.878830</td>\n",
       "      <td>-107.669052</td>\n",
       "      <td>0.201318</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2883</th>\n",
       "      <td>Weston WY</td>\n",
       "      <td>43.846213</td>\n",
       "      <td>-104.570020</td>\n",
       "      <td>0.125633</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>2884 rows × 4 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "             county        lat        long   percent\n",
       "0        Autauga AL  32.532237  -86.646439  0.265878\n",
       "1        Baldwin AL  30.659218  -87.746067  0.215894\n",
       "2        Barbour AL  31.870253  -85.405103  0.513685\n",
       "3           Bibb AL  33.015893  -87.127148  0.261520\n",
       "4         Blount AL  33.977357  -86.566440  0.123719\n",
       "...             ...        ...         ...       ...\n",
       "2879  Sweetwater WY  41.660328 -108.875677  0.282569\n",
       "2880       Teton WY  44.048662 -110.426087  0.541957\n",
       "2881       Uinta WY  41.284726 -110.558947  0.190655\n",
       "2882    Washakie WY  43.878830 -107.669052  0.201318\n",
       "2883      Weston WY  43.846213 -104.570020  0.125633\n",
       "\n",
       "[2884 rows x 4 columns]"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pol.to_csv(\"../data/processed/election-2012.csv\")\n",
    "pol"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Exploratory Data Analysis\n",
    "\n",
    "In the previous section we got a pretty good idea of what the data looks like and managed to condense it to fewer key variables that will be useful when answering the research questions established above. Now, we will compare these variables, attempt to create some plots, and see if we can't uncover any interesting relationships that aren't evident by looking at the numbers alone.\n",
    "\n",
    "Let's start with a scatterplot of the `gb` data set. I have no idea what this is going to look like. Let's see:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkMAAAHTCAYAAAAktxxFAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/P9b71AAAACXBIWXMAAA9hAAAPYQGoP6dpAAB21UlEQVR4nO3dd3wUdeL/8ddsSzaBJARIaAIKbKL0KqIgHRUQscBZ0FOxoqen3p16Rb3fWb6eynnq2RUQRUBAFBSRpmIBpIlCAKVIkVCSEEiy2Ta/P4YsLAklpJJ9Px+PPCAzszuf+WR39r2f+Xw+Y5imaSIiIiISpWxVXQARERGRqqQwJCIiIlFNYUhERESimsKQiIiIRDWFIREREYlqCkMiIiIS1RSGREREJKopDImIiEhUUxgSERGRqKYwdBobNWoUaWlpVV0MKYMHH3yQtLQ0tm/fXtVFOWkHDx7kX//6F3379uWcc84hLS2NdevWlek509LSGDVqVDmVsGY5Xd/np9vf9IUXXiAtLY0lS5ZU6n5Pt3qqqRxVXYCapuik1ahRI+bMmUNMTEyxbfr27cuOHTv46aefcDhq1p9g06ZNvPXWWyxZsoRdu3Zht9tJTk6mefPmdOrUid/97nfUq1evUsrywgsv8OKLLzJhwgTOPffcYusffPBBZsyYwfz582nSpEmllKmi9O3bF4AFCxZU+L6efvppJk+eTJ8+fbj00kux2+0n/JtWZvlOlWmazJ8/n1mzZvHDDz+QlZVFMBikTp06pKWl0bNnT4YOHUqdOnWquqgiEbZv306/fv0YPnw4Tz31VFUX55SNGjWKpUuXsn79+krfd836JK5Gdu7cyfjx47n11lsrbB//93//R0FBQYU9f2l9++233HbbbRQWFtKxY0d69uxJrVq12L17NytXruTrr7+mU6dOlRaGpGIsWrSI5s2b88orr1R1UcpNZmYmf/zjH1m+fDlut5tzzz2X5s2bExMTw549e1i5ciWPP/44Y8eO5fPPP9drWKSGURiqAImJiRiGwWuvvcaVV15JcnJyheynUaNGFfK8p+qRRx6hsLCQp556iuHDhxdbn5GRQWJiYhWUTMrT7t276dq1a1UXo9zk5eUxevRoNmzYwEUXXcSjjz5aYuvPqlWrePbZZ/F6vVVQShGpSApDFSA2NpabbrqJJ598kpdeeom///3vJ/W46dOns3DhQtauXcuePXtwOBx4PB6uvvpqhg0bVmz7o5sUZ8+ezX333ccNN9zAww8/XGx7n8/H+eefT0xMDIsWLYq4RDdr1iwmT57MunXrKCwspEmTJgwdOpTRo0fjcrlOWPZ9+/axdetWateuXWIQAkhPTy9x+a5du3jjjTf48ssv2bVrF7GxsTRt2pQ+ffowZsyY8Hbfffcds2fPZvny5ezatYtAIEDTpk256KKLuOWWWyIuSRZdigS4/vrrI/a3fv36iD4Y/fr1C/+/cePGEZdycnJyePPNN5k3bx47duzA6XTSpk0bbrnlFi644IKI550+fToPPfQQTz75JMnJybzyyitkZGTgdDo577zzuO+++2jevPkJavKwTz75hHfffZeMjAz8fj/NmjVjyJAh3HjjjeG/yZIlSyKO78jjOtkm8927d/Pyyy/zxRdfsHv3bmrVqkWXLl24/fbbadOmTXi7otcbwNKlS8P76tatG++8806Jz13a8mVlZTF27FgWLlxITk4OzZo146abbuKKK64o8fm/+uorJkyYwA8//EBeXh4NGjRgwIAB3HHHHSQkJJzw2AHefvttNmzYQJcuXRg7diw2W8ldKTt06MCECRMIhUIRy+fNm8ecOXNYs2YNmZmZAJx11llcdtllXHfddRHPd9999zF79mzeeecdunXrVmwfn332GX/4wx+49tpr+cc//hGxzufz8dJLL/Hxxx+ze/duGjRowGWXXcatt95a4nv022+/5Y033mDNmjXk5+fTuHFjBgwYwG233Ubt2rUjtj3eZcxjXW5OS0ujW7duPPfcc/znP//hyy+/ZO/evTz++ONcfvnl4e1K8zcNhUJMnjyZDz74gE2bNmGaJi1atOCKK67gd7/7XYl/m9IcJ8CPP/7I2LFjWbFiBYZh0K5dO+65555i2xX5/vvveeONN1i7di1ZWVkkJibSuHFjevXqxV133XXMxx3J5/Px+uuvM2PGDHbt2kVKSgpDhw6NOL8dKTMzk6lTp7J48WK2bdvG/v37SUpK4txzz+WOO+6gZcuW4W2L/j4AM2bMYMaMGeF1Tz75JJdffjk+n48pU6bwxRdf8PPPP7Nnzx7i4uI455xzuPHGG7nwwgtP6jjA6jM4fvx4Pv30U3bu3IlpmtStW5c2bdowevToiHMGwOrVq3nzzTdZvnw5+/fvp27dulx44YWMGTOG1NRU4PBlviJHnieOd34pTwpDFeTaa6/l3XffZfLkyYwaNeqkPgQfffRRWrZsSdeuXalfvz45OTl88cUX/PnPf2bz5s3ce++9x318//79qV27NrNmzeLPf/5zsf5I8+bNIzc3l5tuuili3UMPPcT06dNp0KABAwcOJCEhgVWrVvH888/z7bff8vbbb5+wb1Pt2rVxOBzk5+eze/duUlJSTni8AGvWrGH06NHk5OTQtWtXBgwYgNfr5eeff+bFF1+MOFm8/vrrbN68mY4dO3LhhRfi8/lYsWIFL7zwAkuWLGHcuHHY7XbACkDz589n6dKlDB8+nMaNG0fs96677mLevHlkZGRw/fXXhz84jzx57tixg1GjRrFjxw66dOlCz549KSgoYOHChYwePZp//vOfjBgxotgxzZ07l6+++or+/fvTrVs31q1bx2effcaSJUuYNGkSZ5111gnr5bnnnuPVV1+lTp06DBkyhLi4OL766iuee+45Fi9ezJtvvonL5aJx48bcddddjB8/HoAbbrgh/Bxnn332Cfezbds2rrnmGnbv3k337t0ZPHgwv/32G3PmzGHRokW88MIL9OnTB7DCS7du3XjxxRdp3LhxOPQeXbdHKk35cnNzufrqq3G5XAwaNAifz8ecOXN4+OGHsdlsxUL2iy++yAsvvEBSUhK9e/cmOTmZDRs28NZbb/Hll18yefJkatWqdcI6+OCDDwC48847jxmEihiGEX6NFXnmmWew2Wy0a9eO1NRUDhw4wHfffcfjjz/OmjVr+Pe//x3e9uqrr2b27NlMnjy5xDA0efJkAH73u98VW3fPPfewZs0aLrroIhwOB/Pnz+eFF17gxx9/5OWXX8YwjPC277//Po8++ihut5uLLrqIunXrsnTpUl5//XUWLlzIpEmTTjosHk9OTg4jR44kLi6OgQMHYhgGdevWDa8v7d/0T3/6E7NmzaJhw4ZceeWVGIbBvHnzeOyxx1i+fDnPPvtsxPalPc4VK1Zw44034vf7GTBgAM2aNWPdunWMGjWK7t27Fzu+L7/8kttuu41atWrRt29fUlNTycnJYdOmTbz33nsnFYZM0+Tee+9l/vz5NG3alOuuuw6/38+0adPYsGFDiY/5/vvvef311zn33HMZOHAgcXFxbN26lc8++4wFCxYwadKk8JfLbt26cf311zNhwgTS09Pp379/+HmK3mP79+/n8ccfp2PHjvTo0YPk5GT27NnDwoULufXWW/nXv/7FVVdddVLHMnr0aFauXEnHjh256qqrsNvtZGZmsmTJErp06RIRhj744AP+8Y9/4HK56Nu3Lw0aNGDr1q1MnTqVBQsWMGXKFBo1akRCQgJ33XUXM2bMYMeOHRH1erzzS7kypVx5PB6zZ8+epmma5qeffmp6PB5zzJgxEdv06dPH9Hg8pt/vj1i+devWYs9XWFhoXn/99eY555xj7tq1K2LdddddZ3o8nohlf//7302Px2MuWLCg2HPdcsstpsfjMTMyMsLLpk2bFi5jQUFBxPb//e9/TY/HY44bN+4kjtw07777btPj8Zj9+vUz33jjDXPVqlVmfn7+MbcvLCwM18VHH31UbP1vv/0W8fuvv/5qhkKhYtuNHTvW9Hg85uzZs0ss/3fffVfi/v/yl7+YHo/H3LZtW4n
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plot_naive_scatter1 = sns.scatterplot(data=gb, x=\"long\", y=\"lat\")\n",
    "_ = plot_naive_scatter1.set(\n",
    "    xlabel=\"Longitude\",\n",
    "    ylabel=\"Latitude\",\n",
    "    title=\"Naive Scatteplot of the Gaybourhoods data set\"\n",
    ")\n",
    "\n",
    "plot_naive_scatter1.figure.savefig(\"./images/graphs/1-naive-scatter1.png\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Here the data does appear to trace the outline of the United States, probably because most of the cities the data set covers are along the coast. We can do something similar with the election data we cleaned earlier:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAArEAAAHTCAYAAAAwDkDhAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/P9b71AAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOxdd3xUVdp+7tyZSSYJkAIJSguQTIYSEjpSQzEWmhRFURRcFyzo6vqtrt+u277dddfd1bWX3RUsq4ICimDBAtIElRJqEpQeJZEUAkkmM3PnfH+8c2dun0kIEuU8v58/zNx27jl35j7nPc/7vAJjjIGDg4ODg4ODg4PjBwTb+W4ABwcHBwcHBwcHR1PBSSwHBwcHBwcHB8cPDpzEcnBwcHBwcHBw/ODASSwHBwcHBwcHB8cPDpzEcnBwcHBwcHBw/ODASSwHBwcHBwcHB8cPDpzEcnBwcHBwcHBw/ODASSwHBwcHBwcHB8cPDpzEcnBwcHBwcHBw/ODQqknsnDlzkJOTc76bwXEW+OUvf4mcnBwcP378fDclZpw5cwZ//OMfMW7cOPTu3Rs5OTnYv3//+W7WDxrLly9HTk4Oli9frvp83LhxGDduXEz7ni/s3r0b8+bNw9ChQ5GTk4OpU6ee7yadU/wQv7Mc5wbN+S621ufn+PHjyMnJwS9/+cuzOs/WrVuRk5ODJ554ooVapkdr5D45OTmYM2fO+W6GDvamHiB37MUXX4z3338fcXFxun3GjRuHsrIy7N27F3Z7ky/RqnHw4EG88MIL2Lp1K06cOAFRFJGamorMzEwMGDAA1157Ldq3b/+9tOWJJ57Ak08+iZdeeglDhw7Vbf/lL3+JFStW4OOPP0bnzp2/lzadK8hE55NPPjnn13r44YexZMkSjB07FlOmTIEoiqZjunTpUjz44IO4/PLL8dhjj1me99lnn8Wjjz6KOXPm4Ne//nWT2hRtrFsCc+bMweeffx7z/tOmTcNf/vKXc9KW1oIzZ85gwYIFaGxsxNSpU5GSkvK9fb/PFb6PZ6m5qK6uxkcffYR169ahtLQU5eXlcDgccLvdmD59OmbMmAGbzTj2sn37djzzzDMoKiqC1+tFt27dMGPGDMyZMweiKKr2ra2txdKlS7F//37s378fhw8fhiRJWLRoEYYPH647N2MMGzZswKeffoovv/wS33zzDbxeLy6++GKMHj0aCxYsaNJzIY+BEnFxcbjoooswfPhwLFiwAB07doz5fByEnJwcDBkyBC+//PL5bkoYW7duxY033hh1v5KSku+hNeb4Pt+xLYlmM8xvvvkGL774IubPn9+S7VHhr3/9KxoaGs7Z+ZuKzz77LPxC69+/P0aNGoWkpCRUVFRgx44d2LRpEwYMGPCDf8ld6Fi3bh0yMzPx7LPPRt134sSJeOihh/Dxxx+jqqoKqamphvsxxvDmm28CAGbNmtWi7W0pTJs2DUOGDFF99tFHH6G4uBjjx49Hr169VNu0f7cULr30UuTl5SE9Pf2cnL8p2LVrFyorK3HPPffg1ltvPd/N+dHj/fffx+9+9zt06NABQ4cOxcUXX4yTJ0/iww8/xK9//Wts2LABjz32GARBUB330Ucf4a677kJcXByuuOIKtGvXDmvXrsVDDz2E7du34/HHH1ftX1ZWhr/97W8AgI4dOyIlJQUnT540bZfP58NPf/pTOBwODB48GMOHD4ckSdiyZQteeuklvPvuu/jvf/+LzMzMJt3vkCFDwt+56upqbNq0Ca+++iree+89LF26FF27dm3S+Voarem7eLbIyMjAu+++izZt2pyX63fq1AnTpk07L9duCbz77rtwuVznuxk6NIvEtmvXDoIg4Pnnn8fMmTNNX9xni4svvvicnLe5+O1vf4vGxkb85S9/MXwYi4uL0a5du/PQMo6WREVFBQYPHhzTvomJiZg0aRKWLl2Kt99+G/PmzTPcb8uWLTh27Bj69++P7Ozslmxui2H69Om6z8rKylBcXIwJEyYYbj8XaNOmzXl70WhRUVEBAD+Kl/gPAZmZmXjmmWdQUFCgirj+/Oc/x9VXX40PPvgAa9aswWWXXRbedubMGTz44IOw2Wx46aWXkJubCwC4++67cdNNN+GDDz7A6tWrMXHixPAxF198MRYvXoxevXohOTk5vGplBpvNhrvvvhuzZ89W/cYHg0H87ne/w5IlS/CXv/wlpomvEkOGDMGdd94Z/tvv9+OnP/0pPvvsMzzzzDN46KGHmnS+lkZr+i6eLRwOB3r27Hnert+pUyfVWP/QcD77zgrN0sTGx8fjtttuw+nTp/HUU0/FfNzy5ctx5513Yvz48ejXr194+f3tt9823F+rC1m9ejVycnLw5z//2XB/n8+HwYMHY+TIkQgEAqptq1atwpw5czBo0CDk5ubiiiuuwNNPPw2fzxdT2ysrK3HkyBG0adPGdDbl8Xhw0UUX6T4/ceIE/vjHP6KwsBD9+vXDkCFDMHPmTF3fbdmyBQ8++CCuvPJKDBgwAP369cOkSZPw5JNPorGxUbXvuHHjwstRN954I3JycsL/AbSsIv8ojx8/PrxNqz+sqanBP/7xD1xxxRXo168fBg4ciJtuugkbN27U3YdSH7Vu3Tpce+21yM/Px+DBg3HXXXfh8OHDMfWljHfffRfXX389Bg4ciH79+mHy5Ml47rnnVGMi64/KyspQVlamus9YtU0VFRX4/e9/j3HjxqFv374YNmwYFi5ciD179qj2k583xhg+//zz8HWi6YCuueYaAMAbb7xhus/SpUsBqKOwhw8fxn333YdRo0ahb9++GDlyJO677z5dP0YbaxkNDQ147rnnMHXqVOTn56N///6YNWsWVq1aZd1BzcCePXvwxz/+EVOmTMGQIUOQm5uLwsJC/OUvf8GpU6fO6txWOrwTJ07gD3/4A8aPH4++ffti6NChuPXWW7Fr1y7dvmfOnMFTTz2FSZMmYcCAAejfvz8mTJiAu+++Wzf2Wsj6ufvvvx8A8MADD4T7XG6XrP07duwYXn75ZUyePBn9+vVTPS+xjjFAS8w5OTnYunUrVq1ahenTpyMvLw8jR47EQw89FP5efPbZZ5gzZw4GDBiAwYMH4xe/+AWqq6tj6ttYnyUZr7/+OiZPnozc3FwMHz4cDz74IE6fPm24b1PGxgyXXHIJxo0bp5MMdOjQAddeey0A6CQv77//PqqqqjBx4sQwgQVoef5nP/sZAOC1115THdOuXTtccsklSE5OjqldDocDt912my5IYbPZcMcddwCg36qzhcPhCP+eKPtN/h367rvv8Ktf/QqjRo1Cr169VN+RoqIi3HXXXRgxYgT69u2LMWPG4De/+Q3Ky8t11zl27BgefPBBXHrppeF30uTJk/Gb3/xG9SxZfRc3b96M2bNnIz8/H0OGDMHtt9+Or7/+2vL+mtJG+fc4EAjg2WefRWFhYfiYv/3tb6r3hNxOAKrfbqV21UwTe+jQIfz973/H9OnTMWzYMPTt2xdjx47Fgw8+iBMnTljez/nAhg0b8NOf/hRDhw5F3759MWHCBPz1r39FbW2t4f7R+Ees71izd+Hp06fxj3/8A5dddhlyc3MxePBg/OQnP8HmzZt1+yr1xPv378f8+fMxaNAg5OXl4YYbbsD27dub3B/NlhNcf/31+O9//4slS5Zgzpw5MS2j/O53v0NWVhYGDx6MDh06oKamBp9++inuu+8+HDp0CHfffbfl8RMmTECbNm2watUq3HfffTq97UcffYTa2lrcfPPNqm0PPPAAli9fjo4dO6KwsBBt27bFzp078dhjj+Gzzz7DokWLomp327RpA7vdjvr6elRUVMQcmdm9ezduueUW1NTUYPDgwbj00kvh9Xrx1Vdf4cknnwz/AALAv/71Lxw6dAj9+/fHmDFj4PP5sH37djzxxBPYunUrFi9eHNZ23Xjjjfj444/x+eefY9q0aejUqZPqugsXLgwvBd94441o27Zt+D5klJWVYc6cOSgrK8OgQYMwatQoNDQ0YO3atbjlllvwhz/8IfyDqsSaNWuwYcMGTJgwAUOGDMH+/fvxwQcfYOvWrXj
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plot_naive_scatter2 = sns.scatterplot(data=pol, x=\"long\", y=\"lat\")\n",
    "_ = plot_naive_scatter2.set(\n",
    "    xlabel=\"Longitude\",\n",
    "    ylabel=\"Latitude\",\n",
    "    title=\"Naive Scatteplot of Vote Tallies from the 2012 Presidential Election\"\n",
    ")\n",
    "\n",
    "plot_naive_scatter2.figure.savefig(\"./images/graphs/2-naive-scatter2.png\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "This plot looks a little weird, which probably has to do in part with the fact that we are very thoughtlessly trying to project a portion of a geoid onto a cartesian plane using a measurement that is not cartesean coordinates\n",
    "\n",
    "Another issue with these plots is the fact that the data is very messy. It's not all that meaningful for more reasons than just the fact that it's messy, but being messy doesn't help all that much\n",
    "\n",
    "One thing that stands out is the outliers in the bottom left of the plot, which is clearly Hawaii. Hawaii is not covered by the Gaybourhoods data set at all so we can remove these rows like so:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAArEAAAHTCAYAAAAwDkDhAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/P9b71AAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOxdd3xV5f1+zjl3JiETEpBhZNwEyGIoCAhhGBdgWaIoitSCWrRSq9a2dv3a2tqqdVXpEBxVQUBFcOAAERSskIShAZQNkkgGIePu9/fHe8+9Z7xn3Juwz/P5+JHcs97znvW83/f5Pl+OEEJgwYIFCxYsWLBgwcJZBP50N8CCBQsWLFiwYMGChXhhkVgLFixYsGDBggULZx0sEmvBggULFixYsGDhrINFYi1YsGDBggULFiycdbBIrAULFixYsGDBgoWzDhaJtWDBggULFixYsHDWwSKxFixYsGDBggULFs46WCTWggULFixYsGDBwlkHi8RasGDBggULFixYOOtwRpPYmTNnIi8v73Q3w0Ib8POf/xx5eXk4dOjQ6W6KaTQ1NeEPf/gDxowZg379+iEvLw9ff/316W7WWY3ly5cjLy8Py5cvl/0+ZswYjBkzxtS6pwvbtm3DrbfeiiFDhiAvLw/XXnvt6W7SScXZ+MxaODlI5Fk8U++fQ4cOIS8vDz//+c/btJ9NmzYhLy8PTz31VDu1TI0zkfvk5eVh5syZp7sZKtji3UDs2AsuuADvvfcenE6nap0xY8bg8OHD2LFjB2y2uA9xRmPPnj14/vnnsWnTJhw9ehSCICAzMxO5ubkYOHAgrr/+enTs2PGUtOWpp57C008/jRdffBFDhgxRLf/5z3+ON954Ax999BG6det2Stp0siASnY8//vikH+uRRx7B4sWLMXr0aEycOBGCIGhe0yVLluChhx7ClVdeiSeeeEJ3v8899xwef/xxzJw5E7/61a/iapPRtW4PzJw5E1988YXp9SdNmoQ///nPJ6UtZwqampowd+5c+Hw+XHvttcjIyDhlz/fJwqm4lxJFfX09PvzwQ6xduxa7du1CdXU17HY7PB4PJk+ejClTpoDn2bGXLVu24Nlnn0VlZSW8Xi8uvPBCTJkyBTNnzoQgCLJ1GxsbsWTJEnz99df4+uuvsW/fPoRCISxcuBDDhg1T7ZsQgk8//RSffPIJvvzySxw5cgRerxcXXHABRo4ciblz58Z1X4jXQAqn04kuXbpg2LBhmDt3Ljp37mx6fxYo8vLycMkll+Cll1463U2JYtOmTbj55psN19u5c+cpaI02TuU3tj2RMMM8cuQIXnjhBcyZM6c92yPDX/7yF7S2tp60/ceLzz//PPpBGzBgAC677DKkpKSgpqYG5eXl2LBhAwYOHHjWf+TOd6xduxa5ubl47rnnDNe95ppr8PDDD+Ojjz5CXV0dMjMzmesRQrB06VIAwPTp09u1ve2FSZMm4ZJLLpH99uGHH6Kqqgpjx45F3759ZcuUf7cXLr/8chQXFyM7O/uk7D8ebN26FbW1tZg/fz5uv/32092ccx7vvfcefvvb36JTp04YMmQILrjgAhw7dgwffPABfvWrX+HTTz/FE088AY7jZNt9+OGHuPvuu+F0OnHVVVchLS0Na9aswcMPP4wtW7bgySeflK1/+PBh/PWvfwUAdO7cGRkZGTh27Jhmu/x+P370ox/Bbrfj4osvxrBhwxAKhbBx40a8+OKLeOedd/Df//4Xubm5cZ3vJZdcEn3m6uvrsWHDBrzyyit49913sWTJEvTo0SOu/bU3zqRnsa3IycnBO++8gw4dOpyW43ft2hWTJk06LcduD7zzzjtwu92nuxkqJERi09LSwHEc/vnPf2Lq1KmaH+624oILLjgp+00Uv/nNb+Dz+fDnP/+ZeTNWVVUhLS3tNLTMQnuipqYGF198sal1k5OTMX78eCxZsgRvvfUWbr31VuZ6GzduxMGDBzFgwAD06dOnPZvbbpg8ebLqt8OHD6Oqqgrjxo1jLj8Z6NChw2n70ChRU1MDAOfER/xsQG5uLp599lmUlpbKIq4//elPMW3aNLz//vtYvXo1rrjiiuiypqYmPPTQQ+B5Hi+++CIKCwsBAPfccw9uueUWvP/++1i1ahWuueaa6DYXXHABFi1ahL59+yI9PT06a6UFnudxzz33YMaMGbJ3fDgcxm9/+1ssXrwYf/7zn00NfKW45JJLcNddd0X/DgQC+NGPfoTPP/8czz77LB5++OG49tfeOJOexbbCbrejV69ep+34Xbt2lV3rsw2ns+/0kJAm1uVy4Y477sCJEyfwzDPPmN5u+fLluOuuuzB27FgUFRVFp9/feust5vpKXciqVauQl5eHP/3pT8z1/X4/Lr74YowYMQLBYFC2bOXKlZg5cyYGDx6MwsJCXHXVVfjHP/4Bv99vqu21tbXYv38/OnTooDmays/PR5cuXVS/Hz16FH/4wx9QVlaGoqIiXHLJJZg6daqq7zZu3IiHHnoIV199NQYOHIiioiKMHz8eTz/9NHw+n2zdMWPGRKejbr75ZuTl5UX/A+i0ivhSHjt2bHSZUn/Y0NCARx99FFdddRWKioowaNAg3HLLLVi/fr3qPKT6qLVr1+L6669HSUkJLr74Ytx9993Yt2+fqb4U8c477+DGG2/EoEGDUFRUhAkTJmDBggWyayLqjw4fPozDhw/LztOstqmmpga/+93vMGbMGBQUFGDo0KGYN28etm/fLltPvN8IIfjiiy+ixzHSAV133XUAgNdff11znSVLlgCQR2H37duH+++/H5dddhkKCgowYsQI3H///ap+NLrWIlpbW7FgwQJce+21KCkpwYABAzB9+nSsXLlSv4MSwPbt2/GHP/wBEydOxCWXXILCwkKUlZXhz3/+M44fP96mfevp8I4ePYrf//73GDt2LAoKCjBkyBDcfvvt2Lp1q2rdpqYmPPPMMxg/fjwGDhyIAQMGYNy4cbjnnntU114JUT/3wAMPAAAefPDBaJ+L7RK1fwcPHsRLL72ECRMmoKioSHa/mL3GAJ1izsvLw6ZNm7By5UpMnjwZxcXFGDFiBB5++OHoc/H5559j5syZGDhwIC6++GLcd999qK+vN9W3Zu8lEa+99homTJiAwsJCDBs2DA899BBOnDjBXDeea6OFSy+9FGPGjFFJBjp16oTrr78eAFSSl/feew91dXW45pprogQWoNPzP/nJTwAAr776qmybtLQ0XHrppUhPTzfVLrvdjjvuuEMVpOB5Hj/+8Y8B0HdVW2G326PvE2m/ie+h77//Hr/85S9x2WWXoW/fvrJnpLKyEnfffTeGDx+OgoICjBo1Cr/+9a9RXV2tOs7Bgwfx0EMP4fLLL49+kyZMmIBf//rXsntJ71n87LPPMGPGDJSUlOCSSy7BnXfeiW+//Vb3/OJpo/g+DgaDeO6551BWVhbd5q9//avsOyG2E4Ds3S3VrmppYvfu3Yu//e1vmDx5MoYOHYqCggKMHj0aDz30EI4ePap7PqcDn376KX70ox9hyJAhKCgowLhx4/CXv/wFjY2NzPWN+IfZb6zWt/DEiRN49NFHccUVV6CwsBAXX3wxfvjDH+Kzzz5TrSvVE3/99deYM2cOBg8ejOLiYtx0003YsmVL3P2RsJzgxhtvxH//+18sXrwYM2fONDWN8tvf/ha9e/fGxRdfjE6dOqGhoQGffPIJ7r//fuzduxf33HOP7vbjxo1Dhw4dsHLlStx///0qve2HH36IxsZGzJ49W7bswQcfxPLly9G5c2eUlZUhNTUVFRUVeOKJJ/D5559j4cKFhtrdDh06wGazoaWlBTU1NaYjM9u2bcNtt92GhoYGXHzxxbj88svh9XrxzTff4Omnn46+AAHgX//6F/bu3YsBAwZg1KhR8Pv92LJlC5566ils2rQJixYtimq7br75Znz00Uf44osvMGnSJHTt2lV23Hnz5kWngm+++WakpqZGz0PE4cOHMXPmTBw+fBiDBw/GZZddhtbWVqxZswa33XYbfv/730dfqFKsXr0an376KcaNG4dLLrkEX3/9Nd5//31s2rQ
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "pol = pol[pol[\"long\"] > -140]\n",
    "\n",
    "plot_naive_scatter3 = sns.scatterplot(data=pol, x=\"long\", y=\"lat\")\n",
    "_ = plot_naive_scatter3.set(\n",
    "    xlabel=\"Longitude\",\n",
    "    ylabel=\"Latitude\",\n",
    "    title=\"Naive Scatteplot of Vote Tallies from the 2012 Presidential Election\"\n",
    ")\n",
    "\n",
    "plot_naive_scatter3.figure.savefig(\"./images/graphs/3-naive-scatter3.png\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAlcAAAJKCAYAAAABc2DBAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/P9b71AAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOydeZwcVbn+n1PVe/fsM9mGhADJJGEXiCAiyL4omwpuRAUVF3ABvCIXrteruFzcEQU3lh9wBZSwJiiyKDsJZAFCNkLWmWQy+0zv3VXv749TVdNLdXdVT2eZyfv9fDRDd71d55yqOuepU6feRxARgWEYhmEYhqkJyp4uAMMwDMMwzESCxRXDMAzDMEwNYXHFMAzDMAxTQ1hcMQzDMAzD1BAWVwzDMAzDMDWExRXDMAzDMEwNYXHFMAzDMAxTQ1hcMQzDMAzD1BAWVwzDMAzDMDWExdUE4ZRTTsGcOXPw6quv7umiVM2cOXMwZ84cbNu2bU8XhWFc8Z3vfAdz5szBb37zmz1dlHHBggULMGfOHCxcuHBPF8U1e3PZt23bZvWjzJ7FM5bgRCKBhx56CM899xzWrFmDgYEBCCHQ3NyMQw89FKeeeirOPPNMBAKBWpV3j2F2mp/97GdRX1+/h0uzb7F69Wo89dRTaG9vx0c+8pE9XRxmH2J4eBh33XUXAOBrX/vaHi7N7mNv7e++853v4KGHHsJ73/te3H333TXbthZMhHNlcHAQf/nLX/Dvf/8b7777LmKxGBoaGtDa2opZs2Zh/vz5OOGEEzB9+vS8uKeeegqrV6/Ge9/7Xhx77LE1L9d4HAOqFlfPPPMMvvvd76Knp8f6LBQKQQiBzs5OdHZ24h//+Ad+9rOf4aabbsL73ve+mhR4T3HLLbcAAC688MK9qrOZSBxwwAEAAK/Xm/f56tWrccstt+C9733vuLmwmInB8PCwde2XGzDb2tpwwAEHoKmpaXcVbZfC/V1ppk6digMOOAB1dXV5nzs9V/ZWVq5cia985Svo6+uzPotEIkgkEli7di3Wrl2LRYsW4dRTT8Xvfve7vNinnnoKDz30EK688spdJq7G2xhQlbhauHAhrr/+eui6jgMOOABf+cpXcOKJJ1ody8jICF566SXcc889WLJkCV577bVxL66YXc/f//73PV0EhqmKa665Btdcc82eLgazG7jpppv2dBFqzvDwsCWsZs6ciSuvvBKnnHIKwuEwAKCnpwdLly7FE088ASHEHi7t+MC1uFqzZg3++7//G7qu46STTsLNN99c9Nivrq4OZ555Js4880wsXrwYO3bsqFmBGYZhGIapHYsWLUJfXx98Ph/uuusuTJkyJe/7trY2nHPOOTjnnHOQSqX2UCnHF67F1a9+9Suk02lMnjwZP//5zyuupzrnnHNAREWfp9Np3HvvvVi8eDHeffddZDIZTJ06FR/84AfxhS98AW1tbUUxv/nNb3DLLbfgwgsvxE9+8hPb/ZnP2a+88sq8qdlXX30Vn/nMZ9De3o5nnnkGr7/+Ov7whz9gxYoVSCQSmDlzJi6++GJ8+tOfzlPm5u+ZnHrqqXn7u/DCC/HjH/8YZ5xxBrZs2YL/+q//wiWXXFKyPS655BIsXboUX/rSl3D11VeXbjgbHn30Udxzzz1Yv349vF4v5s6di8suuwwf/OAHK8bGYjHcfffd+Oc//4lNmzYhnU5j6tSpOPHEE/H5z38eU6dOLYpZsGABlixZgh//+Mc455xz8Mc//hGPP/44tm/fjnA4jOOOOw7f+MY3MHPmTNt9PvXUU7j//vuxatUqDA0NIRwOo6mpCYcccghOO+00nHPOOXnbm4swn376aey33355nwHAkiVLihZq/r//9/+gKAouueQSeL1ePP/88yUfzWzduhWnn346iAhPPPEEDjzwwIrtBgAbN27EokWLsGTJEmzbtg09PT3w+/048MADcdZZZ+FTn/pU2esgk8lg4cKFWLx4MdauXYtoNGo9RjrzzDNx7rnnIhQKAZALUs1zbO3atVixYgVuv/12LFu2DH19fbjkkktw/fXXW7/9yiuv4J577sHy5csxNDSEhoYGHHnkkbjkkktKzhZHo1HcddddePrpp7Fx40ZkMhk0NjZi0qRJOPbYY3HhhReio6MjL2bJkiW45557sGLFCvT398Pv96OlpQWzZs3CiSeeiIsvvhiK4u79GDftYlJNv7Fw4UJcd9111tqbZ555BnfccQdWr14NTdMwe/ZsfOYzn8GHP/zhvDjz/DcpPPdy+5hS/U7h8Vy3bh1uvfVWLFmyBMPDw2hvb8e5556LL37xi/D5fCXb6rXXXsO9996L119/Hf39/QiHw5g3bx4+9rGP4UMf+lDRbMKu7O/MvjedTuPpp5/Gs88+izVr1qC7uxvxeBytra046qijcOmll+LQQw8tWac9SW5foygKfve73+G5555Df38/Jk2ahDPPPBNXXHEFIpFIUWxuv2g+onJzrphs27YNd9xxB1544QXs2LEDiqLggAMOwFlnnYVLLrmk6Nw3SaVS+NOf/oTHHnsMnZ2daGhowPz583HFFVdUvb553bp1AIB58+YVCatC/H6/9bd5npnccsst1qNRk7Vr11p/r1q1Ck888QRef/11bN++Hb29vQiHw+jo6MB5552Hj3zkI1BVNS/eyRhQ+CjS7fUCyPHhT3/6E1555RVs377dWj++33774YQTTsBFF12E5ubmsm2Tiytx1d3djX/9618A5MlU+My5FIUV6e/vx+c//3m8/fbbAACfzwev14tNmzbhzjvvxEMPPYQ//OEPOPLII90UzzELFy7EDTfcAF3XEYlEkEqlsHbtWvzgBz/A5s2b8wavSCSC1tZW9Pb2AgCampryDn4kEoEQAh/96Efxy1/+EgsXLiwprrZs2YLXXnsNAFw/N/7+97+Pe++9FwCgKAo8Hg+WLFmCV199Na+8dmzYsAFf/OIX0dnZCQDweDzw+XzYvHkz7r77bjz66KO49dZbcfTRR9vGR6NRfPKTn8Tbb78Nn88HRVHQ39+PxYsX46WXXsJf//pXzJgxIy/ml7/8JW677Tbrv8PhMJLJJDZt2oRNmzbh1VdfLRJXdrS2tiKZTCIajcLr9aKhoSHve6/Xi6OOOgozZ87Epk2b8Pjjj2PBggW2v/Xggw+CiHDUUUc5FlaAfOSzatUqALJjCYVCGBoawsqVK7Fy5UosWrQId911l21H3N3djS996UtYvXo1AHns6uvr0dPTg66uLrz44ouYOXOm7TqFxYsX4z/+4z+QzWZRV1dX1OnktrEQAnV1dejr68NTTz2Fp556CpdffnnRo6qRkRF84hOfwDvvvGOVx4zr6enBqlWroKoqvvWtb1kx999/P7773e9a/x0MBqHrOjZv3ozNmzfj6aefxoUXXpjX6VaimnapRb/x29/+FjfffDMURUE4HEY8HsfKlStxzTXXoLe3F5/73OesbRsaGtDU1ISBgQEA8lzMpdTgV4oXXngBV1xxBZLJJOrq6pDNZrFx40bcfPPNWLVqVdE6FpOf/vSn+NOf/mT9dyQSwdDQEF5++WW8/PLLeOaZZ/Czn/2spLitdX9n8uKLL+Kb3/wmAHn+1dfXQwiBrq4udHV14e9//zt++MMf4oILLnDVTruTNWvW4Prrr8fg4CDC4TCICJ2dnbj99tuxZMkS3HfffUVrQO1we648+eST+Na3vmXNAgWDQWQyGaxatQqrVq3CY489hjvuuKPod2KxGC699FKsXLkSgOz/EokEFi9ejH/961/4wQ9+UHVbAPLxHxE5fvTn9XrR2tqKkZERpFIphEKhstfFZZddhsHBQQCyzsFgEIODg1iyZAmWLFmCf/7zn/jd734Hj2dUmjgZA3Kp5npZtWoVFixYgFgsZv1mMBi0zuUlS5Zg3rx5OPHEEx21CwCAXPDII49QR0cHdXR00DvvvOMmNI/Pf/7z1NHRQfPnz6fFixdTNpslIqI33niDPvzhD1NHRwcdf/zx1NfXlxd38803U0dHB1177bUlf/vaa6+ljo4Ouvnmm/M+f+WVV6ijo4OOOOIIOuS
      "text/plain": [
       "<Figure size 600x600 with 3 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# A prettier way to represent the number of observations in a given region\n",
    "# would be to use a hexbin plot:\n",
    "plot_naive_hexbin = sns.jointplot(data=pol, x=\"long\", y=\"lat\", kind=\"hex\")\n",
    "\n",
    "_ = plot_naive_hexbin.figure.suptitle(\"County density across continental United States\")\n",
    "_ = plot_naive_hexbin.ax_joint.set(xlabel=\"Longitude\", ylabel=\"Latitude\")\n",
    "plot_naive_hexbin.figure.subplots_adjust(top=.925)\n",
    "\n",
    "plot_naive_hexbin.figure.savefig(\"./images/graphs/4-plot-naive-hexbin.png\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "There are quite a few counties concentrated around the east coast, above Florida and below New Hampshire. On the other hand, there are far fewer individual counties along the west coast. This is all pretty self-evident."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Number of data points in the region of Boston: 131\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAlEAAAHBCAYAAABJ8u4AAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/P9b71AAAACXBIWXMAAA9hAAAPYQGoP6dpAACwoElEQVR4nOzdeZxcVZn4/8+5W+3Va7rT6ez7RghhVxAxEBE3RgGRr1FAQAVGxtEZdQTFkR8qw4iMzow4KogzKCiDG4qOICgqsmUhCSQhZO8t6a32utv5/VFJJ52u6vTe1Z3zfr14Affcqjp1q7rqqbM8j5BSShRFURRFUZRB0ca7A4qiKIqiKBORCqIURVEURVGGQAVRiqIoiqIoQ6CCKEVRFEVRlCFQQZSiKIqiKMoQqCBKURRFURRlCFQQpSiKoiiKMgQqiFIURVEURRkCFUQpiqIoiqIMgQqiFEVRFEVRhsAY7w4cj+/7XHHFFWzYsIGzzz6b+++/v6ftySef5IknnmD9+vU0NTURCASYNWsWl112GZdccgmGMfCn5/s+P/7xj3nkkUfYvn07AA0NDbzlLW/hU5/61Eg/LUVRFEVRJriyD6K+//3v9wQ1x7r11lsJhUJccMEFzJs3j2QyyWOPPcbnPvc5fvvb33LvvfcihDjuY7iuy9/+7d/yhz/8gbe97W38zd/8DUII9u/fz/79+0f6KSmKoiiKMgmUdRC1d+9e7rnnHj7xiU9wxx139Gm/6667OOuss3oFSh/60IdYu3YtTz/9NH/4wx8477zzjvs49957L0899RTf/va3Offcc0f0OSiKoiiKMjmV9ZqoW265hfnz57N27dqi7WeffXafkSZd17nooosA2Lp163EfI5vNct9993H++edz7rnnIqUklUohpRz+E1AURVEUZdIq25Gohx9+mBdeeIFHHnkETRtcrNfa2gpATU3Ncc998cUXSSaTrFixgq9+9av8+Mc/JplMEolEeOtb38o//uM/UlVVNaTnACClxPf7D8g0TRz3nBOZuj79U9fn+NQ16p+6PsenrlH/Jtv10TQxoOVAZRlEtba2cuedd3L11VezePHiQd22paWFhx56iIqKClavXn3c83fs2AEU1l7pus7NN99MXV0dv/vd7/jf//1fNm3axE9+8hMCgcCQnguArh8/CNT1479YJzJ1ffqnrs/xqWvUP3V9jk9do/6diNenLIOo2267jaqqKm666aZB3S6dTnPDDTeQSqX4xje+QWVl5YBuA9Dd3c3Pf/5z5s+fD8Bb3/pWAH7+85/z05/+lPe9732DexKH+L4kkciUbNd1jXg8RCKRxfP8IT3GZKauT//U9Tk+dY36p67P8alr1L/JeH3i8dCABkDKLoh67LHHePLJJ7nvvvsIBoMDvl06neb6669ny5Yt3HrrrVx44YUDut3hx1ixYkVPAHXYe9/7Xn7+85/z7LPPDjmIAnDd47+pPM8f0HknKnV9+qeuz/Gpa9Q/dX2OT12j/p2I16esgijbtrn99ts555xzaGxsZPfu3b3ac7kcu3fvJhKJUFtb23M8lUpx3XXXsW7dOm677TauuOKKAT/m1KlTAairq+vTNmXKFKAwSqUoiqIoinK0sgqicrkcHR0dPPPMM6xZs6ZP+7p161izZg0XX3wxd999NwDJZJJrr72WjRs3cvvtt3PppZcO6jFPPvlkAJqbm/u0tbS0AANboK4oiqIoyomlrIKoUCjEPffcU7Tt5ptvZuHChdx44400NDQAhQDqmmuuYfPmzXz5y1/mkksu6ff+29raSCaTTJs2jVAoBEBjYyNnnHEGzz//PBs3bmTFihU95z/44IMAnH/++SPw7BRFURRFmUzKKogyTbMnx1MxNTU1vdqvuuoqNm3axOrVqxFC8LOf/azX+YsWLeq1u+9rX/sajz76KA888ABnnnlmz/Fbb72VK6+8kquvvpoPfOAD1NXV8eSTT/LMM89w7rnn9tsnRVEURVFOTGUVRA3Wpk2bAHjiiSd44okn+rTfdNNNA0qRsHDhQh566CHuuecefvSjH5FOp2lsbOTjH/8411133aDzVCmKoiiTi+d52LY9aXafjSTfF+RyOradx/PKN1eUpulomjag/E8DJaRKzT2qPM+noyNdst0wNKqqInR2pk+4XQ0Doa5P/9T1OT51jfqnrk//ClUsOsnl0iqA6oemafh++V8fTdOJRisJhSL9BlPV1ZGJmeJAURRFUcpFItFBNpumoqIKXbdGdBRjMtF1UdajUIXqIR65XIZEoh3HyVNRMfxNYyqIUhRFUZQifN8jm00RjVYRi1Wokbp+GIY2Ia5PMBgmlTJJpbqJxSrRNH1Y96cW+yiKoihKEZ7nAWBZQy/7pQydEKBroAmJJmCklidbVhCQPa/vcKiRKEVRFEXph5rCG3uaBtJ18NwjgY4QAj0QYLhL00by9VQjUYqiKIqilA0hAM9Dur1HiqSUePnciI1IjYQy6oqiKIqiKCc6IcB3neKNEiijXYAqiFIURVEUpbz0t9GvjDIzqSBKURRFUZTy0t+ypTKaz1MLyxVFURRlCM4557RBnf9P//QFLr74nT3/n0gkeOSRh/jzn//I3r17yedzVFVVs2zZSbzrXZdw+uln9br9pZe+k5aW5gE/3tVXX8eHP/wRvvvde7nvvv/q+f/DDh8H+Pu//zTvec9lfe7jV7/6BXfc8UU++MFruP76G4reFgqLtUOhMLFYjDlz5rJy5Sre+taLmTKlrs99Njc3cdll7zpu/3/84I9pmNpALpfj6uuvZn/Tfv796//Oyaed3mekatOmjdx443XU10/l/vsfJByOHPf+R4IKohRFOS7DKPzymwh5YBRlrFx99XV9jv34xz8klUpx2WXvJxqN9mpbsGBhz3+vX/8St9zyj3R1dTF79hzWrLmIcDjCvn17+MtfnuH3v/8db33r2/j0p2/FsiwALr/8/SSTyV73+etf/5KWlmbe9rZ3MHVqQ6+2U045dcDP5b77/ouLLrp40MHHypWrOOWUU9E0QSaTpb39IC+/vIFnn/0z3/vet7n66utZu/aqoreNRqNcdtn7i7YJAfHKKhAQDAa59bO38rGPf4zbv3I7993/Q8LhcM+5mUyGL33p8wDceus/j1kABSqIUhSlHwFTI6D72J0HkZ5PqKoKTzPJ5FUwpShHj+oc9utf/5JUKsXll7+fhoZpRW+3c+fr/MM/3Ew+n+cTn/gH3vOey3ttu29tbeGzn/0Uv/nNrzEMk89+thAgXH75lX3ua926F3uCqFWrBjcydtj06TPYt28vP/jB/XzkIzcO6rannHIqH/7wR3ol25RS8vTTT3LnnXdw773fBCgaSEWjsaLX8DAhDu3Uk7D8lFV86EPX8L3v/Rf33HNXzzUB+Ld/+1f279/H2rVXc9JJJw+q/8NVPhOLiqKUlaClIboP0vni86Rf30Fm90661r9EfvcOoqHhZflVlBPZ179+F9lsliuv/CDvfe/7+uQtqq+fyp13fp1YLM5jj/2cl1/eMKr9ee9730dt7RQefvhB2tpah31/Qgje/ObV3H77VwG4//7/4uDBg4O+HykLG/H8Q//+4Ac/zJIlS3nssZ/zhz88BcAf//gUv/zlz1i4cHG/AdloUUGUoihFmbikd+7oc9xuP4jX3dEzxacoysA1Ne3nxRefw7IsrrzygyXPq62t5Z3vfDcAP/vZ/45qn4LBINde+1Hy+Tzf/vZ/jNj9rlp1GitWrCSfz/P0008O+/4Mw+DWW/+ZYDDInXf+f7z22na++tX/j0AgwOc//yUMY+wn19R0nqIofZimTr55f8n27L59hJdU4bpj2ClFmQQ2blwPwKJFi4nH4/2ee/rpZ/Lggz8Y9ZEogIsvfic//vEP+e1vf8373nclCxYsGpH7PeWUU9m4cT2vvLK5T1sqleS737236O1qamq45JJL+xyfOXM2H/vYx7n77ju5/voPYds2f/d3n2L27Dkj0t/BUkGUoih9aAI8xy7Z7ruOGsZWlCFoby9Ma9XV1R/33Lq6qQAcPHhgVPsEoGkaN9xwM5/85N/
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Boston is one of the cities covered by the Gaybourhoods data set.\n",
    "# Let's take a closer look at it:\n",
    "\n",
    "gb_ny = gb[gb[\"long\"] > -73]\n",
    "print(f\"Number of data points in the region of Boston: {gb_ny.shape[0]}\")\n",
    "plot_boston_scatter1 = sns.scatterplot(data=gb_ny, x=\"long\", y=\"lat\", hue=\"TOTINDEX\")\n",
    "_ = plot_ny_scatter1.set(\n",
    "    xlabel=\"Longitude\",\n",
    "    ylabel=\"Latitude\",\n",
    "    title=\"Queer Neighbourhoods in Boston\"\n",
    ")\n",
    "\n",
    "plot_boston_scatter1.figure.savefig(\"./images/graphs/5-plot-boston-scatter1.png\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "In the above plot, we colour points darker the higher their \"gaybourhood index\" is. This index is a weighted composite of a few factors. The TOTINDEX boldy weighs the queerness of same-sex married couples more than unmarried same-sex households. While it's subjective, in my opinion, this marginalizes the experiences of queer people who have no desire to marry.\n",
    "\n",
    "If we wanted to look at unmarried same sex households alone, we could use the \"SS_Index\" column:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAlEAAAHTCAYAAAAH7E2QAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/P9b71AAAACXBIWXMAAA9hAAAPYQGoP6dpAADRnklEQVR4nOzdd5hU1f348fctU7cXdmGXXhaQDoqCIFIUazRRozESu8aWxJ+aaIyJydcSSzR2McbejV0UC6iIjQ7Se11gF7ZNn7n3nt8fwy6MO7O9zC7n9Tx5IvfcmXvumdmZz5zyOYoQQiBJkiRJkiQ1itreFZAkSZIkSeqIZBAlSZIkSZLUBDKIkiRJkiRJagIZREmSJEmSJDWBDKIkSZIkSZKaQAZRkiRJkiRJTSCDKEmSJEmSpCaQQZQkSZIkSVITyCBKkiRJkiSpCWQQJUlSvWbMmMHAgQOb/TyPPPIIAwcO5Icffmjza7e1gQMHMmPGjPauRoM15bVpiB9++IGBAwfyyCOPtOjzSlIy0Nu7ApLUWD/++COvvPIKCxYsoLS0FF3X6d69OxMnTuSiiy6iS5cu7V3FNlEdWBQUFDB79mwcDketc6ZMmcKuXbtYtWoVui7/3KXO6ZFHHuHRRx+NOaYoCikpKQwYMIDTTz+dc889t83+Bn744Qd+85vfcO2113Lddde1yTWl9iE/VaUOQwjB/fffz9NPP42u64wfP56TTjqJSCTC0qVLefrpp3nllVd44IEHmDx5cntXt80UFxfz/PPPc8UVV7TaNe655x4CgUCrPb/UeQ0fPpyPPvqIrKysVr/W2LFjGTt2LACGYbBnzx7mzp3LP/7xD5YsWcK//vWvVq+DdHiRQZTUYTz22GM8/fTTFBYWMnPmTAYMGBBT/sknn3DTTTdx3XXX8fLLLzNixIh2qmnbycjIQFEUnnrqKc4++2yys7Nb5ToFBQWt8rxS5+dyuejXr1+bXGvs2LG1en6Ki4s59dRT+fDDD7n++uvp3r17m9RFOjzIOVFSh7Bz506eeOIJbDYbTzzxRK0ACmD69OnccsstRCIR/va3v8WU1TXfY+fOnQwcOJCbb765VlkgEGDmzJmcccYZjBw5klGjRnHuuefy4YcfJqzr119/zeWXX87RRx/N0KFDmTZtGvfccw9VVVW1zp0yZQpTpkzB6/Vy9913M2XKFIYMGdLg+SNOp5OrrroKj8fDY4891qDHVFu+fDm/+93vOPbYYxk6dCiTJk3ir3/9K3v37q11bqJ5SeFwmEceeYSpU6cydOhQpkyZwoMPPkg4HK53TtDs2bM5++yzGTFiBGPHjuX666+Pe+1Dr/Xggw8yZcqUmnZ99NFHCYfDcc//7rvvuPTSSxk7dixDhw5l+vTp3H///Xg8nlrnVr8O8SR671TfX2lpKbfeeisTJ05k8ODBvP322zHnlZWVcdtttzFhwgSGDh3KqaeeyltvvRX3WpZl8eqrr3LWWWcxatQoRo4cyVlnncUrr7yCZVnNvk+AlStXcumllzJq1ChGjx7NRRddxNKlS+OeC7Bo0SJ++9vfctxxxzF06FCOPfZYfvnLX9YaPksk0Zyo6veUYRg8+eSTnHjiiTXvw/vuuy/h69pYBQUF9OnTB4i+Fj+1cuVKrrvuOsaNG8fQoUOZPHkyt99+OyUlJbXO3bdvH/fccw/Tp09n5MiRHHnkkUyfPp2bb76ZHTt2AHDzzTfzm9/8BoBHH32UgQMH1vzv0PdQOBzmqaee4vTTT2fEiBGMHj2a888/n48++qjWdQ/9jNq5cyfXX389Rx99NMOGDeMXv/gFX3zxRYu0ldR4sidK6hDefvttDMPg5JNPrnOS8TnnnMNjjz3GmjVrWLZsGSNHjmzyNauqqrjwwgtZvXo1Q4YM4ayzzsKyLObPn88NN9zAhg0buP7662Me8+ijj/LII4+QmZnJ8ccfT3Z2NuvXr+eZZ55h3rx5vP7666SmpsY8JhwO85vf/IbKykqOPfZYUlNTG/Vr+de//jUvv/wyr7/+OjNmzKB37971PuZ///sff/3rX7Hb7UyZMoWuXbuybds23nzzTebOncsbb7xRb++TEILrrruOL7/8kt69e3PBBRdgGAbvvPMOGzdurPOxr7zyCnPnzmXKlCkcddRRrFixgo8++oi1a9fy3nvvYbfbaz3m97//PT/++CMnnXQSuq4zZ84cHnnkEVauXMkTTzyBoig157722mvcfvvtuFwuTjrpJHJycliwYAH/+c9/+OKLL3j11VdJT0+vt53qU1FRwbnnnovb7ebEE09EURRycnJqyquqqvjVr36F3W5n+vTphMNhZs+ezZ///GdUVeXnP/95zPPddNNNfPjhh3Tr1o2zzz4bRVH4/PPP+fvf/87ixYtrDUc19j6XLFnCxRdfTCQS4YQTTqBXr16sWbOGGTNmcMwxx9S6v3nz5nHllVeSmprKlClTyM/Pp6Kigs2bN/PKK69w7bXXNrsNb7jhBhYvXszEiROZNGkS8+bN4+mnn6asrIy777672c+/e/dutmzZQkpKCn379o0p++KLL2p6rqZPn05BQQGrVq3i1VdfZc6cObzyyiv06NEDiP6g+tWvfsX27ds59thjmTJlCkIIiouLmTNnDtOnT6dHjx5MmzYNgHfeeSdmeBGgsLAQiP7NX3rppSxYsIC+ffty/vnnEwwG+eSTT7j++utZu3Yt/+///b9a97Jr1y7OOeccevTowRlnnEFlZSUfffQRV199Nc8++2zc11BqZUKSOoDf/OY3oqioSLz++uv1nvv//t//E0VFReKpp56qOfbwww+LoqIi8f3339c6f8eOHaKoqEj86U9/ijn+pz/9qdbzCCFEMBgUl1xyiRg4cKBYvXp1zfHvvvtOFBUViXPPPVdUVlbGPOatt94SRUVF4s4774w5PnnyZFFUVCQuvPBC4fP56r23QxUVFYmJEycKIYT4+OOPRVFRkbjmmmviPn8kEqk5tnnzZjFkyBAxbdo0sWfPnpjzv/32WzFo0CBx9dVXxxy/4IILRFFRUcyxd955RxQVFYnzzz9fhEKhmuOVlZVi+vTpoqioSFxwwQUxj6l+HUaNGiXWrl0bU1b9us2aNSvutU888URRUVFRczwYDIpf/vKXoqioSLzzzjs1x3fu3CmGDBkiRo0aJTZu3BjzXH/7299EUVGR+Mtf/lKrnSZPniziSfTeKSoqEkVFReKmm26Kad+flv/5z38WhmHUHN+wYYMYPHiwOPnkk2PO/+CDD0RRUZE488wzhdfrrTnu8/nEz3/+c1FUVCTef//9Jt+nZVk1r8tnn30Wc/5zzz1XU99D7/Paa68VRUVFYs2aNbXub//+/XHb66e+//57UVRUJB5++OGY49Wv689//nNRXl4ec7/Tpk0TgwYNEiUlJQ26RvVrdMEFF4iHH35YPPzww+KBBx4QN998sxg7dqwYO3as+OSTT2Ie4/V6xdixY8WgQYPEwoULY8pmzpwpioqKxMUXX1xzbM6cOXH/hoUQIhQKCY/HU+89V3vyySdFUVGRuOyyy2LeO/v27av5m128eHHN8erPqKKiIvHII4/EPNe8efNqnktqe3I4T+oQSktLAejatWu953br1g2gzqGh+pSXl/P+++8zdOhQLr/88pgyh8PBTTfdhBCCDz74oOb4iy++CMD//d//1erl+MUvfsHgwYNjzj/UzTffjNvtbnJ9TzrpJEaNGsVnn33GokWL6jz31VdfJRKJcOutt5Kfnx9TNm7cOKZMmcIXX3yB1+ut83neffddAP7whz/E9Bylp6dz9dVX1/nYeMOD55xzDhBdfRnPVVddRUZGRs2/HQ5Hza/1Q4fH3n//fSKRCBdccEGtuTjXX389KSkpvPfeey0yXGSz2fjTn/6UcNWXy+XilltuQdO0mmP9+/dn9OjRbNq0CZ/PV3O8+h5uuOEGUlJSao673W5uuukmAN58880m3+eSJUvYsmULRx11VE1vSbULLriAnj17JrzPeCs/W2r+3Y033khmZmbNv91uN6effjqWZbFy5cpGPdeCBQt49NFHefTRR3n
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plot_boston_scatter2 = sns.scatterplot(data=gb_ny, x=\"long\", y=\"lat\", hue=\"SS_Index\")\n",
    "_ = plot_boston_scatter2.set(\n",
    "    xlabel=\"Longitude\",\n",
    "    ylabel=\"Latitude\",\n",
    "    title=\"Queer Neighbourhoods in Boston\"\n",
    ")\n",
    "\n",
    "plot_boston_scatter2.figure.savefig(\"./images/graphs/6-plot-boston-scatter2.png\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "\n",
    "As predicted, the diversity of the data becomes more apparent when we don't hide it. While the difference isn't extremely substantial, it will be worthwhile to continue to analyze multiple factors when exploring this data. This is especially the case when we consider the fact that the data set is already very limited in how it fails to account for gender diversity.\n",
    "\n",
    "Just for fun, let's do something similar for the political dataset:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAApgAAAHTCAYAAACDRmj7AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/P9b71AAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOydd3wU553/31O2q/cu0SSBBAKJXk03NtjG3cQ9js/xOb74ckmc5JK7Sy/nJBfHceJfEtdgGxtjY8A2GIPpvRcBAiTUe98+M78/ll202pVAQGLjzPv14vVCM7MzzzxTns98n28RNE3T0NHR0dHR0dHR0blKiJ91A3R0dHR0dHR0dL5Y6AJTR0dHR0dHR0fnqqILTB0dHR0dHR0dnauKLjB1dHR0dHR0dHSuKrrA1NHR0dHR0dHRuaroAlNHR0dHR0dHR+eqogtMHR0dHR0dHR2dq4ouMHV0dHR0dHR0dK4qusDU0dHR0dHR0dG5qugCU+eKaGxs5Nvf/jbTp09n+PDh5OXl0dHR8Vk3q1/y8vK47777PutmfCGpqqoiLy+Pp59++u96nHfeeYe8vDzeeeedv+tx/pH8o/runwG9L7+YPP300+Tl5VFVVfVZN0XnEpAH+oO8vLygv0VRJCoqiry8PO644w4WLVp01Rr3eefZZ5/l97//Pa+88goTJkz4rJvzmfD000+zdetWbrzxRrKzsxEEAZPJ9Jm2adasWQB88sknn2k7dC6fnTt3cv/99/PEE0/wta997bNujs5V4Omnn2bFihWsX7+ejIyMz7o5QVRVVTF79uygZSaTCZvNRmZmJoWFhdxwww2MHTv2M2rhFwv9Hf3Z8I/u9wELTD9PPPEEAF6vlzNnzrB+/Xp27tzJkSNH+M53vnPVGqjz+cXtdrNt2zYmT57MM88881k355JZs2YNFovls26GzhUwd+5cioqKSEpK+qybovMFIjIykgceeADwjW3t7e2Ulpby+uuv87e//Y2pU6fyi1/8goSEhM+4pTo6n38uW2D2tips376dhx56iJdffpn77rvvc/eFqnP1aWpqQlXVa26QHzJkyGfdBJ0rJDIyksjIyM+6GTpfMKKiosJazCsrK/nud7/Lli1beOSRR3jzzTc/85kaHZ3PO1fNB3PSpEkMHjwYTdM4fPhwYHldXR0//OEPmT17NoWFhUyYMIHHHnuMQ4cOhezj2WefJS8vj507d/L+++9zxx13MGbMmIBZF8DhcPDCCy9w6623MmbMGMaMGcOCBQv48Y9/TFNTU9D+HA4Hf/rTn7j55psZPXo0Y8aM4a677mLVqlUhx965cyd5eXk8++yzHD9+nEcffZSxY8dSVFTEvffey759+4K2nzVrFr///e8BuP/++8nLywv883P27Fn+93//l1tvvZWJEydSWFjIzJkz+f73v09dXV3YfnS73Tz77LOB/po1axa/+c1vcLvdffoOer1e/va3v3HnnXdSXFxMUVERt9xyC6+99hqqqoY9Tl+Ul5fzrW99i2nTplFYWMjUqVP51re+RXl5ecj5z5w5E4AVK1YEzr0/n6f6+nqGDx/OLbfc0uc2jzzyCHl5eZw8eTJo+Zo1a/jSl75ESUkJo0aNYtGiRfzpT3/C7XYHtvFfw+rqaqqrq4OuSc92hevHnvfehx9+yO23305RURHjx4/nqaeeor6+Pmx7Dx06xMMPP8yYMWMoLi7mwQcfZP/+/UH7uxT8vkWVlZW8+OKLXH/99YwcOZLp06fz05/+lK6urrC/O3LkCF/72teYNGlS4P767//+bxoaGq7KMfrzVx2IP9RAnoWnn36a+++/H4Df//73QdfR35/9+WBeTp9UVVXxxhtvsGjRIkaOHMnkyZP5/ve/T2dn50XPDeCNN94gLy+PZcuWBS1fvnw5eXl5FBUVBd2rAHfccQcjR47E6XSG7K+qqoqnnnqKCRMmMHLkSG699VY2bNgQ9thut5sXXniBRYsWUVRURHFxMUuWLGHNmjWX1HaA66+/nsLCQlpaWsKuf+GFF8jLy+O1114LWn6pfZ2Xl8eKFSsAmD17duB69ny3A7S1tfHMM8+wYMECRo0aRUlJCQ888ABbtmwJ266uri5+9rOfMX36dEaOHMn111/Piy++iKZpl3zul0JmZiYvvPACgwcP5vjx47z++ush21zuWLdq1SpuvfVWioqKmDp1Kj/72c8C98r27du57777KC4uZty4cXzzm9+ktbU1bBsHct/DpY+lPd8br776KosWLWLUqFGB94Lb7ea1117jK1/5CjNnzqSwsJDx48fz4IMP8umnnwYd81Lf0RdD07SLvsMURWHGjBkUFxfT3d0ddj8/+tGPyMvL48MPP7yk4w5EfzQ0NPA///M/zJo1i8LCQiZOnMgTTzzBkSNHQvbb33jRlz/xQN5dF+v39vZ2ioqKmDNnTp/PzmOPPUZeXl6QvrsYl23BDIe/YYIgAHD06FEefvhh2tvbmTp1KvPmzaO1tZWPP/6YJUuW8NxzzzFjxoyQ/bz44ots3bqVmTNnMmHChEBHtbe3c//991NaWsqgQYO47bbbMBgMVFZWsnz5cubOnRuYuujo6OCBBx7g2LFjFBQUcNttt6GqKlu2bOEb3/gGp06d4qmnngo59pEjR/jzn//M6NGjueOOO6ipqWHt2rU8+OCDvPvuuwwePBjwicr169eza9cuFi9eTHp6esi+1q1bxxtvvMGECRMoLi7GYDBw6tQp3nrrLTZs2MDy5ctJTk4O6r+vfe1rbNy4kZycHO699168Xi8rVqygrKwsbJ97PB4ee+wxtmzZwqBBg1i4cCEmk4mdO3fyox/9iIMHD/KrX/3qkq7foUOHeOihh+ju7mbWrFkMHTqUM2fOsHLlStavX8+LL77IqFGjAudfXV3NK6+8Qn5+PnPmzAFg+PDhfe4/OTmZyZMns2XLFk6cOBHiz9vQ0MC2bdsoKCggNzc3sPzXv/41f/rTn4iNjWXhwoVYrVY2b97Mr3/9a7Zs2cJf/vIXjEYj6enpPPHEE7z88ssAgamui7WrJ0uXLuWTTz5h1qxZjBs3jkOHDrFmzRpKS0t57733MBqNgW13797Nww8/jKqqzJ07l6ysLE6ePMn999/PxIkTL+l4vfnpT3/Knj17WLBgAZGRkWzZsoWXX36ZPXv28PrrrwdZTTZs2BCwtsyfP5+0tDSOHj3K66+/zvr161m6dCmZmZlXdIyrxUCeBf+9tGLFCsaPH8/48eMD+wn3nPXkcvvkV7/6FVu2bGHmzJlMmTKFnTt3smzZMioqKnjllVcuen6TJk0CfILgzjvvDCzfvn07AE6nk/379wd8tTs7Ozl69CglJSWYzeagfVVXV3PHHXeQmZnJzTffTHt7O2vWrOHxxx/nxRdfDLq33G43X/7yl9m1axeDBw9myZIlOJ1OPvroI5566ilKS0v593//94u2f/Hixfz6179m9erVYT8o3n33XQwGAwsXLgwsG0hfP/HEE3z88ceUlpZy//33ExUVBRBkha6urua+++6jurqasWPHMm3aNBwOBxs2bOCRRx7hhz/8YVDfut1uHnzwQQ4fPkx+fj6LFi2is7OTP/zhD+zateui5zxQLBYLDz/8MP/5n//J+++/z4MPPhhYd7lj3WuvvcamTZuYM2cO48ePZ+vWrbz00ku0t7cze/ZsnnrqKa677jruuusu9u/fz8qVK2ltbeXPf/5z0H4Get8PZCz185Of/IQ9e/YwY8YMZsyYgSRJgX395Cc/YcyYMUyePJm4uDgaGxvZsGEDjz76KD/+8Y+54447AK7KOxou7R0mSRJ33HEHzz77LKtXrw66d8D3TK5cuZLExMQQ/9twDKTPKisrWbJkCQ0NDUycOJEbb7yR2tpaPvzwQzZu3Mizzz4bMNBcKZfy7rpYv0dHR3PDDTfwzjvvsG3bNqZMmRJ0jNraWjZt2kRBQQEjR4689MZpAyQ3N1fLzc0NWb5161YtLy9Py8vL06qqqjSPx6PNmTNHKyws1Hbu3Bm0bV1dnTZ16lRtypQpmsvlCiz/3e9+p+Xm5mpFRUXa0aN
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plot_democrat_scatter = sns.scatterplot(\n",
    "    data=pol,\n",
    "    x=\"long\",\n",
    "    y=\"lat\",\n",
    "    hue=\"percent\"\n",
    ")\n",
    "_ = plot_democrat_scatter.set(\n",
    "    xlabel=\"Longitude\",\n",
    "    ylabel=\"Latitude\",\n",
    "    title=\"Percentage of voting population who voted Democrat by county\"\n",
    ")\n",
    "\n",
    "plot_democrat_scatter.figure.savefig(\"./images/graphs/7-plot-democrat-scatter.png\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "A better way to represent this data would be to use a heatmap, but even this scatterplot reveals something we probably already knew: Democrats tend to be concentrated around the urban parts of the United States"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Renewed Analysis Plan\n",
    "\n",
    "The research questions posed earlier in this analysis still feel sufficiently specific, having clear, measurable goals and room for additional analysis. Here are two very rough step-by-step plans to answer the two research questions I seek to answer in this analysis:\n",
    "\n",
    "**Is there a correlation between political alignment & living in neighbourhoods with large quantities of LGBT people?**\n",
    "1. Unify the political data with the gaybourhoods data set\n",
    "    1. Establish the best way to measure the distance from a given `gb` and a county (Euclidean distance? Some other measurement?)\n",
    "    2. Find the county that is closest to each observation by minimizing the function established in step (a)\n",
    "    3. Merge the two tables. Each `gb` observation should then include a political breakdown of the nearest county during the 2012 presidential election\n",
    "2. Use this information to plot queerness by different metrics against political alignment and measure the correlation\n",
    "\n",
    "**Is there a correlation between geographical stratums & being LGBT?**\n",
    "1. Explore different metrics of queerness and analize more qualitatively how different metrics reveal different information in different places by graphing them\n",
    "2. Graph and measure the clusteredness of neighbourhoods surpassing different threshholds of queerness\n",
    "3. If the previous steps show that this is a relevant line of inquiry: measure and graph the rate of change in queerness radially outward from clusters of queer neighbourhoods\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Analysis"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Visualiztion strategy\n",
    "Seaborn doesn't provide an easy interface for plotting data over maps. To do this, we will take a best-effort approach to rendering maps over an interval of coordinates and overlaying relevant plots on top of them using Matplotlib. The accuracy of these renderings will be limited by my access and knowledge of GIS tools, however, the goal is to represent the data spatially and the physical maps are entirely ornamental.\n",
    "\n",
    "To create a render for a city:\n",
    "1. Draw a square over the region in Google Earth\n",
    "2. Note the coordinates of the corners of the square\n",
    "3. Screen shot the map over the marked region\n",
    "4. Render the graph over the interval of coordinates\n",
    "5. Render the image beneath the data\n",
    "6. Adjust the position of the image to minimize obvious inaccuracies (for example, if a neighbourhood is plotted in a body of water)\n",
    "\n",
    "For example:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkwAAAJKCAYAAADePRreAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/P9b71AAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOy9d7xuRXU+/szs8tbTbgGkqCieqwFEY0GQSCyABaNGbF/FYDQkij9r1BiNJhFr7KiJGrsmorGDJSo2YgNFMQRQQVQQbjn97XvPzO+PmTV77X3ec88FLucCdx4+l3PO++4ye/bMmmeetWaNMMYYBAQEBAQEBAQErAm5rwsQEBAQEBAQEHBrRyBMAQEBAQEBAQHrIBCmgICAgICAgIB1EAhTQEBAQEBAQMA6CIQpICAgICAgIGAdBMIUEBAQEBAQELAOAmEKCAgICAgICFgHgTAFBAQEBAQEBKyDQJgCAgICAgICAtZBIEwBNxsPechDsG3bNvzoRz/a10W53eP2VNd/93d/h23btuGcc87Z10W51eOcc87Btm3b8Hd/93f7uigl3FrLFRBwSyDe1wUI2PdYXFzEueeei+985zv47W9/i6WlJUxOTuJOd7oTTjzxRDzpSU/CzMzMvi7m7RqXX345vvGNb+CQQw7Bn//5n+/r4tyq8aMf/QhPf/rTS58JIdBsNtFut3HooYfiyCOPxIknnogHPvCBEELso5JuDD784Q9jZWUFj3vc43DooYfu6+JsCLZt2+Z/f8UrXrGqPRDyPMeRRx4JAPjoRz+KY489Ftdccw0e85jHYDAY4O///u/xF3/xF7u91+te9zp85CMfwZYtW3DeeeeVbGG/38fnPvc5fPe738UVV1yBhYUFCCGwadMmHHXUUXjoQx+KU045BfV6fS88dcC+RlCY9nN86UtfwkknnYS3vvWt+MlPfoL5+Xk0m00sLCzgpz/9Kd72trfh5JNPxpe+9KV9XdTbNS6//HK8613vwuc+97ndHnfYYYfh8MMPR6PR2KCS3boxMzODLVu2YPPmzZBSYseOHfjJT36Cj370o3jmM5+Jk08+GT/4wQ/2dTFvNmZmZnD44Ydj69atq7776Ec/ine961247rrrblXl2ii8733vQ7/f3+Pj73znO+N5z3seAODtb387fv/736957M9+9jN87GMfAwC86lWvKpGlCy64ACeddBL+6Z/+Cd/61rdw/fXXQ0oJKSWuu+46fO1rX8NLX/pSnHTSSbeLNhgQFKb9Gp/85Cfxj//4jzDG4Mgjj8Tzn/98HHfccUjTFFmW4Yc//CHe8Y534Be/+AVe8pKXoNvt4slPfvK+LvZ+jY985CP7ugi3KvzXf/1XSVXJ8xxXXnklvvWtb+ETn/gEfve73+EZz3gGXvva1+Lxj3/8PizpzcPTnvY0PO1pT9vXxViFW0O5du7ciY997GM488wz9/icM844A1/96ldx6aWX4h/+4R/w4Q9/eNUxo9EIr3jFK6C1ximnnIJTTjnFf/fZz37Wf3f44Yfj2c9+Nh70oAd5QrWysoLvf//7+PjHP44f//jHuPjii3Hcccfd7GcN2LcICtN+iv/7v//D2WefDWMMHvrQh+Lcc8/FiSeeiDRNAQBJkuBP/uRP8MlPfhIPfehDYYzB2Wefjcsvv3wflzwgYG3EcYwjjzwSz33uc3HeeefhmGOOgTEGr371q3HFFVfs6+IF7GU86EEPAgB84AMfQKfT2ePzoijC61//eiRJgh/84Af41Kc+teqY97znPfj1r3+N6elpvPrVr/afX3HFFXj1q18NrTVOPPFEfP7zn8djHvOYkvo0MTGBU045BR/72Mfwtre9Da1W62Y8ZcCtBYEw7ad4+9vfjizLcMABB+BNb3oTkiQZe1wcx3jjG9+IrVu3IssyvOMd79jtdf/whz/gFa94BU488UQcffTReMhDHoI3vvGNWFlZWfOcb3zjG/irv/orHH/88TjyyCNx//vfH6eccgpe9KIX4ctf/vKa511wwQV49rOfjQc+8IE46qijcNxxx+Fv/uZv8L3vfW/s8Z/97Gexbds2nH766QCAL37xi3ja056GY489Ftu2bcM3vvENPOMZz8C2bdvwxje+cbfP+apXvQrbtm3DWWedVfr8sssuw5vf/GY85SlPwZ/+6Z/iqKOOwrHHHovTTz8dn/70p6GUWnWtbdu24eUvfzkA4Mc//jG2bdtW+scDvNcL+t61axfe8IY34OEPfziOOeYY3Oc+98Fpp52GD37wgxiNRmPP4cHXSil8+MMfxp/92Z/hmGOOwf3vf3/89V//NX7xi1+MPVcphe985zt41atehT//8z/H8ccfj6OOOgonnHACzjrrrH3qiti8eTPOOeccTExMIMsyvPvd717z2IsvvhgvfOEL8aAHPci/szPOOAPnnXcejDGrjv/Rj36Ebdu24SEPeQgA4Cc/+Qn++q//Gsceeyzuec974s/+7M/w8Y9/fOy5ADA3N4c3vvGNOPXUU3Gve90LRx99NE488UQ8+clPxjve8Y5V7rVxwdX0GR379Kc/vdRuqJ3T5+u16Ze97GXYtm0bXvziF+/2uPXKRaByXHvttfjDH/6AV77ylb5+yS7cGJIzDk960pNwhzvcAYuLi/jgBz94o8494ogj8JznPAcA8KY3vQnbt2/3311xxRX493//dwA2Rmrz5s3+u7e//e0YjUY48MAD8Za3vGXd+KRHPvKReMYznnGjyhZw60QgTPshbrjhBnz3u98FYCX1dru92+MnJia87P7tb38bN9xww9jjfve73+Hxj388/uu//gvLy8sQQuC6667DBz/4QTz+8Y/Hjh07Vp3ztre9DWeddRa++93vYm5uDrVaDYPBANdccw3OP/98vPa1r111TpZl+Nu//Vs8+9nPxgUXXIBdu3ahVqthfn4e3/rWt/CsZz0L//Iv/7LbZzr77LPxkpe8BD/5yU9gjIGUtiuceuqpAIAvf/nL0FqPPTfLMnzta18rHU/4y7/8S7z//e/HT3/6UywuLqLRaGBxcRE//vGP8cpXvhLPfvazked56ZwtW7b4d5AkCbZs2VL6txaZreLSSy/Fox71KHzoQx/Cb37zG0RRhCzL8Itf/AJvfOMb8YQnPAFzc3Nrnp/nOc4880y8/vWvx9VXXw0pJZaWlvDtb38bT33qU3HJJZesOueqq67CmWeeiXPPPReXXXYZhsMhkiTBzp078Y1vfANnnHEG3vve9+5R+W8JHHjggT6I/oILLkCv11t1zL/8y7/gqU99Kr785S9j+/btqNVqWFpawg9+8AO8+MUvxotf/OI12wJgifjpp5+O73znO1BKYTgc4sorr8RrXvMavO51r1t1/HXXXYfHPOYx+OAHP4hf/epXGI1GaDQa2L59Oy655BK85z3v8f1zd2g2m9iyZYtvu1NTU6V2MzU1BQB4whOeAMBOEKptj9DpdHyb3tuuyyuuuAKPe9zj8OlPfxqdTgfGGG8X/uIv/gJZlt3ka6dp6knPRz7yESwsLNyo888880zc4x73wMrKCl71qlcBsJOAV7ziFciyDA9+8IPxZ3/2Z/747du349vf/jYA4PTTT8fExMQe3ef2vvBgv4EJ2O/whS98wczOzprZ2Vnz61//eo/O+dWvfuXP+eIXv1j67sEPfrCZnZ0197nPfcxJJ51kLrroImOMMUop8/Wvf90ce+yxZnZ21jzjGc8onff73//e3P3udzezs7PmLW95i5mbm/Pfzc3Nma9+9avm5S9/+aqyvPa1rzWzs7PmpJNOMl/+8pdNt9s1xhizsrJiPvGJT5h73/veZnZ21nzpS18qnfeZz3zGzM7Omnvd615m27Zt5pxzzjFLS0v+3F27dpnl5WVz9NFHm9nZWfOjH/1obF1861vfMrOzs+be9763GQwGpe9e9KIXmfPOO8/s2LHDf9btds3nP/9588AHPtDMzs6a97///auuSWV72tOeNvaeBKrrH/7wh6XPFxcX/fVPPfVU8/Of/9wYY0ye5+YrX/mKud/97mdmZ2fNGWecseqaL3vZy8zs7Ky5733va+5
      "text/plain": [
       "<Figure size 600x600 with 3 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "gb_ny = gb.query(\"`long` > -74.2735 & `long` < -73.2336\")\n",
    "\n",
    "ny_map = Image.open(\"images/new-york-adjusted.png\").convert(\"LA\")\n",
    "ny_map.putalpha(int((5/8) * 256))\n",
    "\n",
    "## Colormap stuff\n",
    "color_array = plt.get_cmap('copper_r')(range(256))\n",
    "\n",
    "# Set the alpha level for each color. These will follow the graph of y=x^(1/10)\n",
    "# which will make regions with fewer observations less visible than those with more\n",
    "# observations. Additionally, we cap the alpha level to .5 so the map isn't hidden\n",
    "color_array[:,-1] = np.array([np.power(x, 1/10) for x in np.arange(0, .5, step=(1/256)*.5)])\n",
    "\n",
    "# Construct the reusable colormap\n",
    "copper_alpha = LinearSegmentedColormap.from_list(name='reds_alpha5',colors=color_array)\n",
    "\n",
    "ny_gb_density = sns.jointplot(data=gb_ny, x=\"long\", y=\"lat\", kind=\"hex\", cmap=copper_alpha, color=(1, .745, .4))\n",
    "\n",
    "plt.xlim((-74.2735, -73.4536))\n",
    "plt.ylim((40.3520,41.0608))\n",
    "\n",
    "_ = ny_gb_density.ax_joint.imshow(ny_map, aspect=\"auto\", extent=[-74.2735, -73.4536, 40.3520, 41.0608])\n",
    "\n",
    "ny_gb_density.ax_joint.set(\n",
    "    xlabel=\"Longitude\",\n",
    "    ylabel=\"Latitude\"\n",
    ")\n",
    "\n",
    "ny_gb_density.figure.suptitle(\"Observational Density in NYC\")\n",
    "ny_gb_density.figure.subplots_adjust(top=0.925)\n",
    "\n",
    "ny_gb_density.figure.savefig(\"./images/graphs/8-ny-gb-density.png\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Is there a correlation between geographical stratums and having same-sex partners?\n",
    "\n",
    "### Quantitatively measuring queerness\n",
    "\n",
    "Before we attempt to visualize queerness geographically, it's worth taking the time now to consider what that actually means.\n",
    "\n",
    "To begin, we must acknowledge the role statisticians have played presently and historically in systematically eliminating minorities. One local example of this is the way that Canada's Indian Act works to incrementally strip Indigenous people of their legal recognition of being Indigenous through the malitiously-named process of \"enfranchisement\"[^1], in a process many now refer to as \"statistical genocide.\" Discretely categorizing people enables oppression and marginalizes deviation. That second issue is particularly pertinent in the case of the queer community, which is predicated on \"bending rules,\" so to speak. For that reason, I'm hesitant to even use the phrase \"queer community,\" as it implicitly makes the assumption that the constituents of the so-called \"queer community\" have a universal experience, which is completely untrue.\n",
    "\n",
    "Somewhat depressingly, the article for which this dataset was produced admits the following:\n",
    "\n",
    "> Currently, there’s no comprehensive way to quantitatively measure gayborhoods, or even where LGBTQ Americans live. Most of the existing data sticks to a narrow view (i.e. traditional marriage, the male/female gender binary) of the queer spectrum and “rainbow-washes” any intersectionality of race, ethnicity, class, gender, and sexuality. This project aims to paint a slightly more complete picture, combining several metrics to create a gayborhood index, but even then it admittedly underweights and undercounts areas with non-binary and minority populations. Still, this is some of the most complete data that we have. (More about the limitations in the methodology section.)[^2]\n",
    "\n",
    "Essentially, this dataset fails to represent queerness outside the context of monogamous partnerships between cisgender people (or at least, those who have been statistically represented as such). For this reason, this analysis seeks to be very upfront that it only explores \"same-sex\" partnerships.\n",
    "\n",
    "Data Analysts at The Pudding created the `TOTINDEX` column in the Gaybourhoods dataset in an attempt to cumulatively measure queerness in a given neighbourhood. Before we tear it appart and try to come up with something better, let's take a look at how this looks spatially in New York City, and compare it to how the spread changes when we only factor in the number of same-sex partnerships that reside there:\n",
    "\n",
    "[^1]: https://www.cbc.ca/radio/unreserved/how-the-indian-act-continues-to-impact-the-lives-of-first-nation-people-1.5614187\n",
    "\n",
    "[^2]: https://pudding.cool/2018/06/gayborhoods/"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAABFEAAAH9CAYAAADS/+OVAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/P9b71AAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOx9Z4AUVdb2U1WdJjIkUREFdXtAchYJShDERQRFQQXFAJjQVZdXUFFXTKsgSjBhQlEUBCSIKEkRUEBEQBYBQXIOE7unu6vqvj+qzq1bPT3DzDADrnuf7/NdprrCjec895xzz1UYYwwSEhISEhISEhISEhISEhISEsVCPdMFkJCQkJCQkJCQkJCQkJCQkPhvgDSiSEhISEhISEhISEhISEhISJQA0ogiISEhISEhISEhISEhISEhUQJII4qEhISEhISEhISEhISEhIRECSCNKBISEhISEhISEhISEhISEhIlgDSiSEhISEhISEhISEhISEhISJQA0ogiISEhISEhISEhISEhISEhUQJII4qEhISEhISEhISEhISEhIRECSCNKBISEhISEhISEhISEhISEhIlgDSiSEhIlAqdOnVCZmYmVq1aVW7vHD9+PDIzMzF8+PBye2dJ3p2ZmYnMzEzs3bu33L97KvhvLDPhz16+8saGDRtw9913o3Xr1qhbty4yMzMxfvz4M10sCYky41TncEXoiPJGReqcPwMGDBiAzMxMzJw5M+HvJ5NbeXl5eOGFF9ClSxc0aNAAmZmZ6NSp0+kqvsSfEDNnzkRmZiYGDBhwposiIfGngOdMF0BC4q8KXdcxZ84cfPnll9iyZQuysrKQlJSEatWqoVatWmjRogUuvfRSNGrU6EwXFQCwefNmLFq0CDVr1sR11113povzX4cPPvgAubm56N27N84777wzXZxyx6pVq7B69WrUq1cPXbp0OdPF+VNg586duPXWWxEOh6GqKipXrgxVVZGcnFzqd23YsAEzZ87E6tWrcejQIcRiMVSrVg2NGjVC9+7d0bVrVyiKUgG1+Ovg0KFDmD59OlauXIndu3cjKysLfr8f55xzDho3boy///3vaNOmjWzHMkLqiIrF3r170blz50LXk5OTkZqairPPPhv169fHZZddho4dO8Lr9ZbpOyWRW0OHDsXKlSsBAKmpqahUqRIqV65ctor9l6A8xnd8Hw4YMABPPPFEsc/06tULmzdvBgDUrFkTS5YsKdO3E2HmzJnYt28funTpgnr16pXbeyUkJKQRRUKiQnD8+HEMGjQIv/76K7/m9/vBGMMff/yBHTt24LvvvkNaWhp++umnM1hSB5s3b8aECRPQqlWrYglErVq14PP5kJSUdBpLVzGoU6cOAJSZjIr48MMPsW/fPrRq1eqUjSiVK1dGnTp1UL169VMuV3lh9erVmDBhAnr37l2sEaU82/TPjs8++wzhcBgtWrTAG2+8gfT09FK/IxKJYOTIkZg9eza/5vf74fV6sW/fPuzbtw9fffUVGjdujHHjxuHss88uzyr8JcAYwxtvvIE333wTkUiEX09PT0dBQQG2bduGbdu24fPPP0fDhg0xYcIE2Y7FoKg5/L+oI84UKlWqxNs/Fovh6NGjOHz4MDZs2ICpU6eievXqeOyxx3D11VcnfP6cc85BnTp1kJaWVui3k8mtbdu2YeXKlfB6vZgyZQqaNGlS7vX7M6Kk47s0+PLLL/Hoo48WqQ+3bNnCDSgVgVmzZmH16tWoWbPmKRtR0tLSUKdOHZxzzjnlVDoJif9uSCOKhEQFYNiwYfj111+RkpKCe++9F9deey1fEOfl5WHDhg1YuHAhvvvuuzNc0tJj8uTJZ7oI5YYFCxac6SIkRP/+/dG/f/8zXYwy4c/aphWB33//HQDQvXv3MhlQYrEY7rrrLqxevRqqqqJfv37o378/LrroIgDAsWPHMHfuXEyYMAHr169Hv3798Omnn0oDQBwef/xxzJgxAwDQrl073H777WjRogUCgQAA4MCBA1i6dCnef/99bNy4Ebt27ZJtWAxOdQ7/lXTEmcL48ePRunVr/rdpmti+fTtWrFjBDfYPPfQQtm/fjqFDhxZ6/qWXXiry3SeTW9u2bQMABIPB/xkDSkXg3HPPxf79+7Fs2bKEEUYA8MUXXwCwIlD27dt3GktXelx55ZW48sorz3QxJCT+NJA5USQkyhnbt2/H8uXLAQDPP/887rrrLldEQWpqKi677DI89dRT+Oqrr85UMSUkJE4RBQUFAFCm7TsA8Morr3ADypgxY/DUU09xAwoAVK1aFQMHDsSnn36KKlWq4MCBAxg2bBgYY+VS/r8CPv30U25AGTp0KN599120a9eOG1AAyyt/880346uvvsKQIUOgqpL6SPx3QVVV/O1vf8PAgQMxd+5cdOzYEQAwYcKEUjtjTia3KJorJSXlFEos0bNnTwBwRRmKMAwDc+fOhaIouOaaa05n0SQkJMoBkklISJQztm7dyv9NRKco+P3+QtfEhHemaeKDDz5Az5490aRJE7Ru3Rp33303NmzYUOQ7N23ahNGjR+Omm27CFVdcgQYNGqB169YYMGAApk+fDsMwCj2TmZmJESNGALC2bVBiQfpPTBBYXNLANWvW4Nlnn8UNN9yAdu3aoUGDBmjTpg3uvPPOCo1Q2LFjBx5++GG0adMGjRo1wlVXXYUJEyYgGo0W+1xxCRRXr16NBx54AB06dECDBg3QvHlzdO3aFffeey8+/fRTmKYJwOkv8iLdeuutrrYTk7DF9+2UKVPQp08ftGjRApmZmTyst6RJD7du3YqHHnoIbdu2RcOGDXHVVVdh4sSJRda7JAkf49tk7969yMzMxIQJEwBY4cHx40Nsv5Mlpdy9ezeefPJJdO7cGQ0bNkTLli1xyy23FDk2AXeSxIKCAowfPx7dunVDo0aN0KZNGzz00EPYuXNnsW1VHEzTxPTp09G/f3+0atUKDRs2RKdOnTBy5Ejs2rWr0P3UjqtXrwYAjBgxgte7pMkXDx06hI8++ggAcNNNNxUZlg8AF198MZ588kkA1rj89ttvXb+XJOFfScbUkiVLcM8996Bt27Z87t599934/vvvi61LNBrFlClTcPPNN6NVq1Zo0KABOnbsiBEjRmD79u0Jnxk+fDhPZhmNRvHGG2/gmmuuQdOmTZGZmYmcnJxivwlYi71x48YBsGTt/fffX+z9Ho8HDz/8MFq0aFHot6NHj+LFF1/EVVddhcaNG6N58+bo06cP3nvvvSLnU3wdXn/9dXTv3h2NGzfGFVdcgWeffRbZ2dn8/l9//RX3338/2rZti0aNGuH666/HokWLEr67JP0lfl/EqlWrXGNx7dq1GDJkCFq3bo1GjRqhZ8+emDJlSpHGuKLmdXnoCKBs4wUAFi1ahEGDBuGyyy5D/fr10apVK3Tr1g0PP/ww5s+fX+RzJ0Np9OyECROQmZl50u0eM2bMQGZmJi6//HKuK8oLKSkpGDNmDN82+tprrxW6J1Fi2ZPJrfgxF9/P8f2Zn5+PN998E9dffz2aN2+Ohg0bomvXrnj22Wdx4MCBhGUXy5WTk4OXX36Zz7n4eXmqcsUwDN6vjRs3RqtWrTBkyBBs3Lix0HOlGd8lxZVXXonk5GQsXbo0oTxbsWIFjhw5gubNm5doC/Dx48cxZswYLiebNGmCHj16YOzYscjKynLdSzohUV/H66l4/TFnzhz0798frVu3RmZmJpdRJdEz27dvx5NPPolu3brxPr3mmmvw7LPPura3EypqTktInA7I7TwSEhWIQ4cO4fzzzy/Ts4wxPPDAA1i4cCE8Hg+SkpKQlZWFpUuXYtmyZRg9enTChdcdd9zBFWpSUhJ/bvXq1Vi9ejUWLlyI119/HR6PM/2rVauGgoIC5OXlwev1olKlSq53liS/RX5+vmsLSkpKCvx+P44fP47ly5dj+fLl6Nu3L5555pkytUdRWLNmDQYNGoRwOAzAivTZu3cvxo8fj+XLl6NVq1alfudnn33GF62A1Y6maWLXrl3YtWsXFi9ejN69e8Pv9yM5ORnVqlXD8ePHYZqmay87gEJtCVh
      "text/plain": [
       "<Figure size 1300x500 with 2 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "figure, axes = plt.subplots(1,2, sharex=True, sharey=True, figsize=(13, 5))\n",
    "figure.suptitle(\"Spatial distribution of Queer Comunities by Different Metrics\")\n",
    "axes[0].set_title(\"Queerness by TOTINDEX\")\n",
    "axes[1].set_title(\"Queerness by Same-Sex Partnerships\")\n",
    "\n",
    "sns.kdeplot(data=gb_ny, x=\"long\", y=\"lat\", weights=\"TOTINDEX\", color=\"orange\", fill=True, cmap=copper_alpha, ax=axes[0])\n",
    "sns.kdeplot(data=gb_ny, x=\"long\", y=\"lat\", weights=\"SS_Index\", color=\"orange\", fill=True, cmap=copper_alpha, ax=axes[1])\n",
    "\n",
    "plt.xlim((-74.2735, -73.4536))\n",
    "plt.ylim((40.3520,41.0608))\n",
    "\n",
    "for a in axes:\n",
    "    a.imshow(ny_map, aspect=\"auto\", extent=[-74.2735, -73.4536, 40.3520, 41.0608])\n",
    "    a.set_xlabel(\"Longitude\")\n",
    "    a.set_ylabel(\"Latitude\")\n",
    "    \n",
    "figure.savefig(\"./images/graphs/9-ny-queer-metrics.png\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The distinction is somewhat subtle but there remains a meaningful difference. It's clear that when we factor in gay bars and pride parades, we narrow our scope. This makes quite a bit of sense, as not all queer people frequent bars and not everyone has the privilege of celebrating their queerness. Going forward, we will measure the queerness of a neighbourhood based only on how many same-sex partnerships reside there.\n",
    "\n",
    "### Measuring gaybourhood clusteredness\n",
    "To facilitate answering this research question, we've added an additional column to the gaybourhoods dataframe that discretely classifies neighbourhoods into 7 categories labeled `0` through `6`, with zero indicating there are very few same-sex partnerships and 6 indicating that there are many. The choice to divide the dataframe into seven categories was arbitrary, although inspired by Alfred Kinsey's research into the fluidity of human sexuality[^3]. Similarly to the Kinsey scale, our factor will be linear. The consequence of this is that the vast majority of the observations will be lumped into the lower categories. Nevertheless, this is a fair approximation that will make further analysis significantly easier.\n",
    "\n",
    "As seen below, virtually all observed neighbourhoods have a same-sex partner index of less than 25:\n",
    "\n",
    "[^3]: https://en.wikipedia.org/wiki/Kinsey_scale"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 29,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAowAAAHTCAYAAACtA0lpAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/P9b71AAAACXBIWXMAAA9hAAAPYQGoP6dpAABq5ElEQVR4nO3deXxM1+M+8Gcmm+wrIiG0dBJFRBJBrEkQak+tIWqvokp1oVqttpZaaq0lag8JGrtKEGoXPpZobVVqSYiSRVZZ7++P/Ga+xszcTDKTjef9eqWac869c+7JnZln7j33jkQQBAFERERERBpIK7oDRERERFS5MTASERERkSgGRiIiIiISxcBIRERERKIYGImIiIhIFAMjEREREYliYCQiIiIiUQyMRERERCSKgZGIiIiIRL22gdHf3x/+/v4V3Q2tuLq6IiQkpKK7oRc7d+6Eq6srdu7cWdFdKZGpU6fC1dUV8fHxFd0VrcTHx8PV1RVTp06t6K4QVWqxsbFwdXXFsmXLtF5m2bJlcHV1RWxsbBn2rHxUpfdCqtwM9bkyV1dXAICTkxOioqJgYmKi0sbf3x8JCQm4du0aDA1L//AhISE4f/48bt26Vep1VHX//vsvIiIicP78eSQkJCAzMxPm5uaoW7cuvL290a1bNzRu3Liiu0mVRE5ODsLCwhAVFYW7d+8iJycHNjY2qFGjBjw8PNClSxf4+PhUdDfL1OnTp7F161bExcUhNTUVpqamsLOzg6urK7y9vRESEgKJRFLR3VS8vr3MzMwMb731Fjp37oxhw4ahWrVqZfb4rq6u8PHxwebNm8vsMd508r/xpk2b0KJFi4ruTrmR5wRN5syZg6CgoHLqDZWEXgOj3KNHj7Bx40aMGTOmLFavlQ0bNlTYY5c1QRDwyy+/4JdffkFhYSEaNWqErl27wsbGBpmZmbh16xY2b96MdevWYcaMGRg8eHBFd5kqWGZmJkJCQnDt2jVUr14dnTt3RvXq1ZGVlYWbN29i+/btSEtLe60D46pVq7Bo0SIYGhqibdu2eOuttyCVSvHw4UNcuHAB0dHRCA4O1umDrL716dMHzs7OEAQBiYmJOHz4MBYtWoSYmBhs3boVRkZGFd3F19LgwYPx3nvvwcnJqaK78tqaMGGC2vKGDRuWc09IW3p/ZbS2toZEIkFoaCj69u0LOzs7fT+EVlxcXCrkccvDL7/8gmXLlqFWrVpYuHAhvLy8VNokJSVh48aNSE9Pr4AeUmWzceNGXLt2DW3atMHKlSthbGysVP/8+XPcuXOngnpX9hISErB06VJYWFhg69atKkc5CgsLcfr0aRgYGFRQD9Xr06eP0tGnSZMmoU+fPrh69Sr279+PPn36VGDvXl92dnYV9t71pvj4448rugtUQnqfw1itWjV89NFHSE9Pxy+//KL1cjt37sTHH3+MgIAAuLu7w9PTEwMHDsSePXuU2snnbslP17i6uip+Xp4HqGneRm5uLkJDQ9GjRw80bdoUnp6eCA4Oxu+//67S9uV5YvHx8Zg8eTJatGiBJk2aICgoCMeOHVNZJj09Hb/++iuGDh2Kdu3aoXHjxmjZsiXGjh2Ly5cvaz0emjx8+BArV66EkZER1qxZozYsAoC9vT0+/fRTjBo1Sqn833//xYIFCxAUFISWLVuicePG8PPzwzfffIPExESltidPnoSrqyumTZum9jFyc3PRokULtGjRArm5uSr1f/zxBwYOHAgPDw80b94cEydOxL1799Su67///sPMmTPh7++vGLMJEybgr7/+UmkrNr9I09w++RzFhw8fYvPmzejRowfc3d3Vzh2NiIhAjx490KRJE/j6+uKbb77RGLz/+usvfPzxx2jVqpViLL/77jv8999/Om8nAGRkZGDOnDlo164dmjRpgi5dumD9+vUQBEFte03k+96gQYNUwiJQ9EHP09NTqezJkydYvnw5Bg4ciNatW6Nx48Zo06YNpkyZgn/++UdlHS+P/YMHDzBx4kS0aNECzZo1w4gRI/D3338DAJKTk/HNN9+gTZs2aNKkCd5//32cO3dObb/z8/OxZcsW9O/fH56enmjatCl69+6NsLAwFBYWar39cXFxKCgoQIsWLdSeEpNKpWjbtq3a09FxcXGYOHGiYgzat2+PGTNm4MmTJ0rtxo0bB1dXV2zatEllHYsXL4arqyu++uorrfusTo0aNdCpUycAwNWrVwEU7YM//vgjevbsCR8fHzRp0gSdO3fG3Llz8fz5c5V1vDzP+MSJEwgJCYGXl5eiTD4+58+fV3p9lc8BLM3rotz+/fsREhICb29vNGnSBF27dsWKFSvUvn7IX9OfPn2K6dOno23btmjYsKFifvSzZ8/w008/ITAwEB4eHvD29kZgYCCmTp2Khw8fqn38GzduYMyYMfD29kbTpk0xZMgQXLp0SaWdptcYeZ+ePHmCzz//HK1atYK7uzuCgoKwb98+jdtdEvL3rqysLPz000/o0KEDGjdujE6dOiE0NFTtc18QBISFhaFbt25o0qQJ2rZti++//77YAwba/D2eP3+ueL169XWqsLAQISEhcHV1xe7du/Wy/YDyPvbvv/9i0qRJaNWqFdzc3JT+JidPnsTo0aPRokULNG7cGB07dsRPP/2EtLQ0tes9c+YMgoOD4eHhAR8fH4wbNw537txRO4e9uLmvYnNDS7Ofv/y62LhxY3Tr1g2RkZEax+jUqVMYO3as4r2nffv2+Oijj3DmzBnF2Ojy/q1OmZx7GTx4MLZs2YJt27YhJCQE9erVK3aZ7777Dg0aNEDz5s1RvXp1pKam4vjx4/jiiy8UOwwAWFlZYcKECdi1axcSEhKUDms7OzuLPkZubi5GjhyJ8+fP4+2330ZwcDBevHiB6OhoTJ48GTdv3sSnn36qslxCQgL69euHOnXqoFevXnj+/Dl+//13jBs3DuvXr0fLli0Vbe/cuYPFixfD29sbHTp0gJWVFR4/foyjR4/i5MmTWLlyJdq1a6fdQKoRGRmJ/Px8dO/eHe+8806x7V89vXb48GFERESgRYsW8PT0hJGREW7fvo0dO3bg2LFjiIyMRM2aNQEAbdq0gYuLCw4ePIivvvoKlpaWSuuKjo5GamoqRowYoRJCDh06hJMnT6Jjx47w8fHBjRs3EB0djdjYWISHh+Ptt99WtH348CGCg4Px33//oWXLlujWrRseP36MqKgo/PHHH1i2bBn8/PxKO2RKZs2ahf/9739o37492rdvr3JEaf78+Th16hT8/PzQunVrxMbGYvv27bh//75KEDh27JjiU3JgYCCcnJxw7do1hIeHK04Z1qlTp9TbmZubi2HDhuHPP/+Em5sbevTogfT0dKxYsUJlfltxbGxsABR9YNDW//73P6xZswYtWrRA586dYWZmhvv37yM6OhpHjx5FeHg43NzcVJaTP1/q16+PPn36ICEhAYcPH0ZISAi2bduGUaNGwcLCAl27dlU8l0aPHo3o6GilU4B5eXkYO3YsTp06hbfeegvdu3eHiYkJYmNj8cMPPyAuLg7z58/XaltsbW0BFP0NCgoKtD6S+Ntvv2HGjBkwNjaGv78/HB0dcf/+fezYsQNHjx7F9u3bFX2ePXs2+vTpg/nz58Pb2xvvvvsuAODs2bNYvXo1GjRogG+++UarxxUjDwzycLt9+3YcOXIEzZs3h6+vLwoLC3Ht2jWsX78eJ06cwPbt22FhYaGynujoaJw8eRLt2rXDwIED8ejRIzRs2BATJkzA8uXL4ezsrHQE89XpCiV5XQSAadOmYefOnXB0dETnzp1hZWWFK1euYMmSJTh79izWr1+v8nqVmpqKAQMGwMzMDJ07d4ZEIoG9vT2ys7MxaNAgPHjwAK1bt4a/vz8EQcCjR48QExODwMBApeceUBSsf/31V3h4eKBfv3549OgRDh06hGHDhmH37t1Kr0linj9/jkGDBsHS0hJBQUFIT0/HwYMH8dlnn+HJkycqH9JLIy8vDyNHjsR///2Hdu3awcDAAEeOHMHChQuRm5urcjp31qxZ2Lx5M6pXr44BAwbA0NAQMTExiIuLQ25urto
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "ss_index_hist = sns.histplot(data=gb.query(\"`SS_Index` != 0\"), x=\"SS_Index\")\n",
    "\n",
    "_ = ss_index_hist.set(\n",
    "    title=\"National Gaybourhood Same-Sex Partnership Index Frequency\",\n",
    "    xlabel=\"Same-Sex Partnership Index\",\n",
    "    xlim=(0,30),\n",
    ")\n",
    "\n",
    "ss_index_hist.figure.savefig(\"./images/graphs/10-ss-index-hist.png\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "As such, we will choose 25 as our threshold for the highest level of queerness in a single community, such that anything higher than it will recieve a rating of `6`."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkMAAAHTCAYAAAAktxxFAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/P9b71AAAACXBIWXMAAA9hAAAPYQGoP6dpAABHcElEQVR4nO3deVhUdf//8dcgiLigKFrumTVg7vuSG0qaW+aSlUtqlpao6W33nVZ3dZtad2lmZamZZi6U+5a5pqmZS26ZS7elomBSBojgwjLn90df5ufIoDAMDHCej+vqKs45c+Z93vNheHXOZ85YDMMwBAAAYFJeni4AAADAkwhDAADA1AhDAADA1AhDAADA1AhDAADA1AhDAADA1AhDAADA1AhDAADA1AhDAADA1AhDkCQNGDBAQUFBni7Drc6ePauwsDA9+OCDCgoKUqNGjTxdkgN39fzDDz9UUFCQ9u7dm+vPnduCgoI0YMAAT5eRaa68Nma1d+9eBQUF6cMPP/R0KQVKfvud8RTCkBsFBQUpKChIISEhunHjhtNt2rVrp6CgIKWkpORydeaSmpqqsLAwfffdd2rbtq1GjBihoUOH3vFxvIZA9qUFm4z+CG/fvl316tVTnTp1tGXLllyuLm9Iex+JjIz0dCmQ5O3pAgqiCxcuaP78+Zn644ucERkZqV9//VV9+vTRm2++meXH58Zr+N///lfXrl3Lsf0DedHy5cv12muvqVixYvrkk0/UsGFDSVKdOnW0fv16BQQEeLhCmBFnhtysZMmSKlWqlGbPnq2YmBhPl2Naf/zxhySpXLlyWX5sbr2GFSpUUPXq1XNs/0BeM2vWLL388ssKDAzUokWL7EFIkvz8/FS9enWVLl3agxXCrAhDblakSBE9//zzunLlimbMmJGpx9zpWnm7du3Url07h2UrVqxQUFCQVqxYoe+//159+/ZV/fr11axZM40fP17x8fGSpOPHj2vYsGFq3Lix6tevr+eee+62p2WTkpI0bdo0tWvXTrVq1VJoaKg++ugjJSUlOd3+t99+07hx49SmTRvVqlVLLVq00NixY3X69Ol0244bN05BQUE6f/68FixYoG7duqlOnTqZvp79888/a+TIkWrevLlq1aqlkJAQvfHGG/bgkyYoKEj9+/eXJH300Uf2S1+ZnYvgymuY5siRIxo1apQefPBB1apVS23atNFrr72m6OjodNtmNG8nKSlJH374odq3b69atWqpXbt2mjZtmpKSku54/X/Dhg3q3bu36tatqyZNmmjMmDFOn/vm58rK6/3DDz9oyJAhatKkiWrVqqWOHTtqypQpunLlSrptnY3bNBnNpUk7vj///FOvvPKKWrVqpRo1amjFihUO28XExOjf//63WrZsqVq1aqlLly5avny50+ey2WwKDw9Xr169VL9+fdWrV0+9evXS4sWLZbPZsn2c0t9jc8iQIapfv74aNGigQYMG6dChQ063laQff/xRzz33nFq3bq1atWrpwQcfVJ8+ffTRRx9l+Jib3fyecejQIQ0aNEgNGzZU/fr1NWTIEB09etTp41JSUrRo0SL16dNHDRo0UN26dfXoo49q4cKF6XoRGRmpoKAgjRs3TmfOnNHo0aPVvHlzBQcHZ3kOlGEYmjhxot577z3dd999+vLLL3X//fdneEw3S/s9SUlJ0cyZM9WhQwf779a7777rdKxmpb/Xrl3TrFmz1L17d9WrV0/169fX448/rnXr1jlst3PnTgUFBWn8+PFOjzEpKUlNmzZV06ZNM/z9yYy034GsjPGkpCTNmDFDoaGh6d4zMpLZsXD+/Hk1atRITZo0UVRUlMM+rl69qk6dOqlGjRr5fl4cl8lyQL9+/bRo0SJ99dVXGjBggO65554ce65vv/1W27dvV9u2bfXEE0/o0KFDWrFihSIjIzV27Fj7m2Tv3r31v//9T9u2bVNkZKTWrFkjL6/0WfiFF17Q0aNH9fDDD8vb21tbt27Vhx9+qJ9//lmffPKJLBaLfdsdO3Zo5MiRSklJUUhIiKpUqaLo6Ght2rRJ27dv1xdffKGaNWume45Jkybpxx9/VJs2bdSmTRsVKlTojse5bds2jRw5UpLUsWNHVahQQceOHVN4eLi2bt2qxYsXq3LlypKkESNGKCoqSitXrlSTJk3UpEkTSbL/OzNceQ2XLVum1157TYULF1a7du109913KyIiQkuXLtW3336rJUuWqEKFCrfdh2EYGjlypLZv36577rlH/fv3V0pKilauXKlff/31to9dvHixvv32W7Vr106NGzfWTz/9pPXr1+vkyZNavXq1ChcunO4xWXm9v/zyS73xxhvy8/PTww8/rDJlymjfvn369NNPtW3bNoWHh8vf3/+OfbqTuLg4Pf744ypatKg6dOggi8WiMmXK2NfHx8frySefVOHChdWxY0clJSVpw4YNevnll+Xl5aUePXo47O+f//yn1q1bp/Lly6t3796yWCzasmWL/vOf/+jAgQOaOnWqw/ZZPc6DBw9q8ODBSk5O1kMPPaSqVavqxIkTGjBggJo1a5bu+Hbs2KFhw4apePHiateune666y7FxcXp9OnTWrx4sUaMGJHpXh05ckSzZs1SixYt1K9fP0VERGjz5s3av3+/5s6d6/ChgeTkZD333HPatWuXqlWrpq5du8rX11d79+7Vm2++qSNHjujdd99N9xznzp1Tnz59dM8996hbt266fv26ihcvnukak5KS9NJLL2n9+vWqX7++Zs6cqVKlSmX68WnGjh2rAwcOqFWrVmrTpo127NihOXPmKCYmRm+99ZZ9u6z0Nz4+XgMHDtTx48dVs2ZN9erVSzabTbt27dLYsWN16tQpjRkzRpLUsmVLValSRd98841efvlllShRwqG+jRs3Ki4uTk8//bTT37WsyMoYNwxDo0eP1tatW1WlShX1799fycnJWr58uf73v/853X9WxkLlypU1ceJEvfDCCxo7dqwWLlwob++/o8N//vMfnT59WiNHjlTTpk2zdcweZ8BtrFar0apVK8MwDOObb74xrFarERYW5rBNSEiIYbVajeTkZPuyPXv2GFar1fjggw+c7jckJMQICQlxWLZ8+XLDarUaNWrUMPbu3WtfnpqaagwaNMiwWq1G48aNjdWrVzs8bvz48YbVajU2b97ssLx///6G1Wo1OnToYMTFxdmXX79+3ejTp49htVqNlStX2pfHxcUZjRo1Mpo0aWKcOnXKYV+//PKLUa9ePePRRx91WP7SSy8ZVqvVaNmypXHu3Dmnx+pMQkKC0aRJEyM4ONjYv3+/w7pZs2YZVqvVGDx4sMPyO/U0I66+hqdPnzZq1qxphIaGGhcvXnTYfvfu3UZwcLAxfPhwh+VpPb/ZypUrDavVavTt29e4ceOGffnly5eNjh07Glar1ejfv7/DYz744APDarUa9evXN06ePOmw7h//+IdhtVqNr7/+2ulzZ/b1joyMNGrWrGnUr1/f+PXXXx329frrrxtWq9V49dVX0/Xp1nF7a8179uxxWG61Wg2r1Wr885//dOjvretffvllIyUlxb781KlTRo0aNYxOnTo5bL927VrDarUajz76qJGQkGBfnpiYaPTo0cOwWq3GmjVrXD5Om81mf11u/Z36/PPP7fXefJwjRowwrFarceLEiXTH99dffznt163SxrfVajUWLFjgsG7z5s2G1Wo1HnroISM1NdW+PK3nEyZMcOhdSkqK0/eF8+fP259j6tSpmarr1vp69Ohhfz8aNmyYce3atTs+5tbf2bSx2qNHDyM2Nta+PDEx0QgNDTWCg4ONP/74w748K/1Ne0+aPXu2w/Lr168bTz/9tBEUFGQcP37cvnzOnDlOe35znadPn87wGG+W9j5y/vx5h+VZHeNr1qwxrFar0adPH+P69ev25bGxsUb79u1v+56R2bFgGP9//E+ZMsUwDMNYsWKFYbVajQEDBjiMs/yKy2Q55OGHH1b9+vW1efNm/fjjjzn2PF26dHE44+Hl5aXu3btLku6//3498sgjDts/+uijkqSTJ0863d/zzz+vkiVL2n/29fXVP/7
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "gb_kinsey_bar = sns.barplot(\n",
    "    data=gb.groupby(\"kinsey\").aggregate(len).reset_index().drop(0),\n",
    "    x=\"SS_Index\", y=\"kinsey\", palette=\"Pastel2\", orient=\"h\"\n",
    ")\n",
    "\n",
    "_ = gb_kinsey_bar.set(\n",
    "    title=\"Number of Neighbourhoods per Kinsey Index\",\n",
    "    ylabel=\"Kinsey Index\",\n",
    "    xlabel=\"Count\"\n",
    ")\n",
    "\n",
    "gb_kinsey_bar.figure.savefig(\"./images/graphs/11-gb-kinsey-bar.png\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Excluded from the above graph is all neighbourhoods with a Kinsey index of 0, which won't be relevant to our next line of inquiry. It shouldn't come as a major suprise that as an under-represented minority group, most recorded neighbourhoods would have a Kinsey index of `0`. We now turn our attention to measuring their clusteredness.\n",
    "\n",
    "Without trying to quantify this question any further, we can preliminarily visualize the clusteredness of queer populations as we did before using topological diagrams.  We predict from personal experience that queer communities will tend to cluster together, representing crests in their respective diagrams:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 31,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAABJ4AAApGCAYAAAApQQH1AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/P9b71AAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOx9d5wVRdb2U33T5MAMAgImdBAwIGACFUFExbCiYpZlcXVxV11f9zWuYV/z6n7rymLYgBnRlSDBBCqgiIAERUEEDOSByXnm3ttd3x/Vp7r6zp0EM6LueX7rcufe7urq6uoTnnPqlJBSSjAYDAaDwWAwGAwGg8FgMBjtDGtfd4DBYDAYDAaDwWAwGAwGg/HzBBNPDAaDwWAwGAwGg8FgMBiMDgETTwwGg8FgMBgMBoPBYDAYjA4BE08MBoPBYDAYDAaDwWAwGIwOARNPDAaDwWAwGAwGg8FgMBiMDgETTwwGg8FgMBgMBoPBYDAYjA4BE08MBoPBYDAYDAaDwWAwGIwOARNPDAaDwWAwGAwGg8FgMBiMDgETTwwGg8FgMBgMBoPBYDAYjA4BE08MBoPBYDAYDAaDwWAwGIwOARNPDAaD8SPFjBkz0Lt3b/Tu3Rvbtm3b19350WDHjh245557MGLECBx55JF6jN5777193TUGg5GAZcuW6Xd02bJle9XW7bffjt69e2P48OHt1DsGo2WwLmYwGIy9R3Bfd4DB+DkgGo1i3rx5+PDDD7FmzRqUlZWhuroaGRkZ6N69O4488kicccYZOOGEE2BZzPcyGHuKHTt24IILLkBZWVm7trt+/Xq8+eabWLJkCQoLC1FRUYHs7Gx07doVgwcPxtlnn43DDz+8Xa/JAL744gvMmzcPy5Ytw86dO1FWVoZwOIxOnTqhb9++GDx4MEaNGoWsrKx93VUGo0Mwc+ZM3H777QAAy7KwYMECdO3adR/3ivFTRu/evQEAxx13HF566aUWj1+2bBnGjh0LALj++utxww03JD2urKwMM2bMwEcffYSNGzeioqICQghkZ2eje/fuOPzww3HMMcfg5JNPRl5eXqv6unz5clx11VX676lTp2LAgAGtOpfBYPy0wMQTg7GXmDdvHh555BFs37690W/l5eUoLy/H2rVr8eqrr+Kggw7CHXfcgVNPPfWH7+h/EW6//XbMnDkT3bt3xwcffLCvu+ODaeC9+OKLOP744/dxj35aePrpp1FWVoZgMIibbroJxx57LNLS0gAA+++/f5vbq6ysxAMPPIDZs2dDSun7rbi4GMXFxfjyyy/xr3/9C6NHj8Yf//hHZGRktMu9/Ddj+/btuP/++7FgwYJGv8ViMdTU1GDr1q1499138fDDD2PcuHG47rrrkJKSsg96y+hI/P3vf8ekSZMAAF9//fU+7s0Pj5kzZ+rPjuNg9uzZuPbaa/dhj/47sG3bNpx22mkAgIcffhgXXHDBPu7Rjxvvv/8+7rzzTpSXlzf6raioCEVFRfjss8/w6quv4uijj8Z//vOfVrVrzn8AeOONN5h4YjB+pmDiicHYCzz55JOYOHGi/nvIkCEYPnw4evXqhaysLFRUVOC7777DBx98gCVLluD777/H448/zsQTo1W44IIL2BhOwJIlSwAAp512Gq655pq9amv37t24+uqrsWHDBgBAly5dcOGFF2LgwIHIzc1FeXk5Vq5cienTp6OwsBAzZszA+vXr8e9//7vV0VxGY3zxxReYMGECiouLAQDdu3fHOeecg2OOOQZ5eXmIxWLYtWsXlixZgvnz56O8vBzPPPMMzjzzTPTp02cf957RVhx//PHtRig98sgjeOSRR9qlrR8Ddu7cieXLlwMA0tLSUFtbi1mzZjHxxPhRYcWKFfj973+PWCyGQCCAs88+G8OGDUOPHj0QCARQXFyMdevW4aOPPsLq1atb3W59fT3effddAN78f+edd3DXXXchHA531O0wGIx9BCaeGIw9xPTp0zXplJeXh7/97W847rjjGh03ePBgXHHFFdiwYQMefvhhlJaW/tBdZTB+Nti1axcA4OCDD96rduLxOG688UZNOp199tm4//77kZ6e7jtuyJAhuPrqq3H33XfjzTffxLp163DTTTfhhRde4GWze4CioiIf6XTdddfht7/9bVInY9SoUbj99tvx3HPP4R//+McP3VUGo8Mxa9YsSCkRCoVw22234d5778WmTZvwxRdf4Mgjj9zX3WMwAKiMMCKd/v3vf2Pw4MGNjhk6dCiuu+46bN++HZ988kmr2p0/fz5qamoAAHfddRfuvPNOVFRU4IMPPsCZZ57ZrvfAYDD2PdhqZjD2ALt27cL9998PQEVpXnrppaSkk4mCggJMnjwZ48eP/yG6yGD8LBGLxQAAweDexU2ee+45HZk9+eST8dhjjzUinQjp6el47LHHcNJJJwFQNSlefvnlvbr+fyvuueceTTr9/ve/x0033dRsZDsjIwM33HADnn/+eV7iyPjZYdasWQCUDBo9ejQyMzN93zMY+xq7du3Cl19+CQAYMWJEUtLJRPfu3XHRRRe1qu033ngDgKpJdeGFF+qAEn3PYDB+XuCMJwZjD/D888+jrq4OAHDjjTeiV69erTrPsiz84he/aPL3FStW4LXXXsPKlStRVFSESCSCHj16YOjQofjlL3+JTp06JT0vWd2gt956C6+99hq+/vpr1NbWolu3bjjttNNw7bXXIicnp8W+Llq0CHPmzMHq1atRUlIC27bRuXNn9OvXD8OGDcNZZ52F1NTUpOeuXbsWr732GpYtW4bdu3dDSokuXbrghBNOwLhx45rMVpkxYwbuuOMOAKqewP7774/XX38dM2fOxDfffINYLIaePXti1KhRGDduXKPrm7VCAFVHhopsmjCXfQwfPhzbt2/H6NGj8cgjj+DLL7/Eyy+/jOXLl2P37t2IxWK+47du3Yr58+dj+fLl2LBhg3ai8/LycPTRR+OCCy7AKaec0uiaZj0JAj0zE2aticTx6NGjR9JxKy0txQsvvIBFixZh27ZtaGhoQOfOnTFw4EBccsklGDRoUNLzkt3/t99+i2effRZLlizB7t27kZWVhWOOOQbXXHMN+vfv32Q7rYXjOJgzZw7mzp2LdevWoaKiAhkZGTjssMNw5plnYsyYMY2ICHMcCJMmTfI9a+p/axCNRvHCCy8AAMLhMO677z4EAoFmzwkEArj//vtxxhlnIBqN4tlnn8Xll1/uI8DaUr+L5mVzhVwBYOnSpZg5cyZWrFiB4uJiBAIBdO/eHSeddBLGjRuHLl26tHi/e9pGYu2dqqoqvPjii5g/fz62bduGqqqqNtVG2bBhg6651qdPH/zmN79p1XkAmp3D7Tn/165di2effRYrVqxAWVkZunTpguHDh+M3v/mNT/6uWrUKzz//PNasWYOSkhJ06dIFZ5xxBq677romCbKrrroKy5cv1wV/N2/ejMmTJ+Pjjz9GUVER8vPzccIJJ+B3v/sdunfv7hu35557TsukTp064dRTT8WNN97Y5JLPxGs1hZbqKyXO0zVr1uD555/HihUrUFpaitzcXJxwwgmYMGFCk3qwqfci2XudTF6bsq+19fuqqqrwyiuvYMGCBfj+++9RXV2NnJwcHHHEETj//PNxxhlnQAjR5Pnz58/HzJkz8eWXX6K0tBShUAidOnVCly5dcNxxx2H48OE46qijmjy/NVizZg2+/fZbAMB5552HSCSCM844A9OmTcObb76J2267DaFQqMV2SktL8corr2Dx4sXYvHkzKisrkZaWhoMOOgjHHnsszjvvvEYbI7RV7wFqfrz88stYtmwZdu3aBcuysP/++2PIkCEYO3Zsk/oJAGzbxqxZs/Dmm29i/fr1qKioQCQSQV5eHrp164YTTzwRI0aMwKGHHtro3O+++05fd/v27YjFYsjJyUFeXh769u2Lk08+GSNGjGjT0qzEeXbHHXc0movNyWbHcdpkmyQbj9mzZ+Odd97B2rVrUV5ejvT0dBxyyCEYOXIkLrvssh9NPbudO3fqzwceeGC7tbt7926dGXXeeefpf5944gksXrwYpaWlTdq
      "text/plain": [
       "<Figure size 1500x3000 with 15 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "figure, axes = plt.subplots(5,3, figsize=(15, 30))\n",
    "figure.suptitle(\"Concentration of Queer Communities Across the USA\", fontsize=20)\n",
    "figure.subplots_adjust(top=0.95)\n",
    "\n",
    "# Atlanda\n",
    "for i, city in enumerate(CITIES):\n",
    "    ax = axes[(i - i % 3)//3][i % 3]\n",
    "    sns.kdeplot(\n",
    "        data=gb[(gb[\"long\"] > city[\"pos\"][0]) & (gb[\"long\"] < city[\"pos\"][1]) & (gb[\"lat\"] > city[\"pos\"][2]) & (gb[\"lat\"] < city[\"pos\"][3])],\n",
    "        x=\"long\", y=\"lat\", weights=\"SS_Index\", fill=True, cmap=copper_alpha, alpha=.7, ax=ax\n",
    "    )\n",
    "    ax.imshow(\n",
    "        Image.open(f\"images/{'-'.join(city['name'].lower().split())}.png\").convert(\"LA\"),\n",
    "        aspect=\"auto\", extent=city[\"pos\"], alpha=.5\n",
    "    )\n",
    "    ax.set_title(city[\"name\"], fontsize=18)\n",
    "    \n",
    "    ax.get_xaxis().set_ticks([])\n",
    "    ax.get_yaxis().set_ticks([])\n",
    "    ax.set_xlabel(None)\n",
    "    ax.set_ylabel(None)\n",
    "    \n",
    "figure.figure.savefig(\"./images/graphs/12-queer-concentration-nationally.png\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "As expected, in all 15 cities studied, we see a relatively sharp \"peak\" in gay residents in one area. Further, neighbourhoods tend to get less queer radially outwards of this peak.\n",
    "\n",
    "Another interesting observation is that with the exception of Chicago and Miami, all of the queerest communities in each city tend to be clustered around the geographical city centre. This is in line with conventional wisdom that the inner-city tends to be inhabited primarily by poor people and other marginalized groups, while the more privileged groups tend to live outside the city, commuting in for work.\n",
    "\n",
    "Our use of KDE plots to illustrate this phenomenon works, but isn't ideal. It'd debateable (and outside the scope of this project) whether or not queerness is discrete or continuous in a spatial environment. The observations collected in the gaybourhoods dataset certainly are discrete, however, and by visualizing them using a KDE plot we are effectively \"undoing\" that discretization. While this may be a meaningful action to take, it will inevitably introduce error. We deem this acceptable because our goal is to illustrate the clusteredness of queer communities on a macro level.\n",
    "\n",
    "Another way we can quantify the clusteredness of queer communities more numerically is to assess how likely a given neighbourhood is to be surrounded by neighbourhoods with a Kinsey index that is equal to or higher than its own. To do this, we algorithmically select a small group of observations in the vicinity of a given observation and calculate the mean Kinsey Index of that group, assigning it as a new column to the original observation as its \"Neighbourhood Kinsey Index.\"\n",
    "\n",
    "This approach is loosely inspired by differential calculus, in how the differential of a function is considered the \"instantaneous\" rate of change, or what the rate of change of the function would be about a point.\n",
    "\n",
    "While it has proven to be sufficient, there is one notable limitation. Our data fundamentally discrete, and there is a limit to how small the smallest possible neighbourhood about an observation can be."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 32,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAABBEAAAI2CAYAAAAGghDfAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/P9b71AAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOzdd3xT5f4H8E/SlTSdTJFRClJaLIgXKAqIslEUKMpSQFFEZOi9iojzqter6E+WKCAXRYagjLJXsaAyhMooqxQUSoECZXRndOX8/ijnmDRJm92k/bxfL1+WM59zkrR5vud5vl+ZIAgCiIiIiIiIiIiqIK/uBhARERERERGRd2AQgYiIiIiIiIiswiACEREREREREVmFQQQiIiIiIiIisgqDCERERERERERkFQYRiIiIiIiIiMgqDCIQERERERERkVUYRCAiIiIiIiIiqzCIQERERERERERWYRCBiKgK06dPR+vWrTF9+vTqbgp5iHnz5qF169YYPXp0dTeFiIiI4PjfZv5ttx6DCFQl8QMl/rd169Yq9xk/frzRPleuXHFDSz1DQkKCdN3R0dFITU2tdHtx24SEBJN1Yue1Z8+elR7j66+/lo7zyiuvoLi4GABw5cqVSo9fW4j3YN68eRa3ycnJwdChQ9G6dWvExsZi27ZtbmxhzWD43q9Nn3kiVzL8G3zfffchKyvL4raGv/MPHTpksn706NFVfkHW6/V47733pOP897//hSAIAIBDhw5VevyaKDs7GwsWLMDIkSPRuXNn3HvvvejcuTMGDBiACRMmYNGiRTh69Gh1N9Nj1bT7Z/h3zvC/2NhYdO3aFc8//zzWrFmDkpKSamnfvHnzMG/ePP4NJpfzre4GkPdJSEjAgAEDLK7PysrCvn373NgizyUIAmbOnIlvv/3WZcf/9NNPsXTpUgDAsGHD8OGHH0IuZ3zQFllZWXj++efx119/QalUYt68eXjooYek9fXr10dkZCTq169fja0kotpOp9Ph66+/xkcffeSS4xcXF2PatGnYvn07AOCVV17BpEmTXHIub/D777/jn//8J3Jzc6VlgYGBKC0txV9//YW//voLe/bsAQCcPXu2mlrpuWr6/QsPD4ePjw8AQKPR4NatW7h16xb279+PH3/8Ed999x1CQ0Pd2qavvvoKABAXF4cmTZq49dw1QXh4OCIjI9GoUaPqborHYxCBrBYeHo6ioiIcOHAA169fx1133WV2u40bN6KsrAyNGzdGZmamm1vpefbt24fff/8dDz74oFOPW1ZWhnfffVcaYTBu3Di88cYbTj1HbZCRkYGxY8ciMzMTISEh+Oabb/CPf/zDaJvXX38dr7/+ejW1kIjob+vWrcPYsWMRGRnp1ONqtVpMnjwZ+/btg0wmw3vvvYdnnnnGqefwJlevXsXEiROh0WjQuHFjTJo0CX369EFISAgAoKCgAMePH8fu3buxZcuWam6t56kN92/t2rVGHfXLly9jzpw52LJlC06dOoX33nsPX375ZTW2kGw1atQojBo1qrqb4RX4uJKsFhgYiH79+kGv11c6NH7dunUAgCFDhriraR6rR48eAICZM2dKw0Gdobi4GK+++qr0OkydOpUBBDukpaXhmWeeQWZmJurXr48VK1aYBBCIiDxBo0aN0Lp1a5SWlmL27NlOPXZ+fj7Gjh2Lffv2wc/PD//3f/9XqwMIAPDTTz9Bo9HAz88PK1aswJNPPil1gAEgODgY3bp1w/vvv4/ffvutGlvqmWrj/WvatCm++OILdOzYEQCQmJiImzdvVnOriFyDIxHIJkOGDMH69euxfv16TJw40WT94cOHcfHiRTRt2lT6JVqVX375BevWrUNKSgpycnKgVCoRFRWFAQMG4KmnnoK/v7/JPnl5eUhMTMS+fftw4cIFZGVlQavVol69evjHP/6B0aNHo3379mbPN2/ePHz11VeIi4vD8uXL8fvvv2PJkiU4ceIE1Go1mjRpggEDBuDFF19EQECATfenotdeew2//vorTp48iR07duDRRx916HgAoFarMWnSJPz+++/w8fHBhx9+iKFDhzp8XHN69uyJzMxMfPrpp+jbty+++eYbJCYm4tq1a1AqlejQoQNeeukl3HfffZUe58iRI1i1ahWOHDmCW7duwd/fH5GRkejbty+eeeYZqFQqk32mT5+O9evXIz4+Hp9++inWrl2LhIQEXLhwAbm5ufj0008dClQdPXoUL730EvLz89GkSRMsWbIEzZo1M7utYVtmzJhhtG706NFITk7G5MmTMXnyZKxZswZr1qzB+fPnIQgCoqKi8PTTT2PQoEFmj11aWop169Zhy5YtOHfuHAoLCxEUFISwsDDExMSga9euFl/fK1euYOnSpThw4ACuXr0KvV6PRo0aoVu3bnj++edx9913G23/r3/9C9u2bUP37t3xv//9z+K9ycjIQL9+/SAIApYtW4bOnTtXdiudxpbfBbdv30b37t1RWlqK+fPno1evXhaPO3fuXMyfPx/NmjXDrl27TNbb8/4kcje5XI7XX38d48ePx86dO3HixAm0a9fO4ePevHkTL7zwAs6ePQuFQoEvv/wSDz/8sBNa/Lddu3Zh8uTJ8PPzw969exEeHm5x22eeeQaHDx/Gk08+iU8++URafv36dXz33XfYv38/MjMzUVpairCwMDRo0AAdO3bE448/7pT7ITpz5gwAICYmxuR3aUUKhcLs8pSUFOzatQvHjh3DtWvXcOvWLQQEBKBFixbo3bt3pb9fWrduDQBYtmwZoqKisHDhQiQlJeHGjRuoW7cuevTogcmTJ6NOnToAgMzMTCxatAh79+7FzZs3UbduXfTv3x+TJ09GUFCQxbZnZ2dj6dKl+PXXX3H58mUUFxejQYMG6Ny5M8aOHYtWrVpVea/Mccb9A8pzdGzZsgWbN2/G6dOnkZ+fj6CgILRp0wZDhgzBgAEDIJPJpO0zMjIQHx8PtVqN0aNH49133zU5ZmFhIYYMGYKMjAx069YNixcvNjqGI2QyGQYOHIjDhw9DEAScOnUKPXr0gFarRVJSEn777TecPXsWWVlZKCwsRFhYGNq1a4fhw4db/NwlJCTgrbfeQuPGjbF7924cPHgQy5Ytw4kTJ3D79m3p+8X69eulfcaMGWN0DHFfoDy3ibj+7NmzyMjIwMKFC3HgwAHcvn0bderUQffu3TFlyhQ0bNjQ4rUWFxdjzZo12LFjB86dOwe1Wo3Q0FC0a9cOI0aMsHg9hu/te+65B4sWLcIvv/yC69evQ6fTGU1t2bt3L3766SecOHEC2dnZ8Pf3R3h4OCIiItC1a1c8+eSTCAsLs9hGW7/jV+wjGKr4vfTHH3/EunXrcOHCBem6nn76aTzxxBMW27Nt2zYkJCQgNTUVeXl5UCqVqFOnDlq0aIGHHnoITz31lMN9D3dhEIFs0qlTJzRr1gyXLl3CH3/8gU6dOhmtF5+Mx8fHV/kLWafTYdq0adi5c6e0LCgoCAUFBTh8+DAOHz6MjRs3YtGiRSZzypYtWybN+/Lx8ZH+QF69ehVXr17F1q1b8fbbb5v8Eq1o8eLF+OKLLwCUR8VLSkpw4cIFzJs3D8nJyViyZIk0380eUVFRGDRoENavX485c+agT58+8PW1/2OXm5uLF198ESdOnICfnx9mzpyJfv362X08a+Xn5+Opp55Ceno6/Pz8EBAQgNzcXCQlJWHPnj34z3/+g6eeespkP71ej08++cToF3FgYCC0Wi1OnjyJkydPIiEhAd9++y0aN25s9tyCIODVV1/Fzp07IZfLERwc7HDOh71792LKlCnQarWIiorCt99+iwYNGjh0zLKyMkyaNAlJSUnw9fWFQqGAWq1GSkoKUlJSkJGRgVdeecVkn/Hjx2P//v3SsuDgYGg0GuTm5uLixYvYvn272SDCpk2b8M4770hJNP39/SGXy5Geno709HQkJCTgyy+/RLdu3aR9RowYgW3btmHfvn24evWqxS92a9asgSAIaN68uVsCCPb8Lqhbty66deuGX375BRs3brQYRBAEAZs2bQIAk0COM96fRO708MMPIy4uDsn
      "text/plain": [
       "<Figure size 1200x600 with 2 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "fig, axes = plt.subplots(1, 2, figsize=(12, 6))\n",
    "mean_kinsey = gb.groupby(\"kinsey\").mean().reset_index()\n",
    "\n",
    "mean_nk_bar = sns.barplot(data=mean_kinsey,\n",
    "    x=\"kinsey\", y=\"neighbourhood_kinsey\", ax=axes[0])\n",
    "\n",
    "mean_nk_scatter = sns.scatterplot(data=gb, x=\"SS_Index\", y=\"neighbourhood_kinsey\", marker=\"x\", ax=axes[1])\n",
    "\n",
    "mean_nk_bar.set_title(\"Mean NKI per Kinsey Level\", fontsize=18)\n",
    "mean_nk_bar.set_xlabel(\"Kinsey Index\", fontsize=15)\n",
    "mean_nk_bar.set_ylabel(\"Mean NKI\", fontsize=15)\n",
    "\n",
    "mean_nk_scatter.set(ylabel=None, xlim=(0,30), ylim=(0,3))\n",
    "mean_nk_scatter.set_xlabel(\"Same-Sex Partnership Index\", fontsize=15)\n",
    "_=mean_nk_scatter.set_title(\"NKI vs. Same-Sex Partnerships\", fontsize=18)\n",
    "\n",
    "fig.savefig(\"./images/graphs/13-neighbourhood-kinsey-comparison.png\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The graph on the left shows that neighbourhoods with a higher Kinsey index tend to have a higher Neighbourhood Kinsey Index on average, which strongly suggests that queer neighbourhoods tend to cluster together. Interestingly, this is less the case with neighbourhoods with a Kinsey index of 6. This makes sense when we consider that the number of neighbourhoods that fall into that category is small, and we've defined it such that nothing can be higher.\n",
    "\n",
    "Neighbourhoods with a KI of 1 are about three times more likely to be surrounded by other queer neighbourhoods than neighbourhoods with a KI of 0. This further illustrates the fact that queer communities tend to cluster together.\n",
    "\n",
    "The graph on the right illustrates the same fact by measuring NKI versus the Same-Sex Partnership Index. While the results are consistent, this graph reveals that there's quite a bit of variance.\n",
    "\n",
    "All together, it's very clear that there is a strong coorelation between geological strata and wanting/having same sex partners."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Is there a correlation between political alignment & living in neighbourhoods with large quantities of LGBT people?\n",
    "\n",
    "To tackle this question, we need to unify the gaybourhoods dataset with the election data from the 2012 US presidential election. To accomplish this, we used our `nearest_neighbour` algorithm to locate the county closest to each observation in the gaybourhoods dataset and create an additional column connecting the gaybourhoods to their county's election outcome."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 33,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAABJ4AAApGCAYAAAApQQH1AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/P9b71AAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOydd3gUVdvG79nZvukJCYTQS+gdkY+uNFGQF1BBBLvYe1fAXnjt4GtFRQREkK4gSgdp0kFCDzWk1+07M98fZ8+ZmWRTSSh6ftfFRXZ3dvbMzCnPuZ/nPEdQFEUBh8PhcDgcDofD4XA4HA6HU80YLnUBOBwOh8PhcDgcDofD4XA4/0y48MThcDgcDofD4XA4HA6Hw6kRuPDE4XA4HA6Hw+FwOBwOh8OpEbjwxOFwOBwOh8PhcDgcDofDqRG48MThcDgcDofD4XA4HA6Hw6kRuPDE4XA4HA6Hw+FwOBwOh8OpEbjwxOFwOBwOh8PhcDgcDofDqRG48MThcDgcDofD4XA4HA6Hw6kRuPDE4XA4HA6Hw+FwOBwOh8OpEbjwxOFw/lUkJyezf2fOnLnUxeFUkAULFrDnNm7cuEtdHA6Hw+FwOJwKw+0Yzr8d46UuAKd6kSQJ27Ztw5o1a7Bz505kZWUhOzsbNpsNsbGxqF+/Pnr37o1rr70WtWvXvtTF5YQgOTm51M/MZjPCw8MRHh6O6OhoJCcno02bNujcuTMaN258EUvJ4Vz5jBs3Dtu2bQv5mSiKCA8PR1hYGCIiItC0aVO0bt0a7dq1Q8eOHSEIwkUuLYdTOtdccw3Onj1bred8+OGH8cgjj1TrOTlXHlOnTsW0adPY686dO2P27NmXsEScmuT06dMYMGAAFEVh761cuRINGjS4hKXicDj/BLjw9A9i3bp1mDJlCo4ePVriM5/Ph/z8fBw/fhxr167FW2+9hTFjxuDhhx9GVFTUxS8sp0r4fD5kZ2cjOzsbqamp2LVrF/usc+fOGD16NG644QYYDDyY8ULQGtr/+c9/8M4771ziEnEuNpIkIS8vD3l5eQCAv//+G0uWLAEA1K9fHzfffDNGjx6N8PDwS1hKzpXG1q1bMX78eABA3bp1sXr16ktcIs6VgFYkf/vttzFixIiL8ruKomDRokW693bs2IFTp06hfv36F6UMnIvL4sWLdaITACxatAiPPfbYJSoRh8P5p8CFp38AiqLgzTffxMyZM3XvG41GJCUlIT4+Hk6nE2lpacjJyQEABAIBzJw5EytWrMD06dPLjLLhXDqaN2+O+Ph43XuFhYUoKCjA2bNn4fP52Ps7duzAjh07MG/ePLz77rtITEy82MXlcK5Y4uPj0bx5c917LpcLBQUFSE9PR2FhIXv/1KlTeO+99zBnzhy8++676Nq168UuLoejo2vXrmjUqFGZx2zfvh1erxdA6PpeHC4scP7666+QS9IXLVqERx999BKUiFOThBIaASJGPfroozzSl8PhXBBceLrCURQFzzzzDJYuXcrei4qKwiOPPIIhQ4YgJiZGd+zu3bvxzTffYOXKlQCAzMxM3HbbbZg+fTratWt30cvPKZs777yzVM+mz+fDwYMHsXLlSvz0008oKCgAAGzbtg3Dhw/HnDlz0KRJk4tZ3CuCQ4cOXeoicC5DevToUWpkm6IoOH78ODZu3IgffvgBp06dAgCcPXsW48ePx5QpUzB06NCLWVwOR8e7775b7jHa5Xhl1XcOh7Jw4UL2d+/evbF+/XoARHh65JFHuBDxD2PHjh04ffo0ALCldSdPnsTZs2exbds2dOvW7VIW74pnxIgRFy1akcO5HOHrca5wZsyYoROd2rVrh+XLl+O2227TiU4AIAgCOnbsiKlTp+Ldd9+FKIoAgIKCAjz++OMoKiq6qGXnXBhmsxnt27fHM888gzVr1mDIkCHss/z8fNx///3Izc29hCXkcP4ZCIKAJk2a4Pbbb8dvv/2Gp59+mvWfsizjpZdewu7duy9tITkcDqcacbvdWLFiBQDSB06aNAlJSUkAiOi+ffv2S1k8Tg2wYMEC9vfQoUN1DpVQkVAcDodTGbjwdAVz7NgxvP/+++x148aNMX369BKCUyiGDx+OV155hb0+e/Ys3nzzzZooJuciEBYWhg8//BBjxoxh7506dQpTp069hKXicP55GAwG3HvvvXjvvffYe16vFxMnTiyRF4PD4XCuVFauXAmn0wkA6NixI+rVq4cbbriBfa6NhuJc+WiFRgAYNmwYhg0bxl6vWLECLpfrUhSNw+H8Q+BL7a5gpk+fznL8CIKA119/HRERERX+/s0334zffvsNGzduBEDWcD/88MOoW7duiWOrkmy5Kskwz549i6VLl2LDhg04ffo0cnJy4HA4ULt2bfzf//0fRo4ciaZNm1b4GgGyJG358uVYs2YN/v77b2RlZUGSJMTGxqJ9+/YYNGgQBg0aVG7I+PPPP88MLbrbjyzLWL16NZYsWYKUlBRkZmbC5XJdst2AXnzxRezYsQOHDx8GAMyfPx8PPPAAatWqVe539+7di19//RVbtmxhOW0iIyPRoEED9OnTBzfddFO5omZpiXPPnDmDuXPnYu3atTh37hxkWUZCQgK6d++OO+64o8K7pZw9exYLFizA1q1bcfz4cba80OFwIDExEU2bNkWnTp1wzTXXICEhIeQ5tPnMVq1axTy4QOhdzhYuXFiqgU2/P3r0aJbo/Z577sEzzzxToevxeDzo2bMnyx/04Ycf6iLXqsLZs2exYcMGbN++HYcPH0ZaWhrcbjccDgfi4uLQsWNHDB48GL169arQ+ULdr6KiIixYsAC//vorTp06hcLCQsTExKBTp064+eab0b1790qVefXq1ViwYAH27duHnJwcREdHo2HDhszjarVaK3W+i8GQIUOwZcsWzJ07FwBw+PBhrFq1Cv379y/3u9XVz4V6NtnZ2Zg3bx5WrlyJc+fOoaioCImJiejTpw/uu+++En2Bz+fDwoULsWzZMhw+fBhOpxMxMTHo2rUr7rjjDrRt27ZS92Xr1q345Zdf8NdffyEzMxNerxexsbFo2rQp+vXrh+HDh8Nut1fqnIFAACtXrsS6deuwZ88e5OTkwOl0wuFwoF69emjbti369u2LXr16sUg0LdolZt9//z26deuGoqIiLFu2DCtWrMDJkyeRmZkJv9/PPtdy4MABbNq0CTt37sTRo0eRnZ0Nv9+PiIgIJCYmomvXrhgxYgSaNWtW6jUU3xkMIPWgtPyKFzOBdHFC3a+yOHPmDK699lr2urTlzKFsAp/PhxUrVuCXX37B0aNHWZ3RXv/F6INOnDiBRYsWYePGjUhLS0NBQQGioqKQlJSEXr16YcSIEahTp06556nuexeqfrzwwgt44YUXSrxf3cnqtREuVIC48cYb8fnnnwMAfvvtN0yaNAk2m61S573U7Rmo3n6qqKgIS5Yswdq1a3Ho0CHk5eUhEAjAZrMhPj4eDRs2RPv27dG3b98y86n6fD78+uuv+OOPP3Dw4EHk5OTA6/XCarUiLi4ODRo0QLt27dCrVy+0b9++2pc5/v7770xobN++PbPLOnTogN27d8PlcmHlypUYPnx4hc53OY5POTk5WLJkCdavX48TJ04gOzsbZrMZ8fHx6NatG2688UZ06NCh3PNUtQ4uWLCAtd2rrrqqRG7eUCiKgvXr12PVqlXYtWsXMjMzUVhYCJvNhrp166JVq1Zst3Kz2RzyHH6/H9u2bcOff/6Jffv24cSJE8jPzwcAREZGokmTJrj66qtx8803Vyh4gcOpKlx4ukLJycnRLbHr3bs3unTpUunzPPHEE0x4kiQJM2fOxPPPP19t5awogUAAn3zyCb799ltdwmwAbGeplJQUzJgxA+PGjcOzzz4b0iApzsaNG/Hqq6+ynCxazp49i7Nnz+LXX39FmzZt8NFHH6FevXoVLnNWVhaefvppbN68ucLfqWnMZjPuuecePPvsswBIJMaiRYtw7733lvqdnJwcvPLKK/jtt99KfJaVlYWsrCzs2LEDX3/9NV566aUKGx2URYsW4ZVXXoHb7da9f+LECZw4cQL
      "text/plain": [
       "<Figure size 1500x3000 with 15 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "figure, axes = plt.subplots(5,3, figsize=(15, 30))\n",
    "figure.suptitle(\"Queer Density and Democrat Turnout Across America\", fontsize=25)\n",
    "figure.subplots_adjust(top=0.95)\n",
    "\n",
    "for i, city in enumerate(CITIES):\n",
    "    ax = axes[(i - i % 3)//3][i % 3]\n",
    "    sns.kdeplot(\n",
    "        data=select_region(gb, *city[\"pos\"]),\n",
    "        x=\"long\", y=\"lat\", weights=\"SS_Index\", color=\"orange\", ax=ax\n",
    "    )\n",
    "    sns.kdeplot(\n",
    "        data=select_region(gb, *city[\"pos\"]),\n",
    "        x=\"long\", y=\"lat\", weights=\"percent_democrat\", color=\"blue\", ax=ax\n",
    "    )\n",
    "    \n",
    "    ax.imshow(\n",
    "       Image.open(f\"images/{'-'.join(city['name'].lower().split())}.png\").convert(\"LA\"),\n",
    "       aspect=\"auto\", extent=city[\"pos\"], alpha=.4\n",
    "    )\n",
    "    ax.set_title(city[\"name\"], fontsize=20)\n",
    "    \n",
    "    ax.get_xaxis().set_ticks([])\n",
    "    ax.get_yaxis().set_ticks([])\n",
    "    ax.set_xlabel(None)\n",
    "    ax.set_ylabel(None)\n",
    "\n",
    "orange_patch = mpatches.Patch(color='orange', label='SS Index')\n",
    "blue_patch = mpatches.Patch(color='blue', label='Democrat Turnout')\n",
    "_ = axes[0][0].legend(handles=[orange_patch, blue_patch])\n",
    "\n",
    "figure.savefig(\"./images/graphs/14-queer-vs-democrat-density.png\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "With the exception of Washington DC, Miami and Chicago, the peaks in the respective topological maps almost completely overlap, implying a relationship between queerness in a given community and voter turnout for the Democratic Party. This aligns with conventional wisdom once again in that urban queer people tend to support liberals more often than conservatives in elections.\n",
    "\n",
    "Similar to before, using a weighted KDE plot to illustrate election results can be somewhat misleading in certain cases, as some cities contain fewer counties than others.\n",
    "\n",
    "Alternatively, we can use our Kinsey index once again to assess this relationship more categorically:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 35,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkAAAAHTCAYAAADPgKdGAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/P9b71AAAACXBIWXMAAA9hAAAPYQGoP6dpAABiX0lEQVR4nO3dd1xT9/4/8FfCkC0ouHCPAAqigBsHbouKExc4qlWv29a2or0darW1rrrXFQdKHeBCRdS6FQdatYqjblDRCsgQWTm/P/wlX0OCQggEOK/n4+GjzTknJ+/zycnhlfP5nBOJIAgCiIiIiEREqu8CiIiIiIoaAxARERGJDgMQERERiQ4DEBEREYkOAxARERGJDgMQERERiQ4DEBEREYkOAxARERGJDgMQERERiQ4DEBGJXmhoKBwcHBAaGqrvUvLkwoULcHBwwLJly/RdSqni4OAAf39/fZdBRcRQ3wVQ8eXg4KDy2MjICBYWFqhcuTLq16+Pzp07w9PTEwYGBnqqsPTw9/fHxYsXcefOnTwtHxoaioCAgHy9Rl7XXZJNnz4du3fvxrFjx1C1alV9l5MnFy5cwNChQ9G0aVNs2bJFbf6JEycwZcoUyOVyLFq0CB07dtRDlfrVvn17xMbGlqj3lYo/BiD6pAkTJgAAsrOzkZycjHv37mHv3r3YtWsXnJ2dsWDBAtSqVUvPVYqLk5OT8n1RiI2Nxe7du2Fvb4/evXvrqTLSpZCQEHz//fcwNzfHqlWr4O7uDgBo2LAhDh48CBsbGz1XSFRyMQDRJ02cOFFt2r///ovZs2cjPDwcI0aMQEhICMqXL6+H6sTJyckJTk5OKtMuXLigDECa3jMqWdasWYNFixahUqVKWL9+PerVq6ecZ2pqijp16uixOqKSj2OASCu2trZYvHgxmjZtiufPn2P16tVqyyQmJmLhwoXo1q0bGjZsCHd3dwwbNgxnzpxRW/bDMRhnz57F4MGD0bhxYzRv3hwBAQFISkoCANy6dQtjxoxBkyZN0LhxY4wdOxYxMTEaa3z06BG++eYbtG7dGs7OzvD09MQ333yDR48eaVw+OzsbwcHBGDhwINzd3dGwYUN06tQJM2fOVHnOsmXL4ODggAsXLmD//v3o378/GjdujPbt26tsz8SJE9GhQwc0bNgQbm5uGDhwIPbu3avymjExMXBwcMDFixcBvO92VPzT1ViET40Xad++vUrtivoV78epU6fg7+8Pd3d3Zbfoh/MjIyPh7++Pxo0bw83NDaNHj8b9+/c1vtbLly/x008/oX379nB2dkbz5s0xYcIE/P3332rLftjOOSnabfr06cppDg4O2L17NwCgQ4cOynbMuW2fcuLECQwcOBCNGjVCkyZNMGnSJLV95ssvv1R533I6fPgwHBwcMGvWrHy9NgAIgoA5c+Zg0aJFqFu3Lv744w+V8APk/p76+/vDwcEBWVlZWL16NTp37gxnZ2e0bdsWv/32GzIyMtRe7/Llyxg7dizatGkDZ2dntGrVCr6+vli+fLnasmlpaVizZg18fHzQqFEjNG7cGAMGDEBYWJjKcqdPn4aDg0Ou3bQZGRlo1qwZmjVrprGmvFJ8TuLj4/Hf//4Xnp6ecHZ2hre3N0JCQnJ97RUrVqBjx45wdnZG+/btsXjx4o/WkZWVha1bt8LX1xdubm5wdXVFr169EBQUBLlcrlzu6dOn8PDwQNOmTREbG6uyjrdv36Jbt25wcnLSuE9T0eMZINKaVCrFuHHjcPHiRRw4cAAzZsyARCIB8L47xt/fH7GxsfDw8EDr1q2RlpaG48ePY9SoUZg1axZ8fX3V1vnnn3/ixIkTaNeuHQYOHIirV68iNDQUMTEx+OqrrzB8+HC4u7ujX79+uHv3Lo4fP46YmBjs27cPUun/5fnr169jxIgRSE1NRfv27VG3bl08ePAA+/btw7FjxxAYGIiGDRsql8/IyMDYsWNx9uxZVK5cGd27d4eFhQViY2Nx9OhRuLu7o2bNmiq1BgYG4uzZs/Dy8kKzZs2QnJysnPfjjz+ibt26aNKkCezs7JCYmIiTJ0/im2++wcOHDzFlyhQAgJWVFSZMmIDdu3cjNjZWpVvL3t5eF29TgRw+fBinT59GmzZtMHDgQDx79kxl/okTJ3Ds2DG0bt0aAwcOxP3793Hy5EncuHEDBw4cQLly5ZTLPn36FIMHD8bLly/RvHlzeHt74/nz5wgPD8eJEyewbNkyeHl5aV3rhAkTcPToUdy+fRtDhw6FlZUVAMDS0jLP64iIiMDp06fRsWNHNG3aFNHR0Th8+DAuXLiA4OBg1K5dGwAwaNAgHDhwANu3b0fTpk3V1rN9+3YAwMCBA/O1DRkZGfj2229x8OBBNG7cGKtXr4a1tXW+1gEAX331FaKiotC6dWu0bdsWp06dwvr16xEfH4958+Yplzt16hTGjBkDCwsLtG/fHhUrVkRiYiIePHiAbdu2qeyPSUlJGDZsGG7duoUGDRqgb9++kMvlOHPmDL766ivcu3cPU6dOBQB4enqievXqOHToEGbMmKH2Hhw+fBiJiYn4/PPPYWxsnO/t+1BSUhIGDRoEY2NjdOnSBRkZGQgPD8eMGTMglUpVuoMFQcCUKVNw7NgxVK9eHX5+fsjMzERISAju3r2rcf2ZmZkYO3Yszpw5g1q1aqF79+4oU6YMLly4gNmzZ+PatWv47bffAADVqlXDnDlzMHnyZHz11VcICgqCoeH7P7M//fQTHjx4gIkTJ6JZs2YF2mbSEYEoFzKZTJDJZB9dJj09Xahfv74gk8mEJ0+eKKf7+fkJDg4OQlhYmMryb968EXr27Cm4uLgIr169Uk4PCQkRZDKZ4OTkJFy4cEE5PTs7Wxg+fLggk8mEJk2aCHv37lVZX0BAgCCTyYQjR44op8nlcqFr166CTCZTW/7AgQOCTCYTunTpImRnZyunL1y4UJDJZMKYMWOE9PR0tW18/fq18vHSpUsFmUwmuLq6Cjdv3tTYLo8fP9bYVkOHDhXq168vvHjxQmWen5/fJ9v6UyIjIwWZTCb4+flpnL506VKNz/Py8hK8vLxUpineDwcHB+HkyZNqz/nw/Tp37pzKvAULFggymUxYu3atyvTPP/9ckMlkwsqVK1WmR0VFCU5OTkLTpk2FlJQU5XRFO0dGRqq9/tOnTwWZTCZ8++23KtO//fZbQSaTCU+fPtW4rblRbI9MJhP+/PNPlXkbN24UZDKZMHToUJXp3t7egrOzsxAfH68y/cmTJ4KDg4MwYMCAPL224v3p3bu3cl8fM2aMkJaW9snn5HxPFftR7969hYSEBOX01NRUoWPHjoKjo6Pw8uVL5fQJEyYIMplMiI6OVnuND/d5Qfi/ts35vr579074/PPPBQcHB+HWrVvK6evXrxdkMpmwZcsWtXUr6nzw4EGu2/ghLy8vje+r4j2bMWOGkJWVpZx+7949wcnJSejWrZvK8vv27RNkMpng6+srvHv3Tjk9ISFB6NChg8bPj2I/nDVrlsprZGVlaTz+CIIg/PDDD4JMJhMWLFggCIIghIaGCjKZTPD391c57pB+sQuMCsTY2Fj5DTUhIQEAcPv2bVy8eBGdO3eGt7e3yvJWVlaYOHEi0tPTcfjwYbX1eXt7q3yjlkql8PHxAQDUq1cPPXv2VFm+V69eytdUuHLlCh48eIDGjRurLf/ZZ5/B3d0dDx8+RFRUFID3XV/btm2DiYkJfvrpJ7VvpMbGxipnMhR8fX1Rv359je1SvXp1tWnGxsYYMmQIsrKycP78eY3PK246dOiANm3a5Dr/s88+Q4sWLVSmKc7s3bhxQzntxYsXOHPmDKpUqYJRo0apLO/m5gZvb28kJibiyJEjOqw+/5o3b652FsrPzw/Vq1dHZGSkSrfGoEGDkJGRoex2U9ixYwcEQcj32Z+bN2/i3LlzqFWrFpYvXw4TExOtt2PatGkqZ47MzMzQo0cPyOVyjd2NZcqUUZv24T6fkJCAffv2wdnZGV988YXac7/++msIgoD9+/crp/fp0wdlypRRng1TePDgAS5evIhmzZrp5OIJU1NTBAQEqFyNWrd
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "gb_pol = gb.groupby(\"kinsey\").mean().reset_index()\n",
    "\n",
    "gb_pol[\"percent_democrat\"] = gb_pol[\"percent_democrat\"] * 100\n",
    "\n",
    "kinsey_democrat_plot = sns.barplot(data=gb_pol, orient=\"h\", x=\"percent_democrat\", y=\"kinsey\", palette=\"Pastel2\")\n",
    "_ = kinsey_democrat_plot.set(\n",
    "    xlabel=\"Percentage of the Population who Voted Democrat in 2012\",\n",
    "    ylabel=\"Kinsey Index\",\n",
    "    title=\"Democrat Turnout by Kinsey Index\",\n",
    "    xlim=(50,80)\n",
    ")\n",
    "\n",
    "kinsey_democrat_plot.figure.savefig(\"./images/graphs/15-kinsey-democrat-plot.png\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "While this once again reveals that queerer communities tend to vote blue more often, it is worth noting here that our data only ranges between 60% and 75%. That is a meaningful different (one so meaningful that it ultimately won an election) but it is useful to keep in mind how substantial this data is.\n",
    "\n",
    "Of course, in the 2012 presidential election, the Democrats won, so it's not suprising that across the nation, people tended to vote for the Democrats more than they did for the Republicans. Nevertheless, it's still clear that communities with a larger queer population tend to vote more blue than others.\n",
    "\n",
    "Nonetheless, it's clear that queer communities tend to lean democrat across the United States"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.10.10"
  },
  "vscode": {
   "interpreter": {
    "hash": "b2baa059f790e7ad780c83135aaea020c73a7a7a6921010b599b8b664933698d"
   }
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}