{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "## Chapter 7\n", "### This notebook and the figures below are made by Yeseul Lee, under the guidance of Steven Skiena." ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": true }, "outputs": [], "source": [ "from scipy import stats\n", "import pandas as pd\n", "import matplotlib.pyplot as plt\n", "import numpy as np\n", "from sklearn.metrics import roc_curve\n", "%matplotlib inline" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Figure 7.3 and 7.6 Left\n", "#### Data Source: https://www.statcrunch.com/app/index.php?dataid=1406047" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": false, "scrolled": false }, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYAAAAD7CAYAAABjVUMJAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAGTFJREFUeJzt3X2QHPV95/H3Z0UUxxACCg4ULCsBesCABcK2bM5YzIVg\nlJTIJiGAiCsnEuoQdrhcLjnAplxo11Xc2aok9uUoF3EibJk4p0tkDj1gsMoHE1UKActJrBDoYWNA\n0kqAA4KEhwiE9L0/ulf0jmZ3Zndnd3q2P6+qqe3+9W9mvts729/+PXSPIgIzMyuetmYHYGZmzeEE\nYGZWUE4AZmYF5QRgZlZQTgBmZgXlBGBmVlDHNTsAAEmei2pmNgoRodE+NzctgIho2ceyZcuaHsNk\njL/ez0Ve42/1/T/ZY58M8Y9VbhKAmZlNLCcAM7OCcgJogFKp1OwQxsTxN1crx9/KsUPrxz9WakQ/\n0piDkCIPcVi+SGpIP6fZZJX+j7T+ILCZmU0sJwAzs4JyAjAzKygnADOzgnICMDMrKCcAM7OCqisB\nSFooaYekXZJur7J9qqRVkvokbZLUkdk2V9JjkrZJ6pU0tZG/gJmZjU7NBCCpDbgbuBI4H7he0rkV\n1W4EDkTELOCbwPL0uVOA+4CbIuICoAQcalj0ZmY2avW0AOYDfRGxOyIOAauAzoo6ncDKdHk18Mvp\n8ueA3ojYBhARr/uKLzOzfKgnAZwB7M2s96dlVetExGHgXyRNA2YDSHpY0lOSbh17yGZm1gjj9X0A\nA5cmHwd8BvgEcBD4v5KeiohHx+l9zcysTvUkgH1AR2a9PS3L6gfOBPan/f4nRsQBSf3Axoh4HUDS\nD4GLgWMSQFdX19HlUqlU+Js0mZlVKpfLlMvlhr1ezZvBpQf0ncDlwEvAk8D1EbE9U+eLwAUR8UVJ\ni4HfiIjFkk4CfgxcCrwPPAT8eUQ8VPEeHhqwY/hmcGbDG+vN4Gq2ACLisKRbgA0kYwYrImK7pG6g\nJyLWAyuA+yT1Aa8Bi9PnviHpz4GngCPAg5UHfzMzaw7fDtpyyy0As+H5dtBmZjYqTgBmZgXlBGBm\nVlDjdR2ATUbf/nbtOjfdNP5xmFlDuAVgZlZQTgBmZgXlBGBmVlAeAzBrgGrDIx4OsbxzC8DMrKCc\nAMzMCsoJwMysoJwAzMwKygnAzKygnADMzArKCcDMrKCcAMzMCsoJwMysoJwAzMwKyreCsOrqufWz\nmbU0twDMzArKCcDMrKDcBVQ0vm3lUd4VVnRuAZiZFZQTgJlZQTkBmJkVlBOAmVlB1TUILGkh8E2S\nhLEiIr5esX0q8D3g48CrwHURsUfSdGA7sCOt+nhEfLFRwVuDeM6/WSHVTACS2oC7gcuB/UCPpDUR\nsSNT7UbgQETMknQdsBxYnG77p4i4uMFxmzWVc6ZNBvV0Ac0H+iJid0QcAlYBnRV1OoGV6fJqkmQx\nQGOO0szMGq6eBHAGsDez3p+WVa0TEYeBNyRNS7fNkPT/JD0q6dKxBmxmZo0xXheCDZz1vwR0RMTr\nki4GHpB0XkS8VfmErq6uo8ulUolSqTROoZmZtaZyuUy5XG7Y69WTAPYBHZn19rQsqx84E9gvaQpw\nYkQcSLe9BxARmyX9BJgNbK58k2wCMJsMfKWxNVrlyXF3d/eYXq+eLqAeYKak6elsn8XA2oo664Al\n6fI1wCMAkk5JB5GRdDYwE3h+TBGbmVlD1GwBRMRhSbcAG/hgGuh2Sd1AT0SsB1YA90nqA17jgxlA\nC4CvSnoPOAIsjYg3xuMXsZzwaa9Zy6hrDCAiHgbmVJQtyyy/C1xb5Xn3A/ePMUYzMxsHvhuoFYLn\n7ZsdywnArAYnD5usnAAmMx+5zGwYTgA2KTn3mdXmBGCW4cRhReLbQZuZFZRbANYclafavlbAbMK5\nBWBmVlBOAGZmBeUEYGZWUE4AZmYF5QRgZlZQTgBmZgXlBGBmVlBOAGZmBeUEYGZWUE4AZmYF5QRg\nZlZQTgBmZgXlBGBmVlC+G+hkkteb2Tcwrmov5RuJmo2OWwBmZgXlFoDZBHILxvLELQAzs4JyAjAz\nK6i6EoCkhZJ2SNol6fYq26dKWiWpT9ImSR0V2zskvSnpjxsVuJmZjU3NBCCpDbgbuBI4H7he0rkV\n1W4EDkTELOCbwPKK7X8G/HDs4ZqZWaPUMwg8H+iLiN0AklYBncCOTJ1OYFm6vJokYZDW7wSeB95u\nRMCWyuuUTzNrGfUkgDOAvZn1fpKkULVORByW9IakacC7wG3AFcCtYw+3oHywN7NxMF7TQJX+7AK+\nERHvSMqWH6Orq+vocqlUolQqjVNoZmatqVwuUy6XG/Z69SSAfUB2ULc9LcvqB84E9kuaApwYEQck\nfQq4WtJy4GTgsKR/i4hvVb5JNgGYmdmxKk+Ou7u7x/R69SSAHmCmpOnAS8Bi4PqKOuuAJcATwDXA\nIwARsWCggqRlwJvVDv5mQ3Hvl9n4qZkA0j79W4ANJLOGVkTEdkndQE9ErAdWAPdJ6gNeI0kSZlaH\nyiTnK4NtotQ1BhARDwNzKsqWZZbfBa6t8Rpja6uYmVlD+UpgM7OCcgIwMyso3w3U8sGjvWYTzi0A\nM7OCcgIwMysodwFZrp27cXDX0I4FniNp1ihuAZiZFZQTgJlZQTkBmJkVlBOAmVlBOQGYmRWUE4CZ\nWUF5GqjlwsaNzY7ArHjcAjAzKygnADOzgnIXUN74pmjDqrwyGPA3qJiNklsAZmYF5QRgZlZQTgBm\nZgXlBGBmVlBOAGZmBeUEYGZWUE4AZmYF5esAms3z/s2sSZwAJpIP9maWI3UlAEkLgW+SdBmtiIiv\nV2yfCnwP+DjwKnBdROyR9Ekge9TrjogHGhK5Wcp51Wx0aiYASW3A3cDlwH6gR9KaiNiRqXYjcCAi\nZkm6DlgOLAaeAT4eEUcknQb0SlobEUca/puYTWL1JjnfFcNGop5B4PlAX0TsjohDwCqgs6JOJ7Ay\nXV5NkiyIiIOZg/3PAT7wm5nlRD0J4Axgb2a9Py2rWiciDgNvSJoGIGm+pG1AL3Czz/7NzPJhvAaB\nNbAQEU8CF0iaA3xP0kMR8d44va9VmHHHHew5cKBmvenTpnHDJZfQ/eCDQPIH3HrnnZx/+ulH65x2\n66389M03EXD4nnvGKeL8ueOOGRw4sKdmvWnTpnPJJTfw4IPdg8qlNo4/fhpnnfVprrzyNmbOvHS8\nQjUbkXoSwD6gI7PenpZl9QNnAvslTQFOjIhBR52I2CnpLeACYHPlm3R1dR1dLpVKlEqlOkKzWtqk\nD7JxjXoDBpaWrVvH6qVLB5XX81qTTTIMVvs3T+odXTu6FBG89dZrPPPMerZt+yFLl/6Aiy6q7EU1\nq61cLlMulxv2evUkgB5gpqTpwEskg7vXV9RZBywBngCuAR4BkDQD2BsRh9PnzwFerPYm2QRgjfP8\nXXcNWm+7+WaAqmfx3evWHV0O4IGnn2Zrfz9z29sHlRctCdx11/OD1m++eeBAL+655/CgbevWfXD2\nv2jRMhYtupODB9/i/vtvY+PGe4gIVq/+EycAG5XKk+Pu7u6hK9eh5hhA2qd/C7ABeBZYFRHbJXVL\nWpRWWwGcIqkP+CPgS2n5pSQzfzYDPwC+UNkysHw6rq2NiGBZJinY6HzoQyfQ2TmQiINXX32Bt9/2\nv4E1X11jABHxMMnZe7ZsWWb5XeDaKs/7G+BvxhijNcHn589n5eOPs7a3ly179jCvo6P2k2xInvtg\neeR7AVlVn5wxg6vmziWAO9eubXY4Le3gwTdZs+YrR9c/8pFzOP74aU2MyCzhW0HYkLoWLWLd1q38\ncNs2el58saGvvXFjQ18ul9av72L9+q5BZVIbV1+9vDkBmVVwArAhzevooPPCC1nb2+tWwKgkw+WS\n+PCHT+bssy/hiiv+K7NnLxi3d6x2xbCvDrahOAHYsLqvuoq1vb1seO65QVNFrbaBWUBmeeUxABvW\n3PZ2fnPePAI4EtHscMysgdwCsJq6r7qK/7NlCz78Twzf3dQmilsABTNwNa+G6M6pVnr+6adz7Sc+\nUfO5xZHsiaH3Q9H3j7UKRQ6a9ZIiD3GMO5/aHVXPLKDLvr+Uf/j8X9ast2OBRzmH40HgyUsSETHq\nMw63AMzMCspjAOPFZ/tmlnNOANbyzt14bLJ1t5BZbe4CMjMrKCcAM7OCcheQjbsi3PfHrBW5BWBm\nVlBOAGZmBeUuILNJrnJGsi8MswFuAZiZFZQTgJlZQTkBmJkVlBOAmVlBOQGYmRWUE4CZWUE5AZiZ\nFZQTgJlZQTkBmJkVVF0JQNJCSTsk7ZJ0e5XtUyWtktQnaZOkjrT8VyQ9JalXUo+kf9/oX8DMzEan\n5q0gJLUBdwOXA/uBHklrImJHptqNwIGImCXpOmA5sBj4Z2BRRLws6XzgR0B7o38Js1r8pTFmx6qn\nBTAf6IuI3RFxCFgFdFbU6QRWpsurSZIFEdEbES+ny88CH5L0Mw2J3MzMxqSeBHAGsDez3p+WVa0T\nEYeBNyRNy1aQ9NvA5jSJmJlZk43X3UA1aCXp/vnvwBVDPaGrq+vocqlUolQqjVNoZmatqVwuUy6X\nG/Z69SSAfUBHZr09LcvqB84E9kuaApwYEQcAJLUD9wO/GxEvDvUm2QRgNlbV+vzNWl3lyXF3d/eY\nXq+eLqAeYKak6ZKmkgzurq2osw5Yki5fAzwCIOkkYD1we0Q8PqZIzcysoWomgLRP/xZgA/AssCoi\ntkvqlrQorbYCOEVSH/BHwJfS8j8AzgHulLRF0mZJpzT8tzAzsxGrawwgIh4G5lSULcssvwtcW+V5\ndwF3jTFGy6lqX/a+YMHEx2Fmo+Mrgc3MCsoJwMysoPyl8NZQ1bqFzCyf3AIwMysoJwAzs4JyAjAz\nKygnADOzgnICMDMrKM8Csrp5ho/Z5OIWgJlZQTkBmJkVlBOAmVlBeQygEb7te89ba6v2Eb7JX5k8\n6bkFYGZWUE4AZmYF5QRgZlZQTgBmZgXlBGBmVlCeBWRmVXlm0OTnBGBWMJ61bAOcAEbK/z1mNkl4\nDMDMrKDcArDCOnfjsa25HQvcyW3F4RaAmVlBOQGYmRVUXQlA0kJJOyTtknR7le1TJa2S1Cdpk6SO\ntHyapEckvSnpLxodvJmZjV7NMQBJbcDdwOXAfqBH0pqI2JGpdiNwICJmSboOWA4sBg4CXwEuSB+t\nZ5LN+qn2rV4LFkx8HGbWfPW0AOYDfRGxOyIOAauAzoo6ncDKdHk1SbIgIt6JiMeAdxsUr5mZNUg9\ns4DOAPZm1vtJkkLVOhFxWNIbkqZFxIHGhGl58d2t6/juMw8OuX3h2ZfwpUuW8J2t61iZ1pPg3l+7\nk7NOOv1ovd/8wa28fvBNJHj0d+6ZkDgGtEn8/NTjOe+Us7j+vCuZ+0szR/z+ZpPBeE0D1Uif0NXV\ndXS5VCpRKpUaGI4NZzTdQhriL6x0gyrqfWfrOr66YGldrzESI40DIAj+9b232LTvGZ7Yv42vfnYp\nl5550diDMRtn5XKZcrncsNerJwHsAzoy6+1pWVY/cCawX9IU4MSRnv1nE4C1hiUfW8QNH1tUs14E\n/GP/0/zk9X7OObl9UHkjkkC9cQzUe+fQQe7Zcj9r+zYSBN/avNoJwFpC5clxd3f3mF6vnjGAHmCm\npOmSppIM7q6tqLMOWJIuXwM8UuV1GvCvbq1oSlsbQXDv1nXNDgWAD//Mh/iPFybDWBHw0tuv8q/v\nvt3kqMwmXs0WQNqnfwuwgSRhrIiI7ZK6gZ6IWA+sAO6T1Ae8RpIkAJD0AvDzwFRJncDnKmYQWQ5V\n6xYarV+ZMZ8fPf84j+3rZdeBPcye1lH7SePsCNHsEMyarq4xgIh4GJhTUbYss/wucO0Qzz1rLAFa\nfn1363q+u3X9oLK7LvsCn2m/cFDZub84gzffe4fH+rdy79a1fK10S1PiGPDOoYP8de+ao+unn/AR\nTvzZ4xsak1krKO69gHyz8zGr7L/XML18N3xsEY/1b+WJ/dvY/tqLTYmjWqJok7h53tXDvr7vGfSB\nyn8b/8u0tuImABuzegdfAWZP6+Az7Rfy2L5e7u2tHEKamDgGEoUQJ0z9MOefcjaLP3oFF546u6Hx\nmLUKJwCbML8/9yoe29fLUy8/N2xrYbyMJGGZFYETgE2Yc05u57Pt89i4dwson4Ow1bp7zCYr3w3U\nJtTvzb0KKZl+aWbN5RaAjVi9F29Vq3fWSadT6vgE5T1PJXXG0BU0ljisMTyXorUpcnAqJikmPI56\nP7kFuBtoXl32/aX8w+f/stlhVFXUWUD1cAKYOJKIiFGf4rgFkDXJDvY2cTxV1FqRxwDMzArKCcDM\nrKCcAMzMCsoJwMysoDwIPIm10owfM5t4bgGYmRWUWwBm1lC+OKx1FCMBeH6/NZjvGWSTweRLAD7Y\nm5nVxWMAZmYF5QRgZlZQk68LqMA87dPMRqK1EoCnF1gL8Q3iLO9aKwGY2aTh87nm8xiAmVlBuQUw\ngar10S9YMPFxWPPUc/1AkbuJ6pnF7VZC4zgB5IwHcm0y8uU5+VRXF5CkhZJ2SNol6fYq26dKWiWp\nT9ImSR2ZbV9Oy7dL+lwjgzczs9Gr2QKQ1AbcDVwO7Ad6JK2JiB2ZajcCByJilqTrgOXAYknnAdcC\nHwXagR9LmlXXFwC30ClDeedOSnPmNDuMUdvyyk7mner4m6WV49+5s8ycOaVmhzHqAeVyuUypVGp4\nPK2ini6g+UBfROwGkLQK6ASyCaATWJYurwb+Z7r868CqiHgfeFFSX/p6TzQg9kQOEkV5165RJ4A8\ndPk8/cqulj0AgeNvpl278pEARssJoLYzgL2Z9X6Sg3jVOhFxWNK/SJqWlm/K1NuXlh1r06bB6++/\nD8cNDi+vg6i7dx8bWx7issmj3msK6qlX743smjEYXRlbtRjGcs6Xg/PFXBmvQWCN+BkvvTR4/ciR\nBoViZmZVRcSwD+DTwMOZ9S8Bt1fUeQj4VLo8BfhptbrAwwP1Kp4ffvjhhx9+jPxR6xg+3KOeFkAP\nMFPSdOAlYDFwfUWddcASkr79a4BH0vK1wPclfYOk62cm8GTlG0TEyFsMZmY2JjUTQNqnfwuwgWTa\n6IqI2C6pG+iJiPXACuC+dJD3NZIkQUQ8J+nvgOeAQ8AX65oBZGZm404+HpuZFdOE3AtI0gpJr0ja\nWmXbn0g6ks4aGij7i/TisaclXTQRMQ6nWvySlknql7Q5fSzMbMvVxW9D7X9J/ymN8RlJX8uU5yb+\nIfb9qsx+f0HS5sy23MSexlMt/gvTCya3SHpS0icz21rhsz9X0mOSeiWtkXRCZlve9n+7pEckPZt+\nzv8wLT9Z0gZJOyX9SNIvZJ6Ti7/BMLH/tqRtkg5LurjiOSPb/2MZQKj3AVwKXARsrShvJxkYfgGY\nlpb9KvBguvwp4PGJiHGk8ZNc9/DHVep+FNhC0r02A/gn0pZWzuIvkXTrHZeun5LH+If67GS2/ynw\nlTzGPsy+/xHwuXT5V4FH0+Vfa5HP/pPApenyDcBX0+Xzcrj/TwMuSpdPAHYC5wJfB25Ly28Hvpb5\ne+TibzBM7HOAWSRjrRdn6o/48z8hLYCI+Efg9SqbvgHcWlHWCXwvfd4TwC9IOnV8IxzeMPFXG7zu\nJL34LSJeBAYufmuaIeL/AsmH/v20zqtpea7iH2bfD7gW+Nt0OVexw5DxHwEGzjhPIrk+BpILJ1vh\nsz8rLQf4MXB1unz0ws8c7f+XI+LpdPktYDvJiWcnsDKttjJdhxwdf4aI/YyI2BkRfRx7/Bnx579p\nt4OW9OvA3oh4pmJT5YVnQ1881nx/kDYT/zrThGyV+GcDCyQ9LulRSR9Py1slfiR9Fng5Ip5Pi1ol\n9v8C/KmkPSS3TflyWt4q8T+b/v9CkoDb0+Vcxy9pBklr5nHg1Ih4BZIDLTBwkM/l75CJfbi7KIw4\n9qYkAEk/B9zBB7ePaEXfAs6JiIuAl4E/a3I8I3UccHJEfBq4Dfj7JsczGtcD/6vZQYzCF4D/HBEd\nJMng3ibHM1K/T3Ly0wMcD7zX5HhqSscpVpPs97dI5tBn5XY2TJXYG6ZZLYBzSPqoeiW9QHIGsVnS\nL5FkrTMzddv5oImcGxHxz5F2vAF/xQdNrZaIn+RM4X6AiOgBDkv6RZJYOzL1chm/pCnAbwH/O1Pc\nKvt+SUQ8ABARq4GBQeCWiD8idkXElRHxSWAV8JN0Uy7jl3QcyQH0vohYkxa/MtC1I+k04Kdpea5+\nhyFiH8qIY5/IBKD0QURsi4jTIuLsiDiL5P5C8yLipyQXj/0HAEmfBt4YaKo12dH44eiHZsBvAdvS\n5bUkd0KdKukshrj4rQkGxQ88APwygKTZwNSIeI0k/utyFn9l7ABXANsjYn+mrFX2/T5JlwFIupyk\nrxZa57P/kfRnG/AV4J50U173/73AcxHxPzJla0kGsEl/rsmU5+lvUC32rOznauT7f4JGs/+W5FbS\n7wJ7gN+r2P486SygdP1ukhHsXjKj3M16VIufZKBoK/A0ycH01Ez9L6fxbyed7ZHD+I8D7gOeAZ4C\nLstj/EN9doDvADdVqZ+b2IfZ9/8u3edbSG6WOC9TvxU++39IMiNlB/Dfcr7/PwMcTv9PtwCbgYXA\nNJIB7J0ks+FOytvfYJjYf4OkBf9vJHdneGi0+98XgpmZFZS/FN7MrKCcAMzMCsoJwMysoJwAzMwK\nygnAzKygnADMzArKCcDMrKCcAMzMCur/AwxLFv6XR+0/AAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "path = \"/home/yeseul/practice-notebook/weight.xls\"\n", "data = pd.ExcelFile(path).parse('Sheet1')\n", "men = data[data[\"Gender\"]==1]\n", "women = data[data[\"Gender\"]==0]\n", "\n", "bsize = 75\n", "n, bins, patches = plt.hist(men[\"Height\"], bins=range(140,max(men[\"Height\"]),1), normed=True, \\\n", " color=\"blue\", alpha=0.4, linewidth=0)\n", "\n", "plt.hist(women[\"Height\"], bins=range(140,max(men[\"Height\"]),1), normed=True, \\\n", " color=\"red\", alpha=0.4, linewidth=0)\n", "\n", "plt.ylim(0,0.065)\n", "plt.vlines(168, 0, 0.08)\n", "\n", "plt.text(158, 0.027, \"TN\", fontsize=16, fontweight=\"bold\")\n", "plt.text(162, 0.007, \"FN\", fontsize=16, fontweight=\"bold\")\n", "plt.text(169, 0.007, \"FP\", fontsize=16, fontweight=\"bold\")\n", "plt.text(175, 0.027, \"TP\", fontsize=16, fontweight=\"bold\")" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### For Figure 7.6, need data that's not normalized." ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "collapsed": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX0AAAEACAYAAABfxaZOAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAFIRJREFUeJzt3X+wZGWd3/H3BwkVa80isHEoGRVBQDa1LsuWSkV37KiL\nwCbM/siyq9nIj61yy99Zt1yBbBV3kqoVrLWIZmvLmEUKLYEV3TikYiKhoJclK6Dhp/yYGSMCjmHU\nDVhlkjIOfPNHnxnanjt37nT3ne6+z/tVdYvTzzl9+suh+fTTz3nO6VQVkqQ2HDbrAiRJh46hL0kN\nMfQlqSGGviQ1xNCXpIYY+pLUkAOGfpKrkuxKcv9I+3uTPJzkgSSXD7VfkmRHt+7MtShakjSew1ex\nzdXAvwU+vachSQ/4J8DPVdXuJD/TtZ8KnAecCmwEbk5yUnkxgCTNhQP29KvqduCpkeZ3ApdX1e5u\nm+937ZuB66tqd1V9C9gBvGZ65UqSJjHumP7JwKYkdyS5Nckvdu3HAU8Mbbeza5MkzYHVDO/s73lH\nVdUZSV4N3ACcML2yJElrYdzQfwL4S4Cq+mqSZ5Icw6Bn/9Kh7TZ2bftI4ji/JI2hqjLuc1c7vJPu\nb48vAm8ESHIycERV/S1wI/BbSY5I8nLgFcBd+9tpVS3s32WXXTbzGqx/9nW0WP8i174e6p/UAXv6\nSa4FesAxSR4HLgM+BVyd5AHgR8DbuxB/KMnngIeAHwPvqmlUKUmaigOGflW9bT+r/vl+tv8w8OFJ\nipIkrQ2vyB1Tr9ebdQkTsf7ZWuT6F7l2WPz6J5VZjb4kceRHkg5SEuoQnMiVJK0Dhr4kNcTQl6SG\nGPqS1BBDX5IaYuhLUkMMfUlqiKEvSQ0x9CWpIYa+JDXE0Jekhhj6ktQQQ1+SGjLuzyWqBZ/85Oq2\ne8c71rYOSVNjT1+SGmLoS1JDDH1JasgBQz/JVUl2Jbl/mXV/kOTZJEcPtX08yY4k9yY5bdoFS5LG\nt5qe/tXAW0Ybk2wEfhl4bKjtbODEqjoJ+D3gE1OqU5I0BQcM/aq6HXhqmVVXAh8cadsMfLp73p3A\nkUk2TFqkJGk6xhrTT3Iu8ERVPTCy6jjgiaHHO7s2SdIcOOh5+kmeD1zKYGhHkrRAxrk460TgeOC+\nJAE2AncneQ2Dnv1Lhrbd2LUta2lpae9yr9ej1+uNUY4krV/9fp9+vz+1/aWqDrxRcjzwH6vq55ZZ\n9yhwelU9leQc4N1V9StJzgD+TVWdsZ991mpeWzPkFbkrWu7wNHoodAgloaoy7vNXM2XzWuBvgJOT\nPJ7kwpFNCghAVX0JeDTJN4B/B7xr3MIkSdN3wOGdqnrbAdafMPL4PZMWJUlaG16RK0kNMfQlqSHe\nWlnPWe2JW0kLy56+JDXE0Jekhji80wInlEvq2NOXpIYY+pLUEENfkhpi6EtSQwx9SWqIs3da1eiF\nWE5kUuvs6UtSQwx9SWqIoS9JDTH0JakhnsiVVqHR895ah+zpS1JD7Olrcs6DlBbGan4Y/aoku5Lc\nP9T2kSQPJ7k3yReS/PTQukuS7OjWn7lWhUuSDt5qhneuBt4y0nYT8A+q6jRgB3AJQJKfBc4DTgXO\nBv4sSaZXriRpEgcc3qmq25O8bKTt5qGHdwC/0S2fC1xfVbuBbyXZAbwGuHNK9UpzzZEuzbtpjOlf\nBFzXLR8HfGVo3c6uTYeK00wkrWCi0E/yL4EfV9V1B9x4GUtLS3uXe70evV5vknIkad3p9/v0+/2p\n7W/s0E9yAXAO8Mah5p3AS4Yeb+zaljUc+pKkfY12iLds2TLR/lYb+un+Bg+Ss4APApuq6kdD290I\nfDbJlQyGdV4B3DVRhdKYHOmS9nXA0E9yLdADjknyOHAZcClwBPBfu8k5d1TVu6rqoSSfAx4Cfgy8\nq6pqrYrXghlN4Tk+w+kHhtar1czeedsyzVevsP2HgQ9PUpQkaW14Ra7WBXvm0uoY+mqeHxhqiTdc\nk6SGGPqS1BBDX5IaYuhLUkMMfUlqiKEvSQ0x9CWpIYa+JDXE0Jekhhj6ktQQb8Ow6Ob1HgLzWpfU\nOHv6ktQQQ1+SGmLoS1JDHNPX3Fvu9MAc/+iWNNfs6UtSQwx9SWrIan4Y/SrgHwO7qupVXdtRwF8A\nLwO+BZxXVT/o1n0cOBv438AFVXXv2pTeIKdBLiSHpzRPVtPTvxp4y0jbxcDNVXUKcAtwCUCSs4ET\nq+ok4PeAT0yxVknShA7Y06+q25O8bKR5M/CGbvka4FYGHwSbgU93z7szyZFJNlTVrinW3AZ79ZLW\nwLhj+i/aE+RV9SSwoWs/DnhiaLudXZskaQ5Ma8pmjfOkpaWlvcu9Xo9erzelciRpfej3+/T7/ant\nb9zQ37Vn2CbJscB3u/adwEuGttvYtS1rOPQlSfsa7RBv2bJlov2tdngn3d8eNwIXdMsXAFuH2t8O\nkOQM4GnH8yVpfqxmyua1QA84JsnjwGXA5cANSS4CHgPOA6iqLyU5J8k3GEzZvHCtCpckHbzVzN55\n235WvXk/279noorUNCctSWvLK3IlqSHecE2zs2y3vo1LVUf/1b1CV4eKPX1JaoihL0kNcXhHc+WV\nt+075PPIJsc+pGmxpy9JDTH0Jakhhr4kNcTQl6SGGPqS1BBDX5IaYuhLUkMMfUlqiBdnzQNvLbmi\n5S7Y8mY10njs6UtSQ+zpa2Zuu23WFUjtsacvSQ0x9CWpIYa+JDVkotBP8vtJvp7k/iSfTXJEkuOT\n3JFke5LrknjeQJLmxNihn+TFwHuB06vqVQxOCr8VuAL4aFWdDDwN/O40CpUkTW7S4Z3nAT/V9eaf\nD3wH+EfAF7r11wC/NuFrSJKmZOyhl6r6TpKPAo8D/we4CbgbeLqqnu02+zbw4omrXE+8EEvSDI0d\n+kleCGwGXgb8ALgBOOtg9rG0tLR3udfr0ev1xi1nfhnykibQ7/fp9/tT298kJ1nfDHyzqv4XQJL/\nALwOeGGSw7re/kZg5/52MBz6kqR9jXaIt2zZMtH+JhnTfxw4I8nfTRLgTcCDwK3Ab3bbnA9snahC\nSdLUjB36VXUX8HngHuA+IMAngYuBDyTZDhwNXDWFOiVJUzDRHPqq2gKMftd4FHjtJPuVDsRTJdJ4\nvCJXkhri1bLSnFrttxl/WkAHw56+JDXE0Jekhhj6ktQQQ1+SGmLoS1JDDH1JaoihL0kNcZ6+Dpnb\nbpt1BZLs6UtSQ+zpSwtu9Mpdr9DVSuzpS1JD7OlLc8C7hupQsacvSQ2xpz9NdtckzTlDXwvplbf9\n5AfsI5s8eymthsM7ktQQQ1+SGjJR6Cc5MskNSR5O8mCS1yY5KslNSbYl+XKSI6dVrCRpMpP29D8G\nfKmqTgV+HngEuBi4uapOAW4BLpnwNSRJUzJ26Cf5aeCXqupqgKraXVU/ADYD13SbXQP86sRVSpKm\nYpLZOy8Hvp/kaga9/K8B/wLYUFW7AKrqySQvmrxMLRpvribNp0lC/3DgdODdVfW1JFcyGNqpke1G\nH++1tLS0d7nX69Hr9SYoR5LWn36/T7/fn9r+UrXfTF75ickG4CtVdUL3+PUMQv9EoFdVu5IcC9za\njfmPPr/Gfe255cVZex3qnr7z9FfmTdjWjyRUVcZ9/thj+t0QzhNJTu6a3gQ8CNwIXNC1nQ9sHfc1\nJEnTNekVue8DPpvk7wDfBC4Engd8LslFwGPAeRO+hiRpSiYK/aq6D3j1MqvePMl+JUlrwytyJakh\nhr4kNcTQl6SGGPqS1BBDX5Ia4o+oaN0a/aEV8CIuyZ6+JDXE0Jekhhj6ktQQx/S1Liw3fi9pX/b0\nJakhhr4kNcTQl6SGGPqS1BBDX5Ia4uwdHZTlfgZx06ZDX4ek8djTl6SG2NPXxA71j6BLGp89fUlq\nyMShn+SwJHcnubF7fHySO5JsT3JdEr9NSNKcmEZP//3AQ0OPrwA+WlUnA08DvzuF15AkTcFEoZ9k\nI3AO8OdDzW8EvtAtXwP82iSvIUmankl7+lcCHwQKIMkxwFNV9Wy3/tvAiyd8DUnSlIw93p7kV4Bd\nVXVvkt7wqtXuY2lpae9yr9ej1+vtd1sdes7KkWav3+/T7/entr9U1XhPTP4Y+B1gN/B84O8BXwTO\nBI6tqmeTnAFcVlVnL/P8Gve159Yn19ftfddj6Lf6c4nvaPNfe11KQlWtunM9auzhnaq6tKpeWlUn\nAL8N3FJVvwPcCvxmt9n5wNZxX0OSNF1rMZ3yYuD6JP8auAe4ag1eY/bWWa9eUhumEvpV9VfAX3XL\njwKvncZ+Ja2d5fotDgOtf16RK0kNMfQlqSGGviQ1xNCXpIYY+pLUEO+AKWkvZ/Ssf4b+ajgnXwvO\nt7D2cHhHkhpiT19NeeVt+3Z5W70fj9pkT1+SGmLoS1JDDH1Jaohj+qOc5iBpHTP0G7Dcj6Fs2nTo\n65A0ew7vSFJDDH1JaojDO41yyEdqkz19SWrI2KGfZGOSW5I8mOSBJO/r2o9KclOSbUm+nOTI6ZUr\nSZrEJMM7u4EPVNW9SV4A/PckNwEXAjdX1UeSfAi4hMGPpWvOLTfkI2l9GTv0q+pJ4Mlu+YdJHgY2\nApuBN3SbXQP0mZfQ976xWqXRe/R4fx6tF1MZ009yPHAacAewoap2wd4PhhdN4zUkSZObePZON7Tz\neeD9XY+/RjYZfSxpgYx+QfbL8WKbKPSTHM4g8D9TVVu75l1JNlTVriTHAt/d3/OXlpb2Lvd6PXq9\n3iTlSGNZ7nbL0rzo9/v0+/2p7W/Snv6ngIeq6mNDbTcCFwBXAOcDW5d5HvCToS9J2tdoh3jLli0T\n7W/s0E/yOuCfAQ8kuYfBMM6lDML+c0kuAh4DzpuoQknS1Ewye+e/Ac/bz+o3j7tfSfPNSXCLzdsw\neCtljcmfXtQi8jYMktQQQ1+SGmLoS1JDDH1JaoihL0kNcfbOOuOdMiWtxJ6+JDUkVbO5H1qSWtPX\nbnT+vT392XKe/nO8YGttJKGqMu7z7elLUkPWx5h+o716STpY9vQlqSGGviQ1ZH0M7zTME7eSDsb8\nh773cdUC8c6bmnfzH/qS1g1/b3f2HNOXpIbY019jo2PumzbNpg7Nzmp+eL3VIaDVzrb2G8H0GPpz\nwJOxWo+8fGY+rdnwTpKzkjySZHuSD63V60iSVm9NevpJDgP+FHgT8B3gq0m2VtUjKz5xgboG/W3b\n6J1yyqzLGNs9u7bxCxusf1YWuf5t2/qcckpv1mWMPbGv3+/T6/WmXs+iWKvhndcAO6rqMYAk1wOb\ngZVDf7Xm4MOhv337WKE/L0M59+7avrChA9Y/S9u3z0foj8vQXxvHAU8MPf42gw+Cn/TXf/3ccpa/\nadxyITkPJ0Mfe8yTtFpbq53zP7rdarbZn0N9Qnm1/47j9vPmoH84d2Z7Ivd733tu+TBnj0rSWluT\n++knOQNYqqqzuscXA1VVVwxtM5sb+UvSgpvkfvprFfrPA7YxOJH7P4G7gLdW1cNTfzFJ0qqtyfBO\nVT2T5D3ATQymhV5l4EvS7M3s5xIlSYfeWl6cdVWSXUnuX2bdHyR5NsnRQ20fT7Ijyb1JTlurulZr\nufqTXJbk20nu7v7OGlp3SVf/w0nOnE3Ve2tZ9tgneW9X3wNJLh9qn5vau3qWO/bXDx33R5PcPbRu\nEer/+SRfSXJPkruSvHpo3SK891+V5G+S3Jdka5IXDK2bm+OfZGOSW5I82L3P39e1H5XkpiTbknw5\nyZFDz5mb479C/f80ydeTPJPk9JHnHNzxr6o1+QNeD5wG3D/SvhH4L8CjwNFd29nAf+qWXwvcsVZ1\nTVI/cBnwgWW2PRW4h8Fw2fHAN+i+Rc1R7T0Gw22Hd49/Zh5rX+m9M7T+T4A/WqT6gS8DZ3bLZwO3\ndsvnLMh7/y7g9d3yBcC/6pZ/dp6OP3AscFq3/AIG5xZfCVwB/GHX/iHg8qH/FnNz/Feo/xTgJOAW\n4PSh7Q/6/b9mPf2quh14aplVVwIfHGnbDHy6e96dwJFJNqxVbauxQv3LnTXfDFxfVbur6lvADpa7\nLuEQ2U/t72TwRt/dbfP9rn2uaocVj/0e5wHXdsuLUv+zwJ7e5QuBnd3yuSzGe/+krh3gZuA3uuVz\nmaPjX1VPVtW93fIPgYcZdDQ3A9d0m13TPYY5y5791H9cVW2rqh3smz8H/f4/pJPjk5wLPFFVD4ys\nGr2Ya2fXNo/e3X0N/POhr4iLUP/JwKYkdyS5Nckvdu2LUPteSX4JeLKqvtk1LUr9vw/8SZLHgY8A\nl3Tti1L/g93/vzD40N3YLc9t/UmOZ/CN5Q5gQ1XtgkGwAnuCfRHqv3OFzQ66/kMW+kmeD1zKYIhk\nUf0ZcGJVnQY8CXx0xvUcjMOBo6rqDOAPgRtmXM+43gpcN+sixvBO4P1V9VIGHwCfmnE9B+siBh2e\nrwI/Bfy/Gdezou6cw+cZHPMfAqMzVuZ6Bssy9U/Noezpn8hgzOm+JI8y6CncneRFDD6dXjK07Uae\n+/o7N6rqe9UNpAH/nue+Ri1C/U8AfwlQVV8FnklyDIM6Xzq03TzWDuy9/uPXgb8Yal6EYw9wflV9\nEaCqPg/sOZG7EPVX1faqektVvRq4Hvgf3aq5qz/J4QwC8zNVtbVr3rVn2CbJscB3u/ZFqX9/Drr+\ntQ79dH9U1der6tiqOqGqXs7gfjy/UFXfBW4E3g57r+Z9es9XsRnbWz/sfbPs8evA17vlG4HfTnJE\nkpcDr2Bw4muWfqJ24IvAGwGSnAwcUVV/y6D235qz2mHf+gF+GXi4qr4z1DaPxx72rX9nkjcAJHkT\ng7FXWJz3/t/v/nkY8EfAJ7pV83j8PwU8VFUfG2q7kcEJaLp/bh1qn7fjv1z9w4bfVwd//NfwLPS1\nDG6r/CPgceDCkfXfpJu90z3+UwZnnu9j6Oz0rP6Wq5/BCZ/7gXsZhOiGoe0v6ep/mG6WxpzVfjjw\nGeAB4GvAG+ax9pXeO8DVwDuW2X7u6wf+YXfc7wG+wqDDs2f7RXjvv4/BTJJHgD+e1+MPvA54pvt/\n9B7gbuAs4GgGJ6C3MZjF9sJ5PP4r1P+rDL6t/18Gdzn4z+Mefy/OkqSGeGtLSWqIoS9JDTH0Jakh\nhr4kNcTQl6SGGPqS1BBDX5IaYuhLUkP+P/7AaWEYevYLAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "nM, binM, patchesM = plt.hist(men[\"Height\"], bins=range(140,max(men[\"Height\"]),1), \\\n", " color=\"blue\", alpha=0.4, linewidth=0)\n", "\n", "nW, binW, patchesW = plt.hist(women[\"Height\"], bins=range(140,max(men[\"Height\"]),1), \\\n", " color=\"red\", alpha=0.4, linewidth=0)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Figure 7.6 Right" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAVAAAAFRCAYAAAAitPV8AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XmYFNW9xvHvj11FQBQ3NhFFcd+CJm7jFjEx0bjEBYwG\nbzSJmGiSeyXJncCAxiXGuCaaREUFRa9bFKMRUdxxQ1xBdByGRcQF2UQRmN/9o6qdpqenp6enq5fq\n9/M8/dBdfbrq1DD9zqk6p06ZuyMiIq3XrtgVEBEpVwpQEZEcKUBFRHKkABURyZECVEQkRwpQEZEc\nKUClIMysi5k9aGZLzezOYtcnF2b2hJmNyNO6+ptZg5kV/TtoZn3NbLmZWRZlS6bepUA/hAiY2Vwz\nWxX+Un5gZjeb2YYpZb5lZlPDMp+Z2b/MbHBKmY3N7Eozqw/LvWtmV5hZzwzb/oWZvWFmK81snpnd\naWY7R7WvrXAC0AvYxN1PyscKzaxr+POoM7MV4c/9LjMbko/1F0BeBmGb2Wgzuy3N8gYz27bFSrjP\nd/dunv2g8GbLhXW5Ncv1lD0FaDQc+K67dwP2APYEfpt408y+CfwHuA/YChgAvA48a2bbhGU6Ao8D\ng4Fvh+v6JvApkDYgzOxq4FxgJLAJMAi4H/hua3fAzNq39jMt6A/MacWXNGNdzKwT8ASwM/AdoBvB\nz2oSMLRtVc2vCH6W6aT7ueoqmai5ux55fgB1wKFJry8FHkx6/RRwTZrP/RsYHz7/L2ARsEGW29wO\nWAvsnaHME8CIpNenA08nvW4Afg7MAWqBvwJ/SlnH/cB54fOtgLuBj8Ly5zaz3THAauArYDnwY8CA\n/wXmAh8C44FuYfn+YV1GAPXAtDTr/C9gIdClhZ/LjsCjBH94ZgEnJr13M3AtMDms1/PAgKT3jwg/\n8xlwDTAt5ec3Ang7XPfDQL/mfpZp6pXYx5+E+7EQ+HX43hbA5wSt9UT5vcKfc/s06xoN3JpmeQOw\nbfjcgFHAe8DHBH9oeqTUpV34ehvgSWBZ+LO7FrgtpeyPwv+bj4Dfhe8dGf4/rwZWAK8W+7sY9aPo\nFYjjg6QABfoQtC6vCF9vQBB0B6f53BnAwvD5HcDNrdjm2UBdC2XSBehTSa8bCFrG3YHOwIFAfdL7\nPYBV4RfcgJeB3wPtwy/de8ARzWx7vS95GD5zwi/khsA9ifeTvqTjw59X5zTruwO4qYX93RCYF37Z\nDdg9DI8dw/dvDl/vTXA0NgG4PXxvU4JQ/UG4f+cBaxI/P+CYsP6Dws/+Dni2uZ9lmrol9nEi0AXY\nJQyjxO/NZODspPJXAFdl87NNqUMiQH8JPEfwR68j8Lekfe0PrKMxQJ8j+KPfAdifIEhT/29uADoB\nuwFfAjtkqktcHzqEj879Zrac4Au8mKAVBtCT4Au3KM1nFgGbhc83baZMc1pbvjl/dPdl7r7a3Z8G\n3MwOCN87AXjO3RcTnEbYzN0vcvd17j4X+CdwcpbbOZXgj0q9u68iOMVxclLnhAOj3f0Ld1+d5vOb\nEbRcATCz3cNzycvMbFa4+GiCPyq3euA1gqA+MWk997n7K+6eCLM9wuXfAd509/vC/bsyeXsEf7Au\ndvc54WcvAfYws75JZb7+WWb4OYxx9y/d/U2CQD8lXH4rcFq4b+3C5U3OcyY5ycyWJD0+Y/1D+LOB\n37v7IndfA4wFTkjtDDKzfsA+BD/7te7+LPBAyrY8rPdX7v468BrBH6eKowCNzjEenLc8mOAwMhGM\nnxH8Bd8qzWe2Aj4Jn3/aTJnmtLZ8cxakvL6Txi/1qQQhA9AP6J3yhf0tsHmW29ma4BAwoZ6gxbNF\nhrokW29/3f01d98EOI6g9QxBa2m/lDqemrKN5FBcBXRNqt/8lG0mv+4PXJVYd1gfB3pnWX/C8sll\n6sPtAvwLGGxm/YFvA0vd/eUM67rT3XsmPTYhaHUn1/e+pPq+TdCi3iJlPVsBS9z9y6RlqT8HCBoF\nCck/t4qiAI2OAYStuFuAP4evVxGcazsxzWd+CDwWPn8MONLMNshye1OBPma2V4YynxMc1iZsmaZM\nasfDHQQtlX7AvgQtOAi+VO8nf2Hdvbu7fy/L+n5A8KVO6E/whU7+YmbqBJkKfLuFn898gvOnyXXs\n5u4js6jfIoI/EsmSW5fzCQ6xk9fd1d2nZ1n/dOvsR/BzIWy13kXQCh1O5tZnNuYBR6XUdyN3Tz1q\nWQT0NLMuzdSxJRXVcaUALYwrgSPMbNfw9SjgdDMbGQ7F2cTMLgT2Izi0guALMx+4x8x2sMCmZvZb\nM2vSy+zu7xF0+txhZgebWUcz62xmJ5nZ/4TFZgLHmdkGZrYdcGZLFXf3mQStq38Cj7j78vCtF4EV\nZvY/4RjP9ma2s5ntk+XP5A7gfDPbxsy6AhcBk8LDYVi/9ZTOrQRf9vvC7bYzs87AN5LKTAYGmdlw\nM+sQ/kz2MbMdsqjfQ8BOZnZsuG+/ZP0/ONcDvzOznQDMrLuZnZDFepMZUB3+f+xM0Lk2Ken92wjO\ni3+PtgfoDcAfwz+EmFkvM/t+Sl1w93kE57bHhD+vb4bbT613cxYD25i1PKY0DhSg0Vjvr7C7f0LQ\nCv1D+PpZgh7L4wlCoI7gHNL+7l4blvkKOByYDUwhOJE/neBc5wtpN+r+S4Ie0+sIThW8BxwLPBgW\n+QtBK+9DgvNtEzLVO8ntwGE0Hr4TBt3RBOcM6wg6QP5BMJwoGzcRhMJTBD34q4BfZFGXxPZXA4cQ\nHIo+RPDzmU3QIfTDsMxKgsPfkwladh8QnKvsnGaVqev/lOAo4VKC0yoDgWeS3r8/XNckM1tK0FGY\n/Ictm5aYE/R2v0fwf3yZu09N2sZzBKd7Zrh7usPobNafcBXBaYFHzWwZQUfRkGbKDgO+RbDfYwlC\nfXUzZVNf/x9BwH5qZplOOcSCuVdUi1ukrJjZVGCiu99UxDpMAma5e02x6lCq1AIVKVFm9g2CizAK\neulreJpj2/C00VDg+wTjfyVFh2JXQESaMrPxBGNNf+Hunxd481sC9xIMuVsA/DQcAiYpdAgvIpIj\nHcKLiOSobA7hzUxNZRGJhLvnNOyqrFqgrb1OdfTo0UW/VjaKRxz3S/tUHo847lNblFWAioiUEgWo\niEiOYh2gVVVVxa5CJOK4X9qn8hDHfWqLshnGZGZeLnUVkfJhZngldCKJiJQSBaiISI4iDVAzu9HM\nFpvZ6xnKXG3B3SZnmtkezZUTESk1UbdAbyaYti0tMzsKGOju2xPccuD6iOsjIpI3kV6J5O7PhLck\naM4xBBPj4u4vhJPSbuHBPXdEpELV1dVTXT2ehQsb6N27HePGncGAAf3TvnfWWYfz978/xsKFDXTr\ntgqztSxb1u3r54sXd+LDD9+jR4/NWbr0I3r02JwPP5zHRx+tpaFhTZvqWexLOXuz/v1WFobLFKAi\nJShTsOWrTF1dPUcccQ21tTXARsDnTJ8+milTzgVo8t6dd57D2rUXAIMJ7lpTTXDT183C5z8Pn48G\nLgRuBC4P/z0T2Cn3H0gBLpPqD7zezHsPAt9Kev0YsFczZV1EsvP++3N92LAxXlX1Bx82bIy///7c\nNpd9//25PnDgrx1WOrjDSh848Nfrlc9HmZNPHufwefhe4rHKhw69wYcOvcFhVZP3YGLS6y8cbs/w\nvPHfQVznYbbklm+5fjDrDWQO0OuBk5Jezwa2aKasjx49+uvHE0880ewvhEgc5TPocik7bNiYpHKJ\nx+f+3e9e588/7z5livtBB90RBlRymdW+005P+69+5X7OOe7bbfeyw1cpZdZ5ly7LvXNnd7M1Ke8F\nj65dF3rXrgvTvgfvpLyeleH5Ew7neE/O8fPZsOQDdBvgjWbe+w7wUPh8P2B6hvVk/u0SibG2B90q\nP+ywW/yBB9wnTXK/6Sb3a65x32OPRx1Wp5Rd41tv/Y4fcoj7vvu677KL+7bbunfqtDxteG200Ye+\n777uhx3mvummb6ctM2DAFP/Tn9yvvtp90KAH0pb51rf+5KtWuZ96arr6r/Rhw8Y0s28rHcY08zrd\n85U+iJE+n639DI4p3QAluBnZBwQ3pJpHcNfBs4GzkspcS3BTrdeaO3x3BajETGsOsd0ToZJ6WBu0\n7M480/3YY90POMB98GD3jh1XpA2oHj1q/eij3X/4Q/cf/zhoDfbt+3TasrvuOsGnTnWfPt399dfd\n33vP/bjj/tRssCU0F3CtKZPpj0W69zp0ON3hbW8MyfMd5qZ5/muHt30QP/b5bOFncGT4udwDVJdy\nihRYuk6S/v3/zJVX/oSGhq2YPx8WLAgeiefz5q3FvWmf77bbPsaoUYez6aaw2WbBo7r6cu6992fh\nuhM+Z9iwy5kwYfR6nx8+vIaJE3+TVdl09R44MOjcydQBlGuZ6urxfPBBA1tvnb4XPvFeohf+gw8a\n2HjjoOd9+fJuXz9fvLgjH35YS48evej58fvc/tFzXLLBtly9cksaGtYCT+A5XsoZaQs0nw/UApUY\nWLLE/dBDb01z2LzOe/Zc6Mcc4z5ypPsll7hPmOD+5JPutbWJjpXMLbuEqM6BJsoPGzbGDzmk+ZZz\nvsrk3ezZ7r17u99443qLUQtUpDgyDcf54guYORNefBFeein4d9EiaN9+LsuWbdNkXYccMprHH09/\n5+BsWm3p6pWuBdeWsmXrnXfgsMNg7FgYMWK9t9oymYgCVCRHTUNtFb16Tebww49i9uyNmT0bBg+G\nb3wDhgwJHoMHw+mnZ3/YnLq92AddFBLhOW4c/PjHTd5uS4AW/dA82wc6hJeIZduxs3at+5tvun/z\nm/d60+E4a3zvvf/tzz/vvmpV89tpzWGztEHisP2mm5otQhsO4Yt9JZJISWju6peHHz6X1av788or\nfP147TXYaitYtmw7oGPKmjrQrdt09tvvqGa3NWBAf6ZMOZfq6suTWpPpD8WlDVpoeeaDAlQEqK4e\nnxSeABtRW3spO+20jgEDYO+9g8cPfgB77gk9esDw4fcyceK2pB6Kb711y3P0DBjQP+PhurRRAcIT\ndA5UYqila7Hdob4e3nij8fHAA4tZtWqLJus68MA/8tRTv2t2O63p2JECaWV4tuUcqFqgEivpQu2J\nJ27hrLNO5oMPevLGG/Dmm9CtG+y6a/D4zndg2bL7ePjh00htTfbr1/xsPToUL0EFankmqAUqseEO\nxxxzLQ8++BOgc9I769huu1c577x92HVX2GUX6Nlz/c+qNRkDOYanhjFJxVq5Eh5/HB55JHgsXLiM\nr77q3qRcpjGWCRomVMba0PLUIbzESqZzmO7BIXgiMF98EfbdF4YOhQcegIsvvpLbb286xlIdOzFW\n4MP29eQ6/qnQDzQOtCKkGyO5zTa/92uv/chHjHDfemv3AQPcf/5z9wcfdF+xouXPa4xljGUxzrMl\n6FJOiYv0k1s4W2/9HqNGbc/QobDddmAZDrh0KF4h8tTy1DlQKXvuwSD1449/mnnzDmzyfjbnMKWC\n5PGwXedApSw1NMD06XDPPcGjUyfo2tWAL4ANkkpmdw5TKkQxz3mm0G+lRKKurp7hw2s45JDRDB9e\nQ11dPQBr18ITT8DIkdC3L5x1FnTtCg8+GHwvJk/uy8CB1QQ3B4PEcKJx484o0p5ISSmh8AQdwksE\n0o2p3Gqr2znooJN5/PGN6dsXjj8+eOywQ/rP6xymNBFReOocqJSU9B1B69hrr6ncffe3GTCgWDWT\nshVhy7MtAapDeMmbREfQtGnfZv3wBGhP9+7PKjyl9UrssD2ZOpGkzebPh4kT4bbbglnYe/b8ioUL\nVwEbJpVSR5DkoITDE9QClSyk6xBasQJuuQUOPxz22APq6uCGG6C2Fv71r20YOPAPqCNI2qTEwxN0\nDlRa0LRD6Eu6dn0Lsz2pqmrHj34ERx8NXbo0/Zw6giRnBQxPdSJJZNJ3CK3m+OOv5u67/7tY1ZI4\nK/SUdOpEkii4w4wZO9G0Q6gzS5asLEaVJO7K4LA9mQJU0nr1VTj4YPjwwwMIrgxKpg4hiUCZhSco\nQCXFxx/D2WfDUUfBaafBiy9+pSuDJHplGJ6gc6AVK3XOzT/84Qz+/e/+XHQRDB8Oo0cHN05LLqsO\nIYlEIjzHjoURIwq+eXUiSauk61nv2HEZ++3XjRtu2IDBg4tdQ6kYRQ5PUCeStFLTW/h2Yc2a7vTt\ne5nCUwqnBMKzrXQlUgWaN68dTXvWu7BoUUMxqiOVqEzPeaZSC7TCPPIIvPrqOcDalHfUsy4FEpPw\nBAVoxViwAE44Ac45B66+eg0DB45CPetScDEKT1AnUuytWQNXXw0XXxyE56hRsMEG6lmXIijR8FQv\nvKT17LPws5/BFlvAddfBoEHFrpFUrBINT9A9kSpe6pjO888fwXXX9eXRR+GKK+DEEzPfxVIkUiUc\nnm2lFmiZazqmczXt2q3h9NMbuPLKbnTrVuwaSkUrg/DUONAK1nRMZ2caGtrz1Vd/UXhKcZVBeLaV\nArTMLVzYQNMxnRvwwQca0ylFVAHhCQrQste+/UAg9dSGxnRKEVVIeIICtKzdey/MnDmMrbf+GxrT\nKSWhgsIT1IlUtm68EaqrYfJk2GQTjemUElCm4alxoBXmssvgr3+FRx/V2E4pEWUanqBxoBXDHS64\nAB56CJ55Bvr0KXaNRCjr8GwrBWgJSx4gv9VW7Vm37nzq6zfmqadg002LXTsRKjo8QQFaspoOkF/L\nBhvU8uKLG7Dppv2KXT2Rig9PKEAvvJkNNbPZZjbHzC5I835fM3vczGaY2UwzOyrqOpWDpgPkO/DF\nF3245JKbi1ktkYDCE4g4QM2sHXAtcCSwM3CKme2YUux/gTvdfS/gFOCvUdapXKQfIL+RBshL8Sk8\nvxZ1C3QI8K6717v7GmAScExKmQYgcdFhD2BhxHUqC1tu2RFNeiwlR+G5nqi/jb2B+UmvF4TLktUA\np5nZfGAycG7EdSp5a9fC0qXns9FGb6MB8lIyFJ5NlEIn0inAze7+FzPbD5hAcLjfxJgxY75+XlVV\nRVVVVSHqV1ANDfCTn0BDw0a8/HJ3Lrzw8qQB8udqgLwUR4zCc9q0aUybNi0v64p0IH0YiGPcfWj4\nehTg7n5pUpk3gSPdfWH4uhbY190/SVlX7AfSu8N558Err8B//gMbpZ4CFSmGGIVnOqU8nd1LwHZm\n1t/MOgEnAw+klKkHDgcws8FA59TwrBSjR8NTTwWXZyo8pSTEPDzbKtJDeHdfZ2YjgUcJwvpGd59l\nZjXAS+4+GfgN8A8zO5+gQ+n0KOtUqv78Z7jrriBAe/Qodm1EUHhmQdfCl4B//hMuvBCefhr69i12\nbURoDM+xY2HEiGLXJlK6Fr7MJF+iuWbNrrz33rE8/XQHhaeUhgoKz7ZSgBZY00s0nT59/kyHDicC\n6mGXItNhe6toVHaBNb1E01iw4GdUV48vYq1EUHjmQAFaYLpEU0qSwjMnCtAC22KLTsC6lKW6RFOK\nSOGZM31rC8gd1qz5JV27vo4u0ZSSoPBsE3UiFdDVV0NtbVemT9+Uiy/WJZpSZArPNtM40AJ5/HE4\n9VSYPh222abYtZGKp/D8WilfyinA3LlBeE6cqPCUEqDwzBsFaMRWrYJjj4VRo4LfWZGiUnjmlQ7h\nI+QetDw7doRbbgHL6SBBJE8UnmnpUs4SknyZ5rJlB7BmzUG8+GJnhacUl8IzEgrQPGp6mWYD/ftf\nxIcf/ki97FI8Cs/I6BxoHjW9TLMd9fW/0mWaUjwKz0gpQPNIl2lKSVF4Rk4BmkfB5Zi6k6aUAIVn\nQeibnUebbHIeXbrMQ5dpSlEpPAtGnUh5cuedMHlyd556ajlXXaXLNKVIFJ4FpXGgeTBjBhx5JEyZ\nAnvsUezaSMVSeOZEl3IW0eLF8IMfwPXXKzyliBSeRaFD+BwkBsvPn2+8/fYZnHpqN44/XrfSlCJR\neBaNDuFbqelg+bVsu+1veeyxkTrXKYWn8GwzHcIXUNPB8h14//0xGiwvhafwLDoFaCtpsLyUBIVn\nScgqQM2sk5ltF3VlykHv3u2AL1KWarC8FNCcOQrPEtHit97Mvgu8AUwJX+9hZvdFXbFSNW7cGWy8\n8evAV+ESDZaXApozBw49VOFZIlrsRDKzV4DDgCfcfc9w2RvuvmsB6pdcj5LoRPrsM9hmm3UcccQV\nLFmyMhwsf4Y6kCR6ifAcOxZGjCh2bWIj6vlA17j7Ult/QsviJ1mRjB8P3/teeyZM+O9iV0UqicKz\nJGUToLPM7IdAOzMbAPwCmB5ttUpTQwP87W9BiIoUjDqMSlY2PR8jgb2BBuBeYDXwyygrVaqmToUN\nN4RvfrPYNZGKofAsadmcAz3O3e9taVnUSuEc6LHHwne+A2edVdRqSKVQeBZEW86BZhOgM9x9r5Rl\nr7j73rlsMFfFDtB584Jr3efNg65di1YNqRQKz4KJpBPJzI4EhgK9zeyKpLe6ERzOV4TEde9PP13F\n5pt34eOPt6RrV/W4S4QUnmUjUyfSR8CbwJfAW0nLVwCjoqxUqWh63fsqjjjiD0yZojk+JSIKz7KS\nzSF8F3f/skD1yVSPgh/CDx9ew8SJv2H9Szc/Z9iwy5kwYXRB6yIVQOFZFFGPA+1tZhcBOwFdEgvd\nfVAuGywnuu5dCkbhWZayGcY0HrgZMOAo4C7gzgjrVDJ03bsUhMKzbGV1Kae77518+aaZvezu+xSk\nho31KPghfF1dPbvt9jErV+4GdCJx3bvOgUreKDyLLupD+NVm1g6oNbOfAguBjXPZWLnZYIP+tGvX\nhxNOuJxPP12lm8RJfik8y142LdB9gbeBTYCLgO7Ape7+bPTVW68eBW+B1tTAokXB/Y5E8krhWTIi\nHUjfzAZ7u/vCXDaYq0IH6OrVsM028NhjsPPOBdusVAKFZ0mJ7JYeZvYNMzvWzDYLX+9sZrcCL+Sy\nsXJy112wyy4KT8kzhWesNBugZnYxMBEYBjxiZmOAJ4DXgFgPYXKHq66CX1bklCkSGYVn7GTqRDoG\n2N3dvzCznsB8YFd3f78wVSue55+HpUuDiUNE8kLhGUuZDuG/dPcvANx9CTAnl/A0s6FmNtvM5pjZ\nBc2U+aGZvWVmb5jZhNZuI9+uugrOPRfaabin5IPCM7aa7UQys6XA44mXwCFJr3H341pceTD8aQ7B\nLUE+AF4CTnb32UlltiMYmH+Iuy83s83c/ZM064q0EykxaUhtbVdmzPg5L730Mbvt1i+y7UmFUHiW\nvKjGgR6f8vraHNY/BHjX3esBzGwSwamB2UllfgJc5+7LAdKFZ9SaThqyhuOOu1oD5qVtFJ6x12yA\nuvvUPKy/N8G504QFBKGabBCAmT1DcEqhxt3/k4dtZy1oeSbCE6AjtbU1VFdr0hDJkcKzImRzJVLU\nOgDbAQcB/YCnzGyXRIs02ZgxY75+XlVVRVVVVV4qoElDJK8UniVt2rRpTJs2LS/rijpAFxKEYkKf\ncFmyBcB0d28A5prZHGB74JXUlSUHaD4Fk4Z8Tuq0dZo0RFpN4VnyUhtfNTU1Oa8r64Qws845rP8l\nYDsz629mnYCTgQdSytxP0EFFOGB/e6CgQ6XGjTuDgQNHA+vCJcGkIePGnVHIaki5U3hWnBYD1MyG\nmNkbwLvh693N7JpsVu7u6wju6vkowaz2k9x9lpnVmNnRYZn/AJ+a2VvAVOA37v5ZbruTmwED+jN5\n8rl06LCOAw+8iGHDLlcHkrTOnDkKzwqUzWQi04GTgPvdfc9w2ZvuvksB6pdcj0iHMT33HIwcCTNm\nRLYJias5c+DQQxWeZSqya+ETZRLDkJKsS1uyjD37LBxwQLFrIWUnEZ5jxyo8K1A2ATrfzIYAbmbt\nzew8gsHxsfLMMwpQaaV33mkMzxEjil0bKYJsDuE3B64GDg8XPQaMLPSA9ygP4RsaYPPN4bXXoHfv\nSDYhcaMOo9iIekb6te5+ci4rLxfvvAPduik8JUsKTwllcwj/kpn928xON7NY3srj2Wdh//2LXQsp\nCwpPSdJigLr7QOBCYG/gDTO738xi1SLV+U/JisJTUmQ1kN7dn3P3XwB7AcsJJloue3V19QwfXsOd\nd37Kgw/+lbq61MEGIiGFp6SRTSdSV4IZlE4GBgP/Au5y94Le1iPfnUiNMzCNBTZEtyyWZik8Yy3S\nm8qZ2VzgQYLQfDqXjeRDvgN0+PAaJk78DanXvw8bphmYJInCM/ai7oXfNpzoI1Y0A5O0SOEpLWg2\nQM3sz+7+a+AeM2vS9MtmRvpSphmYJCOFp2Qh0y09hrj7i2Z2WLr38zThctaiOwd6GUFfms6BSkjh\nWVGiPgc60t2vbWlZ1KK4Eumtt+ax++5bccABf6RPn2BaO4VnhVN4VpyoA3SGu++VsuzVxMxMhRJF\ngL7ySnAJ82uv5XW1Uq4UnhUpkk4kMzuJYOjSADO7N+mtjYGluWys1Lz1Fuy8c7FrISVB4Sk5yNQL\n/yLwKcFtOK5LWr4CeDXKShWKAlQAhafkLNNdOeuAOoLZl2LprbfgzDOLXQspKoWntEGmQ/gn3f1g\nM/sMSD75aIC7e8/IaxcxtUArnMJT2ijTMKZ27t5gZu3TvR/e76hg8t2J9Pnn0KsXrFgB7dPuocSa\nwlNCkdzSI+nqo75A+zAwvwmcTdNLeMpKXV09J5zwd9q1W8Tpp9doEpFKo/CUPMlmGNNM4BsE93d/\nBJgMbO/uR0dfvfXqkZcWaOMA+ouAzmgAfYVReEqKqG8q1+Dua4DjgGvc/XygbOdur64eT21tDUF4\nAmxEbW0N1dXji1grKQiFp+RZNgG61sxOBE4jaH0CdIyuStHSJCIVSuEpEcgmQEcAhwCXufv7ZjYA\nuCPaakWncRKRZJpEJNYUnhKRFs+BAphZB2C78OV77r420lqlr0PezoEeeuhNzJ07hmBEls6BxprC\nU1oQ9bXwBwK3AQsJEmdL4DR3fzaXDeYqn8OYJk5czPnnf8kuu9zE1lu30yQicTVnTnDfdoWnZBB1\ngL4M/Mjd3w5fDwZuc/d9ctlgrvIZoH/5C9TWwrUFnU9KCkrhKVmKuhe+UyI8Adx9FtApl42Vitdf\nh913L3aqtmxUAAAUrUlEQVQtJDKJ8Bw7VuEpkcomQGeY2fVmdkD4+BtlPpnI66/DbrsVuxYSiXfe\naQzPESOKXRuJuWwO4bsAvwASd05/mmA86JcR1y21Hnk5hF+7Frp1g48/ho3K+noqaUIdRpKDyG4q\nZ2a7AgOB+9z9slw2UGrmzIE+fRSesaPwlCJo9hDezH4H3A8MA6aYWSyOh3T4HkMKTymSTLMxvQUM\ncffPzawX8G93/0ZBa7d+fdp8CF9XV8/3vjePzz9fw/77P63hS3Gg8JQ2iuoQfrW7fw7g7h+bWVlf\nqtM4icglQAfmzt2X6dM1gL6sKTylyDK1QJcCjydeElzOmXhd8PvCt7UFOnx4DRMn/obU+8APG3Y5\nEyaMbnP9pMAUnpInUbVAj095XdbDzjWJSIwoPKVEZLon0tRCViRqjZOIrN8C1SQiZUbhKSUkq8lE\nSkFbD+Hr6uqpqrqdefMuIBh8oElEyo7CUyIQ6bXwpSIfvfB33/0hP/3pWnbb7R+aRKTcKDwlIpEN\npE/ZSGd3X53LRkpFhw5b8q1vwQMP1BS7KtIaCk8pUS2eADSzIWb2BvBu+Hp3M7sm8ppF4OOPYfPN\ni10LaRWFp5SwbHpQrgaOBj4FcPfXCIY0lZ2PPlKAlhWFp5S4bAK0nbun3ve3oPeEzxcFaBlReEoZ\nyCZA55vZEMDNrL2ZnQfMibhekVCAlgmFp5SJbAL0Z8CvCO4LvxjYL1yWFTMbamazzWyOmV2Qodzx\nZtZgZntlu+7W+ugj6NUrqrVLXig8pYy0GKDu/pG7n+zum4WPk939k2xWHl4/fy1wJLAzcIqZ7Zim\nXFeCOUent6762amrq2f48BpeeGExV155PXV1qWckpCQoPKXMZDOh8j+AJoXc/awWV262HzDa3Y8K\nX48KPuqXppT7C/Ao8D/Ar919Rpp15TQOtHESkRqCq5BWMXDgHzSAvtQoPKVIor4n0mPA1PDxLLA5\nkO140N7A/KTXC8JlXzOzPYE+7v5wlutslerq8UnhCbAhtbU1VFePj2JzkguFp5SpFgfSu/udya/N\n7DbgmXxs3MwMuAI4PXlxPtadoElESpzCU8pY1lciJRkAbJFl2YUEnU8JfcJlCRsTnBudFobplsC/\nzOz76Q7jx4wZ8/XzqqoqqqqqWqyAJhEpYQpPKYJp06Yxbdq0vKwrm3Ogn9F4DrQdsAQY5e53tbhy\ns/bAO8BhwCLgReCU8NbI6co/AfzK3Zvc9bPt50AvBLqgSURKhMJTSkSUN5UzYHcaW40NrUkxd19n\nZiMJOojaATe6+ywzqwFecvfJqR8hz4fwAwb0Z8qUczn++Kf4+OONOfjgRxk3TuFZVApPiYlsWqBv\nuvsuBapPpnq0aTamX/wCtt0Wzjsvj5WS1pszJ7hvu8JTSkTUvfAzw57yslZfD/36tVxOIpQIz7Fj\nFZ4SC80ewptZB3dfC+wJvGRmtQS9MUYwljOyK4aiMG8e9NdRe/Ekh+eIWNwhWyTjOdAXgb2A7xeo\nLpGaN08t0KJJnPNUeErMZApQA3D32gLVJTIrV8IXX8BmmxW7JhVIHUYSY5kCtJeZ/aq5N939igjq\nE4lE69Py2r8vLVJ4SsxlCtD2QFfyPKyoGNSBVAQKT6kAmQJ0kbuPLVhNIqQOpAJTeEqFyDSMqexb\nnhBciXTNNU/z2GOPM3x4jaayi5rCUypIswPpzaynuy8pcH2alctA+sbLOC8haGzrMs5IKTylDEUy\nkL6UwjNXjVPZJc5UbKSp7KKi8JQKFOspiTSVXYEoPKVCxTpAG6eyS6ap7PJK4SkVLNZJMm7cGQwc\nOJrG2fiCc6Djxp1RtDrFisJTKlyLszGVilxnY5ozp54dd+zLwQfX0Lt3O8aNO0MdSPmg8JSYaEsn\nUuwD9NNPYfvtYUnZd4mVEIWnxEjU09mVteXLoVu3YtciRhSeIl+riADt3r3YtYgJhafIemIfoMuW\nqQWaFwpPkSZiH6A6hM8DhadIWgpQyUzhKdKsWAdoXV09V145mSeffEkTieRC4SmSUWyHMTVOJPJH\noBOaSKSVFJ5SITSMKY3GiUQ6hUs0kUjWFJ4iWYltgGoikRwpPEWyFtsA1UQiOVB4irRKbNOkcSKR\nteESTSSSkcJTpNVi24kEQUfS/vt/Ra9eL7DrrrWaSKQ5c+bAoYcqPKUiaTKRDA48EC66CA46KIJK\nxUEiPMeOhREjil0bkYJTL3wGK1fCxhsXuxYlSuEp0iYVEaBduxa7FiXonXcUniJtlOm+8LGgAE1D\nHUYieaEWaKVReIrkTWwDtK6unmHDali5soGzzx6r6+BB4SmSZ7HshW+8Dr6G4GokXQev8BRJT73w\nKRqvg09cylnh18ErPEUiEcsA1XXwSRSeIpGJZYDqOviQwlMkUrFMlMbr4L8Il1TgdfAKT5HIxbIT\nCYKOpBEjnuX11/fhqKPuqKzr4BWeIlnTtfDNmDQJ7rsP7rwzokqVIoWnSKuoF74ZK1ZU2HXwCk+R\nglKAxoXCU6TgYnstfF1dPePH17F8OXz88ZPxPgeq8BQpilieA228EulioCOxvhJJ4SnSJiV9DtTM\nhprZbDObY2YXpHn/fDN7y8xmmtkUM+vb1m02XonUMVwS0yuRFJ4iRRVpgJpZO+Ba4EhgZ+AUM9sx\npdgMYG933wO4B/hTW7dbEVciKTxFii7qFugQ4F13r3f3NcAk4JjkAu7+pLt/Gb6cDvRu60ZjfyWS\nwlOkJESdKL2B+UmvF5A5IM8EHm7rRhuvRFoTLonRlUgKT5GSUTK98GY2HNgbOLit6xowoD9TppzL\nEUe8SYcO89hnn5mMGxeDDiSFp0hJiTpAFwL9kl73CZetx8wOB34LHBQe6qc1ZsyYr59XVVVRVVWV\nr3qWPoWnSF5MmzaNadOm5WVdkQ5jMrP2wDvAYcAi4EXgFHeflVRmT+D/gCPdvTbDuip3GJPCUyQy\nJTuMyd3XASOBR4G3gEnuPsvMaszs6LDYZQRd5v9nZq+a2f1t3W6shjEpPEVKVuTnQN39EWCHlGWj\nk54fke9txmYYk8JTpKTFZFzP+mIxjEnhKVLyyihRstc4jGltuKTMhjElwnPsWIWnSAmL5bXwEHQk\nVVWtoGvXWey559vlM5lIcniOGFHs2ojEniZUbsYJJ8BJJ8GJJ0ZUqXzTYbtIwSlA06irq+egg1bR\nvftM9thjTum3QBWeIkWhAE3ROA70EoKBBiU+DlThKVI0JTsOtFgax4EmRmmV8DhQhadI2YplgJbN\nOFCFp0hZi2WAlsU4UIWnSNkroUTJn5IfB6rwFImFWHYiQdCRtP/+a+jV63l23fX90umFV3iKlBT1\nwjdj6FA477zg35Kg8BQpOeqFLwcKT5HYiW0LNHEIv9lmz7PbbkU+hFd4ipQsHcKnKKmB9ApPkZKm\nQ/gUJTOQXuEpEmuxDNCSGEiv8BSJvVgGaNEH0is8RSpCLAO0cSD9unBJAQfSKzxFKkYsO5Eg6Ega\nMqQd/fr9h8GDFxamF17hKVJ21InUDPd2QIHOeyo8RSpOLFugjcOYLgXaE/kwJoWnSNlSCzRF4zCm\n9uGSCIcxKTxFKlYsA7Rgw5gUniIVLZYBWpBhTApPkYoXywBtHMaUaHHmeRiTwlNEiGknEgQdSXvu\n2YVBg+5m0KBP8jeMSeEpEivqRMogr38fFJ4ikiSWLdDGYUyXEfyNyMMwpkR4jh0LI0bktg4RKTlq\ngaZoHMaU2L02DmNSeIpIGh1aLlJ+gmFMnwBfAhOApcAZuQ1j0mG7iDQjlgHavfty4CrgcuCXBEOa\nqunWrZWnKxSeIpJBLA/h3TsA40g+hIdx4fIsKTxFpAWxDNDlyzck3ZVIK1ZsmN0KFJ4ikoVYBmhw\nCF8NrAUuBOrJ+kokhaeIZCl250Dr6up59dUGYBTB7p0PVNO370rGjft95g8rPEWkFWLXAq2uHs+8\neX+k8RA+OP+5557dM48BVXiKSCvFLkCbm4kp4/lPhaeI5CB2ARrMxDQLuAhYA9QAs5o//6nwFJEc\nxS5AzzrrcDp0uBQ4D+gI/IYOHS7lrLMOb1pY4SkibRC7AP373x9j7doLgCsJWqCXs3btBfz974+t\nX1DhKSJtFLsAfe+9z4AbSW6Bwo3U1n7WWEjhKSJ5ELsAXbx4PnAmyS1QOJMPP5wfFFB4ikiexG4c\naI8emxNcxrkOcGAm8A49evRSeIpIXkXeAjWzoWY228zmmNkFad7vZGaTzOxdM3vezPq1ZXsffDCb\noOW5E2DALsAaNpj3qsJTRPIq0gA1s3bAtcCRwM7AKWa2Y0qxM4El7r49wXH3ZW3Z5pIlq4GvgOeB\nacCzDGIJdy2ZGZvwnDZtWrGrkHfap/IQx31qi6hboEOAd9293t3XAJOAY1LKHAPcEj6/GzisLRtc\nu3YlsABYDTzJIJYylaepplcswhPi+UusfSoPcdyntog6QHsD85NeLwiXpS3j7uuApWbWM/dNriA4\nfH+InixlKoupZh/G0zH3VYqIpFGKvfA53Zuk0VbA9fThU07nVqq5kPFMCZeLiORPpDeVM7P9gDHu\nPjR8PQpwd780qczDYZkXzKw9sMjdN0+zrvK4+52IlJ1cbyoX9TCml4DtzKw/sAg4GTglpcyDwOnA\nC8CJwOPpVpTrDoqIRCXSAHX3dWY2EniU4HTBje4+y8xqgJfcfTLBZUO3mdm7wKcEISsiUvLK5r7w\nIiKlphQ7kVqt0IP1CyGLfTrfzN4ys5lmNsXM+hajnq3R0j4llTvezBrMbK9C1i8X2eyTmf0w/L96\nw8wmFLqOucji96+vmT1uZjPC38GjilHPbJnZjWa22Mxez1Dm6jAjZprZHlmt2N3L+kHwR+A9oD/B\n7CEzgR1TyvwM+Gv4/CRgUrHrnYd9OhjoEj7/aRz2KSzXFXgSeA7Yq9j1zsP/03bAK0C38PVmxa53\nnvbrBuDs8PlgoK7Y9W5hnw4A9gBeb+b9o4CHwuf7AtOzWW8cWqAFH6xfAC3uk7s/6e5fhi+n03R8\nbanJ5v8JgokMLiG4EqLUZbNPPwGuc/flAO7+SYHrmIts9qsB6BY+7wEsLGD9Ws3dnwE+y1DkGODW\nsOwLQHcz26Kl9cYhQIswWD9y2exTsjOBhyOtUdu1uE9mtifQx91LfV8Ssvl/GgTsYGbPmNlzZnZk\nwWqXu2z2qwY4zczmA5OBcwtUt6ik7vNCsmiUxG42pizFZkiUmQ0H9iY4pC9bZmbAFQRD2r5eXKTq\n5FMHgsP4g4B+wFNmtkuiRVrGTgFudve/hOO9JxDMd1FR4tACXUjwi5nQh6aHEwuAvgDhYP1u7r6k\nMNXLSTb7hJkdDvwW+F54qFXKWtqnjQm+gNPMrA7YD/hXiXckZfu794C7N7j7XGAOsH1hqpezbPbr\nTOAuAHefDnQxs80KU71ILCTMiFDa71yqOATo14P1zawTwTjSB1LKJAbrQ4bB+iWkxX0KD3evB77v\n7p8WoY6tlXGf3H25u2/u7tu6+wCC87rfc/cZRapvNrL53bsfOAQgDJjtgfcLWsvWy2a/6oHDAcxs\nMNC5DM7vGs0f1TwA/Ai+voJyqbsvbnGNxe4dy1MP21DgHeBdYFS4rAY4OnzemeCv5bsEX8xtil3n\nPOzTFIKru2YArwL3F7vObd2nlLKPU+K98NnuE/Bn4C3gNeDEYtc5H/tF0PP+DEEP/QzgsGLXuYX9\nuR34gKBzch7wY+Bs4KykMtcSjD54LdvfPQ2kFxHJURwO4UVEikIBKiKSIwWoiEiOFKAiIjlSgIqI\n5EgBKiKSIwWoNMvM1oXTlb0a/tvsNIDhoOs38rDNJ8Jp1Gaa2dNm1uqrdszs7PASV8zsdDPbMum9\nv6e5tXZb6/mCme2WxWd+aWZd2rptKR0KUMnkc3ffy933DP+d10L5fA0qPsXd9yCYHefy1n7Y3W9w\n98S8m2eQNCmEu5/l7rPzUsvGev6N7Op5HrBhnrYtJUABKpk0uewtbGk+ZWYvh4/90pTZKWyVJSbb\nHRguH5a0/G/hBCKZtvsUkPjsYeHnXjOzf5pZx3D5JWb2Zridy8Jlo83s12Z2PLAPMCH8bJew5bhX\n2Eq9LKnOp5vZ1TnW83lg66R1/dXMXgwnUB4dLjs3LPOEmU0Nl307nKHpZTO708wUrmVGASqZbJB0\nCH9PuGwxcLi770NwjfQ1aT73U+BKd9+LIMAWhIfNJwHfCpc3AMNa2P73gTfMrDNwM8FlkLsTTPL7\ns3BKwmPdfZewJXhh0mfd3e8BXgZODVvQXya9fw/wg6TXJwGTcqznUIJr3hN+5+5DgN2BqnD2pWsI\nJqeocvfDzGxT4PcEl0DuQzDp8q9b2I6UmEqdzk6ysyoMkWSdgGvDWx6sI/3MQs8Dv7fgNiP3uvt7\nZnYYsBfwUtii60IQxulMNLMvgLkE80zuALzv7rXh+7cAPweuA74ws38CDxHMS5lOkxaku39iZrVm\nNoTg+ucd3P05MzunlfXsDGxEMNt5wslm9hOC79eWwE7Am6w/mcV+4fJnw+10JPi5SRlRgEprnQ98\n6O67hVMDfpFawN3vMLPpwNHAQ2Z2NkFw3OLuv89iG6e6+6uJF2FrLV0IrgsD8DCCWbZG0rq7DdxJ\n0NqcDdyX2Fxr6xmeCrgWON7MtiFoSe7t7svN7GaCEE5lwKPu3lLrVkqYDuElk3Tn/roTzAIFwfRf\n7Zt8yGyAu9eFh60PALsBU4ETzKxXWGaTDL36qdt9B+hvZtuGr08DngzPGfZw90eAX4XbSbWCxltP\npLqP4FYOJxPctoIc6/kHYF8zGxRuayWwwoJbQiTfbG15Ul2mA/snnR/eMJcRB1JcClDJJF2v+l+B\nM8zsVYLbVXyepswPw46dVwkmSb7V3WcB/ws8amavAY8SHN62uE13X00w/djd4WfXEcyF2g2YHC57\niqB1nGo8cH2iEyl5/e6+FJgF9HP3l8Nlra5neG71z8B/u/vrBFO8zSKYpf2ZpM/8A3jEzKZ6MHfm\nj4E7wu08R3CqQsqIprMTEcmRWqAiIjlSgIqI5EgBKiKSIwWoiEiOFKAiIjlSgIqI5EgBKiKSIwWo\niEiO/h/2QdAdH2qInQAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "tpr = []\n", "fpr = []\n", "#total num people\n", "totalW = len(women[\"Height\"])\n", "totalM = len(men[\"Height\"])\n", "\n", "#total num people checked.\n", "sumW = 0\n", "sumM = 0\n", "total = len(data[\"Height\"])\n", "\n", "for i in range(len(binM)-1):\n", " tprate = (totalM-sumM)/float(totalM)\n", " fprate = (totalW-sumW)/float(totalW)\n", " tpr.append(tprate)\n", " fpr.append(fprate)\n", " sumM += nM[i]\n", " sumW += nW[i]\n", " \n", "plt.figure(figsize=(5,5))\n", "plt.plot(fpr, tpr, 'bo', fpr, tpr, 'b-')\n", "plt.plot(np.linspace(0,1,100), np.linspace(0,1,100), 'r-')\n", "\n", "plt.xlim(-0.05, 1)\n", "plt.ylim(0,1.05)\n", "plt.xlabel(\"False Positive Rate\")\n", "plt.ylabel(\"True Positive Rate\")\n", "plt.title(\"ROC Curve for Gender by Height\")" ] } ], "metadata": { "kernelspec": { "display_name": "Python 2", "language": "python", "name": "python2" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 2 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython2", "version": "2.7.12" } }, "nbformat": 4, "nbformat_minor": 0 }