Skip to content

Instantly share code, notes, and snippets.

@adamgreenhall
Created September 10, 2015 21:53
Show Gist options
  • Save adamgreenhall/266e9027d50637af0fef to your computer and use it in GitHub Desktop.
Save adamgreenhall/266e9027d50637af0fef to your computer and use it in GitHub Desktop.
class-weighted logistic classificiation intercept correction
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 45,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"import matplotlib.pyplot as plt\n",
"from sklearn.linear_model import logistic\n",
"from sklearn.datasets import make_classification\n",
"from sklearn import calibration"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"X, y = make_classification(n_samples=100000, weights={0: 0.8, 1: 0.2})"
]
},
{
"cell_type": "code",
"execution_count": 49,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.text.Text at 0x10b30bed0>"
]
},
"execution_count": 49,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": [
"iVBORw0KGgoAAAANSUhEUgAAAY4AAAEXCAYAAAC6baP3AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\n",
"AAALEgAACxIB0t1+/AAAIABJREFUeJzsnXmcU9X1wL+HYRhgAFEEFbGOoNalbj+N+4IZpFVbbbW9\n",
"GtSCXYRYt2rjlrobrKYuWGtQW8U1elvR1qUVnThWbdW4tm61oqiIAoIIguz398d9MWGYJTPZk/v9\n",
"fPgw77373j058ybnnnPvPUeMMTgcDofDkS29Si2Aw+FwOCoLZzgcDofD0S2c4XA4HA5Ht3CGw+Fw\n",
"OBzdwhkOh8PhcHQLZzgcDofD0S16F6sjpdSmwOXAzlrrPdu53hf4LTAb2Ab4jdb6f8WSz+FwOBzZ\n",
"UUyPYz/gQUA6uH4GMEtr/RvgWuCP3Xm4Ump0TtJVEU4XaZwu0jhdpHG6SNMTXRTNcGit7we+7KTJ\n",
"YcC/vLavA7sopQZ0o4vRPZeu6hhdagHKiNGlFqCMGF1qAcqI0aUWoIwY3d0bymmOYxiwJON4sXfO\n",
"4XA4HGVEORmOecCgjONB3jmHw+FwlBFSzFxVXiwtqrX2eccbAqu11kuUUucAa7XWUaXUTsANWuuD\n",
"unjW6NTxxIkTLyqk7A6Hw1Gt3HTTTZdkHLZqrVs7a180w6GUOhD4MfBtIAZcA1wCLNRaX5mxquoT\n",
"YGsgorV+N9vnt7S0mObm5o4m3h0Oh8PhISJ1wHBjzEc9+e4sqsdRSJzhSCMio40xraWWoxxwukjj\n",
"dJGmFnUhIj8DZhtj/p55viffnUXbx+FwOByOkvJv8jRv7AxHFVJrI6nOcLpI43SRptS6aByx7bSG\n",
"IcOb2p5fsWDOrKWz35mQjz68cNR44A5jzGpjzAv5eC44w+FwOBxFp2HI8KZRx1+43uKfmXddmvOz\n",
"Nxs6ctrgQcOathu5FytXfbVVfe+Gn2w/au/VixbPm/XJ/Pcm5NwBznBUJbUYv+0Ip4s0ThdpqlUX\n",
"ItJ726Y9djhq7Om+jNPfAJg+Y0re+imnfRwOh8NR9fiiic169xu4SYEev9Oq1Su2LNCzv8Z5HFVI\n",
"NY6keorTRRqnizTF1IUvmugPHAgc4v0b0auh/5p8PV9EegO9jDErjTGvbD9q77cocNYNZzgcDocj\n",
"j/iiiV7ArlgjMRbYE3gFeBz4OfDSykVznwA63ODcTS4FZgE3AzT06bdhnp7bIc5wVCHVGr/tCU4X\n",
"aZwu0uRbF75oYgvSHsUYYAEwA7gOaE2G/Jl5+GicMmlWexPhKxbMmdWD7icDy8LBuAChAf0Hb9eD\n",
"Z3QLZzgcDoejm/iiiYFYj2Es1lhsDLRgvYpzkiH/h53dn8uSW2+Z7V+AnxtjPjHGfBkOxvtjS1Fs\n",
"/dEn//3L9BlT1gtVLVo8b1ZP+2yLMxxViBtVpnG6SON0kaa7uvBFE3XA7qQNxf8BL2ANxfHAK8mQ\n",
"f22exWwXY8waEbkImAsQDsabgAewG/wO/GDOm18VWgZnOBwOh6MdfNHEVqTnKfzAx1hDcQXwdDLk\n",
"X1osWUSkCTjSGDMFwBjzEkA4GD8YuAf4DXB9JBYoSg4pZziqEBfLTuN0kcbpIk17uvBFE4OBg0kb\n",
"iwFYQ/FX4LRkyD+n2HJm8AUZ9Yq8+YxTgfOB4yKxQEsxhXGGw+Fw1CS+aKIeu+IpFX7aCfgn1lgc\n",
"BfwnGfKXLAusiGwPrDTGzDTGfA7cChAOxvtiM4z/H7BPJBZ4v9iyOcNRhbhRZRqnizS1rgtfNCHY\n",
"kg1j97iq5RDgQeB97OqnC4FnkiH/8hKK2JYDgU+BmakT4WB8c2A6dvntvpFYoGjhskyc4XA4HFWL\n",
"L5rYCGgm7VX0xnoU9wEnJUP+sqoyKiLDjTFzAIwxN2VeCwfj+wEauB64qljzGe3hDEcV4mLZaZwu\n",
"0tSCLnzRRB9gX9LzFN8E/oE1FtcCbyVDfuPpotyMRi/gARE52hgzO/NaOBg/CbgMmBCJBf5WEgEz\n",
"cIbD4XBULF74aXvShuIA4G2sofgV8K9kyL+ydBJ2jYj0MsasNcasFZF9jTFfpyMJB+N9gCnYPSP7\n",
"R2KB/5VM0Ayc4ahCqn1U2R2cLtJUiy580cQw7O7s1E7t1dh5imnAj5Mh/4KunlEuuhCRbwFXY0tq\n",
"08ZobAL8GbsLfe9ILLC4JEK2gzMcDoejrPFFE32B/UnPU2wFtGKNxWTg3VKufsqRN4AT254MB+M+\n",
"4H7sbvDLIrFAUTYXZoszHFVILcSys8XpIk2l6MILP+1M2qPYF7sr+nHgFOCFZMi/Kpc+SqkLERkP\n",
"fGWM0cYYA6yzPyQcjP8Y+C1wUiQWeLAUMnaFMxwOh6Pk+KKJ4djw01jv/yVYjyIGqGTI/0UJxcs3\n",
"L5OxmS9FOBjvDUSB7wKjI7HAm8UWLFuc4ahCKmFUWSycLtKUky580UQj69ao2Jx0ksALkiF/QTe1\n",
"FVMXIiLAJGCaMeYrY8x/2rYJB+MbY5cIrwT2jMQCnxdLvp7gDIfD4Sg4Xo2K3UivfvJhR94zgJ8C\n",
"LyVD/rwVNyonjDFGRAYDg4D1EhCGg/FdsEkKNRCOxAJlrwdnOKqQSollFwOnizTF1oUvmvgG69ao\n",
"mIf1KK4Bnmpbo6KYFFoXXlW+bxljXgUwxlzRXrtwMH4McANwSiQWuK9Q8uQbZzgcDkde8GpUjCa9\n",
"+mkI8ATwGPCrZMg/u+O7q46RwLkiEvAmwNchHIzXAZcDxwJjIrHAa8UWMBec4ahC3Ag7jdNFmnzr\n",
"whdN9Ab2IO1V7AY8j/UqAsBrxapR0V0K8V54Xka9N4/xDtYorEc4GB+MTYXeF/BFYoHP8i1LoXGG\n",
"w+Fw0Dhi22kNQ4Y3tT2/YsGcWZnV6nzRxEjS8xQHA7Ox8xQRbI2KZUURuDw5B1iGTW3SLuFgfAds\n",
"csVHgV9FYoHVRZItrzjDUYW4uH4ap4s0nemiYcjwplHHX3hQ2/Pv3X15b180cRRpY9Ef61E8CJyS\n",
"DPk/KaDIBaNA78XV2FVRAGw2dOS0wYOGNaWO+zb0H9LYb/A3P//i0+ffePefZ+S576LiDIfD4eiQ\n",
"hqEj9gZOwnoVNwKvV/Au7bziLbN9CDjVGPO+MWadlOyDBw1rOmrs6esZ4+kzppT9qqmucIajCnEj\n",
"7DROF2k604X0ru/X3vnl8z76ZzLk/07BhCoR+XgvvGW25wAf5C5RZeEMh8NRo/iiiUHAj4ATGzbc\n",
"dLf2W5mynNwuFSIyAphgjLkcwBjzRnvtwsG49GtoHFZU4YqIMxxViIvrp3G6SCMio/e4quUf2BTd\n",
"JwJHAE8C0eXzPzoTu5O7JsjhvVgIfCQi0t4yW4BwMD4UiPXrO3DLXGQsZ5zhcDhqAF80sdX2p8Um\n",
"ALdh8yTdht1bMQ+gccqkH8y869L1vghXLJgzq5hyliNe7e86Y8zrxphlwO0dtQ0H40di82vd+fni\n",
"uS9h64NUHc5wVCFuhJ2mlnXh5YM6Gutd7NQ4Ytu4d/xK2wnuzCW3tUA334tdgVXA6x018PZmTMFm\n",
"8v1RJBZ49tah4WnTZ0xZL9S3aPG8Wd2TtvxwhsPhqCK8lOT7YY3FUcCzwO+Bh5Ih/4pSylZJiMiW\n",
"wIfGEu+sbTgYPwRbN+MhYNdILLAU4JP5700ouKAlwhmOKsTF9dPUii580cQWwI+BCdiKeLcBO2Tu\n",
"s6gVXWRDZ7rwltneis1o22Gp1nAwPgC4CpsG/aeRWODxAohaljjD4XBUKL5ooh/wfax3sQc2u+rx\n",
"2EJHbq9FNxGR3saY1d4y2zEdTX4DhIPx/bGlap8Bdo7EAouKJWc5UDTDoZQaA/wAmyHTaK0vbXN9\n",
"K2wRkxeAXYBpWuuaseD5xI0q01SbLrxQ1J5YY6GAF7HexZHJkH+9lN2ZVJsucqGtLkRkFHCniOzn\n",
"hac6WjHVF7gMOA4IRmKBvxRc2DKkKIZDKdUfu9JgB631KqXUn5VSfq11IqPZ2cDTWuspSqldsUnA\n",
"diiGfA5HueOLJjYDTsCGouqxo91dkiH/RyUUq2owxswUkSO78DJ2B+4A3sR6GRWXnDBfFMvj2Af4\n",
"QGudqhP8LHA4kGk4PgVSG2aGARWZA6cccLHsNJWsC1800QB8D+td7AtMx6b/eLYnoahK1kW+EZHR\n",
"wGbAAGPMLQDGmPnttQ0H4/VAGAgCZwD3RmKBmg4FFstwDGPdGruLSRuJFNcC05VSV2PjtacXSTaH\n",
"o2zwQlG7YY3FsdgloLdh624vLaVsVUgSu5CgQ8LB+I5YL2MusFskFphTDMHKnWIZjrnAwIzjDbxz\n",
"mUwDbtFa36eU2hh4VSm1rda63TTNSqnR2KIxAEycOHGdEZU3oqAWj40xreUkjzvu+njQNv935PAx\n",
"Pz5k4MidDwA2WDr7ndZ5zz5w6mcvPnZvqr2cnXt/KUr9eUt1DOyENQQAIzpq36e+rz/w3XN/NGLT\n",
"b/4QOO+Km45/15i126YMR7l8nnwdK6UuJk2r1rqVTpBOQnp5w5vjeA3YUWu9Uin1Z+za8leB1Vrr\n",
"JUqpF4CQ1voppVRvYD6wrda6XfexLS0tLaa5uVkK9RkcjnzjiybqgcOw8xYHA3/FehdPlWsBpEpH\n",
"RE4F7jfGdOg5hIPxrbED2dXAiZFY4P0iiVcSevLdWRSPQ2u9TCkVBK5XSs0HXtNaP6mUuhKb++VK\n",
"4JfAGUqpfYFRwEXZGg3HurhYdppy1IUvmtgJayyOA97FGovxyZB/cSH7LUddFBqvKt8expjnAIwx\n",
"v/POr6eLcDDeC7t341LsyqnfRWIBZ8DboWjLcbXWT2DrD2eeOyfj52exk+YOR9XhiyY2AsZhDcYm\n",
"2HDJAcmQv8MNZo68sBnwSxE5tosVU1tgN/0NAvaPxAJvF0vASsRtAKxCam1U2Rml1IVXk3ss1liM\n",
"Bf4GnA+0JEP+ohfzqZX3wvMy+hljlhhjPgKOadsmpYtwMC7YHfe/xS7QuapSy7kWE2c4HI4844sm\n",
"tsMaixOwNblvA05Khvw1tbu4hJwCDAAu76xROBjfFLgJaALGRGKB1wovWnXgDEcVUoux7I4oli58\n",
"0cQG2JHtidgvojuBQ5Ih/5uF7jtbqvm9aFMf40a6WGb7g0NOu2iHrfcJYpMTqkgs4BJAdgNnOByO\n",
"HuKLJnoBfqyxOBw7h3c58Fgy5HfhjuLyoIicb4x5wxizsqNG4WB8I+CGrbfcbX/g+5FY4LniiVg9\n",
"OMNRhVTrqLIn9FQXjSO2ndYwZHhT2/MrFsyZtcPpUy/DhqLGA59hQ1GnJ0P+sk5BUeXvxa+AmamD\n",
"zYaOnDZ40LCmzAYNffpv1Nhv0KjRex3zhz71fbeLxALt7hFzdI0zHA5HOzQMGd406vgLD2p7fvaj\n",
"t+yK3XtxN/C9ZMjv4uIlQEQ2xS7hP9fLSbjO6rTBg4Y1HTX29PV+fw+33vRaJBZwWSlyxBmOKqSa\n",
"Y9ndJd+6WL1s8Wzg/5Ihf4fhkHKlyt6LBdhNxd1i5crli6DqdFF0epVaAIejHOnVp98G7Z1fs3zp\n",
"Z5VoNKoBEdleRPYEMMasMsbc09HeDJFedcWVrrZwHkcV4kZSabqrC180cSBwcZ9BQ7YrjESlowre\n",
"i5HYDXovdNYoHIz7N9xgkz06a1MFuigpznA4HIAvmjgAuBi7lPay5Z/NrgMOLKVMjq8LLL1vjFlr\n",
"jHmks7bhYHwgXinXpcu++B82oaGjADjDUYW4+G2arnTR1mAAdydD/lWNUyaNnnnXpeuFQVYsmDOr\n",
"MJIWngp9L34LXAj8p7NG4WC8GbsnIwHsNPezWddNnzFlYdt2ixbPmwUVq4uywRkOR03iGYyLgK2w\n",
"ey/uSob8qUJjLJ39zoQSiVbziEifjL0YR3WRY2oQ1ss4HDgpEgv8DSASC0wouKA1jDMcVYgbSaVp\n",
"qwtfNLE/1sMYifUw1jEY1UwlvBciMhx4RER298JTnRmNQ4BbsBsvvxWJBb7Itp9K0EU54wyHoyZo\n",
"YzAuB+6sFYNRSRhj5oiI3xjTYTpzz8v4LfAd4OeRWOCxognoAJzhqEpc/DbNFt+ddMqmB6kjsTVe\n",
"atpglOt7ISJHA03GmKsBjDGfd9Q2HIyPxXoZM4CduuNltOmzLHVRKTjD4ahKfNHEfsDFw/Y5ckfg\n",
"19SwwagAngNe6qxBOBjfAOtljMV6GTOKIZijfZzhqEJqeSSVMhjA1sDlvfo03OEMhqWc3gsROQ24\n",
"zxgz1xjzcWdtw8H4t7Fext+wXkbOlRLLSReViDMcjqrAF03sizUY22BDUs5glDeLgfrOGnhextXA\n",
"GOCnkVjg8WII5ugaZziqkFqK33ZlMGpJF11RSl14Vfn2T/VvjJnWWftwMH4otsjSo8DO+fAy2sjj\n",
"3osccIbDUZF4BuMiYFsggjUYLodU+bIhMFFEnjbGdFg2NxyMDwauwdY5+UkkFniiWAI6sscZjiqk\n",
"mkdSvmhiH6yH8U3SHkaHBqOaddFdiq0Lz8sYZIxZaIyZDwQ6ax8Oxg/DehkPY+cylhRKNvde5IYz\n",
"HI6KoI3BiAC3Ow+j7BmPXQZ9fmeNwsH4hsC1wEHA+EgskCiCbI4ccIajCqmm+K1nMC4CtqMHBqOa\n",
"dJErxdBFm9rf04AOd34DhIPxw7Fexl+wXsaXhZQvhXsvcsMZDkdZ4osm9sZ6GD0yGI6S8ScRudIY\n",
"k+xiLmND4DrgAOCESCzwZNEkdOSMMxxVSCWPpDyDcRGwA9ZgTMvFYFSyLvJNkXRxFvBRZw3Cwfj3\n",
"gBjwAHbFVFG8jEzce5EbznA4yoJ2DMaRzsMof0RkY+BS4FRjzBpjzAcdtQ0H4xthvYz9gOMjsUBr\n",
"caR05BtnOKqQSorf+qKJvbAhqYIYjErSRaEpkC4WAi10PZdxBNbLuB/rZSzNsxzdwr0XueEMh6Pg\n",
"NI7YdlrDkOFNmed69ek7sK5v4ybfOOIXa4DJOA+jYhCR7YFNjDGtXhbb+1PXNhs6ctrgQcOaUse9\n",
"evXqPaD/htvU1fXud/Bex34vEgs8VQKRHXnGGY4qpNxGUg1DhjeNOv7Cg9qe/2D6de8AOydD/hWF\n",
"6rvcdFFK8qiLjYER7V0YPGhY01FjT1/vd/3A49c/XU5Gw70XudGr1AI4apfVyxZ/Ukij4cgfIvJN\n",
"EakHMMY8bYy5qzv3d1Zfw1F5OMNRhYjI6FLLkEld38Yhpeq73HRRSnLUxfnALp01CAfj0q+hcVgO\n",
"fRQN917khgtVOQqGL5poBK6uH7Dh1qWWxdF9RKSfMeYrAGPM+M7ahoPxEUCsX9+BWxZFOEdJcR5H\n",
"FVIO8VtfNLE7tjhP/+UL5rxYKjnKQRflQnd0ISJDgGQqPNUR4WC8VzgYPwl4BXjx88VzS/a77g7u\n",
"vcgN53E48oovmqgDQsCZwGnJkP/eximTps2869L1dhGvWDBnVrHlc2SHMWaBiOxjjOmwpkk4GN8a\n",
"W2CpP3BwJBZ4/dah4a2mz5iy3tLcRYvnzSqctI5i4wxHFVKqNeq+aOIbwB3e4R7JkP9DgKWz35lQ\n",
"bFlSuPX6abrShYh8D9jDGHMRgDGm3ey04WC8N3AGcC527831kVhgDcAn89+bkGexC4J7L3IjJ8Oh\n",
"lNpCa91pegFHbeCLJo4BfoetpRBNhvwd5ilylC3/Av7TWYNwML4T8EfgS2CvSCwwsxiCOcqLDg2H\n",
"UuobXdwrQBzYN5uOlFJjgB8A8wCjtb60nTZnAmuB4cAQrfVPs3m2Y12KOZLyRRODgBuAvYBDkyH/\n",
"S8XqOxvcqDJNe7oQkdOBvxhjZhljPgM+a+/ecDDegF1ZdTJwHvDHSCzQ6W7xcsa9F7nRmccxK1+d\n",
"KKX6Y9MN7KC1XqWU+rNSyq+1TmS0OQFYo7We4h3vlK/+HYXBF03sB9wFPAb8XzLkL2kaCUeP+AQ7\n",
"WOuQcDC+N9bLeBfYNRILfFwMwRzlS2eG4x9a69EASqmfYQvL36e1XqiUGgKMI/tVWfsAH2itUxNt\n",
"zwKHA5kFW44DnlZKnQZsAtya9adwrEOh47e+aKIeuAA4CZiYDPn/Uqi+csXFstN4exeeBsYYYx4D\n",
"MMbojtqHg/FGbJXFY4HTgT9VspeRiXsvcqOzL35/xs/HaK1jWuuFAFrrBVrr3wFHZNnPMCBzom2x\n",
"dy6TbwBDtdbXY43G35VSkuXzHUXCF01sDTwD7AnsWs5Gw9EujcDxItKns0bhYHwMdr5jY+BbkVhA\n",
"V4vRcOROhx6H1jrTfd1UKbWl1vrrlMlKqa2wnkE2zAUGZhxv4J3LZDHwnNf3TKVUI7AF8GGWfTg8\n",
"CjGS8kUTApwIXAlcBtyQDPnLPo2EG1WCiNQBQzJ0cUJHbcPB+GDgauAQYFIkFni08BIWH/de5Ea2\n",
"q6quAd5USv0TWAAMxYafTs7y/ueALZVSfbTWK7ET6r9XSm0IrNZaL8GmZh4JoJQaANQBn3b0QKXU\n",
"aGB06njixInruJ+plALuOPdjXzQxZOUX86fXNTSOqOvb/+BkyP+6iIyWs8tDPnfc5fGPgKNE5MbO\n",
"2h/R/Iv9v7XNfkHgwT/86bxJ8xZ8sCwSC1AG8rvjAh8rpS4mTavWupVOEGOy8z6VUjsCPwQ2A+YA\n",
"92ut38jqZr5eVfVDYD6wUmt9mVLqSmCh1vpKpdQgbMH697Bhqwe01n/P9vktLS2mubnZhbbIb/zW\n",
"F02MwdaO1sD5yZB/eT6eWyxqNZYtIr1SiQVFJBWSPrA9XYSD8U2wS6l3AX4eiQX+UTRBS0Stvhft\n",
"0ZPvzqwNRwql1BCt9YJu3VQEnOFIk48/Cl800YDd3HUscGIy5H88H7IVm1r9ghCRu4FbMj97W12E\n",
"g3HBhq2i2HnFSyOxwFdFFrUk1Op70R49+e7MKlSllBqIDVcdB3yqlPIBDwE/1lq/221JHQUlD0Zj\n",
"R+BurPe3azLkb3dtfyVQw18OIdqEetsYjS2BqcCmwKGRWODlokpXYmr4vcgL2S6nnQp8jq0V/LHn\n",
"cZyIdW8dVYIvmhBfNHEq8BR2U9/RlWw0agkR2UBEbkutljLGzGmvBoaXlPAX2ASUTwN71prRcORO\n",
"tpPjm2utjwNQSq0G0Fr/VynV6ZI+R2noiRvuiyY2xYYrNgb2SYb8/yuEbMWmhkISi7ElXFd31GCP\n",
"b4094dsHnHgSdsC4fyQWeLtYwpUbNfReFIRsPY4Gb9Pf1yilBmOzYjoqHF808T1sWuyXgf2qxWhU\n",
"O15VvsMAjOXhDryM+nAwfu6YfU/4HXAfcEAtGw1H7mTrccSAt5VSDwFbKaVuwO78/nXBJHP0mGxH\n",
"Ur5ooj92zf6hgEqG/E8XUq5SUOWjyv6sv5F2HcLB+G7YdCHz6up67xKJBT7orH2tUOXvRcHJynBo\n",
"re9QSs0EJgBvYdOPjNNa/6uAsjnyROOIbac1DBnelHmuV33DgPqBG20z4rCfPwTskgz5vyiNdI7u\n",
"ICLbAR8aY5YZY14BXtls6MhpgwcNa2rTrlef+r6DDx990qbYifI73M5vR77IdlXVblrrZ7E5plLn\n",
"fqKUWqK1fr1g0jl6RNv4bcOQ4U2jjr/woLbtZt1/7VvJkP/4ogpXZKowln0KcC827QsAgwcNazpq\n",
"7Onr/X5b/nX3fGDnSCzwKVSlLnqM00VuZDvHcU07597Arv92VChrvloyr9QyOLpGRAakfjbGnGKM\n",
"eaaz9im+WPLZmymj4XDkk049DqXUgdi6G4O9nzPpx7r5pxxlghtJpal0XYjIQOAFEfk/Y0y7m/NE\n",
"yGrzVqXrIp84XeRGV6GqOwCD3SR0e5trS4CbCiGUI7/0qm8Y0HUrRzlijFkiIru3ZzS8nd9HbDho\n",
"U18JRHPUMJ0aDq11E4BS6jat9YlFkciRM5nxW180MbrP4GE7l1ikklGJsWwRGQt8xxhzJoAxZlnb\n",
"NuFgfHfsiriNl371xbtAl4XPKlEXhcLpIjeyXY4b9ErJfqa1Xubt4ajXWs8voGyOHPFFE0cCf1j2\n",
"yczEzLsu7dv2+ooFc2YVXypHFvwLaHcvTTgY3wKYDDQDFwG3zf3sgz9MnzFlYdu2ixbPm1VIIR21\n",
"S7aGYzJwEPAT4DVgI+A+pdQVWuvphRLO0TO8VOjjsbUzDlv0xj+TpZapVFTKqFJETgGeNMa8YYxZ\n",
"wrqFzwgH4wOBc4FJwI3ANyOxwBKASCwwIZs+KkUXxcDpIjeyNRx7A3unSr9qrd9TSu2PraHhDEeZ\n",
"4Ysmfgn8Ejg4GfK/VWp5HFnxPrBezfZwMN4b+BnWu5gB7BKJBWYXWTaHYx2yNRyrMuqFA6C1XqGU\n",
"WlMAmRw9xKvSd9malV/9uK5Pv/2TIX/NV08s11i2V5Xve8aYBwGMMY9kXvcmvg/FLnmfCxyeazLC\n",
"ctVFKXC6yI1sDcc8pdQFwG3YQkzDsLvI25Z/dZQIXzRRh81WvNfM2y867Yt3Xqx5o1Hm1APfFZEZ\n",
"bSe/w8H4LsBvsaWTQ8DDbte3o5zI1nD8ArgLuCTj3ONAVe86rhR80UQf7NLpTYGDv3jnxcUlFqls\n",
"KKdRpedlbGqM+dgYsxwbgvqacDA+HLgcOAy4FLglEgusWv9JPaOcdFFqnC5yI9tcVfOAsUqpzYER\n",
"wEda6zkFlcyRFb5oohH4M7ASODQZ8tdEBbcK5TvAd4Fg5slwMN6I9SxOBW7BTny73GGOsiVbjwMA\n",
"rfXHwMepY6XUfVrrY/IulSMrfNHEhsDDwLvAT5Mh/2pw8dtMSq0Lz8tYa2yN5ke9fwCEg/E6YDxw\n",
"GbZ41u6RWGBWAWVx74WH00VudGg4lFK3AKd7+zbe76DZJoURy9EVvmhiM+Ax7Mq2s5Ih/3p1GBxl\n",
"wVTgr8BDnvEAIByMH4Kdx1gCHBWJBZ4vkXwOR7fpzON4Elju/fwFcDqslxPn2kII5egcXzQxEjvH\n",
"dBsQSYb860ycupFUmjLQxfnAgtRBOBjfEbtSalvgbOCBYk18l4Euygani9zo0HBore/JODxBa/2f\n",
"tm2UUuMLIpWjQ3zRxE7A34DJyZD/xlLL41gXEemPNeg/McYsNcbMBwgH45tgJ7x/gN1Q+/1ILLCy\n",
"dJI6HD0n28nx9YyGx3WAP3/iODrDF03sCzwAnJEM+eMdtXPx2zTF1oUxZpmI3Ax8BRAOxvtjN2P+\n",
"ErvybbtdldDpAAAgAElEQVRILLBeepBi4N6LNE4XudHZHIeLmZcRvmji29gl0T9Ohvx/K7U8jjQi\n",
"sg2wqzHmTwDGmJZwMN4rHIyfAESA54G9IrHAzFLK6XDki848jheAY7DzGkcAm2ML3S8EhgDjsJsB\n",
"HQXGF00cg93c9/1kyP9sV+3dSCpNkXRhgEGpg3AwPhqbuXYVEIjEAl3+zoqBey/SOF3kRmeG46jU\n",
"Xg2l1KFa60Mzrs0CXlJKzQB+U0D5ah5fNDEJuAAYkwz5/11qeRwWEdke+NQY87kx5l3g3XAw/k3g\n",
"KmBnbEJC7XZ8O6qRzibHMzf4baGUGqy1XpQ6oZTaCOuFOAqAl3fqPOzu4gOTIX/WYQ4Xv01TQF0c\n",
"DySAlnAwPhSbhPBYbEbiYyKxwPLObi4F7r1I43SRG9luALwHeFMp9QjwGTZX1WHA9YUSrJbxRRO9\n",
"sEs2xwL7J0N+t0u/DBCRwcaYRQDGmHA4GO8bDsbPxi6rjWMnvj8rqZAORxHIdlXVZKXUq8CPgN2A\n",
"OcCJWuu/F1K4WqFxxLbTGoYMb/IOpc8GQ7eVut79v3z/P48smfV6t42GG0mlyYcuNhs6ctrgQUNH\n",
"bjl8x923G7nXyyKypl9D47C+fQdscfBex7YA+0Vigf/mLm1hce9FGqeL3Mg65YjWep10CY780TBk\n",
"eNOo4y88qO35mXddOrwU8jjSiIhsN3KvpqPGnn6Ad2r/1LWHW296LRILfL9EojkcJSMrw6GUqse6\n",
"4ycAa4ADgZuASVpr55qXGS5+myYXXYjIgdhcUu2ycuXyRR1dK0fce5HG6SI3emXZ7mrsSpHzgMVa\n",
"6wXA74GbCyWYw1EGPL/5JttcNaD/BluWWhCHo5zI1nDsqrU+Rmv9AF7+Kq31k8AGBZOstmibAywn\n",
"3EgqTXd1ISJBEdk7HIzL+ZPuOWz8Dy75e11d/YACiVdU3HuRxukiN7I1HPVKqb6ZJ5RSDUBj/kWq\n",
"PXr3H7RpqWVwfM3bhx74s0Zsfe/LgZ99seSzN0osk8NRVmQ7OT4deFkpdS+wmVLqLEBhd5I7csAX\n",
"TTTW9R844v37rnx57aoVSzKvrVgwZ1ZPnunit2m60oVXL0MB954/6Z5B50+657vAj7FG48ZILLDq\n",
"1qHhE6bPmLLevYsWz5tVGKkLg3sv0jhd5Ea2y3GjSqk52DrjBru/4DqtdYeJ9hxZc/qI7/z04WTI\n",
"7wpilQYR6bXPaSf8fjB2h/6jwI6RWGBeqsEn89+bUCrhHI5yJNtVVRcBs7XWhxRYnprCF01sDJwJ\n",
"7J3P57qRVJr2dOF5Gd8wxrx//qR7dgP2BPYCjozEAskii1g03HuRxukiN7INVf0EyMloKKXGYGsR\n",
"zAOM1vrSDtodB9wJDNBaL8ulzwrgfODeZMj/bqkFqTH2HTxw6MnhYHwZNgPCecAdkVjAZYR2OLIg\n",
"28nxF4H/tT2plDojm5uVUv2BGHCG1voSYGel1Hp1PJRS2wPbZylTReOLJppI15vOKyIyOt/PrFRS\n",
"uhCROhGpCwfj9edPumf34LjrxgCLsGlCptWC0XDvRRqni9zI1uN4A/ibl6vqC++cACdhizl1xT7A\n",
"B1rrVd7xs8Dh2CRxwNfGJQRMxI7Eq53LgN8lQ/65pRakRrjmgD2OXoL1eueIyIGRWOCtUgvlcFQi\n",
"2RqOU4FXgaMyzgmwWZb3DwMyVwwt9s5lEgEu0VqvUkqlnl+V+KKJXYExwMmFeL6L36YxxrSGg/Et\n",
"z/7Z7U11db13xs4pPViL6c7de5HG6SI3sg1VxbXWBwPN2OJOzVrr0cAtWd4/FxiYcbyBdw4ApdQI\n",
"YDBwrFLqHO/0L5VSu2f5/ErjCuDyZMi/pMuWjh4hIn0a+vT7y1k/+cMVwMu9e9e/JCI7RGKBB2rR\n",
"aDgc+USM6fpvSCk1BJtC/YdAPbay2Z+A07TWXdZP9sJQrwE7aq1XKqX+jE1Z8iqwWmu9pE37tXQx\n",
"Oa6UGg2MTh1PnDjxojFjxhycGkmkYpjldrzHVS29gJtfvfToSau/XLS6EP1lxm9L/XlLcRwOxuWN\n",
"//3z16O+sevpK1d99fqgAUPGT546bqtyka+Ex7saY64rI3lKeXwG8GoZyVOy45aWFnPTTTddQppW\n",
"rXUrnZCt4XgY+Bhb83oBsDF2k9SmWuvvdvkAvl5V9UNsudmVWuvLlFJXAgu11ld6bTYGJgGXYOcA\n",
"bm5TUKpDWlpaTHNzc1mHt7w6Gy8A0WTIX7DNk7W6uUlERm0/au9jf3DIaaOxYdTTJk8dt7YWddEe\n",
"tfpetIfTRZqefHdmazhe0Frvme35UlAhhkNhswzvmQz5q34VTzEJB+MbLFu+5Kr63g3j6nv3CQOx\n",
"SCywqssbHY4apyffndnOcXyslFonL5VSagDwQcbxid3puNbwRRP1wGTgHGc08kef+r47nDzuujOA\n",
"t/v3HdirvnefUZFY4HpnNByOwpHtqqrPgFe85biLgI2waUdalFIXYldAjQduK4iU1cHPgXeTIX9L\n",
"oTuqFTc8HIzvOfHY3/65d+8+S4Ej2tv1XSu6yAanizROF7mRreH4LvB37Mqnwd6557HZcbfCGo5+\n",
"eZeuSvBFEwOxeZAOLbUs1cAh+52w/Z47H/Yr4NBBA4acB9xZCxv4HI5yIVvDEesoRUgKz/NwtM+Z\n",
"QEsy5H+1GJ1V60gqHIzXr1mz6tSdtj3wquUrlt7Ut6Fx+0gs8EVn91SrLnqC00Uap4vcyDY7bqdG\n",
"I9s2tYgvmtgEOA3Yo9SyVDLnTby7uVevXtfX1dV/vGbF0t2m3B78T6llcjhqlWw9Dkc3aByx7bSG\n",
"IcObAOoHbbw1mK9WLV5wW+OUSbOWzn5nQqH7r9T47WZDR04bPGhYU+a5urrefRv69NvygN2PHtTY\n",
"f4PjgL9MuT2Y9Qa+StVFIXC6SON0kRvOcBSAhiHDm0Ydf+FBbU5vPvMu55R1xuBBw5qOGnt6W73x\n",
"2DPTZt0/47oDP/rk7fUSbTocjuKT7XJcRwVRbSOppcu++KCnRqPadJELThdpnC5ywxkOR9nQu67e\n",
"rcxzOCoAF6oqBNKrrqTdV1j8NhyMN65duza8wcCN857UstJ0UUicLtI4XeSGMxwFoM+gIaNKLUMl\n",
"EA7GBVsf41oReWbego8+BTYvsVgOh6MLnOHIM75oYqz0rh808+7Ln8GsXZN5bcWCObOKIUMljKTC\n",
"wfi2xpjr16xZNap37z7jJ08d13rb0F9Pmz5jSlPbtosWz5vV034qQRfFwukijdNFbjjDkUd80cQg\n",
"4JYRh/7s6GTI/1ip5SlHwsF4I7bC48S5n826/fYHL1q4evWqVoBP5r83oZSyORyO7HCT4/nlKuDx\n",
"UhuNcqynHA7GJRyMH2WMeXPt2jUjgZ3/+Ofzz1q9etW4QvZbjrooFU4XaZwucsN5HHnCF000A4cB\n",
"O5ValnIjHIxviy0EtsWjT93yzGtvt7YaY7Kqs+JwOMoPZzjygJfE8A/AxGTI32nupGJQLvHbzLAU\n",
"NqX87157u7UB6LCyY74pF12UA04XaZwucsMZjvzwG6A1GfL/rdSClAOZq6WMWftsLH7Wy4sWz73H\n",
"GLPK1clwOCofN8eRI75o4mDgSGwG3LKglPFbLyz1N2zp3/GTpx43btHiuecB80ohj4tlp3G6SON0\n",
"kRvO48gBXzTRiA1RTUqG/J+XWp5SkhmW+vyLubE//Oncv69ctbwVwBjzUkmFczgcecUZjty4Ang2\n",
"GfI/XGpBMilm/DYzLAU8C+wci/9yuXeu5FRTLLulpeUnQBPQo6JVTzzxBC0tLaPzKVOlUkO6EKAO\n",
"uLm5ufnDfD3UGY4e4osmDgSOpoZXUWWulnr7vRcumD7jumeMMXMisQDAH0srXXXR0tKyH7C2ubnZ\n",
"FUxzdIuWlpYG4KqWlpar82U83BxHD/BFE/2BW4GTkyH/wlLL05ZCx2/DwXhjOBiPAP8EZgC7Tp9x\n",
"XX/K0IhWUSz7MOD2XB5QRbrImVrSRXNz8wrgbOCkfD3TeRw9IwI8nwz5/1JqQYpJ27DUv155aEzi\n",
"uXteBYjEAlNLKlz1s7q5uTnrAlYORybNzc0rWlpa1nTdMjuc4egmvmhiP+AYynB0naIQcf3MsBQw\n",
"fvLUcf8A/iUSP9oYMzvf/eWLKprjyNloVJEucqZGdZG3gYcLVXUDXzTRD7gNOCUZ8i8otTzFoE1Y\n",
"6vGZH762WyQWaDXGrAX2LWej4ahMRORAEXlZRJ4XkdtKLY9jfZzh6B6XAS8nQ/7ppRakM/IRv03l\n",
"lgLeBLYCdp48ddyM+x698q+pNsaYvLm+haKWYtldUUG6iABXGWP2Ap7v6UNEpFVExndwbXRPn+tw\n",
"oaqs8UUT+wDHATuXWpZC0zYsFYkFWgEmTx33CfCTEormqA02Bz4BMMbkMnfm5oQKhPM4siAjRHVa\n",
"MuSfX2p5uqKn8du2YSlg18lTx20pIj/ynmsqLTlhjcay26UQuhCRC0XkExG5RUQeEJEXRORhERni\n",
"Xf+2iPzTG/3/VUQ2a3Pf9SJyt4i8LSK3icj9wGbAdSLyiNd2mIhMF5GnRORZETkyo/8BXt9Pe//+\n",
"ICIbi8gVwK7AuSLypIgcVmhd1BLO4+iAxhHbTmsYMrwJoH7gRiOlrnfflYvm/aJxyqTDl85+Z0Jp\n",
"pcuNzYaOnDZ40LCmzHN9Gxo37t93QNPovY79K7BzJBaYAzB56riXgSUlENNRARhjLhWRrYCDgN2M\n",
"MUtF5CbgehH5NfAnYHdjzP9E5GTgDuCQNvftBfQBTjHGTBaR94HTjTH/8Lq5G3jWGHOxiGwKvC4i\n",
"rxpjPgCuAcQYc4CICPAgsIMx5jwR2Ru4zRhzR/E0Uhs4w9EBDUOGN406/sKD2pw+aOZdl5ZEnu7Q\n",
"VT3lwYOGNR019vS2n42HW296bfLUcccBEydPHTfNGLPcGPOfQspaaFxt6TQF1sUjxpil3s93Ai3A\n",
"G8CLxpj/eefjwA0isokxZq537gljzHJgOTaDcluZNweagZ8CGGM+FZFngICIXAWcABziXTMichbQ\n",
"ZYZq917khjMcjq9ZuXL5Iu/HIcAG2D9mRwUgIhcDGGMu7uS4CWjNpn0PWJTx80KgHhgB7CAiT2Zc\n",
"mwVsAqQMx+IunjvC+/8OEUnNWQwB/g1sDDQAX4ePjTHv9kR4R/dwhqMK6Wok1ae+YWAn9xrsqpaq\n",
"oFZGlW2/8HM97gEbZfy8MbAKmI31OL6buiAig1nXWHQ1gf2R9//RxpgF3jMasN9dXwErgGHAf71r\n",
"w4E1GR5Nu9TKe1Eo3OR4DREOxvuFg/HowMYhZbt50VGRCDBWRAZ4xz/Gzm3cA+wlIt8AO8mN9Xgk\n",
"4z6hfQTAW4zxuPfMFDcBo729RHcAE7zn98Jmq97Ea7cE6Cci3xCRcG4f0ZGJMxwd0Ku+obHUMvSU\n",
"9taoh4PxA7Hu/YjPF89NFl2oEuHW66cpoC4Mdk7jNhF5AdgUOM0YMwsYB9zjhaviwEnGmDUicibw\n",
"bWC8iPwmQ8bp2C/+67z5CoDjgX1F5B8i8g/gbWPMI961M4G1IvI08DTwkDHm3961W7FzI1OB5zIF\n",
"du9FbrhQVQf06tN38IcP/u79VV9+vk42yRUL5swqkUg9IhyMD8RWKPw+cHIkFvjLrUPD0+595Mrh\n",
"IrKmvnfD1zu/Fy2eN6tUcjoqnnnGmFPanjTGPI71GNqevwa7Iqrt+aPaOTcf+FF7nXoT8j/v4NqD\n",
"2FVWjjzjDEc7+KKJTbb8wekbA9skQ/7PSi1Pd0nFb8PB+LeBm4EngG9FYoHPAT6Z/94EEekHrPDc\n",
"/arFxbLTFFAXnYWcyhL3XuSGMxztczKgK9FoAISD8Q2xo7mDgZ9HYoEZIiKTp457GLtWfpYx5qvS\n",
"SumoBkTkQmzI6SsR+dAY43JL1QBFNRxKqTHYtNzzAKO1vrTN9XOw8c1PgD2AC7XW/y2mjL5ooi8w\n",
"CRhdzH7zRTgY//7qNatu6V1Xfx+wUyQWWAJfr3E/F8hbFbBKwK3XT1MIXRhjLgXKf3NTG9x7kRtF\n",
"mxxXSvUHYsAZWutLgJ2VUv42zRq11mdqraPA/UC0WPJlcBzwUjLkf6sEffeYcDA+LByM3wdc9dLr\n",
"j18aiQVOmTx13ODM1STGmNerPTTlcDgKTzE9jn2AD7TWq7zjZ4HDgUSqgdY6syxmHUVOdeGLJgT4\n",
"pfevIvCKKwWwoak7gAlP/PPOVBhqIfCxiIi3P6PmcKPKNE4XaZwucqOYhmMY6xqCxd659VBK9cGu\n",
"2z65CHJlMga7tPCJIvfbI8LB+ObYpYZNwPcisUBSRLafPHVcL2PMG96Kk2mllNHhcFQfxTQcc4HM\n",
"HcsbkE478DWe0bgROF9r/X5HD1NKjSZjHmLixInrxC1T67S7c7zLBX+6rH7gRtclQ37Tk/uLdRwO\n",
"xmXmh69d1bT5jj+vq+t9LXD05Knj9p08ddxoYDiwjYgMLRd5S3mcOlcu8vT0+J133mkaM2ZMTu83\n",
"sKsx5rpy+DxlcHwG8GoZyVPw4xtvvLGpubmZ9q4rpS4mTavWupVOkGJFMLw5jteAHbXWK5VSfwZ+\n",
"D7wKrNZaL/Ha/B6Iaq3fVEodrbW+P5vnt7S0mObm5h4vCfRFE9tjd7VumQz5yzZHUzgY3wq7xHYw\n",
"8JNILPAfb2fuR6lwlJv4S1Mtumhpabm4ubn54lyeUS26yAe1qIuO3qGefHcWbXJca70MCALXK6Uu\n",
"A17TWj8JnOudB7gLOxfye6XUk8DZxZIPOAOYWq5GIxyM9woH46cBSWAGsI9nNAQbjhqValtrfxCd\n",
"4XSRxunCIiIxbJGyQvcTEpFZIrJIbE2Q1L9XROSiTu47Qmx9kic7alNqirocV2v9BG3mD7TW52T8\n",
"vN6u0WLgiyY2BhSwXSn674pwMP5N4I/Y+Zd9I7HAOyLSOxILpJbZNtfq5LfD0RkiMgEYb4w5OOP0\n",
"WdjkiAXFGBMVkf5Ac2b/InIQnSz3N8b8VUQ2xMvBVY64XFWWicADyZC/04yaxSYcjPcOB+PnAs8A\n",
"9wIHeUZja6DV8zZoazRcHp40ThdpnC4sxphlwAFF6q69XfX/xkZXurqvbKn5neO+aKIP8Avs7tey\n",
"IRyM74xN0rYQ8EVigVmpa8aYd0XkKOdl1DaZVSozWbFgzqxsq1Tm4xkiMgq4AeiHreQ3w6TrfHwH\n",
"uABYjU2Dfor3/p4EnIdNPrgE8GFreryJTYz4O2BH7EbgW40xl4hICDgKm7L9VeAsY8wqr5/jgFOA\n",
"ZZ5YEa/Pc4BNvbDPf4AXgMuBt0jXJvFh94wJ1qv/lTHmRRE5ArgSu4jnOWy1wjXAD4zNn9VtPA/I\n",
"GGNu93RwDOnU8qcaY9rdP+aFtr6NrZEzDzjD2KJWA4DrgW2wjsAdxpibeiJbd6h5w4H9xb2ZDPnL\n",
"otJdOBhvAMLY3evnANMisYARkXFAP2PMHwGMMfM6eoaLZaepZl10UKWSjqpUtqeL7j6jLSJSBzwE\n",
"XGGMuVNEBmEr/10sIiOx6dV384zFccBDIrKjMeZmsfXHJwHfwg6QJhtb8nUHYE/gUOwX4gHevScC\n",
"/2eMWS4i92HnQCMisi/wW+BbxpgFInI0cIIx5kSxmXcntAkVjcIuYUdENgD+BhxljPmHiOwP/E1E\n",
"ts4IGf0eG+46V2wd9J9iE4ciIn8BWo0x13aipp0y5is2Ba7IuHaoMWal5w3eBBzYjo53AJQxZkfv\n",
"+BpgW+BT4Fr4unTuAOA1EXndGPNsJ/LkTE0bjowNfxeUWhaAcDC+J9bLeBfYNVX32+MF7AjK4Sgn\n",
"9sYuzLgbwBizWESUdy0APG/SVfnuBW4B9sWGXwX4p/EKNGE9kBSPeB71O8A7IvI4cK+xZWZTz7oC\n",
"61mc6LVPPedB0ql12gv5ZJ77LvCF8eqbG2OeEZHPgSOwJXAF+K+x9c3Bhpm2yrj/XSDz77Q9/p0y\n",
"XCIyPuP8W8DDItIXWzFx5w7uX4L1mo7CGulzsKnke2FTzn/bk/1LEXkYW07XGY4CchDQHzviKBnh\n",
"YLw/cAn2F346oCdPHcfkqeNOwbqei003SmLW4lLDjqhFXfQdtuVBvmgiqzBm32Fb5trdCOBzk5HK\n",
"xhjzr4xrn2WcX+N9KY/IuL+j0rFt64aPAMaJSMpz6IsNG6WuvZrZD3b1YWdsmnFv27DTfGDzDmRc\n",
"gQ3Hpfo6i6752lAZY26Hrz2dh4ETjTHTRWRLoN19a8aYj0TkcKzB+D3WSF8ADMKWzr1KRFLZIgYD\n",
"r2QhU07UnOHIjOn22XCTb61dsXzh6mVfJBqnTMo6pptPvAJLf8S+6DtFYoH5AJFYIFXRbABd12V2\n",
"OL5m+bwPnkqG/KPbnm/PiG5016Wt2AFUT/kI2FBEeqWMh4hsj/0S/BAbUkn1XwdsiC0pC12Xjc3k\n",
"Q+zcydUZzxuScW1Yxvk6YEeTLujUFpPR90fA0DbXh2XI2N69+eCb2A3Rf/eO+7S5/nU/YksgvGmM\n",
"+YGIbAJMxxqRS7CG7BfGmJe8tr2xg+GCUnOGo52Y7hBgm2xjuj1ls6Ejpw0eNKwpdSzSq25A/8Ej\n",
"63v3GTx6r2MCkVjgryLSe/LUcXsaY14AMMZc35O+am2E3RlOF2kKpIvnsOGa44A7RWQj4D5gN2w4\n",
"6VwRGWWMmYmdT5wF/NO7t7OysW2vTQPOEJEbjDErPM/jF8APvWvTRWSIF646FtgdWx1wCdY7QURu\n",
"xE6gC+msFQ9jqw0eYIx5WkT2w2a1+GsnsqUP7HzDC8aYezto3xEfYEPPe2Pz9X2nk372Ao4EfmmM\n",
"mSsi/wV6eUvx78CmZ3rJa3sB1mO6oZvydIuaMxylYvCgYU1HjT19vZHdg0/87tlILJB6SYcDZ4nI\n",
"sW7FlKMrViyYM6u9AU93qlTm+gxjzFoR+R5wg4j8FLuy5xdeuOh9EfkhcLuIrMGuePqed884YDzQ\n",
"ICLTjDETAETkKmAX4BwR2dB4KVKMMXFvMv0pEVmKDWWd5F37l4j8CjvxvgL7xfkTT8QW4GwRuRdr\n",
"LDL7Pc8Yc4W38utqz8Nfi52wXiwifuzIfhNvVdPLGfee4cm2NR2UKvBWgY0HBotIAjvB/Zkn81wR\n",
"ORX4o4i8DvzPu+fv2HBUqt8pwGRgoveMPthJ8d963ZyJNXzPYlebvQxcnM3vLhec4Sgxa9euWSMi\n",
"A40xS4wxH2JHZTlRi3H9jqhmXXQ3tNqeLvIRnjXGvAcc1sG1jkrH3gPc0875s+kgY4TpoNysd+1u\n",
"vAn6NucXYUfsmdyVqQtjzMvYomdt700A27c5/VCbNke0J493LUonpSG8ZbOZS2fP7Kgf7EKD9p7R\n",
"YencQuI2ABaJ3nX1/do7v2r1ys2BU4ssjsPhcPQYZziKQDgY/+EGA4fu1t61+t59Pmbddd05U60j\n",
"7J7gdJHG6SKN00Vu1FyoauWiebM/fmzamhULZj9n1q79el9Ed+LC2RIOxuuBq4AjF3+54N/YCbs2\n",
"iHHzGQ6Ho5KoOcOx/Sk3PAIMSYb8hxayn3AwPgK7umQhsPv8hR9dO33GlC8B1q5d069Xr7qvABYt\n",
"njcr331Xc1y/uzhdpHG6SON0kRs1ZziwWXDvK2QH4WB8DHbX6fXAlZOnjtsE+PTTz94/0XkXDoej\n",
"0qkpw+GLJgYBfmyKgrwTDsZ7AedjS96Oi8QCTwJMnjpuATbBWlFwI6k0ThdpnC7SOF3kRq1Njh8J\n",
"PJUM+Rfl+8HhYHwIdjPRt4E9Jk8d96mI7AFgjFlpjLnbeRsOh6MaqDXDcQyg8/3QcDDuw+7cfAPw\n",
"e8kJR2E3BxUdV3chjdNFGqeLNE4XuVEzoSpfNLEhtnjLuHw9MxyMCzYt9CXAxMlTx/0bWONV5ns4\n",
"X/04HI78IbZ07EK8ehwF6qMJm44lVV9jH2y6lU+wKVDeNsYUJGReDGrJ4/g+0JIM+fOSMDAcjDdi\n",
"J8AnAftFYoEHsLtad8jH83PBxW/TOF2kqUVdiMgEWb9291nAhQXu2gCPGmMO9lKqfwJc5f18bIH7\n",
"Ljg143Fgw1S35eNB4WB8O+B+4PlX33rygEdab04Zo++7eQyHo7wxtnRsoZkLXNfJtc4KP5U9NWE4\n",
"fNHExlhX8Ye5PiscjB+DzTx57uSp4x4DWkVu2cMYs7ZcjIZbo56mmnXRNuNyikWL5836ZP57E9qe\n",
"b08X3X1Ge0iFlo41xhzq3Zv30rFewal207qnrnnZegNtPutz2FT0g40xW4nIPtgB7ycZxaBKUi42\n",
"k5owHMAPgMeSIf+XPX1AOBjvg81IeRgwNhILvOLVzGjOLGLjcBSLjjIuT58xpWjPkMotHbu393Mx\n",
"Sse2izHmZK92SeZn3R+4EZsqPpX59wpgQsatJSkXm0mtzHEcQw6b/sLB+BbAP4Bv3Hxf6LLJU8eN\n",
"Tl0zxnyeu3j5pVpH2D3B6SJNgXSxXulY7CZbaL907JbY0rGQUTrWWNYrHWuMeccY80fsF2fb0rEn\n",
"eD+3Vzr2xow+2iLY1OTQTulYIFU6NtU219KxXZH5WW9tR+avjzPKxd7myfsldhvACRSRqvc4fNHE\n",
"JlgX8NGe3B8OxscCd2AnvqOfff7xCKAufxI6HPll4w03PygcjGcVNt14w827btQ5lVo6NvO5hS4d\n",
"2xVtP2tnDKVE5WIzqXrDgY2JPpIM+b/qsmUG4WC8Dvg1cNKMZ26/78XXH5tmjDGRWOCjgkiZR6o5\n",
"rt9dalEXn33+8VORWGB02/Pt6WL7UVNaqc3SsZtkyF+K0rGdsRJrHFIMzvh5PiUqF5tJLYSqur3p\n",
"LxyMb4z1UPzAHi++/thrrF8T2OFwrFs6lozSsauw4aQ9vDkFyL10rBKRBq+fg0kXQZoGHJZhSI4l\n",
"PSewTulYL9ST+fyHgYEicoDXptulY0Ukm+W17X2mjs5/AGwsIkM9eb+duuAZ51S52BQXtDkuOFXt\n",
"cfiiieHAzqQLwndJOBjfyxij5y+c/cywIVuMj8QCqyOxwK2FkzL/1NoIuzOqWReLFs+b1d4kdkcZ\n",
"l8Iwk/MAAAqSSURBVNvTRXef0c4zXenYDkrHphCRGVgv5lcistR49ck7+awfi8jVwDPYOZUkcKqI\n",
"TDHGnE6JysWu85nKZAVpzrS0tJjm5uZ1LLcvmjgN2D0Z8o/v6n5vF/gpwAWfff7xr26+L3QYcJz3\n",
"B+BwlJSWlpaLm5ubLy61HI7KpaN3qL3vzq6o9lBVVinUw8H4wLVmbXzt2jU/B/a56d5f3WGMObZS\n",
"jYbLw5PG6SKN00Uap4vcqNpQlS+a2AJbaP6JztqFg/EdgT/PnT9rfvzhyY9+tWLpzKII6HA4HBVK\n",
"VRmOxhHbTmsYMrwJoHfj4BG9etevXPnF/BmNUybNWjr7nQlt258/6e5xIr2mAKHbpv/6Tmx8s+Kp\n",
"5rh+d3G6SON0kcbpIjeqynA0DBneNOr4C9suLdx05l2XrnMiHIw3ANd8uWzR+P/89+mfP/n8vfFI\n",
"LFA0OR0Oh6OSqSrD0R4LH72VjZat2HX7UXu3AtTV9W4YNGDIjqtWrfh86Vdf7P3WzOfesKsGq4da\n",
"3LvQEVWki25NXrb7gOrRRc7UqC5yfodSVL3h6LdqDUeOOW0D2mxymj5jystvzXzu9RKJ5XB0l94t\n",
"LS3S3NxcHcsgHUWlpaWlgTxmvKj2VVU1SQ2OpDqkinTxKHYPQY+pIl3kTC3pwjMaVwE35+uZVe9x\n",
"OBzVQHNz87MtLS3fbGlpuZQqWcThKAqC9TSubm5u7nSjYncomuFQSo3BpjefBxit9aVtrvfFpkae\n",
"jU0v/But9f+608eKBXNmtZ0I32jZ8t3ZiQG5yF5p1Gj8tl2qSRfNzc05ZTCoJl3kitNFbhQlVKWU\n",
"6g/EgDO01pcAOyul/G2anQHM0lr/Bptv/o/d7Wfp7HcmLHytdfTn/37qpM///dQhC19rHV2/lpdy\n",
"/gAOh8Ph+JpieRz7AB9orVd5x88ChwOJjDaHYSuCobV+XSm1i1JqgNa6J8WXLgCmAC/mmounEnEj\n",
"qTROF2mcLtI4XeRGsQzHMGyWyhSLyUiD3EWbbhsOY8zXRU2yLX/pcDgcjuwo1qqqucDAjOMNvHOZ\n",
"zAMGZRwP8s5ljYjU90i6KsPl4UnjdJHG6SKN00VuFCU7rjfH8Rqwo9Z6pVLqz9g6vq8Cq7XWS5RS\n",
"5wBrtdZRpdROwA1a6w4LzCilRgOjU8cTJ068qJCfweFwOKqVm2666ZKMw1atdWtn7YuWVt1bVfVD\n",
"bK78lVrry5RSVwILtdZXZqyq+gSb4z6itX634yeu9/yLtdYXF0D0isPpIo3TRRqnizROF2l6ooui\n",
"LcfVWj9Bm0y1WutzMn5ejq2H4XA4HI4yxu0cdzgcDke3qCbD0VpqAcqI1lILUEa0llqAMqK11AKU\n",
"Ea2lFqCMaO3uDVVTOtbhcDgcxaGaPA6Hw+FwFAFnOBwOh8PRLZzhcDgcDke3cIbD4XA4HN2i4upx\n",
"FCM9e6WQhS7OATbBbqrcA7hQa/3fogtaBLrSRUa744A7gQFa62VFFLFoZKMLpdSZ2Loew4EhWuuf\n",
"FlfK4pDF38hWQBR4AdgFmKa1frzoghYYpdSmwOXAzlrrPdu53q3vzYryOIqVnr0SyFIXjVrrM7XW\n",
"UeB+7B9I1ZGlLlBKbQ9sX2z5ikk2ulBKnQCs0Vpfp7U+G7iuBKIWnCzfi7OBp7XWV2H/PtZPpV0d\n",
"7Ac8SMd1x7v1vVlRhoOO07NnchjwL7Dp2YFdlFLVWMipS11orS/MOKxj3ezD1USXuvC+RELAJVQ3\n",
"2fyNHAcMUEqdppSKAFXpeZGdLj4lnal7GNY7rzq01vfTeabxbn1vVprhyCU9e7WR9edUSvUBfgz8\n",
"ughylYJsdBEBLsn4Eulo5FXpZKOLbwBDtdbXA7cCf1dKVaM+stHFtcBeSqmrgTBwVpFkKze69b1Z\n",
"aYajKOnZK4RsdJEyGjcC52ut3y+SbMWmU10opUYAg4FjvXkfgF8qpXYvnohFI5v3YjHwHIDWeibQ\n",
"CGxRFOmKSza6mAbcorU+CzgaeNjzTmuNbn1vVprheA7Y0vsyBNgXeEQptaFSKvWCPIJ1UfHSs7/a\n",
"wyqC5U6XuvD+AG4CrtFav6KUOrpEshaaTnWhtZ6ttT5Ra32l1vpKr801WutqLCuczd9ICzASwAtH\n",
"1GFDNtVGNroYQfqzL8Ia0cbiilkacvnerCjD4a2CCQLXK6UuA17TWj8JnAuc7DWbgn1ZwsCZQFWu\n",
"FulCF0Gv2V3Yl+H3Sqkn+f/27jdEqiqM4/jXRbKgwnRdUyPRyGAxehFoBFqLRUpS0OJDRdCLwCAC\n",
"KSJRioV6kdofihKiiGChZH+CkJJuBWUhpCJBUfSHiMhii61dCiqwP/biPNMOq7R7Y2YHd36fV8O9\n",
"99x59rIzz5xz7n1OmQicdib5f0FEdEbEQ8BJ4MGIWNiSgJtoktdiO3BJfkaeBO6UdKIlATfRJK/F\n",
"fcC9EbEFeB7okzTckoCbKCJWA3cAF0bE1ryL6n9/b7pWlZmZVXJG9TjMzKz1nDjMzKwSJw4zM6vE\n",
"icPMzCpx4jAzs0qcOMzMrJIzrjquWTNFxE3ADmBIUk+r45mMiNgNrAfWSnq31fHY9Oceh1kdSXuB\n",
"x1odRxWSNlCefvZDWTYlnDjMTjUdC/6ZNYyHqqxtRcSlwLPAOZTPwqOSBnP3jIh4GrgSOBvYmPW+\n",
"OoCdwHLKQkhfAJsk/RYRXZSyFXPzfDskvRYRK4AXKEX2dlJKe68GPgaWAa9KuisXV9oK9Eu6PyJu\n",
"APqAE5TChHdLGsrY1wHbgBHgzeZdJbNTucdhbSkiZgL7KF/a11DKzg9ExNI8ZCWlauoq4ClgT7ZZ\n",
"CyyWtCrbzQU6s80rwEe5vRd4KSIWSzoKbKKstvdBzp08AfRQhpcezvbPAO9l0lgC7KbUkboWGAT6\n",
"M/ZOQJRk1kNJQIsaf5XMTs+Jw9rVSmAJpRBkrbz4EcoiRyeBzyV9kseK8qV/FeUX/uURcV32Pm4D\n",
"jkfEImANZX0LJH0PHMr9UIa/fpX0du7fLGkEeIOStADWAQfy9e3AsbrlO3cBa3IJ0BuBHyQdyXPt\n",
"A35v1IUxm4gTh7Wri4BRSX/XbRvO7QCjtY2S/qKU3F4g6TCwEdgMfA08QEkKtXb9EfFOViNeCtSv\n",
"7fDzaeLoZyxxbAAG6uLrrjvXnny/+cAC4Mdx5xmZ+E82awzPcVi7Og5cEBEddcmjC/g0X8+pHZhD\n",
"VLOBoYg4Hzgo6UAOaw0C31F6DgC9kn7KdrOY+DP2OvBiRFwPIOmX3P4Npcexvi6O2ZS5jiuAeePO\n",
"MwezKeIeh7Wrw8CXlKEpMgmsoMxTzKD82l+ex95KSQ7vA7dQehxI+gr4FujISeu3GOs9QFlE6z+f\n",
"Bcl1MAaAlylDYjW7KEuaXpzxdQEHM7b9QFdEXJ37bgbOxXeD2RTxehzWtjJZPEdZ8W0m8AhwFuUB\n",
"wGHgM6AbmMXYXVXLgMezzXnAh8A9kv6MiHmUZXrn51vsl7QtIropCekySsLqlfTvUFjedbUXWFg/\n",
"dJa9kD7gD8odXFtyor12V9V2ypDaIcqcyGjGeayhF8psHCcOMzOrxENVZmZWiROHmZlV4sRhZmaV\n",
"OHGYmVklThxmZlaJE4eZmVXixGFmZpU4cZiZWSVOHGZmVsk/v+y5vtZQCDsAAAAASUVORK5CYII=\n"
],
"text/plain": [
"<matplotlib.figure.Figure at 0x10bf58410>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"plt.plot([0, 1], [0, 1], \"k:\", label=\"perfect\")\n",
"\n",
"for intercept_correction in [False, True]:\n",
" clf = logistic.LogisticRegression(class_weight='balanced', intercept_correction=intercept_correction)\n",
" clf.fit(X=X, y=y)\n",
" fraction_of_positives, mean_predicted_value = calibration.calibration_curve(\n",
" y_true=y,\n",
" y_prob=clf.predict_proba(X)[:, 1],\n",
" n_bins=10\n",
" )\n",
" plt.plot(fraction_of_positives, mean_predicted_value, \"s-\", label='correction: {}'.format(intercept_correction))\n",
"plt.legend(loc='lower right')\n",
"plt.xlabel('observed')\n",
"plt.ylabel('predicted')"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": []
}
],
"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.10"
}
},
"nbformat": 4,
"nbformat_minor": 0
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment