Skip to content

Instantly share code, notes, and snippets.

@NaimKabir
Last active November 14, 2020 13:59
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save NaimKabir/a796ec6d99efbfca865a5e5e465d44eb to your computer and use it in GitHub Desktop.
Save NaimKabir/a796ec6d99efbfca865a5e5e465d44eb to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Simple process: let's say you have two populations. Split-ticket voters, and voters who vote using the straight ticket.\n",
"\n",
"We assume Republican straight-ticket voters always vote for Trump.\n",
"\n",
"Split-ticket voters instead have some small chance of voting for Trump, $p({trump})$."
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": [
"p_trump = 0.30"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Great now let's simulate a bunch of precincts with different fractions of R voters.\n",
"\n",
"We'll assume that all R-voters voted Trump."
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": [
"import numpy as np\n",
"num_precincts = 1000 # let's deal with 1000 precincts.\n",
"# modeling the fraction of R voters in a precint, randomly from a normal distribution close to 0-mean,\n",
"# to model Wayne county.\n",
"r_percentages = abs(np.random.normal(0.04, 0.07, num_precincts))"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Now let's simulate some per-precinct votes and collect % of Trump votes from non-Republican candidates.\n",
"\n",
"This can be sampled a a bunch of weighted coinflips with the probability $p({trump})$--i.e, sampling from a binomial distribution."
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [],
"source": [
"split_trump_percentages = np.zeros(num_precincts)\n",
"# simulating\n",
"for i, percent in enumerate(r_percentages):\n",
" # to closely model Wayne county, let's make pop size variable\n",
" non_r_pop_size = abs(np.random.normal(10, 20))\n",
" split_votes = np.random.binomial(non_r_pop_size, p_trump)\n",
" split_trump_percentages[i] = split_votes / non_r_pop_size\n"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAEdCAYAAADn46tbAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nOydeVxU9f7/X2cWdghBiDStRBOlzCUXcgHNVBTcsvyppWVZlqmZmdV1a3W9Xzeu3muZpmU39RJduGq5kiBWam6IYmbuCKKIwDDMcn5/0Blnzpx1Zs4szOf5ePh4yJmZc97nM2c+78/nvVI0TdMgEAgEAsEKlacFIBAIBIL3QZQDgUAgEOwgyoFAIBAIdhDlQCAQCAQ7iHIgEAgEgh1EORAIBALBDqIcCAQCgWCHJOVQVVUFADh69CiysrJgMBgUFYpAIBAInoUSS4Jbvnw5Ll68iOnTp2PEiBFo2bIlmjVrhk8++cRdMhIIBALBzYjuHHJzc/Hxxx/jxx9/xKBBg7BhwwacPn3aHbIRCAQCwUNIMisFBwfjwIED6NatGwCgrq5OUaEIBAKB4FlElUOjRo0wb948nDx5Ek888QSWLFmC2NhYd8hGIBAIBA8hqhwWLlyI2NhY/Otf/0JwcDAoisLChQvdIRvBz3n55Zexfv16y9/nz59H69at8fe//91yrLy8HI888gju3LnjAQmB559/Hn369MGQIUMwZMgQpKeno3///sjKynLqvH369MGJEydw4sQJTJkyxeHzZGRkYNeuXQDq/YfOykXwHzRib2jcuDGGDx+OM2fOwGQyYfTo0WjcuLE7ZCP4Ob169cLBgwfxwgsvAAD27t2L3r17Y8+ePZg+fToA4ODBg+jQoQPCw8M9Juc777yDAQMGWP4+ceIERo0ahb59+yIsLMypcz/66KNYsWKFw5//+eef0bJlSwDA1KlTnZKF4F+IKod9+/Zh3rx5UKlU+Pe//42BAwdi8eLF6Nu3rzvkI/gxvXr1QkZGBsxmM1QqFfbu3Ytp06bhrbfewqVLl9CsWTMUFBQgJSUFNTU1mDdvHv7880/cvn0boaGhWLJkCYKCgjBo0CD89NNPCA8PB03TGDBgAJYvX45PPvkE7du3x5EjR3Dt2jV06tQJCxcuhEqlwpEjR7BkyRLodDpQFIXJkyejd+/ekuS+dOkSQkJCEBAQAADYs2cPVq9eDYPBgKCgIMycORMdOnTAypUrcfbsWdy4cQPl5eVISEjAJ598YqNQfv75Z3z00UfIyclBdXU1Pv74Yxw5cgRqtRp9+/bFtGnT8Oeff+LDDz9ETU0NSktLkZCQgGXLlmHr1q04efIkFi1aBLVajd27d6NVq1Z46aWXcOjQISxatAg6nQ5arRZvvvkmevXqhczMTOzcuRMqlQoXLlyAVqvFwoUL8fDDDyvyHRO8GFqEp59+mr5+/To9ZMgQmqZp+tSpU3R6errYxwgEl9C3b1/61KlTdEVFBd29e3faZDLRs2fPptetW0fTNE336dOH/v333+nt27fTH330keVzs2fPpj/88EOapmn6tddeo7/66iuapmn6wIED9LPPPkvTNE0/99xz9JQpU2iTyUTfuXOH7tGjB11QUEBXVFTQ/fr1oy9dukTTNE2XlJTQvXr1oq9cuWIn33PPPUf37t2bHjx4MN27d286KSmJnjZtGl1YWEjTNE2fP3+eTktLo2/evEnTNE0XFxfT3bt3p6urq+kVK1bQvXr1osvKymiTyUS/9dZb9IIFC2iapunevXvTx48fpw8ePEgPGjSIpmma/vTTT+lp06bRRqOR1uv19JgxY+iDBw/SCxYsoLOysmiapum6ujo6LS2N3rFjh0W+7du30zRN0zNnzqQ///xz+ubNm3RSUhJ99OhRi0xdunShL168SP/nP/+hO3XqRF+7do2maZr+8MMP6Xfeecfp75Hge4juHMxms40Duk2bNqAoSlGFRSAw9OrVCz///DOio6PxxBNPQKVSoXfv3vj6668tu9f4+HjEx8ejWbNm2LhxIy5cuIBffvkFHTp0AACMGTMGixcvxpgxY/Dtt99i1KhRlvP37t0bKpUKYWFheOCBB3D79m0cPXoUZWVlmDRpkuV9FEXhzJkzaNKkiZ2MjFnp5s2bmDBhAqKiotC2bVsAQH5+PkpLSy2mMeZcFy9eBAAMGDDAYqYdMWIEPv30U8ycOZNzLA4cOID33nsParUaarUaX331FQCgc+fOyM/Px2effYY///wTpaWlqKmp4R3T48ePo3nz5njssccAAK1atULHjh3xyy+/gKIoJCYmIi4uDgDQtm1b7Ny5U+AbIjRURJVDcHAwrl69alEIhw4dQmBgoOKCEQhAvXLYsmULAgMD8eSTTwIAunXrhtmzZ1tMSgCwadMmbN68GWPGjEF6ejoiIyNx+fJlAMATTzwBnU6HgoICHDp0yCagIigoyPJ/iqJA0zRMJhPi4+OxZcsWy2vXr19HVFSUoKxRUVFYtmwZ0tLS0KlTJ6SmpsJsNiMpKQnLli2zvO/atWuIjY3Fzp07oVarLccZ8xkfGo3GZmF27do1BAUF4YMPPoDJZEJqaipSUlJw7do10AK5rWaz2e4YTdMwGo3QarWcY0LwP0SjlaZPn47x48fj4sWLGDlyJCZNmoS3337bHbIRCOjatSuKiorwyy+/oGfPngDqFyxt27bFV199heTkZABAXl4ehg0bhmeeeQYPPfQQ9uzZA5PJBKB+ghs9ejT+9re/IS0tTXRx0759e1y4cAG//vorAKCoqAj9+/dHaWmpqLzNmjXDxIkTMX/+fNTU1KBbt27Iz8/HuXPnANQnlQ4ePBh6vR4AsHv3bty5cwdmsxmbN28W9GskJSXhu+++g9lsRl1dHaZMmYJff/0VeXl5mDRpEgYOHAiKonDs2DHLvavVahiNRpvzPPbYYzh//jyOHz8OADh79ix+/fVXdOnSRfT+CP6D6M6hY8eO2Lx5M3777TeYzWY89thjoisoAsFVBAUF4cEHH4TBYLCJSEpOTsbixYvRtWtXAMD48eMxZ84cZGZmQq1WIzExEcXFxZb3Dxs2DAsXLsTIkSNFrxkVFYUVK1Zg0aJF0Ov1oGkaixYtQtOmTSXJ/NJLLyErKwurVq3C22+/jQ8//BBvvfUWaJqGRqPB6tWrERISAqA+GnDChAm4desWOnfujIkTJ/Ke94033sAnn3yCIUOGwGQyYeDAgejXr5/FBHbPPfcgODgYnTt3tpitevfujYULF9rUQ4uKisLy5cvx0Ucfoba2FhRFYf78+XjooYfw22+/SbpHQsNHtLbSyy+/jM8//9zm2LPPPovNmzcrKhiB4Er+97//4bvvvrN7lj3JypUrcevWLcyZM8fTohAIdvDuHKZMmYLz58/j0qVLSE9Ptxw3Go2CdlECwdt4/vnncePGDaxcudLTohAIPgPvzuHy5cu4cuUKZs+ejY8//thyXK1Wo1WrVoiIiHCbkAQCgUBwL6JmJbPZDJ1Oh8LCQhiNRrRr187prE8CgUAgeDeiDunCwkK89tpraNy4MUwmE65fv45//vOf6NixozvkIxAIBIIHEN05jBkzBpMnT7aU6y4oKMDSpUuJQ5pAIBAaMKKe5aqqKotiAOpjrXU6naJCEQgEAsGziJqVVCoVrly5Yonxvnz5sk1Wp7dw61Y1zGbHMjmjo8NQXl7lYol8CzIGZAwYyDj4xxioVBQaNQrlfV1UOUyaNAkjR45EUlISgPpaMXPnznWdhC7CbKYdVg7M5/0dMgZkDBjIOJAxEFUOffv2RYsWLXDw4EHQNI2JEyciPj7eHbIRCAQCwUOIKodnn30WI0eOxLBhwxAcHOwOmQgEAoHgYUQd0pMnT0ZeXh769OmDOXPm4MSJE+6Qi0AgEAgeRHTn0LNnT/Ts2ROVlZXIzs7GnDlzQNM06UVLIBAIDRhR5QDU11M6ePAg8vLyUF5ejoEDByotF0FhCgpLkJl7DuWVekRHBOKFtEQkNo/0tFgEAsFLEFUOH330EbZv346EhASMGDECy5cvt/TGJfgmBYUl+HL7adQZ65u+lFfqkbHlGMYOaI2kxDgPS0cgELwBUeUQFhaGb7/9Fs2aNXOHPAQ3kJl7zqIYGPQGEzJzzxHlQCAQAEhQDtOmTXOHHAQ3Ul6pl3WcQCD4H6Qxgx+iouQdJxAI/gdRDn4IX+KnnyeEEggEK0SVw6ZNm+yOrVmzRhFhCO4hOiJQ1nECgeB/8PocvvnmG9TW1mL9+vXQ6+/aog0GAzZu3IhXXnnFLQISXM/w5HibaCUACNSqMTyZlEUhEAj18CoHjUaD4uJi1NbWori42HJcrVZj9uzZbhGOoAxMRBLJcyAQCHyINvvZtWsX+vbt6y55HKa8vMrhKooxMeEoK7vjYol8CzIGZAwYyDj4xxioVBSio/lbPov6HLp164YPPvgA48aNQ0VFBebMmYPq6mqXCkkgEAgE70JUOXzyySeIiIhAeXk5AgMDUVVVhTlz5jh10e+//x6DBg3CoEGDsHDhQgBAUVERnn76afTv3x9/+9vfYDQanboGgUAgEBxHVDkUFRVh2rRp0Gg0CA4OxpIlS1BUVOTwBXU6HT755BNs3LgR33//PQ4dOoQDBw5gxowZmD17Nn744QfQNE16VBMIBIIHkdQm1BqTyWR3TA4mkwlmsxk6nQ4hISEwGo3QaDSora1F+/btAQDDhw/HihUrMHr0aIevQyA0FNhFEocnx5MyJwTFEVUOnTt3xuLFi1FbW4v9+/fj66+/RpcuXRy+YFhYGKZOnYrU1FQEBQWhS5cu0Gq1iImJsbwnJiYG169fl3VeIceKFGJiwp36fEOAjIH3jcG+w5ewYccZ6A0mAPUlTjbsOIOI8CCkdFKu3pm3jYMn8PcxEFUOb7/9NtasWYPw8HAsXboUPXv2xKRJkxy+4OnTp/Gf//wHe/fuRXh4ON5++23k5+fbvY+i5NVyINFKzkHGwDvHYH1OoUUxMOgNJqzPKVQs9Ngbx8Hd+MMYiEUriSoHo9GISZMm2SiEo0ePWkxAcsnLy0NSUhKio6MB1JuQ1q5dixs3bljeU1ZWhtjYWIfO760Q0wDBEUiRRIKnEHUevP7665bIIZPJhKVLl+LVV191+IIJCQk4cOAAampqQNM09uzZgy5duiAwMBCHDx8GAGRlZaFXr14OX8PbYPonMD/o8ko9vtx+GgWFJR6WjODtkFInBE8hqhwefvhhvPXWWzh79iyeeeYZHD9+HN99953DF+zRowcGDRqE4cOHY/DgwTAajXjllVewZMkSzJ8/H6mpqdDpdBg7dqzD1/A2uPon1BnNyMw95yGJCL7C8OR4BGhsf6YBGhUpdUJQHNEMaQD48MMP8e233+L999/HmDFj3CGXbLzZ5zB+wR7e1754t49i15WDP9hYxfDWMXC3SdJbx8Gd+MMYOOxzWLduneX/TZs2RaNGjXDkyBHU1dUBAF588UUXitmwiY4I5LQRO2IaIL4L/yMpMY58xwS3w6scrIvtAUDPnj05jxPE4aqC6ohpgKv385fbTwMAmTwIBIJL4VUO8+fPt/z/119/RefOnVFRUYFDhw75RCE+b4KrCqojK34h3wVRDgQCwZWIhrIuXboUR44cwcaNG1FbW4s1a9aguLgYr7/+ujvkazC4wjRAwhqVo6CwBFl5BSi7pSPmOgIBEqKVdu/ejS+++AIAEBcXh6+++grbtm1TXDCCPSSsURkYc13ZLR0AEmpMIAASlIPBYIBWq7X8rdVqZWcvE1wDCWtUBhJqTCDYI2pW6tixI6ZPn44RI0aAoihkZWXhsccec4dsBBau8l0QbCHmOgLBHlHlMHv2bCxfvhzz58+HRqNBUlIS3njjDXfIRuCAhDW6HleGGhMIDQVR5RASEoL33nvPHbIQCB7BVaHGBEJDglc5TJ06FcuXL0d6ejrn69nZ2YoJRSC4E2YnlpV3nkQrEQh/wascJkyYAKDerEQgNHSSEuMwOKVVgy+ZQCBIhVc5PPLIIwDqK6R++umnNq9NnjzZqYY/BAKBQPBueJXD3Llzcf36dRw+fBg3b960HDcajfjjjz/cIhxBGbjqMw1O8e+uVwQCwRZe5TBixAicPXsWZ86cQf/+/S3H1Wo1OnTo4BbhCK6Hrz5TRHiQYp3FCASC78GrHB599FE8+uijeOKJJxAXRxxzDQW+hK8N24uw8NUkD0lFIBC8DdEMaaIYGhZ8iV03/iodQSAQCICEPAeCa/F0Pwa+hK/GjYLdJgOBQPB+RJUDTdN2tZRu376Ne+65RzGhvBlnJndv6MfAl/A1NrWNW65PIBB8A1Gz0vDhw+2OjRo1ShFhvB1mcmdW3nKrd3pDgbekxDiMS02wlIaIjgjEuNQEpHRq5jYZCASC98O7cxg3bhxOnDiB2tpadOzY0XLcbDajTRv/XGU622zHWwq8kfpMBAJBDF7l8I9//AMVFRV4//33bbrCaTQaxMTEuEU4b8PZyZ0UeCMQCL4Cr3IICwtDWFgYNmzYAJPJhKqqKtA0DQCorKxEZKT/xcQ7O7nLKfDmacc1gUDwb0Qd0ps2bcKCBQtgMBgsyoGiKBQVFSkunLfhbPXOpMQ4/H65ArlHr8JMAyoK6P6ovYnHGxzXBALBvxFVDmvXrsW3337rt34GNloNhTpj/f/DgjUY1fdhWdFK+SdKYK7XsTDTQP6JErS8P9LmHHy+jbU5pwAQBUEgEJRHNFrpnnvuIYoBd1fz1bUmy7E6g1ngE/ZIjVbi82GYaZDexgQCwS3w7hwqKioAAO3bt8f69euRlpYGjebu2/3N5yAlUknMTyDVoc3n2+C6pq/SkHwqDeleCAQGXuXQrVs3UBRl8TMsWLDA8po/+hzEJnYpfoLQILXNzoMhNEht8zeXb0OKLL5CQ/KpNKR7IRCs4VUOp0+fdqccXk9YsAZVOqPdcSZSSYqfgJ1pzsA+zkwqa3NOWfwTXNeUgzetbp3NF/EmGtK9EAjWiDqkMzIybP6mKArBwcFo1aoVevbs6dBF9+zZg4yMDNTU1KBHjx6YNWsWDhw4gPnz50Ov1yM1NRXTpk1z6NxKUFBYAl2tvWLQqClLpJKYnwAAp3LhO85MLK7obextq1tvSQZ0BQ3pXggEa0Qd0sXFxfj2229RUVGBO3fu4D//+Q/27t2LFStW4B//+IfsC166dAlz587FqlWrkJ2djVOnTiE3Nxfvv/8+Vq1ahW3btuHkyZPIzc116IaUIDP3HEwcK/hArQpJiXGiDmJmJcm34uc7zlfqQu6E7g1lO6yROw7eTEO6FwLBGtGdQ3l5OTIzMy1Z0RMnTsTUqVPx9ddf4+mnn8akSZNkXXDnzp0YOHCgpRT40qVLceHCBTzwwANo1qy+vk96ejp27NiB5ORkuffjMqzNMHww/gMpk2x5pR4T0tvK3gm4otSFEqtbZ8xUzuaLeBMN6V4IBGtElUNFRYVNuYxGjRqhoqICAQEBNtFLUrlw4QK0Wi1eeukllJWVoXfv3mjVqpXNNWJjY3H9+nVZ542ODpMtizUxMXfbZO47fAkbdpyB3mDvPLb5TKNgxMSE46aESTamUTAGp7RCRHgQNmwvwo1bOjRuFIyxqW0UL3oX0ygYZRz9Ghj5LX/HSGsVyh6f8ko9Nuw4g4jwIEn3Mjgl3CPjIAWpY8DgzffiDHLHoSHi72MgOrs3a9YMf//73/Hss88CALZu3YrmzZvj2LFjUKlErVJ2mEwmHDp0CBs3bkRISAhef/11BAfb9xLgc97yUV5eBTOX91YCMTHhKCu7Y/l7fU6hqGII0KgwtMdDKCu7gyiB0FP2exObR9p1XGOurZTTeGiPhzhXt4xMgP0YCME1PnqDCetzCiW3GhUaB08hZwys8cZ7cQZHx6Eh4Q9joFJRgotq0dn9008/xZUrVzBs2DCMGDEC169fx8cff4zCwkLMnDlTtkCNGzdGUlISoqKiEBQUhCeffBL5+fm4ceOG5T2lpaWIjY2VfW5XIWZuYdv+hyfHI0DDPZRS/QTOlgMXwlW+CwbihCUQGj6iO4eoqCj83//9n93x0aNHO3TB3r17Y+bMmaisrERoaCj279+PAQMGYM2aNbhw4QLuv/9+5OTk4Omnn3bo/K5AqMDe4te72x1nJlmhVb/YrkDpkEhXlukm1WUJhIYPr3KYOnUqli9fjvT0dM7Xs7OzHbrgY489hpdffhmjR4+GwWBA9+7dMWrUKLRo0QKTJ0+GXq9HcnIyBgwY4ND5XYGYk5FvouebfKWEkvrSapw4YQmEhg+vcpgwYQIAYPbs2S6/6IgRIzBixAibY0lJSfjvf//r8ms5gtBOwJGcASm7Al9ajUvZKREIBN+GVzk88sgjAIAuXbqgoqICOp0ONE3DZDLh4sWLbhPQU/DtBBwx/0jZFXhiNW69A4ppFIyhPR6SPMGTbnIEQsNG1OewfPlyrFmzBgCgVqthMBjQsmVLh81Kvo4c8w8z+fJhvStw92qcvQMqu6UjNYEIPoc3lYVpaIhGK33//ffYu3cv+vfvjx9//BELFixAy5Yt3SGbVyI1I5YdfcSGa1eQlBiHxa93x4T0tgCAz7JPYcaqfEVKdHtb1jSBIBclI/wIEpRDVFQUYmNj0aJFC5w+fRpDhgzBhQsX3CGbV8IVtso10XNNvgxCoaTueuB9yQFOIHBBFjjKIqocNBoNLl68iBYtWuDQoUMwGo2orKx0h2xeidScAaFJdvHr3R1yXrsSUhOI4OuQBY6yiPocJk6ciNmzZ2P16tVYvnw5srKykJKS4gbRvBcpzlhHo4/c9cDLdYAT2y7B2/ClCD9fRFQ5GI1GfPnllwCArKwsXLhwAa1bt1ZcMF/H0egjdz3wbAe4ULSSu0p+EwVEkAPJt1EWUeWwdOlS9O3bFwAQHByMhIQExYVqCDgafeTOB956ByRUS8YdDW28recEwfsh+TbKIqocHn74YaxevRqPP/44QkJCLMcTExMVFawh4EgugDc+8O4wdZGOagRHIPk2yiGqHI4dO4Zjx45hy5YtlmMURWH37t2KCubPuOOBZ5twXkhL5K2o6g5TF3EuEgjehahy2LRpk6UxD8PZs2cVE4igPFwmnIwtxzB2QGtOpeQOUxdxLhII3gVvKGtFRQUqKirwyiuv4Pbt26ioqMDt27dx48YN2d3fCN4FlwlHbzDxhsu6uuQ3F1LzRwgEgnvg3TlMnz4d+fn5AICuXbtajqvVajz11FPKS0ZQDEdMOEqburzR10Ig+DO8ymHt2rUAgPfeew/z5893m0CeYN/hS1ifU+jQpMQVfgl49yTnrSYcTzoXCwpLkJVXgLJbOq/8zggEd0PRNC3YW5OmaVy+fBnNmjXDjz/+iOLiYowbNw7h4d7VX9XRNqEFhSV2/aIDNCpZ3dusTTRqCqBUFIymu7JIPZ+74JI7UKvm9Tk0dLjGw9u+M3fiDy0yxfCHMRBrE6qeN2/ePKETzJ07F0eOHEGzZs0wbdo0NGrUCDt27EBqaqqrZXUKna4OwmqOmxVbj+GOzmBzzGSmcaS4DHnHryI8JADNYrkHcMXWY6jSGW2O0QDYOspkpnGhpBL9OjeXL6ACNIsNQ/Q9QbhQUgmd3oToiED0bN8Eu369hH/v/l30vhsaXN+jt31nSlBQWIIVW4/ZfeehoYGoqanztHgexR/GgKIohIQE8L4uGq108uRJbN26FWvWrMGwYcMwffp0j7bwdDVCdnaxRCw5YZbllXrMWJWveBluMXMW+z1MBVjr3ZO/JaD5YxitUNLh4BTvsgoQPIOocqBpGiqVCvn5+Zg4cSIAQKfTKS6Yu+CzvzMIJWKJfZaNkpOulAxjrvd8ln2K83yuTEDz9rIYzvpgvP3+uBBKOhyc0spDUhG8CVHl0Lx5c0yYMAGXL19Gly5dMH369AZVQmN4crydz4ENnwLgiv/n8jlYI2XStZ5sQoPUoCgKVTqj4MQjJcNYqIw4F65YOftCWQxn8jh84f648MfdEkEeosph/vz52LlzJzp16gStVovHH38cQ4cOdYdsbiNAqxJUDmHBGsxYlW+3MuQLvwSAb3YV29mxGYR+gOzJprr2rlxCE4/Yj72gsMShH35BYYlTk5wvlMVg5MjKOy87WskX7o8Lb41YI3gPosohJCQEQ4YMsfw9atQoRQVyJ1xRKlzoao2WiZ49QQv1muZTDkI/QLHVPd/EI2Ti2vjDaeSfcKxZkCOTnPXOhw9vW6EmJcZhcEor2REqvroCJxVNCWKIKoeGjFQzC9tCJGVlKDQ5tIuPduhzQu8ZnhzP6z/IPXrVLoJKKowjXao9XarCZd7rzatrKfjqClzppENf9MMQbPFr5eDM6k7ss0Ir+dyjV9Hy/void+wfkBQnN9fEk5QYx6scHFUMDOyWpcz1uJDj15CzK/HWycaXV+BSkw7ljr2v+mEItvi1chCbiAM0Kmg1lI3d3/qzQgit5M008EXOKZgBS24GEznU5oFIVNzR2+1WrGXim3j47kdFiSuI8BAt7tQYhN+E+l3Tpp1neCcLudFbUvDmyaahl/1wZOx91Q9DsEW0hzRDUVERxowZg6effhr79+9XUia3wVXsjSE6IhDdH+V+kF2xMjTR4EzaK7pQwasYKNQ7zz/LPoUZq/JRUGjrR+ArXpfcvgnvfQZoVJiQ3habPhooWfbqWpPdboKRRY45Rep7vb2RfFJiHBa/3h1fvNtHsD+4L+LI2PuqH4Zgi+Sdw6pVq7B06VIAwIQJE9CzZ0/FhHIXSYlx+P1yBfb+dtXmeIBGhXbx0cg/UWL3wwgL1mBU34dFt9Vrc7h3Dc5AAzaO8c+yT+H3yxV4vn99aLHQKrbl/ZGW48xOgr3KlZu3wWC9KpQa3mutYMXMFmSy8RyOjL2v+mEItvAqh/Hjx+Pdd9/Fww8/DOBujSWNRgORckw+xfFz5XbH6oxmXiduoFYtyd7qrJ1fKnt/q/dfMDLx2ZGl2Je5JnapMJOBUHgvlwKQYrYgk43ncGTsfdkPQ7gLr3KYN28eFi5ciLi4OLz55puYMWMGVqxYAYPBgA8++MAlF1+4cCFu3bqFBQsWoKioCJZHSqUAACAASURBVLNmzUJVVRUef/xxfPDBB9BolHeJ8K2A+CZ3rvdbr3yl2PddDbNLcSbkNKZRMIb2eAjjUhNsJvF28dE4fq7c8rfeYOIM0bWeLIQUFBsp9mklJxuuMWhIZiFncWTsG7ofxl8Qrcq6b98+rF69GiNGjMAzzzzjsgsXFBRg2rRpSElJwYIFC5CWloaPP/4Y7du3x/vvv49HHnkEo0ePlnw+R6uyMmGaUlFRwEtpbXlLUngSxuQFiP8wHa1EKvY568mWou76VUKD1OjS5l4bRSPktAeAL97tY3NdV082pBorN+yKpN4aKaYkpCqriHL4/fffodVqcd9992Ht2rU4cOAAZsyYgXbt2jklFNNhbuDAgTh9+jQmT56McePGYdeuXQCAQ4cOYcWKFdiwYYPkczpTsltoguLCegKRq1yURqOmQJtpO6d27w5NLL4JgF8pRkcEYvHr3QWvwbVTYnYZXH4aPsSiwbjkcOVE5cwYNGT8YWIUwx/GQEw58Npt5s6diytXrkCv16Nr16544403MGTIECxatAghISH49NNPHRZqzpw5mDZtGq5duwYAKC0tRUxMjOX1mJgYXL9+3eHzuxIuM1Gd0YzPsk9hbc4pl5qQAjQqdH80jtffEaChUGcUviBfTae9v13Fr6dLLc50Z5y8zGTM9hWwHfti1BnNCNBqEKBRSTJbuDqklTi6CQR+eJXDsWPHkJWVBQAYNmwY3njjDTRp0gTLli1DQUGBwxfcsmUL7rvvPiQlJSEzMxMAOB3cFEXJOq+QBuRj3+FLWPu/It7XYxoF48Yt/gq0rlAMgVoV6gxmNG4UjLGpbZDSqRk6trmE5d/+ZjPRa9QUJj/bAQCw9N+/ObRLqtIZsWHHGUSEByGmUTDKOO4tplEwYmLulmzed/gSNmwvwo1bOhsZs/IKXGJKq9YZ8dbojpzXYMN1zTqjGVl55xERHiTpHOx7lTIG/oi/3z9AxoBXOcTFxWHevHmora1FYmKizWtJSUkOX3Dbtm0oKyvDkCFDcPv2bdTU1ICiKNy4ccPynrKyMsTGxso6ryNmpYwtx3g/E6BRYWiPh0RrBDlLWLAWi6ffNWGUld1B5Z1a0Cy5jCYaf990BAAQqKVgovgrvwqhN5iwPqeQ19E4tMdDlu00e6VedkuHlZuPovJOLeek6ghREYFIbB6Jha/aPlNcW3q+azJycckptKMY2uMh0THwR/zBpCKGP4yBw2allStXYv/+/QgICED37q6zv65bt87y/8zMTPzyyy+YP38+0tLScPjwYXTq1AlZWVno1auXy67Jh1Al1gBtfdKYM+GdUuBSPJm553gT4QBAb6BBUfVKQm+QryDKK/V2ESVckTpCkUSO5kRYwzYfifkThDLAHcnIlTIGBIK/wqsctFot+vTpw/eyy1myZAlmzZqF6upqtG3bFmPHjnXbtbmo0hltHNXWkTeuJCzY/iuQMunSNBAWHICxA+IFy4NzwYSdWoecFl6swPqcQnyWfcoyMQvZ5Cekt5WlNPmilayjvtZtK7Lshsor9Vi3rcgiJ1BfsJDLryEn7JiN9RgotVr0x2gfgu/j0dpKw4cPx/DhwwEACQkJ2Lp1qyfFEYSmARVFgYJ9JJAz1HHsXqSuypkdgFB5cL7PWbcsLSgs4WwTGhas4c1pYK+6Q4PU0BvMdlnQ1mGhBYUlnEmHQH3/C7aZzGii8c2uYsvn+T7Ll1siliTHnrRfSEtEYvNIwc/IxZvrQhEIQkiurdQQCdDIc3qbaRrBQa7Vp3VGWlKNJC4CtZTDobTWNZEyc8/ZmdjqjGbQNM0pB6NcAFhqCq18Mxk9290H2ziCuzN2QWEJvsg5ZVOT6YucU5Z751Nu1seFEhaF5GSPLyPPl9tP28iTseUY53udwdvrQhEIfMhSDmVlZUrJ4RHGpbaBPPXAP4k5A3uiSEqMw7jUBNGVr95AO2X3ZyYpvnNU15p45WAX3CsoLMH+49dsTG91RhqfZ5/ClOU/4bPsU3Y7LhMNbNp5RrK8fOMRHREoWU4GrklbbzC5fNL21XDZgsISzFiVj/EL9vAqWELDRpZyeOWVV5SSwyMkJcbh5fS2iGkU7FE5mBWu9Q/RutKn0tcWmnQZObjeY70Czsw9xxk9ZV0skIvqWpNlF8JFaJDa8n++qrNiZTS4VurumrSFxtZb2Xf4kt2uikvBEho2spRDQyq4x5CUGIexqW0kmXGUhO+HWFBYApXc7Y0MwoI1vJ3prM0yYpOpEo2T1BQw+qnWlr/ZOypmxwDAZjKTcg0pk7YrVs+OKjRPsmF7ETGFEfy72Q8D14/Bk9QZzfhmVzE27TzDWVrClVTpjNh//Brv64yyCg1SCzY9ckVoK/u8zATKblPKLm0xY1W+6PfHVgZcIcqBWrVNGXG2I5ldIl0KvliEji/x09tNYQTXIks5rFmzRik5PIpQFrSnUMK3wYdYMp2UMhfDk+NtQlGdZfHr3QUjfQBITlDkWqlzTdrW0Up87U7ZJdLZ8IWtyn2/J2nMkzkeGqSW1U+c4NvIUg5ys5Z9Bb4fQ0MhLFiDQK3aqZUfo6z4GgUBQFCA2kapBWgomGlu5RMdEYjaOqPgboQv0uebXcWoM5gl7faEdiDsSds6z0ForNjJddYTvDViYatcyu+LnFPYsKPIktwopbmUq+mcEIttBRfsjuv0ZlTXSu8nTvBtiFkJwNjUNpbSFA2N+uQ956KarDHT9XWeauvqkwSFKtpqNSrRxDeuRDrGB8Ins5RdFbuMuNxcAyEzmfVxsZLtQpnaXMrPRAMmq6z3Kp3RLhlQLnJ3J7+eLuU8bmb5HElf6IaNX+c5MKR0aobeHZp4WgyHUFH15bj5HKw0DZf7LYwmWtI5q2tN2PvbVbSLj+bsr5yUGMfZpzv/RAkKCkscjuhhHNXWpiO5DlYhh7G1XJt2nhHdwTgbGWU00Q47g7nyOcQij+SYWYkfouEiunPQ6XQ4c+YM2rdvjy+//NLSf6FJE9+cTPl4vn8CWt4f6fIy3EpjpsFb4ttb4LLT85ligLsTt9y6Vnx9GIQm55cX7kFy+yZ2Tmah/uLWTmspSjI0SM25eudz8vPJKgTf7kBKpz02csys3hySS3AO0Z3De++9h927d+P48ePYsGEDmjRpgtmzZ7tDNo/gzZMsH74g8ze7ii3/Z69muWCigyhK2s2pKfvVPhOKKoSZrldei7+xNys+3z8BE9LbWibAsGANtBoKn2WfwoxV+Tb3JESN3oR124psVu/rthWhRi99Ryc0CQvtDhzZtXCFdqupenOiNd4ekktwDlHlcOnSJUyfPh179+7FsGHDMHnyZFRUVLhDNrfC/MAIylClM2LK8p8sK1ypuwGpVWcpVjKIFAVkTdGFCuw7fMnuOJMEOCG9LeoMZstKv7xSLzmijOZwyhtNtKxCjkKTsFj1XC6ElE1Kp2Z2+STj09rixYFt7HJMiL+h4SJqVjIYDACAvLw8zJw5EyaTCTU1NYoL5m7kTFiOolZRMPnCMl8hqnRGxcqfW9vlHe3BsWF7kV1fCQZ3PB98BGgowUlYbvVcKSt+vvBbogz8B1Hl0LFjRwwcOBBqtRodO3bEuHHj8MQTT7hDNsVhVrE3K/WQM2X37tAELe+PlNx7WqOm8OLANvhy+2m/Vg5A/YqWr4qqszDmFEcncSE7uycdr2Ya2PjDad6oL77IKq7quSQ/gSAVihapiWEymfDbb7+hdevWCA8Px759+9CrVy+oVN4V6CS3E5xYCKIQARoV4ptGoOiCdPOamoJLS30T7HFW6ahUFD5/pzfna3zVb5naT3IjwjRqCrTZtvx7gEYFiqIlmdKEQnXZr8vF1X0tvDHRTwzSCU6Cz0GtVuPGjRv4/PPPodPpUFVV5XWKwRGcMRPUGc2yFAPgHYohNEgNma25fQpndyNCiwu+Gkmjn2qNoAB56UKBWjUCtSqYaFjqZjE2fKk+FutQXL6aU66cgB2tM+VIKC3BOxB9qtesWYP8/HyUlJTghRdeQEZGBi5cuIBJkya5Qz7F8Lf47NAgNQxGeU5Qf0OoOq+QeUaqeZFBbzBBX+/Ks/SisA49lfpsWr9PqESHszjTsMiRUFqCdyCqHP73v/9hy5YtePbZZ9GoUSNs3rwZI0eO9Hnl4OpCcd6O0gX8PAGzUnbF9xigUWFsahu74wWFJTZtWEOD1JiQ3tZmYuPrmCcV612AUF9zLqw7+imFMxO8r/azIEhQDhqNBgEBAZa/IyIioNH4ftWN4cnxsld8BO+BWW274jusN7dR+L9NRxBltSNg97UG6pXsFzn11VkZB7Er4HOmM72380+UcJpB3VHjyJkJXshZTvBuRJ0H9913H/bt2weKolBXV4fVq1ejadOm7pBNUciW1j0o1Yui+6P1ZhQpPhQhGcKCNdAbzKjSGUHjbvLd+AV7sDbnFGfRQNNfiXNSJkfrZkVicE3+QQEaPN8/QbAzoNK9FpxpWOSL/SwI9Ygqh9mzZ2PdunWWEho//fQT5syZ4w7ZFIesXpTli3f7iDqJA7XSJ09rfj1dioLCEkk+FCEZqnRG3jLjzjq4w4I1GP1Ua6caSTEKiEnGE3ufEjgzwbvDWU5QBkn2oS+//BI6nQ4mkwlhYWH4/ffflZbLLRDTkvP2cj6so3CEJq7V05MF6yzxUaUzSi5fAQAUICuXxVk0asqm1La1I1vOfbIXMJ4w0zibK6Gks5ygHLzKgSmRMWHCBGzcuBE0TYOiKNy4cQOvv/46fvzxR7cJqRRJiXH434HzuFrumV4OnnaKf/FuHxQUliiiIJlVt5ACZiY068mDHROvN5h4lZccpebuIK1ArcpyT+zJkS9nggt2C1euYoQBGhXaxUcr2oiHTPD+B69ymD59OvLz64uWde3a1XJcrVbjqaeeUl4yN1BQWOJRxcDV7cydyJmkHGHK8p94J3A+swR7ElJKeSmNUHSYnGqz+SdKbCracq3i28VH2zisSSMegivgVQ5r164FUF+Vdf78+W4TyJ14qmG6Rk0htlEwXl64B2a63uThCZxVDMwKlW/y5lMMcrqbJSXG8fbSllPy2t0wuyKh7GAppjSukFGunQjJJSC4GtHyGQBw7Ngx7N+/HwaDAT169EDnzp3dIZss5JbPAIDxC/YoJM1d2GUSKAq4LyrYYzsWb4Cv7wIfQrurQK1adm4AGyF/hEZNIVCrkq2EJqS3BQDRshZCuytrvni3D+9rQs9xaJAao59qLUtJiJWO8MVyGHIh5TMkRCtlZWVhypQpuH37Nqqrq/HWW29h8+bNLhXSU4QFK5+v0er+e2zKSdM0/FoxAPU7FqllGKxLfHOFpDqrGKIjAhHK8xyoKKBnu/scOq9Qo521Oacs9y1hbSbqbBZ6ncnLYK7naBkMBlIOw38QnR3Xr1+PLVu2IDY2FkC9g/qll17Cs88+q7hwSlJQWAJdreujdNjIrcHkT4jZxtk7BldXcmV2MHwrb6YREBuxCC+xzG0zDazbVsRrLrNGSsiomA/DRN81oTpaBoOBlMPwH0SVg9lstigGALj33nudLryXkZGB7du3AwCSk5Pxzjvv4MCBA5g/fz70ej1SU1Mxbdo0p64hRmbuOa8ohufv1BnN2LTzjOQWl66kts5o6VUt1/8iVAGWmcyFzms00TCahBWDVJMN87qQ4768Uu+SiZ2Uw/AfRGf5yMhI7Nq1y/L3rl27cM899zh8wQMHDiAvLw/fffcdsrKyUFhYiJycHLz//vtYtWoVtm3bhpMnTyI3N9fha0iBPMzKw5jtxBzu1bUmWS0upRIdESiYZFdda8Jn2acQ2yhYVqJalc4ouIthSm/U1jm2M1VT9WPHTOhSTDZJiXGC5iUhReVM3gVDWLDGKXMVwftQz5s3b57QGx555BG8//77WLt2LdavX4+CggIsXLgQjRs3duiCtbW16Nq1K5o2bQq1Wo0TJ05Ar9ejuroa48aNg0qlAk3TyM3NRd++fSWfV6erk1VxNO/4Vehk9PBlw1Q5JfATHqLFc/1a48Qf5bKaHJnMNC6UVCIkUOPwd6SmgNH9WuOXolLR9964XYtej92Har0RNRJNjdERgZyyRUcEIjwkAF9uP43aOvm7ntAgNUw0LJ/V6U04+Uc5ou8JQrNYfuchAISHBOD4uRt2iosZiwsllbwy9+vc/K4MoYGoqanjvcZJ1vepUVMwGMyo/uvccmT2VoTGwNUUFJZgxdZj+Pfu35F3/CrCQwIkjZujn2OgKAohIQG8r4ualeLi4rBjxw78+eefMJvNeOihh5wqvNeqVSvL///8809s27YNzz//PGJiYizHY2Njcf36dVnnFfK6c/FCWiIythxzyKEZ0ygYY1PbOPx5OQQFqFFb553hmmLcrNQjK++8Q6YhvtVsoLa+J4XYmJhoICvvPGIaBQt2eGM4+ectfDGrH8Z//KPo+8NDtJzPT6BWjRfSErFhe5FD9xzTKBi1eqOduanOaEZW3nkMTrn729l3+BI2bC/CjVs6NP7reRyc0goR4UFYk3UCd2oMFllfGfooUjo1Q0R4EK/MMTHhtrKw/ra+XliIFoEBalTVGND4L5mZ6wnJ7Guwx0AJ9h2+hA07zli+k/JKPTbsOIOI8CCkdGrm8s/JQXSWf/LJJ5GSkoJnnnnGpSGsZ8+exauvvoqZM2dCo9Hg/PnzNq9TMrvSyA1lTWweibEDWlvahIYGa6CrNYr6IQI0Kgzt8RAq79RCo4alLr8S9O7QBPknfHd7HhURKGlilgqT8PVLkbSFQ9ktHSakt5WURMfIObTHQ6IJav/vyVZIbB6JJx65F7lHr8L8V9OeJx65F4nNIx265wCNCo882IjTAc7Ix4RWsh31Zbd0WLn5KCrv1FrOBdz1WSQ2j0RZ2R2bZ97av8O8zsAO42Rf706NAQEaFV7+q3Q5n0PfWmZfw12hrOtzCu0WmHqDCetzCpHYPNLln7NGLJRVVDns3r0bOTk5WLRoEe7cuYNnnnkGw4YNQ1RUlCQBuDh8+DCmTJmC999/H4MGDcIvv/yCGzduWF4vLS21cYIrBZNMxDwIUmr8jEtNAGAfv64EfBOFt8GVJ8BE2azNOeWSKCMVBbtMYKHrA3d7KEtRDhQFjJ69DXdqDAgNUkPIXcD4FPJPlFjuzUzfzWZ2pF1pgFYlqPSsq7vyOZY37TwDg5G2i0ayLi/uSF6CmCOblOV2HEf9QO4IDBD1woWHh2PUqFHYsmULli1bhh9++AHJyckOX/DatWuYNGkSlixZgkGDBgEAHnvsMZw/fx4XLlyAyWRCTk4OevXq5fA1HIWpfMlX4llF8cev+yuhQWqo1fYDxpTUdlX4KRNWyjXuocEap8tC0zQsphGh8FJmwhOaMB255yqdUfC6Or0Zk5flYvyCPbwTQHWtiVMm6/Li5ZV6rNtWJMthLDYRkbLcjuNoOXRnyqhLRVKIRmFhIT7++GO89NJLiIqKwvLlyx2+4Nq1a6HX67FgwQIMGTIEQ4YMQWZmJhYsWIDJkydj4MCBaNGiBQYMGODwNZyF78fNHHeVdmZWg77Y2zlAo8KE9LYICtBwlrw+fq4cgHtWj1U6o11Z6O6P1itxV2bBW094QhOmEvdspmmXlQoxmmhZFW2FIpQAUpbbGRxVrO5QyKJmpfT0dOh0Ojz99NPIzMzEvffe69QFZ82ahVmzZnG+9t///tepc7sKsW2yq6qpWpc1cEcpD2dhzDdM5zSxuPqCwhJRhz1Fwem+1oz5yLqyq6vMfsx3zTbH8JmOGPOXkEnQEbOTGAEaFQK0KsmVauVUtB2eHG/XEQ8AdLX1eSLM2LtCGXh7aQ5Xy+doOXRny6hLQbS2Un5+Prp3l14Hx1M4UluJQcwBB9T/+Lo/GufS1pChQWoEBdTHsysxYSjBhPS2kibesGAN6gxmbjOQSOvL0CA1avQmyUqDXcjPVdVmrWtAWU8KYgX/nCkIGBasQU2tcB4F1/VGP9UagDxfGF+9Ji5n7ORluZz3JLdOlhB8vztP7EK4xsCb5HMFTjukfUExOENBYQmy8gpQdktn0b4AoNVQFqdkWLAGnRNieSczR6muNVl+cL6gGEKD1JL8LQGa+lwVvveNfqo173kcKWVepTPalIEQUgxSmxtp1JTlWWDLopRiYEhu3wT7j1/j7VDHhsm3kVPtld2+1Fr5xTQKxtAeD9lMeHz35EoHqJKlOVyx4ve30iHKV57zYtg/eqZ/MNvUUWcw45ei65wTlTtW/CqKgtlZ24sL0BvMqK4VngyYVbyQyUlo8rJuiyn2Xmus+ygLfSdSzSm01QmkBiAEaFSyQ7DZVOmM2H/8GkwyartYT1DMPyHlqqZg2WkA3KGx7JpL7ohIUioCh+t37ki/C38rHeJckSQfh+9Hz56H64xm3pWTmYZTPYKl4A2KAah3ZPJFcjHUGcz4/bJwsUEhp210RKClciijYKRWz2V+9K5Q1tbF6qT++MelJrik5arRRMvuXMclY4DW/rmMjgjE+LS2NpOi0IqYwR0OUKUicKTcnxDM88hHQw3ZlfSrO3v2LPLz86FSqdCrVy88+OCDCovlHlxllx6eHC+7B7IrcWd/ZEYZ8q2k64xm7BPJz+DLWWDaXbJXeRyRsrzndaXZj/k+pQQgME5xpZ8Dvl2R9QQl1zYuZUWstAOUrxaVKxSQMyt+MfNmQw7ZFVUO33zzDVasWIGePXvCbDbjn//8J2bNmoWBAwe6Qz5FkRN1FKiloDfY/yrbxUfbRGp4IurInfsKKcpQTB4mYczawW99XvYPUaqFxdXmPRVVPzmIlcS29k8IdcZzBVzKmd1DmmuxIGQbl2oyUqqPNN8ELKdjoBDOmMSETIreGE3lSiT1c8jKyrKEsF69ehWvvPJKg1AOUn/IARoVtBoV9Ab7lQ0Tz884vNxNaJAaOr3JLQ5tZhJiFIMzYah1RjOOnyu3OJ8zc895Xa9oM11fBjssWCMYqWbtn5CakQ3U2/65FJ+Yj4mtVNm7ML5PWsvOjr7SqCkbB7iaqi/HMH7BHsUnQb4JOFCrdsk1uZS71BW/0OLRVVFa3oqosTwsLMwmt6FJkyYICOCv5OdLJCXGoc0D3HVIArX1tgwmoYfPlszE81t3x3IXjGNRCcUQHRGI3h2a2CWW7T9+zXKfzrpC3D12gVoKARr5DmPGSTw8OZ5ztWntnwDqa2JJYXxaW8tzZg0Fmve5BOoXJIws5ZV65B7lzhxnY93X2nrMq2tNoM20xbcTHqIFpaIsz7zS3d5c7ehld7sD4HCSnjsykb0V0Z3D448/jrlz52LMmDFQq9X4/vvv8eCDD6KwsBAAkJiYqLiQSlLKUyQtLDgAq6ffXRnwmVGiIwLdUk6DnWehou5OSlLDM61hrxTZcK2Kpiz/SXJ4pRTcNXYALBnodQ6WWWeyioUWCQzP90+QXBeLy1RpovmfS+ZaX+Scsuw6pC4OmJUyn+kuUKvGiqm9MPNfBZxVVpUK2RQz+8gJQ+WLTBqXmuDQSl/KrsPbE/ccRVQ57Ny5EwCwf/9+m+OTJ08GRVHYvXu3MpK5AaGGMuzjQg+JVNNU90fjLFU85WBte+V6+KU6bK15cWAbQYXHIKUYoSPIGTtnCQ1SQ28wc07EcqjSGXknMnZElZg/KyxYI2iGZCYavnPI1dGBWsoyYYk98zd4FJNSuzuh35bcMFRX5yKIOeJdFSbrjYgqhz17vL+sgyMwXyofXM44gPshkTJ5dn80Ds/3T0DL+yMlJ3dxrUKccdhan5c5p9CqyJVlKIC7kTZyxw6on0wDtWrZE5TmL83pqh2PlFISzPv4xk6jpkRzQQI0lGgZDjlo1HctyGIr9cY8PTCcNaXwrbCFflszVuXLmuyVyEUQcsTzKaO1Oacsn/VVRJVDWVkZvvvuO1RU2Mauv/POO4oJ5Q6EzBl8ziq+h0QsmgW467hm/xCEQhO5tsHOrt6s781alpuVekRZ/SgLCktkl9sWC6k10/YlG6SMXYBGZdk5yY0Gk9KrWSoBmvrV96adZ+zOyZj4mDHl+56ts/CFqDPS2HfUdSXbrfN0xEwlY1PbYOXmo5IcuFJNKmIrbL7fltzJXmi3Za28XQXftcw0fH4HIaocXnvtNcTFxaFZM9d0F/IWxHo28D3gUlY+YteTEvoq9+HnUzLW9Zu4fryMLIUXK7A+pxCfZZ/Cpp1noDeYZSkG5ty/X67gXe1yrTy5Vo3t4qPtQlyB+vo+nqTOSAsqJ/b3IrTiFEqqYpDr8BfaXVmPPXvMrQspZuaewwtpiRiXmiA66csxqThq7pEbhipkqlTCZyKkjHy9tIaocjAYDMjIyHCHLG5F6KFzZuXDV/RNKOpB7sPPVxSQK6nMuvIrHwWFJTYtB8VqA3GZhxiSEuN4TWfllXrMWJXPOfkLOQudNW/xFQF0RUVY9nWkooT9nnGWs0NkuVb99mU27rabzNhyDGMHtBZ14MqZ8B0198gNQxUKJVZizMV2vr5cWkP0aU5MTERxcTEefvhhd8jjNuQ+dHw/hG92Fdv8EOSe15GHn5GHvapreX+kQ1ETmbnnZPXCfolVfkFMRmvKK/U2Owuu1SZ7h1ZbZ3TK79E5IdZubFxpz2cQKXBsg6vKvrOp0hmhUVMIVN2NhNIKhO9yPdd6g0nSilfOhO9oIpojmdlK14FiP59CgSa+HPIqqhw6duyIoUOHIiYmBhrN3bf7cpQSIG5vZz+MfD+EKp2tI1Luw+zIw89nrnA0g1XuJCXlGsx7pEQjWa82uXZozsK077Quvy0UjOAocqqxSvG1OEq9n8VWLj5zjzMOXDmTsDOJaHKfa2euJQbX85l/ogTJ7Ztw7tx9ubSGqHLIyMjAkiVL0Lx5c3fI41aYh866djuf+Ugol2DTzjN2kzs78zcz9xzvpK9UWQKpyFnFSl0JMQ5tqTDXVyLvga18XNXXYGB9cwAAIABJREFUmg3X2LAXGtYmtbBgDbQaCtW1JgRqKdQZ+AvuqSggONDxcuBcu1xAuMS4mAOXb/fVLj7a7pjStZncdS0+C8Lxc+WS/DS+hKhyuOeeexpEqQyp8H35Qltz674MjDL5/XKFzUrCm+OfhyfH2/gc+JC6EmIUrFyHNiBvp8D02ZDS+8A6G1upUiOMT8V6B8peaFhPplU6o6XdKl/cPHC3aB4gr5kPG/YuF4BgiXGx55WJwJN63J2LIGufipQFmlSEdlqeXuS5GlHlkJKSgoULF6Jfv342ZTN8PTOaD74vv7rWJLmRS53RzGmD9NbohaTEOESEB2F9TqFlRaurNdrlT3R/1P7h5zLByV39WysdqRVQrX/kLe+PFMxeZj7jzK7EOupLCOtFgJTrsZ8JKate69fkmt3Yz5/QmIk9r0rkFLgSJRLUhNrDNjRElUN2djYA4IcffrAc8/XMaCGE7Khy7MR8q1Pm3N6Wcp/SqRkSm9fX85mxKp9z0mCvCPl+fGLj0+aBSJTe0nHeu9gYc+V/WK/YNv5wmtfU4agD2jrqS0qeBTOpSp0k+UJg+Va91s+J3Jao7PeKKRih19zRAMgZNu084/LObXy/a1/o5CgXv82Q5kPImcW1qtMbTJwTqVBym7en3EtdEfKZ4PhCRCkKeJkn0oldJZSi7EukSzFrCZk6HI0Q6v5onGWSltr5T87Kns9XIeUZ4XteAZqzjhT7WlKUMR9KOn6lwrfIKigsUaS1qbcrRFciqhzWrVvHefzFF190uTDegNi2nr1y47MR8+Ud8JldlDY5ydmpSP0B8P3I+CI6QwK5SzBz9WgO0KjQu8N9dvkQzoRXTkhvK9tmHxqktvkepa4Qpe40+SZTse5l1t9lfNMInLlYATNdvyjp/ih3rglfvgMAXrMc24/C9VmhukNK7o6FFKhQ3SpnJnJvUIjuQlQ5FBcXW/5fV1eHw4cPo2vXrooK5WnkOJYcyTtwZ5IOIN/2KvUHIHclzreSE4oAEUuOY+82+LBObmR/JwA46yWpqXoTap3RXm5mB8EU9bP+rNBOs118NH4pum4ZiwCtCr9frpAcOs023ZVX6m3eyzRSanl/pOToGbYZiys3he954futFBSW2FSOLa+sryTLdQ5HEVKgQs+lMxO5O6OuPA1Fy8ncAXDz5k288847+Pzzz5WSySHKy6tgdtDwZx3K6g6EsqiVaCAi5XrsMZCy6uPbNTHhmWysyztYn1PIjs+uxSR0bT6EWmRan+/fu89aSlWHBqkx+qnWgnkajGzOlJTmk5dvDKWatdjPkhwZZ/6rgLfwntTnc/KyXE75Q4PUWPlmsqRziCH03AhVz10xtZfoud09J3gClYpCdHQY7+vS8/3/IioqCleuXHFKKH/H3VtTR6JKpOyehFbi7PvTqCnoao12DWQAx+y4UiOPpK7skhLjMDilld2EIKWsuZydptQIpgCthrNXt9T1j7XMcneOrijZzbdLdDRPg8FayQn59fh+Y6P6NqxKD0oiy+dA0zROnjyJ6Gj7JBeCdNy9NVXSiSY0MbJLYLAnBsYE4IiylDpRObMTU6LpvVS5q3RGTEhv63AvDea75Uv6E/JxKVWy21nYSo5LMQiZ9Bqq+UcpZPkcAOC+++7z+XLd3oA7E2Y84URj359Q9VlHfshSzCvOTGZKNb2XE8HEjKHccFUANo1yxMKq2cgp2c0HX0UBOcUJ2fDtuvgKQTa0pDR3I/pNNW7cGNOnT3eHLMjOzsbq1athMBjwwgsvYMyYMW65bkPHG1ZRYrsXuT9kMcXgrPLjm4icbXrvSASTXMUQGqS2KBWh6/Apz5ROzVB5p9ap52VU34ftnPxMkyNH4RsHrj4hBOcRVQ779u1zi3K4fv06li5diszMTAQEBOD//b//h65du6Jly5aKXrfeEbnfrmdugIbCuNQ2DWbl4epVlNwwRVfvXoRW4K5Qfkpl//JFMAmF7MqJCmMS9sRkDdCo0C4+2rIrcfWqW4kFiT/lGHgDosrh/vvvx/jx49GxY0eEhoZajrs6z+HAgQPo1q0bIiPrs3T79++PHTt24I033nDpdawpKCzhDGEE6hu7fJZ9ijNaxd/tl44k8fFNikxymVBEFNcEw6dsxKKSpOIpPw0XUjPz2WPIdw9MLgRf7a/BKeGSZRPC1QsSf8ox8AZElQMzWSsdoVRaWoqYmBjL37GxsTh+/Lii18zMPedQb+HySj0+yz6Fz3NOgaaZZjImS0aqs3Zpb8fRJD7ryUKKgpHyHqVMZVwTkZqq73UwfsEeu2xcZ+QQ+zz7XpmcjiqdUfB6QgpU6DscnNJKsuzuxBvMo/6EqHJ48skn0bdvX5tjWVlZLheEK91CKKmJjVC8Lh83nTQRMCKzHW9VOiM+yz6FL/5XBBOHcfyxltH4+LUeTl1bCWJipK0Y+cbtZqVe8jmy8go4J6esvPOWyUnsPYNTwkUnsn2HL2HD9iLcuKVDWIgWAFBVY0DjRsEYm9oGKZ1s298y8g9OCUdEeJDNZ3V621DcDTvO4MqNauw+dNlS0ZY5HhEeZHduXvmsKuLyfV7KvbJh34P1PX/Ok7/BfLdSv0d348g4OIq3joG74FUOe/bsgdFoxKJFi0DTtGXyNhqNWLp0KYYOHepSQe69914cOnTI8ndpaSliY2Mlf96RJLgohbpxMXApBgA49ns50qd/b/mbAkDDs+YqOUk/fOMWFREo+RxcoZLMceYcUt4jBHvnYe1XKrulw8rNR1F5p9Yy3kwfbetV6cJXkwDUJxLeqbH9PvUGE3YcvGDnHNcbTFifU2gpZCjE+pxCu1Lpcj7PhdhOpKzsDkJ4KgxH/WU284UEMCXLc5AkOAHlUFRUhIMHD6K8vBwbNmy4+wGNBi+99JJrpQTwxBNPYOXKlbh58yaCg4Px448/4qOPPnL5dawZnhzP63NwJ8zVrc0mQH1VSeYH7E2mKlfYfqXY9J21+4slnLGbALFX8NYmLKFIGS4crcgq9/NspJrr9Ab7cVFTzpWWcCfeXryyIcCrHCZNmoRJkybh66+/dktI6b333otp06Zh7NixMBgMGDFiBNq1a6foNZmHyLpsgqdhOnax+ykwpiprBzlffLfSuML2K0XBOKuEpEywzHu4+mhbKw8h564zvYPFFKDc1bEUfxCfry04SOMzE6snilf6G6I+B3fmGqSnpyM9Pd1t1wPsyyYUFJbYrNg9gVADFmuYSYlxkH+WfQqBWvtS1wDQu0MTPN8/wWUyuiPU0VklJLVxECC+gufrB9G6eSTOXal0WIEJKUBHVsdSdiJ875H63HkDSoUaE+7ieLpiA0Vs0uOq7RIWrEG1zsjb/9edcCkGANj721X8dOwqTGbvCcWVWr/JUTnFQkCldKBjlAdfn4jSWzqnegcLKUCuJDax1bHQDodpEdoQ8gWcvQdva7bljRDlIBOhyWrjD6dt2oNq1JRsf4ZQRU5nMf01z5RX6rE2p0iw4ijgPUrEUeSEgHL10bZWHkIrVVfsorg+78jqmE8hmum7PaEbQr6AM/dA/BXS4FUOFRUVgh9k8h8Id3m+f4Kg6YatPBjY0UoAbGrhK4FZQqV2a3MVmybRwfh4QpISorkUqRN3UqJtH22pCWVKrrYdvaZWQ4GjXqBl18EUI/SGlbOjK3hnTI7EXyENXuXQrVs3UBTFm39QVFSkqGANETHlwebL7UWcrR69gavlOry6eA9YPlyf3m1Y99Fm44nVttxrSukVwSgbbyhK5+wK3tF7IP4KafAqh9OnT/O9RHADzINfUFjCm0znadiKARDebQD1hevGDmjt8YlJLp7IzpV7TSm9IrzJr+CpFXxD8Lm4A1GfQ11dHXJzc1FdXQ0AMJlMuHjxIqZNm6a4cATx1ZG1qYoCEMATreQt6A0mn61Z5Y7VNpeZRWpPCrGVr7f5FTy1gm8IPhd3IKocpk2bhkuXLqGsrAxt27bFsWPH0KVLF3fIRpCAkKmKq8dylc6IQK3aLqbfGxDadTBtO71VcbgCZ80sSleqdTWeWsGTGk3SEFUORUVF+PHHHzFv3jy8+OKLoGkaH3zwgTtkIziJWGQVV9y+t1Jda7vjUKsoG1NboJbC2AG+XWLdWTOL0pVqXY0nV/De4HPxdkSVQ2xsLDQaDR588EEUFxcjNTUVOh13zRuC78DecbBryVjvOrjQqrl9Du6C7YPRG/hLrLs6AVApnDWz+NqK2Nfk9TdElUNISAiys7ORkJCAzZs3o0WLFqJhrgTfR8rKavE3R1B0wfufhb2/XbXbJXlTrSoGV5hZfG1F7Gvy+hOiymH27NnYsmULZsyYga1bt+K5554jzmgCAGDGqI42f0spPUJRd0ude5IqnRHrthXh98sV+PV0qU3pCE8pDuIoJXgTFM2VyOCDOFKym8EfyvOK4c4x8AV/h4qiEByoQnWtyWLuyDt+1Wan1OaBSDsF6SzeUNaB/B78YwzESnaLKoeJEydyHv/nP//pnGQuhigH5/D0GNRH6nhv0p8YDcle7ulnwRvwhzFwuJ8DQ//+/S3/NxgM2LNnD1q3bu0a6QiEvxCyPXP5NtjRSp6GCTv9/XIFjp8rJw5Wgs8j26xkNBrx3HPP4d///rdSMjkE2Tk4h6+OAV+9Km+BL5TU2nwUFqwBTdM2JixPKhRffRZciT+MgdM7BzZmsxmlpaVOCUUguAqxelV8IblhwRp0TojF/uPXFO0EyJWnwE52s3aGcyUC+kooLqFhIaoc2D6H4uJikiFN8BnEQiVb3h+Jb3YV20zQgVo1jEaTy6rishWTlBpI1jDO++f7J3iFw5rgH8jyOVAUhVGjRqFHjx6KCkUguAtr5WFtSuCahNnRSlJh5yk4Ujso9+hVtLw/kvQhILgNUeUwbNgwd8hBIHgVXDsO5m+uboDREYFoFx+N/BMlonkKUtqXsjHTjpfX4JOX7DoIQvAqh4SEBFAUxftB0s+B4K8Imapa3h8pavYRa1/KhYpyrLwG279h3XecvetglMjNSj2iiPLwe3iVQ0FBAWiaxvLly9G0aVOMHDkSarUamZmZuHrVuxOYCARPIbUvNnC3ppCUKrnJ7ZtYQmTZCJXXEPJvWO86SOtMAhte5dCoUSMAwMmTJ22qsI4dOxbDhw9XXjICoQHDViIFhSV2jnEGJlqJq9ObWHkNMfMV87ojJiu2X6ZdfDTJ8WhAiPocdDod/vjjD7Ro0QIAcObMGRgMBsUFIxD8CUd2HFImYDH/BrPrkGuy4tppWJdEYe88SJSV7yGqHN58802MHDkSrVu3htlsxrlz57BkyRJ3yEYgEFjIrWIq5N+w3nXIrQgrJRyX2XkAICYrH0RUOfTr1w+dOnXC4cOHQVEUOnXqhKioKHfIRiAQnIS92+CLVpJbEVZqtFV5pd5jvaIJzsGrHL7//nsMGTIE69atszl++fJlAMCLL76orGQEAsElyDVZSYlWkhqOK/Q+pXtFE5yDVzlcuHABQH1GNIFAaPgwSkRKXSEp4bjMzoOvo6BSvaKJf8M18CqHKVOmAAAGDhyInj17uk0gAoHg/XA5x4WildzVxIiE5LoOUZ9DRkYG5s6di+HDh+OZZ57Bvffe69QFDx8+jE8//RRGoxGRkZH49NNP0bRpU1RWVuLtt9/GpUuXEBUVhWXLliEmJsapaxEIBOWQ6hx3Z69oV/g3CgpLkJVXgLJbOr/eeUgq2X3u3DlkZmYiJycHCQkJeOaZZ9C3b1+HLtinTx+sWrUKCQkJ2Lp1K3bv3o3Vq1fjww8/RFxcHF555RVkZWVh3759WLZsmeTzkpLdzkHGgIwBgy+Pw/gFe3hf++LdPqKf58ol0agpBGpVXlNS3VWIlexWSTlJfHw8ZsyYgZUrV+LWrVt46623HBKmrq4OU6dORUJCffnh1q1b49q1awCAffv2IT09HQCQlpaGn376ieRTEAgEWfD5MaT6N7h2HkYTbemLzpipCgpLnBPUBxA1K5WXl+O///0vvvvuO5hMJowYMQL/+te/HLpYQEAAhgwZAqC+L0RGRoZlB1JaWmoxI2k0GoSFheHmzZuSzVhCGlAKMTHhTn2+IUDGgIwBg6+OwwtpicjYcsymHEmgVo0X0hIl3dNNCRFUdUYzsvLOY3BKK8ly7Tt8CRu2F+HGLR0aNwrG2NQ2SOnUTPLnPYGkPId+/fph7ty56NSpk+QTb9++HfPnz7c51qJFC6xfvx51dXV49913YTQa8eqrr/KeQ6WStLEBQMxKzkLGgIwBgy+PQ2LzSIwd0NrOv5HYPFLSPUVJDNEtu6WTPEZsU1XZLR1Wbj6Kyju1HjVPOd0JLjc3F2Fh8lflqampSE1NtTteXV2N1157DZGRkVi9ejW0Wi0AIDY2Fjdu3EBcXByMRiOqqqoQGRkp+7oEAsG/kZtFbo3UirlywnB9NQmQVzlMnToVy5cvx6hRozhfz87OduiCM2bMwAMPPIAPP/zQpiR4cnIysrKyMHHiRGzbtg2PP/64RXEQCASCO2Am66y88yi7pUNYsAa6WqNNV0C5Ybi+mgTIqxwmTJgAAJg9e7bLLnbq1Cns3r0bLVu2xNChQwHU7xg+++wzTJ06Fe+++y4GDRqE8PBwUr+JQCB4hKTEOAxOaSXYFVDOil9u3SpvgTeUVaxnQ5MmTRQRyFGIz8E5yBiQMWAg4+DaMeArtT4uNcE3fQ6DBg0CRVGgaRq1tbUIDQ2FWq1GZWUloqOjkZeXp4jABAKB0JBwZxKgK+FVDr/99hsAYM6cOejatSsGDRoEANi9ezd27drlHukIBAKhAeCMk9xTiMaKnjx50qIYAODJJ5/E6dOnFRWKQCAQCJ5FVDmYzWb8/PPPlr9/+uknmygjAoFAIDQ8RPMcZs2ahTfffBNarRZmc71DJSMjQ3HBCAQCgeA5RJXD448/jr1796K4uBgUReHhhx+GRiP6MQKBQCD4MKKzfFlZGb799ltUVFTYHJ81a5ZiQjmCSuWcqcvZzzcEyBiQMWAg49Dwx0Ds/kRLdj/33HMIDw9H27ZtbXwNb7zxhmskJBAIBILXIakq61dffeUOWQgEAoHgJYhGKzVp0gQ1NTXukIVAIBAIXoLoziE2NhZDhw5Fly5dEBQUZDnubT4HAoFAILgOUeXQtGlTNG3a1B2yEAgEAsFLkNRDmkAgEAj+Be/OYdSoUfjmm2/QoUMHzozoI0eOKCoYgUAgEDwH786htLQUsbGxuHLlCucHiamJQCAQGi7ErEQgEAgEO0RDWQkEAoHgf/i1csjOzsbAgQPx1FNP4euvv/a0OB5h7NixGDRoEIYMGYIhQ4bg2LFjnhbJbVRVVSEtLQ2XL18GABw4cADp6eno168fli5d6mHp3Ad7HN577z3069fP8kzs3LnTwxIqS0ZGBgYNGoRBgwZh0aJFAPz3WbCB9lNKSkro3r1707du3aKrq6vp9PR0+uzZs54Wy62YzWa6e/futMFg8LQobufo0aN0WloanZiYSF+6dInW6XR0cnIyffHiRdpgMNDjx4+n9+3b52kxFYc9DjRN02lpafT169c9LJl7yM/Pp0eOHEnr9Xq6rq6OHjt2LJ2dne2XzwIbyTuHAwcOIC0tDX379kVmZqaS+sotHDhwAN26dUNkZCRCQkLQv39/7Nixw9NiuZU//vgDFEVhwoQJGDz4/7d370FVVW8Dx7/IRWTIcRB1GsfSUNSmRK1AxMEwJBSOXAQFCgRREYgxUgYcIZRRwksqGKmUlg1Nmop5ycK8j3pwIC9QCSgBohzFdOKiCHhY7x+87AEPoPkLaGR9/jr7up69zpqzztqXZ8/oUWlSvv/+e+Lj4xk4cCAAubm5vPzyywwZMgQDAwNUKlWPaA+P18ODBw8oLy8nLi4OlUpFSkqKkqr/eTRgwABiYmIwMjLC0NAQCwsLSkpKemRbeFy7t7I2NDRgaGioTKenpyudwsyZM/H09Oz86DpRRUUFAwYMUKYHDhxIbm5uN0bU9aqqqrC1tWX58uU8fPiQgIAAhg0bhp2dXXeH1ulWrVrVarqt9nD79u2uDqvLPV4Pd+/eZcKECSQkJGBiYkJISAh79uxh1qxZ3RRh5xoxYoTyuaSkhMOHD+Pv798j28Lj2h05+Pv7c+rUKWXa2NiYn3/+maNHj2JkZNQlwXUm0cZNWj3tDXfjxo1jzZo1mJiYYGZmhpeXV6vvvCeR7aHJkCFDSE1NpX///vTp00fnd+B5dfXqVebOnUt0dDQvvfSSzvKe2Bba7RzS0tI4ffo04eHhlJWVERsbS1FREXl5eXz66addGWOnGDRoEH/99Zcy3fxcR0+Sk5ODWq1WpoUQPfZFTrI9NCkoKCAzM1OZ7glt4tdffyUwMJDFixfj4eEh28L/a7dz6Nu3L3FxcURERBAfH8+OHTsIDQ0lOjqaoUOHdmGInWPixImo1Wru3btHbW0tR44cwd7evrvD6lLV1dWsWbOGuro6ampq2LdvH1OnTu3usLqFlZUVxcXFlJaWotVqOXToUI9rD9DUGSQmJlJZWUlDQwO7du16rtuERqMhPDycdevW4eLiAsi20KzdvwRarZYzZ85gaGjItm3bOHDgAHPmzCEoKAhnZ+eujLFTDBo0iMjISAICAmhoaMDLy4sxY8Z0d1hdysHBgcuXL+Pu7k5jYyN+fn6MGzeuu8PqFr179yYpKYmIiAjq6uqYPHnyc9HO/6lRo0axYMECfH19efToEU5OTri6unZ3WJ1m27Zt1NXVkZSUpMzz8fGRbYEOnpD+4IMP6N+/P7W1tZiYmLB8+XJqampITU3l6tWrfPnll10dqyRJktRF2h05lJWV8dlnnwHg7u4OgKmpKdHR0RQVFXVNdJIkSVK3aLdzePXVV5k/fz51dXW8/fbbrZZZWFh0dlySJElSN+ow8V5BQQFGRkYMGzasK2OSJEmSupnMyipJkiTp6NGJ9yRJkqS2yc5BApoycc6YMYOwsDAaGhoAqKysZPbs2dTX13dbXCNHjuTevXs68zMyMggJCQFg2bJlnDt3rqtDA5pu3IiIiPjH23333XekpaV1uM758+fbvY00NzeXjz/+uN1t586dq9Tb/PnzuXbt2jOV055nPe72xMbG8ttvv/1r+5P+d7JzkMjPz6eiooIDBw5gbm7OmTNnAEhOTmbhwoX/+XQpq1atYuLEid1Sdnl5OcXFxf94O19fXxYsWPDM5V67dq3DfD9nz55VPn/xxRcMHz78mctqy7Med3vOnTvXZgoTqfvIzkHCyMiI+vp6hBDU1tZiaGhIfn4+Go0GBweHdre7ceMGkydPZu7cubz77rtUVFRw4cIF/Pz88PDwwNPTkxMnTgBN//SDg4MJCgpi+vTpBAUFKT9u/v7+rbJePj69ceNGPDw8cHNzU/bXUsv1T5w4gZubGyqVitmzZ5Ofnw/Ali1b8PLyQqVS4ejoqLyjYNOmTcTExBAcHIyzszN+fn5t/ugWFRXh4+ODp6cnHh4efPvtt2i1WmJjY7l+/TrBwcFt1kdH5SYkJABNowBPT09UKhXh4eF4eHhw/vx5oClLamRkJG5ubjg7O5OTk4NGoyElJYWcnByWLl2qE2vzvDlz5qDRaJgyZQp5eXkA7NmzBxcXF1QqFQEBAWg0mlbb5uTk4ODgoLwj/vjx43h7e+Pu7o6Pjw8XL17UOe6WiouLsbGxUUabWq0We3t7rl27xq1bt1i4cCEqlQpXV1flWakNGzZQUVHBkiVLuHz5MtXV1cTExCh1kpiYyKNHjwBISUlBpVLh6elJcHAwFRUVOscv/Uu6JVG49J+zfv16MWPGDBEXFye0Wq0IDAwUJSUlHW5TVlYmLC0tRXZ2thBCiL///ls4OTkp7wW4deuWsLe3Fzdv3hR79+4VY8eOFX/++acQQoi1a9eKiIgIIYQQ77//vvjpp5+U/bactrS0FFu3bhVCCFFQUCCsra3F3bt3xd69e8WCBQtarX/nzh3xxhtviD/++EMIIURmZqYIDg4WN27cEP7+/qK2tlYIIcShQ4eEq6urEEKIlJQU8c4774jq6mohhBAhISEiOTlZ51iXLl2qxFFRUSE+/PBDodVqRVZWlnBxcWmzPp5U7ooVK0RDQ4Owt7dX3hegVqvFyJEjRVZWlsjKyhKjR48Wly5dEkII8dVXX4mAgAAhhGh1/G2xtLQUd+/eFUII4eDgIHJzc8WVK1eEjY2NKC8vV/YXFxenHINarRaOjo7iypUrQgghiouLhaurq7h3754QQojCwkJhZ2cn7t+/3+q4H/fee+8p39/JkyeFj4+PMn/79u1CCCGqqqqESqUShw4dahWjEELExMSIb775RgghxKNHj8SSJUtEWlqaKC8vF+PHjxd1dXVCCCG2bdsmfvnll3brQPrfPN8ZtaSnFhkZSWRkJAD79+/n9ddfx9TUlMjISGpra/H3928zlbeBgQFjx44F4NKlS9y5c4fw8HBluZ6eHgUFBQDY2dkpt0XPmjULNze3p4rN19cXAEtLSywsLLh48WKb6124cIERI0YwevRoAJycnHBycgJg9erVHDx4kNLSUi5fvsz9+/eV7aytrTE1NQWanu+prKzU2ffUqVOJjo4mNzcXW1tbYmNj6dVLd+Ddsj4GDx7cYbkAhYWFAEyePBmACRMmtEojPWTIEKysrICm1BZ79+59UnW1S61WM2nSJF588UUAAgMDgaZrDs3/6n19fRk1ahTQdGqqoqJCWQ+avs/r1693WI63tzf79u3D2dmZjIwMvL29efDgARcuXGD79u0AvPDCC3h6enL69Gklp1GzkydPkpeXx549ewB4+PAh0JTyZtSoUXh4eGBvb4+9vT22trbPXB9Sx2TnILVSU1NDeno6O3bsYMuWLUpemZkzZ/Ljjz/qrG9kZKRk7dRqtVhYWLB7925l+e3btzEzM+PgwYPo6+sr8xsbG1t2XxMWAAAESUlEQVRNixbnm5sviDdr+SMsOsgSqq+v3yq1shCCgoICtFotYWFhBAYGYmdnx1tvvcWKFSuU9YyNjZXPenp6bZ77dnBwIDMzk3PnzqFWq0lNTWXnzp0d1sfvv//eYbnNMT9eXst6aflOlfZiS05O5vjx4wBMmTKFRYsWPVX9PHz4kJs3byrL0tLSCAsLY9q0aYwZM4bGxkZsbW3ZuHGjso1Go2HgwIHk5OS0WQaAs7Mzn3zyCUVFRWRnZ5OUlIRWq9WJvbGxUTld9Pj85ORk5WHbqqoq9PT06NWrF+np6eTl5aFWq0lMTMTGxobY2Nh2Y5GenbzmILWSmppKUFAQJiYm1NfXY2hoSK9evaitrX3itmPHjqW0tJTs7GwArly5opx7B8jKylLO5+/cuVO5nmFmZqbcqXL9+nVlpNFs3759QNOPbWlpqfJP+nFWVlYUFRVx9epVAI4dO0ZUVBTZ2dm89tprBAUFYW1tzbFjx9Bqtf+oXhYvXszhw4dxcXEhPj4eU1NTNBoN+vr6Op1Zs6cp18LCAiMjI06fPg00XX8oLCx84vsD9PX1lR/WRYsWsX//fvbv3690DC2XN7OxsUGtVivfx86dO1m7di3Q9Ea08ePHEx0dTVRUFLW1tUyYMIGzZ88q6XJOnTrFjBkzqKur6/C4e/fujYuLCzExMTg5OdGnTx9MTU2xsrJS3tVeXV3NDz/8oNxI0DLeSZMm8fXXXyOEoL6+ntDQUNLT08nPz8fV1RULCwtCQkIIDAzUaSvSv0eOHCRFUVERhYWFREdHAzB79mwWL17M5s2bCQ0NfeL2ZmZmpKSkKGnAhRCsWbOGwYMHA02nBaKiorhz5w7Dhw9XLsiGhoYSExPDqVOneOWVV3jzzTdb7besrAx3d3f09PRYv349/fr1a7N8c3Nz1q1bR3R0NFqtFlNTUzZs2EC/fv04cuQI06dPx9DQEFtbWyorK6mpqXnqugkLC2PZsmXs2rULfX19HB0dsba2pqqqCn19fby8vHReRO/q6vrEcg0MDNi0aRPx8fGsX7+eoUOHYm5ujrGxcYcd8rhx49i4cSPh4eGkpqbqLJ86dSp+fn58/vnnyryRI0cSFRXFvHnzgKYOITExkZKSEmUdDw8PMjMzSUpKYsWKFSQkJPDRRx8pI7bNmzdjYmLCiBEjlOPevXu3Tmfm7e1Neno6y5cvV+atW7eOhIQEMjIyqK+vVy4sAzg6OhIZGcnKlStZtmwZq1atQqVS0dDQwMSJE5k3bx6GhoZMmzaNmTNnYmJigrGxsRw1dCL5hLTUJTIyMsjMzGTr1q3dHcp/zurVqwkODsbc3ByNRoObmxtHjx6lb9++3R2a1IPJkYMkdbPBgwcTGBiIgYEBQghWrlwpOwap28mRgyRJkqRDXpCWJEmSdMjOQZIkSdIhOwdJkiRJh+wcJEmSJB2yc5AkSZJ0yM5BkiRJ0vF/ljN/0xpF/GQAAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"import matplotlib.pyplot as plt\n",
"import seaborn as sns\n",
"sns.set()\n",
"\n",
"diffs = split_trump_percentages - r_percentages\n",
"plt.scatter(r_percentages*100, diffs*100)\n",
"plt.title('Wayne Replication')\n",
"plt.xlabel('% republican straight-ticket votes')\n",
"plt.ylabel('% individual trump votes - % straight ticket votes')\n",
"\n",
"plt.savefig('wayne')\n"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"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.7.7"
}
},
"nbformat": 4,
"nbformat_minor": 4
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment