Skip to content

Instantly share code, notes, and snippets.

@econpy
Created March 13, 2013 23:10
Show Gist options
  • Save econpy/5157353 to your computer and use it in GitHub Desktop.
Save econpy/5157353 to your computer and use it in GitHub Desktop.
Logistic Regression in Python
{
"metadata": {
"name": "Logistic Regression"
},
"nbformat": 3,
"nbformat_minor": 0,
"worksheets": [
{
"cells": [
{
"cell_type": "code",
"collapsed": false,
"input": [
"from pandas import read_csv,get_dummies,DataFrame,pivot_table\n",
"import numpy as np\n",
"import statsmodels.api as sm\n",
"import pylab"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 1
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"def cartesian(arrays, out=None):\n",
" arrays = [np.asarray(x) for x in arrays]\n",
" dtype = arrays[0].dtype\n",
"\n",
" n = np.prod([x.size for x in arrays])\n",
" if out is None:\n",
" out = np.zeros([n, len(arrays)], dtype=dtype)\n",
"\n",
" m = n / arrays[0].size\n",
" out[:,0] = np.repeat(arrays[0], m)\n",
" if arrays[1:]:\n",
" cartesian(arrays[1:], out=out[0:m,1:])\n",
" for j in xrange(1, arrays[0].size):\n",
" out[j*m:(j+1)*m,1:] = out[0:m,1:]\n",
" return out"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 2
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"def isolate_and_plot(variable):\n",
" # isolate gre and class rank\n",
" grouped = pivot_table(combos, values=['admit_pred'], rows=[variable, 'prestige'],\n",
" aggfunc=np.mean)\n",
" \n",
" # make a plot\n",
" colors = 'rbgyrbgy'\n",
" for col in combos.prestige.unique():\n",
" plt_data = grouped.ix[grouped.index.get_level_values(1)==col]\n",
" pylab.plot(plt_data.index.get_level_values(0), plt_data['admit_pred'],\n",
" color=colors[int(col)])\n",
"\n",
" pylab.xlabel(variable)\n",
" pylab.ylabel(\"P(admit=1)\")\n",
" pylab.legend(['1', '2', '3', '4'], loc='upper left', title='Prestige')\n",
" pylab.title(\"Prob(admit=1) isolating \" + variable + \" and presitge\")\n",
" pylab.show()"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 3
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"# read the data in\n",
"df = read_csv(\"http://www.ats.ucla.edu/stat/data/binary.csv\")\n",
" \n",
"# rename the 'rank' column because there is also a DataFrame method called 'rank'\n",
"df.columns = [\"admit\", \"gre\", \"gpa\", \"prestige\"]"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 4
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"# plot all of the columns\n",
"df.hist()\n",
"pylab.show()"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAXoAAAEICAYAAABRSj9aAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztnXtcVHX+/1+D4HYREBRBAYPlKheHUZLNzRVDMFMsL6FY\niKK1u3ZZ3X6mWW12UbDWb3ip3bbUTCvabmqlpK4MuplhKWlY4iooCGhyUcALl3n//kBGBpjbucyF\neT8fj3nAmZn3Zc75nA/D+7zP66MgIgLDMAzTY3GydgIMwzCMvPBEzzAM08PhiZ5hGKaHwxM9wzBM\nD4cneoZhmB4OT/QMwzA9HJ7orUhpaSmcnJyg0WgE2bu6uqK0tFTapBiG6XE4WzsBRjj19fXa32fP\nng1/f3+8/PLLVsyIYRhbhL/RMwzD9HB4opeBrKwsBAcHw83NDZGRkdi6dSsAoLW1Ff/v//0/eHl5\nISgoCF999ZWOXXx8PJ5//nn8/ve/h6urKyZNmoSLFy/ioYcegru7O0aMGIEzZ85o3+/k5IRTp07h\nX//6Fz744AO8+uqrcHV1xf3332/Rz8swhjh8+DBUKhXc3NyQkpKC6dOn4/nnn4darYafnx8yMzPh\n5eWFwMBAfPDBB1q7r776CiqVCu7u7hg8eDBefPFFK34KO4cYyfn444+psrKSiIg++ugjuv3226my\nspL+8Y9/UHh4OJWXl1NNTQ3Fx8eTk5MTtba2EhHR6NGjKSQkhE6fPk2XLl2iiIgICg4Opv/85z/U\n0tJCs2bNojlz5mjjKBQKOnXqFBERzZ49m55//nnLf1iGMcD169dp8ODBtGbNGmppaaHPPvuMevfu\nTc8//zyp1Wpydnamp556ipqamig/P59uv/12OnHiBBERqdVq+umnn4iI6OjRo+Tt7U1bt2615sex\nW/gbvQxMmzYNPj4+AICUlBSEhISgoKAAH3/8MRYuXAhfX194eHhg6dKloA5SQwqFAnPmzEFgYCDc\n3Nwwfvx4hIaG4p577kGvXr3w4IMP4siRI3rjEssWMTbGwYMH0draiieeeAK9evXC5MmTMWLECAA3\nx+vLL78MFxcX/OEPf8CECRPw73//GwAwevRoREZGAgCio6MxY8YM5OfnW+eD2Dk80cvAe++9B5VK\nBQ8PD3h4eOCnn37CxYsXUVFRAX9/f+37Bg8e3MXW29tb+/stt9yCAQMG6Gw3NDTImzzDSEhFRQV8\nfX11nvP39wcRQaFQwMPDA7feeqv2tTvuuAMVFRUAgO+++w5jxozBgAED0LdvX7z11luorq62aP49\nBZ7oJebMmTN49NFH8cYbb6Cmpga1tbWIiooCEWHgwIE4e/as9r0df+8OhUJhclxz3sswlmLgwIE4\nd+6cznNnz56FQqEAEaG2thZXrlzRvnbmzBntH4aZM2figQceQHl5Oerq6vCnP/1JcCuyo8MTvcQ0\nNjZCoVCgf//+0Gg02LhxI3766ScAbWWcNWvW4Ny5c6itrUVWVlYX+47lF3NKMd7e3jh9+rT4D8Aw\nEjJy5Ej06tUL69atQ0tLC7Zt24ZDhw4BuPnl5IUXXkBzczP279+Pr776Cg8++CAAoKGhAR4eHujd\nuzcKCgrwwQcf8BcagfBELzERERF46qmncNddd8HHxwc//fQT7r77bigUCjzyyCMYN24clEolYmNj\nMXXq1C4Dt+O2QqEw+no7c+fOxfHjx+Hh4YEpU6bI9OkYxjxcXFzw2WefYf369fDw8MD777+PiRMn\nonfv3gAAHx8feHh4YNCgQUhLS8Nbb72F0NBQAMCbb76Jv/3tb3Bzc8PLL7+M6dOnW/Oj2DUK4it4\nDMNYkLi4OMyfPx933HEH0tLSUFZWZu2UejwGv9Ffu3YNcXFxiImJQUREBJ555hkAQE1NDRITExEa\nGoqkpCTU1dVpbTIzMxESEoLw8HDs2rVL3uwZxoJkZmYiMjIS0dHRmDlzJq5fv27wXGDa2LdvH6qq\nqtDS0oJNmzbhp59+wr333mvttBwKgxP9Lbfcgry8PBQWFuLo0aPIy8vDf//7X2RlZSExMRHFxcVI\nSEjQ1pqPHz+Ojz76CMePH0dubi7mz5/PF0+YHkFpaSnefvttHD58GMeOHUNraytycnL0ngvMTU6c\nOIGYmBh4eHjg9ddfxyeffKLtLuOau4UwteG+sbGRYmNj6aeffqKwsDCqqqoiIqLKykoKCwsjIqIV\nK1ZQVlaW1mbcuHH07bffStLwzzDWpLq6mkJDQ6mmpoaam5tp4sSJtGvXLr3nAsPYEkZFzTQaDYYN\nG4ZTp07hz3/+MyIjI3H+/HntX2Rvb2+cP38eQFvP7O9+9zutrZ+fX5fWKv4LzkgBWfjSkqenJ556\n6ikMHjwYt956K8aNG4fExES950JHeMwzUiF03BvtunFyckJhYSHKy8uxb98+5OXl6bzeXWdI59c7\nQ0QWfaSnp3PMHhTTGpw6dQrZ2dkoLS1FRUUFGhoasGXLFp33GDoXrLmPpThG1s7B0e2JxI17k9sr\n3d3dMWHCBPzwww/w9vZGVVUVAKCyslJ796avr6/OFfTy8vIud8VZg4CAAI7Zg2Jag++//x4jR45E\nv3794OzsjClTpuDbb7+Fj49Pt+eClIjdx1IcI2vn4Oj2YjE40V+8eFHbRXD16lXs3r0bKpUKkyZN\nwqZNmwAAmzZtwgMPPAAAmDRpEnJyctDU1ISSkhKcPHlSq2vBMPZMeHg4Dh48iKtXr4KIsGfPHkRE\nRCA5Obnbc4FhbAmDNfrKykqkp6dDo9FAo9EgLS0NCQkJUKlUSElJwfr16xEQEKAVIYqIiEBKSgoi\nIiLg7OyMN9980ybqk3379uWYPSimNVAqlZg1axZiY2Ph5OSEYcOG4dFHH0V9fX2354KUiN3HUhwj\na+dgzN7NzRP19bUG39OdzLGrqwcuX64RHV9ue9GQhbFCSMrLy+OYPSimNcaQGMTmK3YfS3GMrJ2D\nMXsABJCBR56e5007Ntb+/ETixpHF74xtFzNiGKHY2xiyt3ztkbbKgZB9bD/HRsw4Yq0bhmGYHo5V\nJvqioiJBj3adanNRq9XSfgCOadWYjobYfSzFMbJ2DuI/g3XjW/s8MXrDlBzceec4uLi4m2Wj0TTj\nttuu4/z5M8bfzDASc+LECcyYMUO7ffr0abz88st4+OGHMX36dJw5c0Z7MdbqF94YphNWqdEDOQDM\nlRytgLt7LOrqhH2rZ3oO1q55azQa+Pr6oqCgAGvXrkX//v3x9NNPY+XKld2uM2DtfB0BrtEbhmv0\nDGMme/bsQXBwMPz9/bF9+3akp6cDANLT07F161YrZ8cwXbFK6cbSqNVqxMfHc8weEtPa5OTkIDU1\nFQBM0roBgNmzZ2vvjuzbty9iYmK0+629fqtvOzs726z3S22vVqtRWFiIBQsW2LT9Tdq34ztsFwJY\n0O3rtpJ/d59HrVajtLQUojHUe3n27FmKj4+niIgIioyMpNWrVxMR0QsvvEC+vr4UExNDMTExtGPH\nDq3NihUrKDg4mMLCwujrr7/uthcUyDHS89rd4xy5uw8U1EPqKP3ljhLTyLCVlevXr1P//v3pwoUL\nRETUt29fndc9PDy62IjN1xZ6uK2dA/fRixtHBi0rKyvpyJEjRERUX19PoaGhdPz4cVq2bBmtWrWq\ny/uLiopIqVRSU1MTlZSUUFBQELW2tnZJ1tITPdOzsOZEv3XrVho3bpx2OywsjCorK4mIqKKioluZ\nYmvm6ygYn+j1Pezn2IjJ1WCN3sfHBzExMQCAPn36YMiQIVrZYermosC2bduQmpoKFxcXBAQEIDg4\nGAUFBeL/7WAYG+HDDz/Ulm0A6NV9YhhbwuSLsaWlpThy5IhWb37t2rVQKpWYO3euVvisoqICfn5+\nWpvu9Ojb+AeAZTce2dDtcVUb3Far1V1qWMa2s7OzRdkL2W5/zlLxuottifjZ2dmyx8vOzsayZcuw\nbNkyzJ49G9aisbERe/bs0Vl8fcmSJdi9ezdCQ0Oxd+9eLFmyRPK4HfeHNextIQfxn8G68aU4BqIw\n5Wt/fX09DR8+nD7//HMiIjp//jxpNBrSaDT07LPPUkZGBhERPf7447Rlyxat3dy5c+nTTz/t8u8H\n1+g5phhMHLY2g9h8baE+bO0c5KvRO9+wNe/h6tr1WoyY/E1BzDgy2kff3NyMiRMnYvz48dqrxh0p\nLS1FcnIyjh07pu0fbv9Wc++99+LFF19EXFyc9v3cR8+Ixd760u0tX3tETB+9vfTfy9ZHT0SYO3cu\nIiIidCb5yspK7e+ff/45oqOjAbAePcMwjC1icKL/5ptvsGXLFuTl5UGlUkGlUmHnzp1YvHgxhg4d\nCqVSifz8fLz++usAdPXox48fbzN69Naoj3FMRipsoT5s7RysXaOXssbv5uapXXbSnIcYDN4wdffd\nd0Oj0XR5fvz48Xptli5diqVLl4pKimFskbq6OsybNw9FRUVQKBTYuHEjQkJCWOuGMYu2BVKElpmE\nwVo3jN1hrZp3eno6Ro8ejYyMDLS0tKCxsRHLly9nrRsbwJ5q9NbQ5eGJnrE7rDFxXrp0CSqVCqdP\nn9Z5Pjw8HPn5+fD29kZVVRXi4+Pxyy+/6LyHJ3r54YneMA4hauYotWtHiWkNSkpK4OXlhTlz5mDY\nsGF45JFH0NjYaJbWTfu9AObeeyD2XgUp7nUQey+K3Pa6qNH1XpxsI68b2+5s32HLzPy7+tAXX422\ne41m33iIQHBjpkDAffQcUyRWGLZ06NAhcnZ2poKCAiIi+stf/kLPPfcca93YiD1EaN2YNv90tjfv\nmHbM3/SY0sk1cOmGsTusUQqpqqrCXXfdhZKSEgDAf//7X2RmZuL06dPIy8uDj48PKisrMWbMGC7d\nWAEu3RjGIUo3DCMWHx8f+Pv7o7i4GECbJn1kZCSSk5NZ60ZCrNF66Ag4xETvKLVrR4lpLdauXYuH\nHnoISqUSR48exbPPPstaNxLb32w97PzI0/N8+8NoBmZmLK29tc8Tg330ZWVlmDVrFi5cuACFQoFH\nH30UTz75JGpqavT2DmdmZmLDhg3o1asX1qxZg6SkJIt8EIaRG6VSiUOHDnV5fs+ePVbIhmFMx2CN\nvqqqClVVVYiJiUFDQwOGDx+OrVu3YuPGjd32Dh8/fhwzZ87EoUOHcO7cOYwdOxbFxcVwcrr5jwPX\n6Bmx2FvN297ytSbWqLU7fI1enx69vnUyWY+eYRjG9jB5zdh2Pfq4uDi9vcMVFRVavXrAmB79zzd+\n7wsgBrrrO6Kb7dC2rRu1LkuuVylku/05S8XrGMtS8QBp1iM1tl1YWKhd80CS9TPtDLVa3Lq8Yu1t\nIwc1bs4J9mcvxTEQhSk9mPX19TRs2DCtHr2+3mHWo+eYlsDEYSs5d9xxB0VHR1NMTAzdeeedRERU\nXV1NY8eOpZCQEEpMTKTa2toudmLztXYPuyVzgN4ec3198Kb2wzt2H73Rrpvm5mZMnToVaWlp2tax\n9tu9gTbJ4gEDBgAAfH19UVZWprUtLy+Hr6+vhH+WhGGNv6Qcs+ehUCigVqtx5MgRbUkyKysLiYmJ\nKC4uRkJCQhedGykQu4+lOEbWz8G+7a19ngjSo9e3Tibr0TM9Hep0MUzf9SqGsSXM1qPPzc3V2zvM\nevQcsyejUCgwduxYxMbG4u233wYA1rqRwV6YVg2MvM5aNxYFXKPnmCKxwrAlIqKKigoiIrpw4QIp\nlUrat28fa91IbA+u0ctSo2etG8busIW+9BdffBF9+vTB22+/DbVazVo3EsF99PLEdAgJBIYRy5Ur\nV1BfXw8AaGxsxK5duxAdHa33ehXD2BIOMdE7Su3aUWJag/Pnz2PUqFGIiYlBXFwcJk6ciKSkJNa6\nsVgO9m1v7fPE5BumGMaRCQwMRGFhYZfnPT09WeuGsXm4Rs/YHfZW87a3fK0J1+jliekQpRuGYRhH\nxiEmekepXTtKTEfD+vVxW8jBvu2tfZ4YnOgzMjLg7e2N6Oho7XPLli2Dn5+f9gaqnTt3al/LzMxE\nSEgIwsPDsWvXLvmyZhgr0draCpVKheTkZABATU0NEhMTERoaiqSkJK34GsPYEgZr9Pv370efPn0w\na9YsHDt2DEBb/7Crqyv++te/6rzXFC16gGv0jHisWfP+v//7P/zwww+or6/H9u3b8fTTT3e7NoOt\n5Gst3Nw8b6wWJQSu0Usd0+A3+lGjRsHDw6PL890FYy16pqdTXl6OHTt2YN68edpzgLVuukf/koDG\nHowcCGqvXLt2Ld577z3ExsZi1apV6Nu3rxla9ADr0bMevT3q0S9cuBCvvfYaLl++rH3OHK2bgIAA\nAEDfvn3N2mdi97EUx0jIOXQTNYBCAAs6bAPGz/mO26bYw8DrcsS/sWXm/rvpw1j89t9LIRpjGgkl\nJSUUFRWl3T5//jxpNBrSaDT07LPPUkZGBhGZpkV/o0zEWjccUxQmDFvJ+eKLL2j+/PlE1PaZJ06c\nSET612boiNh87VHrBl30XIxp1RjTnrE1rRtnIf+udHhYVuvG7G/07drzADBv3jztRSlb1aIHHEen\n3VFiWoMDBw5g+/bt2LFjB65du4bLly8jLS1NuzZDu9ZNx/NDKqyvBc968l3tWyC81GR5RV+z2ysr\nKyu1v3/++efajhzWomd6MitWrEBZWRlKSkqQk5ODe+65B5s3b2atG8YuMDjRp6amYuTIkThx4gT8\n/f2xYcMGLF68GEOHDoVSqUR+fj5ef/11ALarRQ84Tn+5o8S0BdrHNmvdmOyB7a2IwdLNhx9+2OW5\njIwMve9funQpli5dKj4rhrFhRo8ejdGjRwNgrRvGPmCtG8busLe+dHvLVwrsSbPGsnbiYgodRw4h\ngcAwjDDc3DyhUCjMfjC2hUNM9I5Su7a3mDyJmIY1a/Q3b3zKg7gbn4TnwPbicYiJnrFN7OnuyWvX\nriEuLg4xMTGIiIjAM888A4C1bhj7gGv0jNWwhuaHGK5cuYLbbrsNLS0tuPvuu/H3v/8d27dv79Fa\nN1xrl9pOXEyu0TOMzNx2220AgKamJrS2tsLDw4O1bhi7wCEmenurXXNM20Sj0SAmJgbe3t4YM2YM\nIiMjzdK6WbZsGZYtW4bs7Gyd/aZWqw1um/t+qe3b6svZnbY7v955G522zbU3N37neObamxvfWDxD\n9p196LNXA1gGYPaNhwgM6SPMmTOHBgwYoKN1U11dTWPHjqWQkBBKTEyk2tpa7WsrVqyg4OBgCgsL\no6+//tqABgZr3XBMc3RGpNP8kIK6ujqKi4ujvXv39nitm5vHyFStGn0aMj1N68bcsZsnwlb8uDdo\nuW/fPjp8+LDORL9o0SJauXIlERFlZWXR4sWLiYioqKiIlEolNTU1UUlJCQUFBVFra6uegWPZiZ6x\nTex1oicieumll+i1116jsLAwqqysJCKiiooKCgsL6/JeW8hXKGKOEdtJH1MoZuvR66tJsh4905O5\nePGitqPm6tWr2L17N1QqFWvdMHaB2eqV+mqSrEevu93+nKXi2asevela4O1ti6WwBpWVlUhPT4dG\no4FGo0FaWhoSEhKgUqmQkpKC9evXIyAgAP/+978lj61Wq0WpR4q1v+EF4hQg2V68gqYIjH3l76xH\nr68myXr0HNNcYIV/Ya2B2Hy5Rm+KPdfoDWF21027/jYAHf1t1qPnmIw82IIeve3pwTuavTjMnuj1\n1SRZj55hGMY2MUuPfuPGjXr1t1mPnmMy8mBNrZsOXtjeqvbiMFuPHoBe/W3Wo2d6KmVlZZg1axYu\nXLgAhUKBRx99FE8++SRqamowffp0nDlzRnsxtm/fvtZOl2F0YK0bxmrYk9ZNVVUVqqqqEBMTg4aG\nBgwfPhxbt27Fxo0bWeume0u2kyGm0HHkEBIIDCMWHx8fxMTEAAD69OmDIUOG4Ny5c6x1w9gFDjHR\nO0rt2lFiWpvS0lIcOXIEcXFxrHWjdxudtlnrpquNMXs1LKJ1IwcA99FzzDZghX5isdTX19OwYcPo\n888/JyL995V0RGy+3EfPffRixz3X6BmrYU81egBobm7GxIkTMX78eO2d1uHh4VCr1fDx8UFlZSXG\njBmDX375RTdbrtGznUQxhY4jhyjdMIxYiAhz585FRESEdpIH9N9XwjC2hENM9I5Su3aUmNbgm2++\nwZYtW5CXlweVSgWVSoXc3Fy995VICffRs71YzBY1YxhH5O6774ZGo+n2NX33ldgSbm6eN9boZRwR\nwTX6gIAAuLm5oVevXnBxcUFBQYFJN49wjZ5px95q9EKxhRo919ptxU5cTIvX6BUKBdRqNY4cOaLV\nnc/KykJiYiKKi4uRkJDQ5cYRhmEYxvKIqtF3/uti+s0j/0Bbf+gytPWXqju8pja4bW4PsFqtRnZ2\ntih7Idvtz1kqXnexLRFfmh5tY9vZuDleZsPR4Bo924tGaF9mYGAgxcTE0PDhw+lf//oXEen2FGs0\nmi49xjfKRMR99ByTyDr9xNZAbL5S9NGL2dem9bEb60831Z776OUY94ItKyoqiIjowoULpFQqad++\nfSbfPMJrxjJE1hnwYpgzZw4NGDBAZyGe6upqGjt2LIWEhFBiYiLV1tZ2sbOFP0ziJ3q2k8bOOuNe\ncOlm4MCBAAAvLy9MnjwZBQUFehclYZiewJw5c5Cbm6vzHF+XYuwBQRP9lStXUF9fDwBobGzErl27\nEB0dbbM3j1ij15tj9jxGjRoFDw8PnecsIWrGNXq2F4ugPvrz589j8uTJAICWlhY89NBDSEpKQmxs\nrOwLJTOMLWGOqFlAQAAAoG/fvmYtqF5YWGjwdVPtb6K+8TPejO1CM9/fOZ4Qe3Pjw8DrcsQ3FM+Q\nfftzxuK3/14KsbDWDWM17LGPvrS0FMnJyTh27BgAwMPDA7W1N29E8vT0RE1NjY4N99GznVQxhY4j\nh5BAYBi54OtSjD3gEBO9o9SuHSWmLWGJ61Jco2d7sTjERM8wUpCamoqRI0fixIkT8Pf3x8aNGy0i\nasYwYuEaPWM17LFGLwSu0bOdVDG5Rs8wDMN0i0NM9I5Su3aUmI5Gx33s5uYJhUJh9kOCLNjeqvbi\ncAg9+sLCQm1PMce0/5g9gYKCAuzcudOk9x48eBD5+fkAcENT3tx/37MBLDTTpjOd+8DZ3rL24pB8\nos/NzcWCBQvQ2tqKefPmYfHixVKHMJu6ujqO2YNi2iLmjvu3396Md945CSDOBO9XkZvb/aInpiHF\nMRLrg+2tiaQTfWtrKx5//HHs2bMHvr6+uPPOOzFp0iQMGTJEyjAMY1MIH/cTADxhQoRlNx4A8JKI\nTBlHRdIafUFBAYKDgxEQEAAXFxfMmDED27ZtkzKEIEpLSzlmD4ppa8g/7kutbG8LOTi6vUhE6592\n4OOPP6Z58+Zptzdv3kyPP/64znvQVmDkBz9EPWwJY+Pe2vuKHz3nIRRJSzemXN0nO+p/ZhhTMDbu\necwz1kbS0o2vry/Kysq022VlZfDz85MyBMPYHDzuGVtH0ok+NjYWJ0+eRGlpKZqamvDRRx9h0qRJ\nUoZgGJuDxz1j60haunF2dsa6deswbtw4tLa2Yu7cudxxw/R4eNwzto7kd8aOHz8eJ06cwLp167Bp\n0yaEhIRg5cqV3b73ySefREhICJRKJY4cOSI6dm5uLsLDw/XGfP/996FUKjF06FD8/ve/x9GjR2WP\n2c6hQ4fg7OyMzz77zCIx1Wo1VCoVoqKiJLmhyVjMixcv4t5770VMTAyioqLw7rvvio6ZkZEBb29v\nREdH632P1GNIKFFRURg0aBB+85vf4P3338eaNWu6vMfQ+CsrK8OYMWMQGRmJqKiobu3b0TeWTPWh\nb2yYYm/oOF+7dg1xcXGIiYlBREQEnnnmmW7j6ztmptgb2oemxge634em2hs6t0zxYcq50traCpVK\nheTk5G5zMHvci2g20EtLSwsFBQVRSUkJNTU1kVKppOPHj+u856uvvqLx48cTEdHBgwcpLi5O9pgH\nDhyguro6IiLauXOnRWK2v2/MmDE0YcIE+uSTT2SPWVtbSxEREVRWVkZERL/++qvsMV944QVasmSJ\nNp6npyc1NzeLirtv3z46fPiwzmLcHZF6DImhsrKSjhw5QkRE9fX1FBoaatb4M8WeyPBYMsWHobFh\nir2x49zY2EhERM3NzRQXF0f79+/XsTd2zIzZGzuHjdkTGd6HxuxNObeM+TDlXFm1ahXNnDmTkpOT\nu/gXMu5l0boxpa+441qbcXFxqKur07sMm1Qx77rrLri7u2tjlpeXC45nakwAWLt2LaZNmwYvLy9R\n8UyN+cEHH2Dq1KnaC4L9+/eXPebAgQNx+fJlAMDly5fRr18/ODuLqwx2t0ZrR6QeQ2Lw8fFBTEwM\nAKBPnz4YMmQIKip0lVYNjT9T7AHDY8kUH4bGhin2xo7zbbfdBgBoampCa2srPD09deyNHTNj9sbO\nYWP2gOF9aMzelHPLmA9j+7C8vBw7duzAvHnzuu3YEjLuZZnoz507B39/f+22n58fzp07Z/Q9YiZe\nU2J2ZP369bjvvvsExzM15rlz57Bt2zb8+c9/BmBaC6rYmCdPnkRNTQ3GjBmD2NhYbN68WfaYjzzy\nCIqKijBo0CAolUqsXr1aVEyheYn94y0FpaWlOHLkCOLi9MsbGBp/+uzNGUv6fJg6NvTZGzvOGo0G\nMTEx8Pb2xpgxYxAREdHlMxg6ZsbsO9LdPjQlvqF9aMzelP1nzIexfbhw4UK89tprcHLqfnoWMu5l\nmehNncw6/7USMwmaY5uXl4cNGzYYrKlLFXPBggXIysrSapJ39xda6pjNzc04fPgwduzYga+//hov\nv/wyTp48KWvMFStWICYmBhUVFSgsLMRjjz2G+vp6wTFNRcoxJAUNDQ2YNm0aVq9ejT59+nT7HkPj\nz5C9qWPJkA9TxoYhe2PH2cnJCYWFhSgvL8e+ffu6VTc1dMxMsQf070Nj9sb2oTF7U/afMR+G9uGX\nX36JAQMGQKVSGZwrzB33skz0pvQVd35PeXk5fH19ZY0JAEePHsUjjzyC7du3GywLSBXzhx9+wIwZ\nMxAYGIgOIdn4AAAgAElEQVRPP/0U8+fPx/bt22WN6e/vj6SkJNx6663o168f/vCHP+DHH3+UNeaB\nAwfw4IMPAgCCgoIQGBiIEydOCI4pJC+xY0gszc3NmDp1Kh5++GG9SwoaGn/G7E0ZS8Z8GBsbxuxN\nPc7u7u6YMGECvv/+e53nTT1m+uwB085hffamno/67M05t/T5MLQPDxw4gO3btyMwMBCpqanYu3cv\nZs2apWMvaNwbreILoLm5mX77299SSUkJXb9+3ejF2G+//Vb0hTRTYp45c4aCgoLo22+/FRXLnJgd\nmT17Nn366aeyx/z5558pISGBWlpaqLGxkaKioqioqEjWmAsXLqRly5YREVFVVRX5+vpSdXW14Jjt\nlJSUmHQxVooxJAaNRkNpaWm0YMECve8xNP5Mse9Id2PJFB+GxoYp9oaO86+//kq1tbVERHTlyhUa\nNWoU7dmzR8fe0DEzxd7QPjTFviOd96Ep9sbOLVN8mHquqNVqmjhxYpfnhYx7WfTo9fUVv/XWWwCA\nP/7xj7jvvvuwY8cOBAcH4/bbb8fGjRtlj/nSSy+htrZWW59zcXFBQUGBrDGlxpSY4eHhuPfeezF0\n6FA4OTnhkUceMVjrlCLm0qVLMWfOHCiVSmg0Grz66qvdXggzh9TUVOTn5+PixYvw9/fHiy++iObm\nZm1MqceQGL755hts2bIFQ4cOhUqlAtD2L/rZs2e1+Roaf6bYS5GDobFhir2h41xZWYn09HRoNBpo\nNBqkpaUhISHB5PPeFHtD+9AUe0OYYm/s3DLFhznnSntJRuzcafE1YxmGYRjL4hBLCTIMwzgyPNEz\nDGPTuLq68roHIuGJnmEYmyE+Ph7r16/Xea6+vh4BAQHWSaiHwBO9DdPS0mLtFBjGbMSMW2vfB9FT\n4YneChw+fBgqlQpubm5ISUnB9OnT8fzzz0OtVsPPzw+vvvoqBg4ciLlz54KIkJWVheDgYPTv3x/T\np09HbW2ttT8C44AEBAQgKysLkZGR8PT0REZGBq5fv272uL127Roefvhh9O/fHx4eHhgxYgQuXLiA\nZ599Fvv378fjjz8OV1dXPPnkkwDabkA6ffo0AKC6uhrJyclwd3fHiBEj8Nxzz2HUqFHaHH/55Rck\nJiaiX79+CA8Px8cff2z5HWWD8ERvYZqamjB58mRkZGSgtrYWqamp2Lp1KxQKBRQKBc6fP4/a2lqc\nPXsWb731FtasWYPt27dj3759qKyshIeHBx577DFrfwzGQfnggw+wa9cunDp1CsXFxXjllVfMHreb\nNm3C5cuXUV5ejpqaGrz11lu49dZbsXz5cowaNQpvvPEG6uvru1XPfOyxx+Dq6orz589j06ZNeO+9\n97T/BTQ2NiIxMREPP/wwfv31V+Tk5GD+/Pn4+eefLbqPbBKjnfaMpOTn55Ovr6/Oc3fffTc9//zz\npFarqXfv3nT9+nXta0OGDKH//Oc/2u2KigpycXGh1tZWi+XMMEREAQEB9NZbb2m3d+zYQUFBQWaN\n25aWFtqwYQONHDmSjh492iVGfHw8vfPOOzrPKRQKOnXqFLW0tJCLiwsVFxdrX3vuuefo7rvvJiKi\nnJwcGjVqlI7to48+Si+++KK4D94DkOWGKUY/FRUVXW5X9vf312pXeHl5oXfv3trXSktLMXnyZB2B\nI2dnZ5w/fx4DBw60TNIMc4OOYlqDBw/WqluaOm4vXLiAtLQ0lJWVYcaMGairq8PDDz+M5cuXaxUc\n9dXpf/31V7S0tHQR9GrnzJkz+O6773RkEVpaWrpICDgiXLqxMAMHDuyi/Hj27Fnt4O48yAcPHozc\n3FzU1tZqH1euXOFJnrEK7XfJtv8+aNAgAOaNW2dnZ/ztb39DUVERDhw4gC+//BLvvfdet3464uXl\nBWdn5y66Sx1jjh49WidmfX093njjDUk+uz3DE72FGTlyJHr16oV169ahpaUF27Ztw6FDhwB0VaQD\ngD/96U9YunSp9gT79ddfRYmiMYxQiAhvvvkmzp07h5qaGixfvhwzZszo9r2Gxq1arcaxY8fQ2toK\nV1dXuLi4oFevXgAAb29vnDp1qlufvXr1wpQpU7Bs2TJcvXoVv/zyCzZv3qz94zBhwgQUFxdjy5Yt\naG5uRnNzMw4dOoRffvlF6l1hd/BEb2FcXFzw2WefYf369fDw8MD777+PiRMnonfv3toLsh35y1/+\ngkmTJiEpKQlubm646667ROnzMIxQFAoFZs6ciaSkJAQFBSEkJATPPfcciMiscVtVVYUHH3wQ7u7u\niIiIQHx8PNLS0rR2n3zyCTw9PbFgwYIuOaxbtw6XLl2Cj48P0tPTkZqaqi0Zubq6YteuXcjJyYGv\nry8GDhyIZ555Bk1NTTLvGduHtW5sgLi4OMyfP1+7agzD2CKBgYFYv3497rnnHmunomXx4sW4cOGC\nVQXt7AHB3+hXr16N6OhoREVFaVdIqampQWJiIkJDQ5GUlIS6ujrJEu1J7Nu3D1VVVWhpacGmTZvw\n008/4d5777V2Wgz0L0i+du1aDBkyBFFRUVi8eLH2+czMTISEhCA8PBy7du2ydLoOx4kTJ3D06FEQ\nEQoKCrBhwwZMnjzZ2mnZPkJadY4dO0ZRUVF09epVamlpobFjx9L//vc/WrRoEa1cuZKIiLKysmjx\n4sXS9Ab1MP71r3+Rt7c39enTh5RKJe3YscPaKTE36G5B8r1799LYsWOpqamJiIguXLhARERFRUWk\nVCqpqamJSkpKKCgoqEe3vQYEBOi0TFqDQ4cOUXBwMN12220UGBhIWVlZVs3HXhBUuvnkk0+Qm5uL\nd955BwDwyiuvoHfv3tiwYQPy8/Ph7e2NqqoqxMfH84UQxu4oLS1FcnIyjh07BgBISUnBn/70py4l\ni8zMTDg5OWm/4d97771YtmwZfve731k8Z4YxhKA++qioKDz77LOoqanBLbfcgh07diA2Nhbnz5+H\nt7c3gLar592tTM5aFowUCPh+IpiTJ09i3759WLp0KW655Rb8/e9/R2xsLCoqKnQmdX0L0vOYZ6RC\n6LgXVKMPDw/H4sWLkZSUhPHjxyMmJkbbHtVOdx0k7dCNRXnleKSnp9u1/57wGeT2b2laWlpQW1uL\ngwcP4rXXXkNKSore91pqzMuxj6X2aQ852pNPMQi+GJuRkYHvv/8e+fn58PDwQGhoqLZkA7QtqTVg\nwABRyQlBbjlTS8il2vtn6GmSsn5+fpgyZQoA4M4774STkxMuXrxo1cXJ5djHUvu0hxztyacYBE/0\nFy5cANB2d9xnn32GmTNnYtKkSdi0aROANuGi7laRZxh744EHHsDevXsBAMXFxWhqakL//v0xadIk\n5OTkoKmpCSUlJTh58iRGjBhh5WwZpiuCtW6mTZuG6upquLi44M0334S7uzuWLFmClJQUrF+/HgEB\nAfj3v/8tZa4m0bdvX7v2b4kY9u5fTtoXJK+uroa/vz9eeuklZGRkICMjA9HR0ejdu7f2dv2IiAik\npKQgIiICzs7OePPNNy1Wj5djH0vt0x5ytCefoiALI3fIvLw8u/ZviRj27t8Kw1YUcuQrxz4W69PV\n1YMAmPVwdfWwaI727FPMOLL4nbEKhUL0hQXGsbG3MWRv+Qql7b8Zcz+nY+wbKRAzjljrhmEYpocj\neKLPzMxEZGQkoqOjMXPmTFy/ft0mJBDUarVd+5cihpubp7a91ZyHm5unTeTPGEeOfSy9T6n92cvn\ntr1zQNBEX1pairfffhuHDx/Wyo3m5OQgKysLiYmJKC4uRkJCArKysqTOlzGB+vpaGC6N5nX7fJud\nY6NP6wYAVq1aBScnJ9TU1GifY60bxi4QUtivrq6m0NBQqqmpoebmZpo4cSLt2rWLwsLCqKqqioiI\nKisrKSwsrIutwJCMGQAggAQ87OPYyJlnd1o3RERnz56lcePGUUBAAFVXVxOR6Vo39rJfxSJs3DnG\nvpECMftKUHulp6cnnnrqKQwePBi33norxo0bh8TERJMkEABg9uzZ2hsK+vbti5iYGMTHxwO4+S8P\nb4vbvkn7drxJ27aSf8ftwsJCbRmwtLQUcjJq1KhuY/z1r3/Fq6++ivvvv1/73LZt25CamgoXFxcE\nBAQgODgYBQUFrHXD2ByCum5OnTqF5ORk7N+/H+7u7njwwQcxdepUPPHEE6itvfnvv6enp86/uYD8\nHQhqtVo7SdijfyliGO9+UOPmRK9jKcmxkXsfyT2GOouabdu2DWq1Gq+//joCAwPxww8/wNPTE088\n8QR+97vf4aGHHgIAzJs3D+PHj8fUqVO75Jueni7pl5vCwkLtwhxSfjmIj48XbD9mzBi0jbub/trG\nmbrD7+i0rUBeXp7J8TrnKubztm9nZ2dL/mVTiuPT/nv7F49NmzYJH/dC/g3IycmhuXPnarffe+89\nmj9/PoWHh1NlZSURta36bo3Sjb33iEsRA0b/hc6T9d9oe++jLykp0ZZuGhsbacSIEXTp0iUiapPq\nvXjxIhERPf7447Rlyxat3dy5c+nTTz+1SL622PvdddzpG2fCx5wtfm5L+RQzjgSLmh08eBBXr14F\nEWHPnj2IiIhAcnKy1SUQ5P62Lbd/y8SQ178l9pGlOHXqFEpLS6FUKhEYGIjy8nIMHz4c58+ft6rW\njRz7WHqfUvuzl89te+eA4BumXn31VWzatAlOTk4YNmwY3nnnHdTX1yMlJQVnz57VSiB0vhXYUW4e\nsSbCblwB7OXmFUuXbjrSsXRz/PhxzJw5EwUFBTh37hzGjh2L//3vf11kEBxlzPMNU/IiahyJ/4fC\nPOQO2ZNKN0JuKb/54NKNEGbMmEEDBw6k3r17k5+fH23YsEHn9cDAQG3XDRHR8uXLKSgoiMLCwig3\nN9di+dpiuaHruOPSjZSIGUeCRc0Y+bnZD28uvNCFUD788EODr58+fVpne+nSpVi6dKmcKTGMaASV\nbk6cOIEZM2Zot0+fPo2XX34ZDz/8MKZPn44zZ85w6UYCxJRguHRjO9hbvkLh0o28iBlHokXNNBoN\nfH19UVBQgLVr16J///54+umnsXLlStTW1na5O9ZRBr0U8ETfPfY2huwtX6HwRC8vVhU127NnD4KD\ng+Hv74/t27cjPT0dAJCeno6tW7eKdW82XW8Ysi//lokhr39L7CNHR459LL1Pqf3Zy+e2vXNAdI0+\nJycHqampAGATd8YWFhaKsre2//abLW62Z6lv/DR3GwZeL9RrL3X+Uvmz1J2xGRkZ+OqrrzBgwABt\n182iRYvw5Zdfonfv3ggKCsLGjRvh7u4OoE3rZsOGDejVqxfWrFmDpKQkWfNjGCGIKt00NTXB19cX\nx48fh5eXFzw8PKx+Z2xPgks33SPnGNq/fz/69OmDWbNmaSf63bt3IyEhAU5OTliyZAkAICsrS9te\neejQIW17ZXFxMZycdP9RdpQxz6UbebFa6Wbnzp0YPnw4vLy8AMAmFgdnGDGMGjUKHh4eOs8lJiZq\nJ++4uDiUl5cD0K91w9gW5sp2SyXXbUuIKt18+OGH2rINAO3i4IsXL7banbFq1roxJQLkvDvWEvvI\nWmzYsEE75isqKnQEzPz8/HDu3Llu7aQuV9qi1s1NOm7Ho/vy4c1tc+J1ztUU+7Y25Ty98YFsADHa\n7fp6hc4Yttbxaf9dknKl0Ab8hoYG6tevH12+fFn7XHV1NSUkJFBISAglJiZSbW1tFzsRIU2iJ90w\nBRFyw3zDlHA6at105JVXXqEpU6Zot1nrRpeu4842bpgy/3wQf7xs7YYpXjPWhuEaffdYQwLh3Xff\nxdtvv43//Oc/uOWWWwBA2zrcXre/99578eKLLyIuLs6i+doKtlqjNz8v2zxevGYsIxHOVl2C0FbJ\nzc3Fa6+9hm3btmkneaCtVJmTk4OmpiaUlJTg5MmTGDFihBUzZZju6XETvdz9q5boj5U/hj7/LRAi\nq9N5CUJb6yE2h9TUVIwcORInTpyAv78/NmzYgCeeeAINDQ1ITEyESqXC/PnzAQARERFISUlBREQE\nxo8fjzfffLOLoJlcyLGPpfcptT+5xpb0Pm3tHBB8Mbaurg7z5s1DUVERFAoFNm7ciJCQEKMSCAxj\ny3SndZORkaH3/ax1w9gDgmv06enpGD16NDIyMtDS0oLGxkYsX76cJRAkxBo1enuo7dvbGLK3fIXC\nNXp5sbjWzaVLl6BSqboo+YWHhyM/P1/bTx8fH49ffvlFsmQdDZ7o9USzszFkb/kKhSd6eREzjgSV\nbkpKSuDl5YU5c+bgxx9/xPDhw5GdnW0TEghyrP9oSf+de3Dlk0AQ67/7eHL0eFtSAsFekONeBel9\nqiH1/Rry3KOhhn3kKQIhPZmHDh0iZ2dnKigoICKiv/zlL/Tcc89R3759dd7n4eHRxVZgSJPhPnpx\nffRC41lyH8k5hubMmUMDBgzQ6aOvrq6msWPHdnt/yIoVKyg4OJjCwsLo66+/tli+3EcvJi/H66MX\nZFlZWUkBAQHa7f3799N9991nE4uD9yTkm+ilt7P0fpGLffv20eHDh3Um+kWLFtHKlSuJiCgrK4sW\nL15MRERFRUWkVCqpqamJSkpKKCgoiFpbWy2ary0hbPzIv2/Mz8s2j5eYvAS1V/r4+MDf3x/FxcUA\n2qSKIyMjbWJxcIYRQ3daN/rkt1nrhrEXBLdXrl27Fg899BCampq00q2tra1ISUnB+vXrte2Vlkbu\n2pglam/yx1CDtW5MR9+1J9a66WrfRsft+A7b8Z1eNz/fzrmaYn8zZvfxO2vdtPuw9vFp/92qWjdC\nkTsk1+iF1CSlLd3Yc42eqKvWjb5rT6x1o0vX8aNvnAkvk3CNXhg97s5Yub9JWuKbqvwx5PXfk77N\nA/rlt319fVFWVqZ9X3l5OXx9fS2Skxz7WHqfUvuTa2xJ79PWzoEeN9EzjNS0y28DuteeWOuGsRcE\nT/QBAQEYOnQoVCqVdnDX1NQgMTERoaGhSEpK0vY+WxK5NSYsoWEhfwx5/duazoc5dNa62bhxI5Ys\nWYLdu3cjNDQUe/fu1apVstaNUY8S+2OtG6EIvhirULSJ83t63lQuzMrKQmJiolYCISsrq4sEAsPY\nMt1p3QBtnWXdwVo3jD0gWOsmMDAQ33//Pfr166d9jiUQusfNzbOLwqPpCNlXLIFgS9hbvkJhCQR5\nsbgEQnvQsWPHolevXvjjH/+IRx55xCYkEGxxu22Sbz9A6hs/403YVpj5/o7bMPK61Ns3tmSShGAJ\nBIYRgdB2nYqKCiIiunDhAimVStq3bx9LIOgBZrct5mnbvLi9svv9aU/IkS+3V4rJi9srTWbgwIEA\nAC8vL0yePBkFBQV629AYpieQmZmJyMhIREdHY+bMmbh+/bpNNCAwjDEE1eivXLmC1tZWuLq6orGx\nEUlJSXjhhRewZ88e9OvXD4sXL0ZWVhbq6upYjx72JTfMNfruKS0txT333IOff/4Zv/nNbzB9+nTc\nd999KCoq4jUYbsA1enmxeI3+/PnzmDx5MgCgpaUFDz30EJKSkhAbG2t1CQSGkQM3Nze4uLjgypUr\n6NWrF65cuYJBgwYhMzMT+fn5ANp0cOLj47nTjLE5BE30gYGBKCws7PK8p6en3jY0S9ETtG7k1qJh\nrRvz8fT0xFNPPYXBgwfj1ltvxbhx45CYmGi1BgTWumGtG7MQf4nAPOQOyRdjrXEx1vmGrXkPV9eu\nF+tN3Z+W5n//+x8NGTKELl68SM3NzfTAAw/Q5s2brdaAwBdjxeTleBdjBffRC8VR6pUdcYQavSVr\n+9YYQx999BF2796Nd955BwCwefNmHDx4EHv37kVeXh58fHxQWVmJMWPGOOy9I1yjlxcx40iU1k1r\naytUKhWSk5MB2IYEAsPIQXh4OA4ePIirV6+CiLBnzx5ERETwGgyMXSBqol+9ejUiIiK0+h7tEgjF\nxcVISEiwykWpnqB1I7cWjf37tzxKpRKzZs1CbGwshg4dCgB49NFH9ergyA1r3UjqVXqPNqZ1I3ii\nLy8vx44dOzBv3jztvxP6VuJhmJ7A008/jaKiIhw7dgybNm2Ci4uLtgGhuLgYu3btQt++fa2dJsN0\nQbAEwsKFC/Haa6/h8uXL2udsQQKh/Tm5JA2E+r9J+3a8kW1z32+uvVj/xvyZ5p8lEIQhR1eT9D6l\n9sd69EIRdDH2yy+/xM6dO/HGG29ArVZj1apV+OKLL+Dh4YHa2pviXZ6enqipqdEN6CAXpjrCF2P1\n29nLxVgx2Fu+QuGLsfJi8YuxBw4cwPbt2xEYGIjU1FTs3bsXaWlpNiGBwDV6R/DPcI1eUq/SezSS\np5ubJxQKhVkPMQia6FesWIGysjKUlJQgJycH99xzDzZv3qx3JR6GYRjmJjcVbc15CEd0H31+fj5W\nrVqF7du3o6amBikpKTh79qxWAqHzxSlH+Te2I1y60W9nT6Wburo6zJs3D0VFRVAoFNi4cSNCQkIw\nffp0nDlzxuHHPJduzIhg4X3FN0xZAJ7o9eECoEWAHawyhtLT0zF69GhkZGSgpaUFjY2NWL58OYua\n3YAnejMi8EQvDlvUujH/oKpxc+EROSbedv/m2pkaT59/6eJZeuK8dOkSVCoVTp8+rfO8qauqpaen\nS9pp9s9/vo2KiraOt7q66ht++xncTkwcg1WrlsumdTNmzBi0Hc+b/oxr3SiQl5dncjwhWjdt51+e\nnvhAV62btpzk1Lq5ua+6y+fmZ2z7vfTG75t4om+HJ3qYYNfu31w7U+Pp8y9dPEtP9IWFhfjjH/+I\niIgI/Pjjjxg+fDiys7Ph5+en7TQjInh6eup0ngHyjPk//GEc9u8fAiDBRItyDByYjYqKE3rfIfbc\n6TrO1TA+DszbN/Kcf2ro5in+eBnL0y6+0V+7dg2jR4/G9evX0dTUhPvvvx+ZmZmoqanhemU3cOlG\nejtLj6Hvv/8ed911Fw4cOIA777wTCxYsgKurK9atW2eVluKpU9Px2Wf3AEg30eIEBg6cZHCiFwuX\nbsyIYOF9Jajr5pZbbkFeXh4KCwtx9OhR5OXl4b///a9NSCAwjBz4+fnBz88Pd955JwBg2rRpOHz4\nMHx8fKzeUswwxhAsgXDbbbcBAJqamtDa2goPDw+TJRBmz56NZcuWYdmyZcjOzu6ivyxmW2p/UvjX\nRY2uNbjO29lmvr/ztrF4Yv0biye1/2wAy248ZsMa+Pj4wN/fH8XFxQCAPXv2IDIy0mqiZr/+WiW5\nT+6jl9CjRe63MQOh+satra2kVCqpT58+tGjRIiIiHW1ujUbTRav7RplIaEiTYD16IfrbptqZGk+f\nf+niWYPCwkKKjY2loUOH0uTJk6muro6qq6spISGBQkJCKDExkWpra7vYyZHvqFFJBLxrxj77hQYO\nDDXok/XopRtfxvIUNvaF5yX6YuylS5cwbtw4ZGZmYsqUKSyB0A1co5fezp7GENfoDVrJfiy5Ri9S\nphgA3N3dMWHCBPzwww82IYEgJ0JuWxZ76zLDMIxYBE30Fy9e1KoJXr16Fbt374ZKpbIJCQQ5a2Nt\nty3nAbLfuqyWINue7J/hGr2kXqX3aGM1ekEyxZWVlUhPT4dGo4FGo0FaWhoSEhKgUqmQkpKC9evX\na9srGYZhGOvS426YkhOutduOnbXGUGtrK2JjY+Hn54cvvvjCaveOcI3ejAhcoxdfo2cYR8IWl89k\nGGMImujLysowZswYREZGIioqCmvWrAFgG4uDy18bk9u/JWLYu3/rYEvLZ3KNXlKv0nvsCTV6FxcX\nvP7664iJiUFDQwOGDx+OxMREbNy4EYmJiVolv6ysLP6Gw/QYbGn5zLq6jm3L6hs/4w1sl918twFR\nM6H56KJvu/v8hMYzZ/lPXT2bzvkUdtrW1aoREr+wsFCC5UXbfy+FaAR34Hfg/vvvp927d1NYWBhV\nVVUREVFlZSWFhYV1ea9EIa0CJLuhiO3E2lmaL774gubPn09EbTfDTJw4kYioy02BHh4eXWzlyHfK\nlFkk9Q1TYhF2POU/lubnZYs5ictL8OLg7ZSWluLIkSOIi4uzicXB5f92oL7x09xtGHld6u2eFK8Q\nQHsZsBTWoH35zB07duDatWu4fPmyzvKZPj4+PfLeEaaHIOavUn19PQ0bNow+//xzIrLet5uOyCmB\nAIBMv71fzDfXPIF2psbT9xlYAsEU1Gq19hv9okWLKCsri4iIMjMzafHixV3eL0e+LIEgJi9j54P4\n42VrEgiCu26am5sxdepUpKWlaW+M6ul3xjJMO+1dN0uWLMHu3bsRGhqKvXv3YsmSJVbOjGG6IqiP\nnoiQnp6Ofv364fXXX9c+//TTT6Nfv35YvHgxsrKyUFdX16OWVeM+etuxs6cxxH30Bq1kP5bcRy+w\n6+abb77Bli1bMHToUKhUKgBAZmYmlixZwnfGMgzD2BqCiz4CkTukKTU8V1ePGzUyIQ9T689iatF5\nAu1MjafvM3CNXg7kyJdr9GLyMnY+iD9etlajF911Y4+0iZORAEtWomQYxv7ocRIIci4MfiOCzP4t\nEcPe/VseW7sb3MvLR3Kf0p87UvuT6/yW3qf885B5CJroMzIy4O3tjejoaO1ztiB/wDBy0X43eFFR\nEQ4ePIg33ngDP//8M2vdMHaBoIl+zpw5yM3N1XnOVgY8a904gn/L4+Pjg5iYGABAnz59MGTIEJw7\nd461bgx7lNgfa90IRVCNftSoUSgtLdV5bvv27cjPzwfQNuDj4+P52w3TI7GFu8FZ64a1bsxBsB59\naWkpkpOTcezYMQCAh4eHdr1YIoKnp6fO+rHagAoF0tPTrSqBMGbMGNy8GKu+8TPehG0F2laYMvX9\n7dscT1y8zhIImyBw2IqmoaEBo0ePxvPPP48HHnhAZ9wDllsnmfvozYjAffTC+3VKSkooKipKu22K\n/MGNPypCQ0oGJGsjZDtr2VmDpqYmSkpKotdff137XFhYGFVWVhIRUUVFhcWE/FjUTM68bDEncXlJ\n1lcJ0ecAAAnUSURBVHVjK/IHXKN3BP+Wh4gwd+5cREREYMGCBdrnrbVOMtfoJfUqvUcbq9FLNtHb\nwsLgDCMX7XeD5+XlQaVSQaVSITc3l7VuGLtAUI0+NTUV+fn5uHjxIry9vfHSSy/h/vvvR0pKCs6e\nPat37UzANrRuWLPG/u2sPYbMgWv0Bq1kP5ZcoxfYdfPhhx92+/yePXsEJcEwDMPIR4+7M5Zr9I7g\nn+EavaRepffYU2v0tkJhYaHxN4mLILN/S8Swd/+Mbh+9NEh/7kg/DuQ5v+0lT+FIPtHn5uYiPDwc\nISEhWLlypdTudXBz84RCodB5LFy4sMtznR/isIS0g9wx7N2/7WHJcQ8ALS3NkvuUXrZE+nEgj7SK\nveQpHEkn+tbWVjz++OPIzc3F8ePH8eGHH+Lnn3+WMoQON1UoOz5e6Oa5zg+GkQ5Lj3uGMRdJZYoL\nCgoQHBysvet1xowZ2LZtG4YMGWLQ7k9/Woj8/IMSZVEqkR9r+bdEDHv3b1sIHfdiaGxskNxnZ1kT\nCTxK7E+OHAH7yVM4kk70586dg7+/v3bbz88P3333XZf3iS+f6Hjr5rlNAu1MYZOJ/sXGa48hNE9j\ndvo+g1TxTN1H9q/xb8q4l3bMtzP7xsM0KiuN59F+L4xwzB8H5u4bYTmadz5IcbyM52m5sS/pRG/K\nzrGn/meGMQVj457HPGNtJK3R+/r6oqzspkpeWVkZ/Pz8pAzBMDYHj3vG1pF0oo+NjcXJkydRWlqK\npqYmfPTRR5g0aZKUIRjG5uBxz9g6kpZunJ2dsW7dOowbNw6tra2YO3eurBekGMYW4HHP2DqS99FH\nRUVh0KBB+M1vfoP333/fpLU1MzMzERISgvDwcOzatcug/2vXriEuLg4xMTGIiIjAM888I6n/dlpb\nW6FSqZCcnCy5/4CAAAwdOhQqlQojRoyQJf+6ujpMmzYNQ4YMQUREBL777jvJYpw4cUIr7KVSqeDu\n7o41a9ZI+hkyMzMRGRmJ6OhozJw5E9evX5d8H0nJ+PHjsXr1ajg7O2PDhg3d9tKr1Wq4u7tr99sr\nr7xi0Gd3S3Z25sknn0RISAiUSiWOHDliNE9jPs3NEdC/nq6YXE3xaU6u+uYNMTma4lPI/gS6zj9i\n8gQgvfByZWUlHTlyhIiI6uvrKTQ0lI4fP06LFi2ilStXEhFRVlYWLV68mIiIioqKSKlUUlNTE5WU\nlFBQUBC1trYajNHY2EhERM3NzRQXF0f79++X1D8R0apVq2jmzJmUnJxMRCSp/4CAAKqurtZ5Tur8\nZ82aRevXr9fup7q6OsljEBG1traSj48PnT17VjL/JSUlFBgYSNeuXSMiopSUFHr33XdlyV8qWlpa\nKCgoiEpKSqipqYmUSiUdP35c5z15eXna8WQK+/bto8OHD+us+9CRr776isaPH09ERAcPHqS4uDjR\nPs3NkUj/OS8mV1N8mptrd/OGmBxN8SlkfxJ1nX/E5in5N3pz19bctm0bUlNT4eLigoCAAAQHB6Og\noMBgjNtuuw0A0NTUhNbWVnh4eEjqv7y8HDt27MC8efO0HRNS+ge6dmJI6f/SpUvYv38/MjIyALSV\nFtzd3SX/DECbkF1wcDD8/f0l8+/m5gYXFxdcuXIFLS0tuHLlCgYNGiRL/lLRsZfexcVF20vfmc7H\n3RCjRo2Ch4eH3tc77o+4uDjU1dXpXcrQVJ/m5gh0f85XVFSIytUUn+bm2nne8PT0FJWjKT7NzRHo\nfv4Rm6esWjemrK1ZUVGh06Hg5+eHc+fOGfSr0WgQExMDb29v7b93UvpfuHAhXnvtNTg53dw9UvpX\nKBQYO3YsYmNj8fbbb0vuv6SkBF5eXpgzZw6GDRuGRx55BI2NjZLGaCcnJwepqamSfgZPT0889dRT\nGDx4MAYNGoS+ffsiMTFRlvylorte+s45KBQKHDhwAEqlEvfddx+OHz8ueczy8nJRPsXm2PGclypX\nfT7NzbXzvBERESE6R2M+hezP7uYfsXnKNtE3NDRg6tSpWL16NVxdXXVeM6Y5Y6wv2cnJCYWFhSgv\nL8e+ffuQl5cnmf8vv/wSAwYMgEql0vuXWGz+33zzDY4cOYKdO3fijTfewP79+yX139LSgsOHD2P+\n/Pk4fPgwbr/99i4LtYuNAbR9i/niiy/w4IMPdmsv1P+pU6eQnZ2N0tJSVFRUoKGhAVu2bJE8fykx\nJd6wYcNQVlaGH3/8EU888YQki/N0HqNiP7eYHBsaGjBt2jSsXr0affr0kSRXQz7NzbXzvNGdwqS5\nORrzaW6Opsw/gvI0+KpAmpubMXXqVKSlpWk/mL6lBjv3IJeXl8PX19ekOO7u7pgwYQJ++OEHyfwf\nOHAA27dvR2BgIFJTU7F3716kpaVJmv/AgQMBAF5eXpg8eTIKCgok9e/n5wc/Pz/ceeedAIBp06bh\n8OHD8PHxkfQY7Ny5E8OHD4eXlxcA6Y7x999/j5EjR6Jfv35wdnbGlClT8O2330qev5SY0kvv6uqq\n/Vd//PjxaG5u7rKQuJiYUnxuoTm2n/MPP/xwt5OZkFyN+RSaa/u88f3334vO0ZhPc3Psbv6ZNWuW\n+DzNvkpgBI1GQ2lpabRgwQKd5xctWkRZWVlERJSZmdnlQtr169fp9OnT9Nvf/pY0Go1e/7/++ivV\n1tYSEdGVK1do1KhRtGfPHsn8d0StVtPEiRMlzb+xsZEuX75MREQNDQ00cuRI+vrrryXPf9SoUXTi\nxAkiInrhhRdo0aJFkseYPn06vfvuu9ptqfwXFhZSZGQkXblyhTQaDc2aNYvWrVsnyzGWiubmZvrt\nb39LJSUldP369W4vxlZVVWnz+u677+iOO+4w6rekpMSki7HffvutSRfljPkUkqO+c15Mrqb4NCdX\nffOGmBxN8Slkf7bTcf4RkycRkeQT/f79+0mhUJBSqaSYmBiKiYmhnTt3UnV1NSUkJFBISAglJiZq\ndxAR0fLlyykoKIjCwsIoNzfXoP+jR4+SSqUipVJJ0dHR9OqrrxIRSea/I2q1WnvVWyr/p0+fJqVS\nSUqlkiIjI2nFihWy5F9YWEixsbE0dOhQmjx5MtXV1Ukao6Ghgfr166f9oyX1Z1i5ciVFRERQVFQU\nzZo1i5qammQ5xlKyY8cOCg0NpaCgIO1x/ec//0n//Oc/iYho3bp1FBkZSUqlku666y769ttvDfqb\nMWMGDRw4kFxcXMjPz4/Wr1+v44+I6LHHHqOgoCAaOnQo/fDDD0ZzNObT3ByJuj/nd+zYISpXU3ya\nk6u+eUNMjqb4FLI/2+k4/4g97oLWjGUYhmHshx63whTDMAyjC0/0DMMwPRye6BmGYXo4PNEzDMP0\ncHiiZxiG6eHwRM8wDNPD+f+4IC6CwPLHdQAAAABJRU5ErkJggg==\n",
"text": [
"<matplotlib.figure.Figure at 0xbbe330c>"
]
}
],
"prompt_number": 5
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"# dummify rank\n",
"dummy_ranks = get_dummies(df['prestige'], prefix='prestige')\n",
"\n",
"# create a clean data frame for the regression\n",
"cols_to_keep = ['admit', 'gre', 'gpa']\n",
"data = df[cols_to_keep].join(dummy_ranks.ix[:, 'prestige_2':])\n",
"\n",
"# manually add the intercept\n",
"data['intercept'] = 1.0"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 6
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"# build logit model\n",
"train_cols = data.columns[1:]\n",
"logit = sm.Logit(data['admit'], data[train_cols])\n",
"\n",
"# fit the model\n",
"result = logit.fit()"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"Optimization terminated successfully.\n",
" Current function value: 0.573147\n",
" Iterations 6\n"
]
}
],
"prompt_number": 7
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"# odds ratios and 95% CI\n",
"params = result.params\n",
"conf = result.conf_int()\n",
"conf['OR'] = params\n",
"conf.columns = ['2.5%', '97.5%', 'OR']"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 8
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"# instead of generating all possible values of GRE and GPA, we're going\n",
"# to use an evenly spaced range of 10 values from the min to the max\n",
"gres = np.linspace(data['gre'].min(), data['gre'].max(), 10)\n",
"gpas = np.linspace(data['gpa'].min(), data['gpa'].max(), 10)\n",
"\n",
"# enumerate all possibilities\n",
"combos = DataFrame(cartesian([gres, gpas, [1, 2, 3, 4], [1.]]))"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 9
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"# recreate the dummy variables\n",
"combos.columns = ['gre', 'gpa', 'prestige', 'intercept']\n",
"dummy_ranks = get_dummies(combos['prestige'], prefix='prestige')\n",
"dummy_ranks.columns = ['prestige_1', 'prestige_2', 'prestige_3', 'prestige_4']"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 10
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"# keep only what we need for making predictions\n",
"cols_to_keep = ['gre', 'gpa', 'prestige', 'intercept']\n",
"combos = combos[cols_to_keep].join(dummy_ranks.ix[:, 'prestige_2':])"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 11
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"# make predictions on the enumerated dataset\n",
"combos['admit_pred'] = result.predict(combos[train_cols])"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 12
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"# plot Prob(admit=1) isolating GRE and Prestige\n",
"isolate_and_plot('gre')"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAEVCAYAAAALsCk2AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzsnXdYFOfa/79LkbrAAiLSRRQpFhQQC00TJfYSiS0nJnpi\n2nviSd4Tk5Nmej/vLyeeJKacFAtYiKIiqCiLYAFrFESa0i1IZ2lbnt8fIwNLUQSWLdyf65qLHWaY\nuWeA+c5dnvsRMMYYCIIgCKIdeuo2gCAIgtA8SBwIgiCITpA4EARBEJ0gcSAIgiA6QeJAEARBdILE\ngSAIgugEiYMWoKenh+vXr3e7vby8HF5eXmhublbJ8e+Hr68vTpw40auf7Sm3b9+Gt7c3Wlpaut3n\n+eefx4cfftin86xZswZvv/12r39eKBSioKCgTzZoO5s2bcKTTz6pVhuKioogFApBVfp9g8RBRbi5\nucHU1BRCoRD29vZ4+umnIZFIVHKuTz/9FE8//TSMjIxUcvz7kZGRgZCQEAB9ezDs2rULU6dOhZmZ\nGcLDw5W2DRs2DOHh4fjhhx+6/fnvvvsOb731Vq/O3YpAIIBAIOjRvmFhYfj555+VvldXVwc3N7c+\n2aDt9PT+qRIXFxfU1dXxtnT1uyIeDImDihAIBDh48CDq6upw4cIFnDt3rss3W5lM1qfzNDc34/ff\nf8fq1av7dBx1Y2Njg1deeQWvv/56l9tXrVqFLVu2qNyOnr5tasJDsCv6+vekSSgUin45jqb+rjQd\nEocBwMHBAREREcjMzATAhXG+/fZbjBo1Cp6engCAH3/8EaNGjYKNjQ0WLlyImzdvKh0jLi4OI0eO\nxNChQ/Haa6/xD7G0tDRYWVnBwcGB3/eXX36Bt7c3LCwsMHLkyE5v3F988QUcHBzg5OSE//73v0rb\n1qxZgxdeeAFz5syBUChEcHAwbt26hZdffhkikQheXl64dOkSv7+bmxuOHTuGhIQEfPLJJ9i5cyeE\nQiH8/Pwe6h7NnDkTjz/+OIYPH97l9sDAQFy/fh3FxcVdbm8fErp79y7mzZsHkUgEGxsbhISE8Pcr\nKysLYWFhEIlE8PX1xYEDB7o8XlVVFebNmwc7OztYW1tj/vz5KC0tBQC8+eabSElJwUsvvQShUIi/\n/e1vAJTDc2vWrMGLL76IefPmwcLCAkFBQUqhuyNHjsDT0xNWVlZ48cUXERoa2u3bbWNjI5566ilY\nW1vD29sbn3/+OZydnfntbm5u+PzzzzFu3DgIhUIoFAqcOXMGU6dOhUgkwoQJE5CcnNztvf/000/h\n4eEBCwsL+Pj4YN++ffy2X3/9FdOnT8c//vEPWFtbw93dHQkJCfz2GzduIDQ0FBYWFpg1axbu3r3b\n7XnEYjGcnJzwySefYOjQoRgxYgR27NjBb1+zZg2ef/55zJkzB+bm5hCLxSgrK8PSpUthZ2cHd3d3\nfPPNN/z+6enp8Pf3h6WlJezt7fHqq68CAAoKCqCnpwe5XN7t7+pB9/+///0vvL29YW1tjYiICBQV\nFXV7XToLI1SCm5sbS0xMZIwxVlRUxHx8fNg777zDGGNMIBCwWbNmsaqqKtbU1MSOHTvGbG1t2cWL\nF1lzczP7n//5HxYSEsIfSyAQsBkzZrCqqipWVFTERo8ezX766SfGGGObN29mc+fOVTp3XFwcu379\nOmOMseTkZGZqasouXLjAGGMsPj6eDRs2jGVmZjKJRMJWrFjBBAIBy8/PZ4wx9tRTTzFbW1t24cIF\n1tTUxGbMmMFcXV3Z1q1bmUKhYG+99RYLDw9Xus5jx44xxhjbtGkTe/LJJ5Vsef7555mVlVWXy/jx\n4zvdtx9//JGFhYV1eU/HjRvH9u/f3+W2NWvWsLfffpsxxtjrr7/OnnvuOSaTyZhMJmOpqamMMcZa\nWlrYyJEj2SeffMKkUik7fvw4EwqFLDs7mz/GW2+9xRhjrKKigv3xxx+ssbGR1dXVsWXLlrFFixbx\n5wsLC2M///yzkg0d76ONjQ07e/Ysk8lkbNWqVWz58uWMMcbKy8uZhYUF27t3L5PL5ezrr79mhoaG\nnY7XysaNG1lYWBirrq5mJSUlbOzYsczZ2Znf7urqyvz8/FhJSQlrampiJSUlzMbGhsXHxzPGGDt6\n9CizsbFh5eXlXR5/9+7d7ObNm4wxxnbu3MnMzMzYrVu3GGOM/fLLL8zQ0JD99NNPTKFQsO+++445\nODjwPxsUFMReffVV1tLSwk6cOMGEQmGnv4FWkpKSmIGBAb9/cnIyMzMz4+//U089xSwtLdmpU6cY\nY4w1NDSwiRMnsg8++IBJpVJ2/fp15u7uzg4fPsyfe9u2bYwxxiQSCTtz5gxjjLEbN24wgUDA5HJ5\nl7+rB93/ffv2MQ8PD3bt2jUml8vZhx9+yKZOndrlNekyJA4qwtXVlZmbmzMrKyvm6urKXnzxRdbU\n1MQY4x4iSUlJ/L7PPPMM27hxI79eX1/PDA0NWWFhIb9/6z8EY4x9++23bObMmYwxxj788EP+odMd\nixYtYl9//TVjjLGnn36avfHGG/y2nJwcpYfamjVr2LPPPstv/+abb5i3tze/fvnyZWZlZcWvtxeH\nd999l61evboHd6d77icO06ZNY1u3bu1yW3txeOedd9jChQtZXl6e0j4nTpxg9vb2St9bsWIF27Rp\nE3+MVnHoyMWLF5lIJOLXw8LCeIFupeN9/Otf/8pvO3ToEBszZgxjjLHffvut08PG2dm5W3Fwd3dn\nR44c4dd/+ukn5uTkxK+7ubmxX375hV//9NNPOz2gZ8+ezX777bcuj9+RCRMmsNjYWMYYJw4eHh78\nNolEwgQCAbt9+zYrLCxkBgYGrKGhgd++cuXKbv8GWsWh/f6RkZHsgw8+YIxx4vDUU0/x286cOcNc\nXFyUjvHxxx+zp59+mjHGWEhICHv33Xc7iV5X4tD+d/Wg+x8REaH0u5DL5czU1JQVFRV1eV26CoWV\nVIRAIEBsbCyqqqpQUFCAzZs3KyWM24cFbt68CVdXV37dzMwMNjY2fBij4/4uLi4oKysDAIhEItTV\n1SmdOz4+HkFBQbCxsYFIJMKhQ4dQUVHBn6vjsTpiZ2fHfzY2NlZaNzExQX19fc9vRD9SV1cHKyur\nbreze6Gjf/zjH/Dw8MCsWbMwcuRIfPbZZwCAsrIypWsHAFdXV/5etqehoQHr16+Hm5sbLC0tERoa\nipqaGqWcxINi2cOGDeM/t79vZWVlcHJyUtq343p7Otrd1b7ttxcWFmL37t0QiUT8cvLkSdy6davL\n4//+++/w8/Pj983IyOD/XgDA3t6e/2xqagoAqK+vR1lZGUQiEUxMTPjt7f+Ou6Kr/VtDqAKBQOna\nCgsL+XO0Lp988gnu3LkDAPj555+Rk5MDLy8vBAYGIi4urtvztv9dPej+FxYW8mHU1tAkAKX/x8EA\niYOaaP/H6uDgoFQCKZFIUFFRAUdHR/577WOeRUVF/LZx48YhJyeH39bc3IylS5fitddew507d1BV\nVYU5c+bwD7Xhw4d3OpYqrqmV5557DkKhsMtl7NixPToGwCVa8/LyMH78+AfaYW5uji+//BL5+fnY\nv38//vWvf+H48eNwdHREcXGx0gO+sLBQ6T63nv+rr75CTk4O0tPTUVNTg+TkZDDO076vnT3BwcEB\nJSUl/DpjTGm9I8OHD1fKtXSVd2lvj4uLC5588klUVVXxS11dHV577bVOP1dYWIhnn30W//nPf1BZ\nWYmqqir4+vr2KDE/fPhwVFVVoaGhQel497s3Xe3fPl/W8TpGjBihdB21tbU4ePAgAMDDwwM7duxA\neXk5Nm7ciMcffxyNjY33vTfAg++/i4sLfvjhB6XzSiQSBAUFPfCe6BIkDhrAihUr8Msvv+DPP/9E\nc3Mz/vnPfyIoKEjprf7LL79EdXU1iouL8e9//xtPPPEEAC5RW11dzb/9trS0oKWlBba2ttDT00N8\nfDyOHDnCHycyMhK//vorsrKy0NDQgPfee0/Jlp48FLrD3t4eBQUFSsf4/vvvUVdX1+Vy5coVfj+F\nQoGmpiZIpVIoFAo0NzdDKpXy29PT0+Hm5tbpzb8ruw8ePIi8vDwwxmBhYQF9fX3o6+tj8uTJMDU1\nxeeffw6pVAqxWIyDBw9i+fLl/DFaj1NfXw8TExNYWlqisrKy030aNmwY8vPzu70X97uPc+bMwZUr\nVxAbGwuZTIb//Oc/3b7VA9zv7JNPPkF1dTVKS0uxefPm+z6AV69ejQMHDuDIkSOQy+VoamqCWCzu\n8s1XIpFAIBDA1tYWCoUCv/zyCzIyMro9dntcXV3h7++Pd999F1KpFKmpqfyD+3607p+SkoK4uDgs\nW7YMQOd7FhgYCKFQiM8//xyNjY2Qy+XIyMjAuXPnAADbtm1DeXk5AMDS0hICgQB6ep0faR1/V3Pn\nzr3v/X/uuefw8ccf4+rVqwCAmpoa7N69u0f3RJcgcVADHf+xZ86ciQ8++ABLly6Fg4MDbty4gejo\naKV9Fi5ciEmTJsHPzw/z5s3DM888AwAYMmQI1qxZg23btgHgBmL9+9//RmRkJKytrREVFYWFCxfy\nx4mIiMCGDRswY8YMjB49GjNnzlSyp2Otf1e1/909mFr/yW1sbODv7/9Q9+T333+HqakpXnjhBaSk\npMDExATr16/nt2/fvh3PP/98tz/f3s68vDw8+uijEAqFmDp1Kl+NYmhoiAMHDiA+Ph5Dhw7FSy+9\nhK1bt2L06NGdjrFhwwY0NjbC1tYWU6dOxWOPPaZ03S+//DL27NkDa2trbNiw4b72tP8eANja2mL3\n7t147bXXYGtri6ysLPj7+3c7TuWdd96Bk5MTRowYgVmzZmHZsmUYMmRIt/fCyckJsbGx+Pjjj2Fn\nZwcXFxd89dVXXZaGent749VXX8WUKVNgb2+PjIwMTJ8+vUfXAQA7duxAWloarK2t8f777+Opp57q\n1i6Ae4EQiURwcHDAk08+iS1btnR5/wGu+uvgwYO4dOkS3N3dMXToUDz77LOora0FABw+fBi+vr4Q\nCoX4+9//jujoaP4e3u93ZWNjc9/7v2jRImzcuBHLly+HpaUlxo4di8OHD9/3unQRAevLqyKhEdy9\nexfBwcG4dOmSWgbCqZo7d+4gLCwMly5duu9DUVtRKBRwdnbGjh07EBoa+sD9v/vuO+zatQtJSUkD\nYF3/IRaL8eSTT3ZbjqwuHvb+DxZU6jkkJCRgzJgxGDVqFJ8U7IhYLIafnx98fX0RFhamSnN0lta3\nH10UBoBLkF+9elWnhOHIkSOorq5Gc3MzPv74YwDoNqZ969YtnDx5EgqFAtnZ2fjXv/6FxYsXD6S5\nOsfD3P/BioGqDiyXy/HSSy8hMTERjo6OCAgIwIIFC+Dl5cXvU11djRdffBGHDx+Gk5PTfQfQEIQu\ncfr0aaxcuRItLS38wLPuxL2lpQXPPfccbty4ASsrK6xYsQIvvPDCAFvcP2jKaOWHuf+DFZWFlU6f\nPo333nuPH0356aefAoBSe4Rvv/0Wt27dwvvvv68KEwiCIIheojLPobS0tFNtdlpamtI+ubm5kEql\nCA8PR11dHV5++eVOjds05U2DIAhC2+jLu7/Kcg49eahLpVJcuHABhw4dwuHDh/HBBx8gNze3036t\nJYa6uLz77rtqt4Guj66Nrk/3lr6iMs+hdcBRK8XFxZ1GJTo7O8PW1hYmJiYwMTFBSEgI/vzzT4wa\nNUpVZhEEQRA9QGWeg7+/P3Jzc1FQUICWlhbs3LkTCxYsUNpn4cKFSE1NhVwuR0NDA9LS0uDt7a0q\nkwiCIIgeojLPwcDAAJs3b8bs2bMhl8uxdu1aeHl58T35169fjzFjxiAiIgLjxo2Dnp4e/vrXvw46\ncdD18l1dvj5dvjaArm+wo/GD4AQCQb/EzwiCIAYTfX12qsxzUDXW1taoqqpStxl9QiQSobKyUt1m\nEARBdEJrxaGqqkrrPQoq0yUIQhX0R7NlarxHEAShA1RUAN9/D4SEABMn9v14JA4EQRBaSkMDEB0N\nLFgAjBwJiMXA//4v0MX8VQ+N1oaVCIIgBiMyGZCYCOzYARw4AEyeDKxcCWzfDgiF/XcenfMc9PX1\n4efnh7FjxyIyMrLLmaEeltauja1Mmzatz8ckCILoKYwBZ84Af/sb4OgIbNoE+PsD164BCQnAX/7S\nv8IAaHEpa3ffFwqF/JzKq1evxqRJk/D3v/+d3y6TyWBg8HAOU/tj9idUpksQxP24do3zEHbsAPT1\ngVWrOC/Bw+PBP9vX54vOeQ7tCQ4ORl5eHpKTkxEcHIyFCxfC19cXCoUC//jHPxAYGIjx48fjhx9+\nAADcvHkTISEhvOeRmpqK119/HY2NjfDz8+ObApqbmwPgJgl54YUX4OXlhVmzZmHu3LmIiYkBAJw/\nfx5hYWHw9/dHRETEfaeBJAiCaKW0FPjXv4BJk4AZMwCJBNi5kxOKd97pmTD0C0zD6c7E7r5vbm7O\nGGNMKpWyBQsWsO+//56JxWJmZmbGCgoKGGOMbdmyhX344YeMMcaampqYv78/u3HjBvvqq6/YRx99\nxBhjTC6Xs7q6OqVjdjzH7t272Zw5cxhjjN26dYuJRCIWExPDWlpa2JQpU9jdu3cZY4xFR0ezZ555\npsfXQBDE4KKqirGffmJsxgzGRCLGnnmGscRExmSy3h+zr88XnUtIt77lA0BISAieeeYZnDx5EoGB\ngXB1dQXAzQJ15coV7NmzBwBQW1uLvLw8BAQE4JlnnoFUKsWiRYswfvz4+54rNTUVkZGRALhJzMPD\nwwEA2dnZyMzMxCOPPAKAm/jIwcFBJddLEIR20tQEHDrEJZITE4GZM4EXXgDmzgWMjdVtnQ5WK5mY\nmODixYudvm9mZqa0vnnzZjz66KOd9ktJScHBgwexZs0avPLKK53ml2jP/WJ6Pj4+OHXq1ENaTxCE\nLiOXA8nJnCDs3QtMmMDlEX7+GbCyUrd1yuh0zqE7Zs+ejW+//RYymQwAkJOTg4aGBhQVFWHo0KFY\nt24d1q5dy4uMoaEhv297pk2bhpiYGDDGcPv2bYjFYgCAp6cnysvLcebMGQDcvBVXr14dmIsjCEKj\nYAy4cAF49VXAxYUbh+DtDVy5Ahw/Dqxdq3nCAOig59BVSwqBQKD0/XXr1qGgoAATJ04EYwx2dnbY\nu3cvxGIxvvjiCxgaGkIoFOL3338HADz77LMYN24cJk2ahK1bt/LHWrp0KY4dOwZvb284Oztj4sSJ\nsLS0hKGhIfbs2YO//e1vqKmpgUwmw9///vdB13GWIAYz+flcldH27YBUylUZJSYCXl7qtqxn6Fwp\n60AjkUhgZmaGiooKTJ48GadOnYKdnV2PflZTroEgiP7h9m1g1y5OEG7cACIjubDR5MnAQLdSG7Rd\nWTWFefPmobq6Gi0tLXjnnXd6LAwEQegG9fXAvn3Atm3cQLX587lBao88AjzkkCqNgjwHNaIL10AQ\ngxGZDDh6lPMQDh4Epk/nPIQFC4AOtS9qo6/PFxIHNaIL10AQgwXGgLNnOUGIjgbc3TlBeOIJYOhQ\ndVvXGQorEQRBqJD8fE4Qtm3jBGL1auDkyQEcqawmSBwIgiA6UF7OJZa3bQOuX+e8g23bgICAgU8s\nqwsKK6kRXbgGgtAVGhqA/fs5LyElhRupvHo1l1g2NFS3dQ8PhZUIgiB6iVwOJCVxXkFsLBAYyAnC\njh393wJb2yDPQY3owjUQhLbBGHDpEicIUVGAgwMnCMuXA/b26rau/6CW3RrI5s2b4e/vD2NjYzz9\n9NPqNocgCAAFBcDHHwM+PsCSJYCpKde+4tw5YMMG3RKG/oDCSirA0dERb7/9Ng4fPtwvM9ERBNE7\nKiuBPXs4L+HqVWDZMuDHH4GpUwdPYrm3kDiogMWLFwMAzp07h5KSEjVbQxCDi6YmIC6OE4Tjx4HZ\ns7lmdxERwJAh6rZOeyBxUCGUTyCIgUGhAE6c4AThjz8APz8uj/Drr4Clpbqt0050Vhz6y2Xsy/O9\nqw6xBEH0H1eucKWn27cDNjacIFy+DDg5qdsy7UdnxUETXtrJcyCI/qekhKsy2raNyymsWsXNqDZ2\nrLot0y1UWq2UkJCAMWPGYNSoUfjss886bReLxbC0tISfnx/8/Pzw4YcfqtKcAYc8B4LoH2prgV9+\n4abSHDcOyMkBvv4aKCwEPv2UhEEVqMxzkMvleOmll5CYmAhHR0cEBARgwYIF8Oow00VoaCj279+v\nKjPUglwuh1QqhUwmg1wuR3NzMwwMDKCvr69u0whCa2hpAQ4f5jyEhARgxgzNmmNZ11GZ55Ceng4P\nDw+4ubnB0NAQy5cvR2xsbKf9dDH08sEHH8DU1BSfffYZtm3bBhMTE3z00UfqNosgNB7GgNOngRdf\nBBwdgc8/50Thxg1uzuWlS0kYBgqVeQ6lpaVwdnbm152cnJCWlqa0j0AgwKlTpzB+/Hg4Ojriyy+/\n7HIqzU2bNvGfw8LCEBYWpiqz+4VNmzYp2UwQxP3JzeU8hO3buQlyVq8G0tK4tthEzxCLxfw89v2B\nysShJ/H2iRMnori4GKampoiPj8eiRYuQk5PTaT960BKE7nHnDrBzJycKhYVc+4roaGDSJBqg1hs6\nvji/9957fTqeysJKjo6OKC4u5teLi4vh1KG+TCgUwtTUFADw2GOPQSqVorKyUlUmEQShZhoauEqj\nuXOB0aOB9HTg/fe5CqT/9/8Af38SBk1BZZ6Dv78/cnNzUVBQAAcHB+zcuRNRUVFK+9y+fRt2dnYQ\nCARIT08HYwzW1taqMokgCDXQ2vl061au82lQEBc22rkTMDdXt3VEd6hMHAwMDLB582bMnj0bcrkc\na9euhZeXF7Zs2QIAWL9+Pfbs2YPvvvsOBgYGMDU1RXR0tKrMIQhiAGEM+PNPLmS0Y0db59PPPqMG\nd9oCtexWI7pwDQTRnqIiTgy2bQPq6zlBWLUK6FDBTgwANNkPQRBqpbq6rfNpRgbw+OPA999znU/1\naFIArYU8BzWiC9dADE6am4H4eE4Qjh4FHn2U8xIeewwwMlK3dQTQ9+cLiYMa0YVrIAYPCgVw6hQn\nCHv2AL6+nCA8/jhgZaVu64iOUFiJIAiVcvUqNzhtxw5u9rQnnwQuXABcXNRtGaFKKCKoAlpaWrB2\n7Vq4ubnBwsICfn5+SEhIULdZBNFjysqAr74CJk7kQkYtLVz7iowM4PXXSRgGA+Q5qACZTAYXFxec\nOHECLi4uiIuLQ2RkJK5cuQJXV1d1m0cQXVJby02Us20b5xksWgR8+SUQGgpQz8jBB+UcBojx48dj\n06ZN/BSigPZdA6F7tLRwHU+3b+e+hodzpafz5gEmJuq2jugLlHPQAm7fvo2cnBz4+Pio2xSC4BPL\n27cDu3cD3t6cIHz3HUANCohWdFYcBO/1T4MW9m7f3uylUilWrVqFNWvWYPTo0f1iE0H0ho6J5dWr\ngXPnADc3dVtGaCIUVlIhCoUCK1euRH19PWJjYztN9qMN10BoN2VlXKO77duB27eBFSs4URg/nhrc\n6ToUVtJQGGNYu3YtysvLcejQIZoFjhgw2ieWz58HFi8GvvgCCAujxDLRc0gcVMTzzz+Pa9euITEx\nEUY0ZJRQMR0Ty2FhwPr1lFgmeg+FlVRAYWEhRowYAWNjYyWP4YcffsCKFSv4dU2+BkLz6ZhY9vJq\nG7FsY6Nu6wh1Q2ElDcTV1RUKhULdZhA6SvvEsokJJZYJ1UDiQBBaQPvE8q1bwMqVXF5hwgRKLBOq\ngcJKakQXroFQHR0Ty4sWcV4CJZaJnkBhJYLQIZqauFbYO3YAR45wrSuefRaYP58Sy8TAQp6DGtGF\nayD6jlwOiMWcIOzdC4wbx41YXrqURiwTvYfmc9BidOEaiN7BGJdE3rED2LmTm1d55UrgiScAZ2d1\nW0foAhRWIggtIjubSyzv2MGVoq5cCRw7RnMsE5oHiQNBqJjSUs472LGD+/zEE1ySOSCAKo0IzYXC\nSmpEF66B6JqqKiAmhhOES5e4SqOVK7mW2FRpRAwElHPQYnThGog2GhqAgwc5QUhK4mZQW7kSmDMH\nMDZWt3XEYKOvzxeaJlRFrF69GsOHD4eFhQXc3d3x0UcfqdskQgXIZFwvo7/8BXBwAH76ifMSioqA\nPXuAJUtIGAjthDwHFZGZmYmRI0fC2NgY2dnZCA0Nxa+//oqIiAh+H02/BqJrGAPOnOE8hF27gBEj\nOA8hMpKrOiIIdVHRUIE9V/cgOjMa4jViqlbSRDrO+mZgYAA7Ozs1WUP0B5mZnCDs2MF5A6tWcY3v\nRo5Ut2XEYKauuQ6x2bGIyohCalEqIjwi8PLklyGGuE/HJXFQIS+88AJ+++03NDc3Y/PmzZg4caK6\nTSIeksJCIDqaE4TKSm6ynL17abIcQr00yZpwKPcQojKicCT/CEJcQ7Bq7CrsfHwnzIeY98s5dDas\nJBb3z39uWFjfbg9jDMnJyXj88cdx6NAhBAYG8tsorKSZ3L3LtcDesQPIyuJaYK9cCUyfDuhRlo5Q\nE1K5FMduHENURhT2Z+/HxOETsdxnOZZ6L4W1Seeh9BpdrZSQkIANGzZALpdj3bp12LhxY5f7nT17\nFlOmTMGuXbuwZMkSZQO1NOfQkeeffx7Gxsb4v//7P/572nYNukxdHRAbyw1QS03lKoxWrgRmzwaG\nDFG3dcRgRcEUSC1KRXRGNPZc3QN3kTtW+K5ApE8khguH3/dnNXaEtFwux0svvYTExEQ4OjoiICAA\nCxYsgFeHoaByuRwbN25ERESETj8opVIpbGgGFo2isRGIi+PCRkePAiEhXNho507AvH88c4J4aBhj\nuHDzAqIyorAzcydExiKs8F2BM+vOwF3kPmB2qEwc0tPT4eHhAbd7M5AsX74csbGxncThm2++weOP\nP46zZ8+qypQBp7y8HMeOHcP8+fNhbGyMxMRE7N69G4mJieo2bdDT0sJ1O42O5sYkBAQAy5cDP/xA\nTe4I9ZJVnoWojChEZ0RDzuRY4bsCCasS4GPn8+AfVgEqE4fS0lI4t+sg5uTkhLS0tE77xMbG4vjx\n4zh79ixe1Am0AAAgAElEQVQE3WT4Nm3axH8OCwtDWFiYKkzuNwQCAb7//ns8//zzYIxh9OjR2Lp1\nKwICAtRt2qBEJuMGpUVHA/v2AT4+XAuLr74Chg1Tt3XEYKagugDRGdGIzohGeUM5nvB5AtuXbIe/\ng3+3z8PuEIvFEIvF/WabysShJxe2YcMGfPrpp3xsrLuwUntx0AZsbW379ZdEPDwKBXDyJCcIe/YA\nrq6ch7BpE3U9JdTL7frb2JW5C1EZUcitzMVSr6X4OuJrTHeZDn293vdW6fji/N577/XJTpWJg6Oj\nI4qLi/n14uJiODk5Ke1z/vx5LF++HABw9+5dxMfHw9DQEAsWLFCVWYQO09oGOzqayxtYW3OCQGMR\nCHVT1ViFP7L+QFRGFM6VncN8z/l4K+QtPOr+KAz1DdVtXpeorFpJJpPB09MTx44dg4ODAwIDAxEV\nFdUp59DK008/jfnz5+tstVJX6MI1qBvGgCtX2gRBX59LKj/xBODtrW7riMGMpEWCAzkHEJURBXGB\nGI+4P4IVviswd9RcmBiqflo/ja1WMjAwwObNmzF79mzI5XKsXbsWXl5e2LJlCwBg/fr1qjo1MQjI\nyeEEIToakEg4D2HPHmDCBBqcRqiPZlkzDucfRlRGFOJz4zHFeQpW+K7A1sVbYWFkoW7zHgqdHQSn\nDejCNQwkBQVcL6PoaODWLWDZMk4UgoJIEAj1IZVLkVSQhJ2ZO7Hv2j742vliuc9yPO79OIaaDVWb\nXRo9CK4/IHEY3JSVcaOVo6OBvDxuXuXly4HgYJoXgVAfMoUM4gIxdmXuwt5rezFSNBKRPpGI9ImE\nk4XTgw8wAJA4aDG6cA2q4O5dbqKc6GhuopyFCzlBmDkTMNTM3B0xCJAr5EguTMauzF34I+sPuFm5\nIdInEo97Pw43Kzd1m9cJEgctRheuob+oqeHGIERHA6dPA489xiWVIyJoPgRCfcgVcqQWpWJn5k7E\nZMXA2cIZkT6RWOa9DCNEI9Rt3n0hcdBidOEa+oJEAhw4wAlCUhIwYwbnIcybB5iZqds6YrCiYAqc\nLDqJXVd3Yc/VPbA3t0ekdySW+SyDh7WHus3rMSQOWowuXMPDUl/P9TPavZvrZzR1KicIixYBlpbq\nto4YrCiYAmdKzmBn5k7suboHtqa2vCCMthmtbvN6hcrFITMzEydOnEBBQQEEAgHc3NwQHBzcaTIb\nVaHt4pCbm4uxY8di2bJl2Lp1q9I2bbmGvlJfz/Ux2r0bSEzkBGHZMk4QqJ8RoS4YY0grTcOuzF3Y\nfXU3LI0s+aTyGNsx6javz6hMHLZu3YpvvvkGNjY2CAwMhIODAxhjuHnzJtLT03H37l28/PLLWL16\nda9P3iMDtVwcZs2ahaamJri5ueH3339X2qYt19Ab6uraBOHYMWDaNE4QFi4kQSDUB2MMZ8vO8oJg\namiKJ3yewDLvZWprcKcqVDYIrqqqCseOHYNQKOxye21tLX799dden3gwEB0dDZFIBG9vb+Tl5anb\nHJVTW8vlEHbv5nII06dzgvDzz4BIpG7riMFKawvsXVd3YVfmLgzRH4InfJ5A3Mo4+Az1eegGd4MF\nyjmoiNraWgQEBCApKQk//PAD8vPzdTKsVFsL7N/fJgghIZwgLFhAgkCoD8YYLt26xAuCnkCPCxl5\nR2LcsHGDQhDU0j7j/fffxzvvvNPrkw4I/fXL7+XNffvtt7Fu3To4ODjo3B9iTU2bIIjFQGgoJwi/\n/QZYWanbOmKwwhjD5duXeUFQMAUifSKxZ9keTLCfoHP/h6qmV+Lw448/ar44qPGN/NKlSzh27Bgu\nXrx4zxTt9g4AoLq6TRCSk4GwME4Qtm6lKiNCfTDGkFmeiV2ZnCA0y5sR6ROJ6KXRmDh8IglCH+hW\nHLrLNQBAY2OjSozRFZKTk1FQUAAXFxcAQH19PeRyObKysnDu3Dk1W9dzqqu5eZV37wZOnADCw7mB\nadu2kSAQ6oMxhit3riAmKwa7M3ejvqUekT6R+H3x7whwCCBB6Ce6zTm4uLggPT0d9vb2nbY5Ozsr\nzdWgSrQx59DY2Ii6ujoA3B/yl19+iYKCAnz//fdK80hr4jVUVbUJQkoKNzBt2TJg/nzAQruaShI6\nBGMM58rOISYrBjFZMZDKpVjitQSRPpEIdAyEnkBP3SZqHCrLOTz55JMoKirqUhxWrFjR6xMOBkxM\nTGBi0tav3dzcHCYmJkrCoElUVrYJQmoq18No1SogKooEgVAfcoUcp4pPISYrBn9k/QETQxMs9VpK\nIaMBgqqV1Ig6r6GykutltHs3N53mI49wHsK8ecB9IooEoVKkcimSC5MRkxWDfdf2wc7MDku9lmKJ\n1xIqO31IBqR9Rm1tLSwsLFBTUwPLAQ42kzj0HxUVbYJw+nSbIMydS4JAqI9mWTOOXj+KmKwYHMg+\nAHeROy8Io2xGqds8rWVAxGHixIm4cOEC/Pz8+AqcgYLEoW+UlHCC8McfwPnzwKOPtgmCublKT00Q\n3SJpkSAhLwExWTGIz4vHWLuxWOq1FIu9FsPF0kXd5ukEAzLOQdsfwoONnBxODPbu5SbImTcP+Nvf\ngFmzAFNTdVtHDFZqmmpwMOcg/rj2BxKvJyLQMRBLvZbiX7P/BXvzzrlNQr2obA5pVSMSibQ+/ijq\npyHEjAEXL3Ji8McfXMXRokXARx9xA9RoghxCXVQ0VCA2OxYxWTFIKUxBqFsolnotxQ/zfoCNqWYW\naBAcWisOlZWV6jZBrcjlXCJ5715uMTAAliwBfvoJmDwZ0KPKPkJN3Ky7iX3X9iEmKwZny87iUfdH\nsXrsakQtjYKFEZW/aQtaKw6DkeZm4PhxzjvYvx8YPhxYvJhrdufr238dQwjiYSmsLsQfWX8gJisG\nmeWZmDtqLl4MeBGzPWbD1JBimdoIiYOGU18PJCRwghAfD3h7cx7CG28A7u7qto4YzORW5PKD0gqq\nC7DAcwH+GfxPzBwxE0YGRuo2j+gjPapWysnJwejRo5GdnQ1PT8+BsItHF6qSHpaKCs4b2LuX63Q6\nZQrnISxcyHkLBKEOGGPIuJPBC0JFQwUWey3GkjFLEOoWCgM9etfUJFReyvr111/j5ZdffuD3VMVg\nEYfS0raS03PnuFHKS5ZwJafU+ppQF62jlGOzYxGbHYsWeQuWei3FUq+lmOI8hdpWaDAqF4euxjZM\nmDABly5d6vVJHwZdFoecnLaEck4OV3K6eDEwezaVnBLqo0HagCP5RxCbHYu4nDg4CB2wcMxCLPRc\nCD97P62vEhwsqEwcoqKisGPHDqSkpCA4OJj/fl1dHfT19XHs2LFen/ShDNQhcWAMuHSpbQxCRQVX\ncrpkCdcCm0pOCXVxR3IHB7IPIDY7FuICMQIcA7DQcyEWeC6Am5Wbus0jeoHKBsFNnToVw4cPR3l5\nOf73f/+XP4lQKMT48eN7fcLBhlwOnDrV5iHo6XFi8MMPQFAQlZwS6iP7bjYfLsq8k4lZI2fhCZ8n\n8Nui3yAyoVjmYEdrG+9pMg0NQGIil1Tevx+wt+fCRUuWAGPHUskpoR4UTIG0kjTsy96H2GuxqGup\nwwLPBVjouRDhbuFUYaRjqCysNG3aNJw8eRLm5uadYowCgQC1tbUPPHhCQgI2bNgAuVyOdevWYePG\njUrbY2Nj8c4770BPTw96enr44osvMGPGjE7n0gZxuHkTOHiQE4PkZGDSJG4e5fnzAQ8PdVtHDFYa\npY1IvJ6I2OxYHMw5iKFmQ7HQk8sfTHKYRAllHWZAGu/1BrlcDk9PTyQmJsLR0REBAQGIioqCl5cX\nv49EIoGZmRkA4MqVK1i8eDHy8vKUDdRQcWAMuHyZE4MDB4DcXCAighODxx6jCiNCfdxtuIu4nDjs\ny96H4zeOw8/ej88fjLQeqW7ziAFiQBrvVVVVobi4GDKZjP/exIkT7/sz6enp8PDwgJubGwBg+fLl\niI2NVRKHVmEAuKk0bW1tH8b2Aae5GRCL28JFhoacGHz6KRAcTAllQn3kV+bz+YNLty7hEfdHsHjM\nYvw0/yfqYUT0igeKw9tvv41ff/0V7u7u0GuXPU1KSrrvz5WWlsLZ2Zlfd3JyQlpaWqf99u3bhzfe\neAM3b97EkSNHujzWpk2b+M9hYWEICwt7kNn9xt27QFwcJwhHj3JtKubP50Yte3lR/oBQDwqmwLmy\nc5wgXIvF3Ya7WOC5AK9NfQ0z3WfC2MBY3SYSA4xYLIZYLO634z0wrDR69GhkZGRgyJAhD3XgmJgY\nJCQk4McffwQAbNu2DWlpafjmm2+63D8lJQXr1q1Ddna2soEDHFZiDMjO5jyD/fuBK1e4AWkLFgBz\n5gB2dgNmCkEo0SxrxvEbxxGbHYv92fthaWyJhZ4LsWjMIppHmeiEysNKPj4+qKqqwrBhwx7qwI6O\njiguLubXi4uL4eTk1O3+wcHBkMlkqKioGPC5lmUybu7k1nBRYyMnBm++CYSHA8b0EkaoiarGKsTl\nxiE2OxZH84/C184XCz0XQrxGjNE2o9VtHqHDPFAc/vnPf8LPzw++vr4wMuJK3QQCAfbv33/fn/P3\n90dubi4KCgrg4OCAnTt3IioqSmmf/Px8uLu7QyAQ4MKFCwAwYMJQXc2Fhg4c4Braubtz4aJdu4AJ\nEyhcRKgHxhiy7mbhUO4hxOXG4XzZeYSPCMdCz4X4z5z/wM6MXFdiYHigOPzlL3/B66+/Dl9fXz7n\n0JPh8wYGBti8eTNmz54NuVyOtWvXwsvLC1u2bAEArF+/HjExMfj9999haGgIc3NzREdH9/Fy7s/1\n623ewdmzQEgIJwiffw44Oqr01ATRLY3SRiQVJPGCIFfIMXf0XLw65VXMGDGDWl4TauGBOYeAgACc\nPXt2oOzpRF/iZnI5kJ7eVm5aXs71L1qwAHjkEaBdsRRBDCiF1YW8GJwoPIEJ9hMwd9RczB09Fz5D\nfah/EdFnVD7O4ZVXXoGRkREWLFjAh5WAB5ey9hcPe4H19VxV0YED3KC0YcPaBqMFBlK7CkI9SOVS\nnC45jbjcOMTlxOG25DYe83gMc0fNxayRs6hdBdHvqFwcwsLCunyLeVApa3/xMBe4Zw/wzDPcNJnz\n53PLiBEqNpAguqFcUo74vHjE5cbhaP5RjBCNwNxRczFn1BwEOARAX09f3SYSOozGjpDuLx7mAqur\nuUSypaWKjSKILlAwBS7evMh5B7lxyL6bjZnuMzHHYw4eG/UYHIQO6jaRGESorJT1q6++4k/QFa+8\n8kqvT6oqrKzUbQEx2KhtrsXR/KOIy41DfF48LI0sMXf0XHw842MEuwZjiP7DjQ8iiL4gk9WgquoY\nKisT+nysbsWhrq4OAoEA2dnZOHv2LBYsWADGGA4ePIjAwMA+n5ggtBHGGLIrshGXE4dDeYdwtvQs\npjpPxZxRc/Bm8JvUu4gYUBhToL7+EiorE1BZmYD6+ouwtJwGa+uIPh/7gWGl4OBgHDp0CEKhEAAn\nGnPmzEFKSkqfT94jAzW08R4xeGiSNUFcIEZcbhwO5R5Ci7wFc0bNwdxRczFzxEyYDaGyN2LgkErv\norLyKCorE1BVdRj6+pawsXkM1tYRsLQMgb4+V/qs8hHSd+7cgWG7jnKGhoa4c+dOr09IENpAcU0x\nLwbiAjHG24/HHI852PvEXoy1G0ulpsSAwZgcdXVnee9AIsmClVUYrK0j4Ob2LkxM3FVy3h4NggsM\nDMSSJUvAGMO+ffvw1FNPqcQYglAXLfIWnCk5w1UX5cShrK4MER4RWOG7Ar8u+hXWJtbqNpEYRLS0\n3EJl5WFUVsajsvIojIwcYW0dgREjPoal5TTo6al+YqYeVSudP38eKSkpEAgECAkJgZ+fn8oNa4XC\nSoQqYIwhrzIPR/KP4Mj1IxAXiDHKehQiPCIwd9RcBDoGUqkpMWAoFFLU1p7ivYOmpgKIRI/A2joC\n1tazYWTUfV+67hiwUtbbt2+jqamJd6ddXFx6fdKHgcSB6C+qGqtw/MZxHLl+BEfyj0Aql2LWyFmY\nNXIWZo6YiaFmQ9VtIjGIaGoq4sWguvo4TEw8YG3N5Q4sLCZDIOjRdDvdonJx2L9/P1599VWUlZXB\nzs4OhYWF8PLyQmZmZq9P+lAGkjgQvUSmkCG9NJ3zDvKPIONOBqa7TOcFwcvWi3IHxIChUDShujqF\nFwSp9A6srWfD2joCItEsDBnSv00VVS4O48aNw/Hjx/Hoo4/i4sWLSEpKwtatW/Hf//631yd9KANJ\nHIiH4EbVDRzOP4wj+UeQVJAENys3TgzcZ2GayzSaBIcYUBob81BZmYCKinjU1KTAzGzsvVBRBITC\nSRCocA4OlVcrGRoawtbWFgqFAnK5HOHh4Xj55Zd7fUKC6E9qm2uRdCOJDxXVNddh1shZWOq1FN/N\n/Q7DzB9uHhKC6AtyuQTV1Um8dyCXN8DaOgL29mvg5bUNhoba00PrgeIgEolQV1eH4OBgrFq1CnZ2\ndjA3Nx8I2wiiE3KFHOdvnseR/CM4nH8Yl25dwhSnKZg1chZiImOozJQYULgy0wuorj6GqqpE1Nam\nQSgMgLV1BHx8/oCZmfb+PT4wrCSRSGBsbAyFQoHt27ejtrYWq1atGrBJeSisRBTVFPF5g2M3jsFB\n6IBZI2dh9sjZCHYJhomhibpNJAYJjDE0NuagqooTg+pqMYyMHGBlNRMi0UxYWYXDwECobjMBqDDn\nwBh7oOL1ZJ++QuIw+KhvqUdyQTKfO6hsrMSjIx/FLPdZeHTko9TAjhhQmpvLUFV17J53cAyAACLR\nTIhEj8DKagaMjIar28QuUZk4hIaGYt68eVi4cCFGj1aeqzY7Oxv79u1DXFwcTpw40euT98hAEged\np7WbaeuYg3Nl5xDgEMB7B+Ptx0NPhYk7gmiPTFaD6mox7x20tNyClVU4RKJHIBLNhInJKK0IFalM\nHJqbm7F9+3ZERUUhIyMDQqEQjDHU19fD19cXq1atwsqVKzFkiGq7TpI46CY3qm4gqSAJidcTcfT6\nUQw1HcqXmIa6hlK/ImLAUCiaUVNzis8bSCSZsLAI4sXA3NwPAoEWDYhkDLh2DQJvb9UPgpPL5bh7\n9y4AwNbWFvr6A3ejSBx0g5LaEiTdSEJSQRKO3ziOZnkzwt3CMXPETMwaOQvOls7qNpEYJDAmR339\nJd4zqK09DVNTb14MLC2nQk9Pi0qeGQPy84GkJOD4cUAsBoyMICgsVI04NDY24vvvv0deXh7GjRuH\ntWvXwsCgbyP2egOJg3Zyu/42kgo4MUi6kYTKxkqEjwhHuBu3jLEdoxWuOaH9MMbQ1JSPqqrEe4Jw\nHEOG2PFiYGUVBgMDLZsMprCQE4NWQVAogBkzgPBwbhkxQnVhpcjISAwZMgTTp09HfHw83Nzc8PXX\nX/f6RL2FxEE7qGioQHJhMo7fOI6kgiSU1ZUhxDUEM9xmIHxEOHztfClvQAwYLS237wkB5x0wJuPF\nQCSaCSMjR3Wb+HCUlbWJQVISUF8PhIW1CcKoUdw0mO1QmTiMHTsWV65cAQDIZDIEBATg4sWLvT5R\nbyFx0ExqmmpwovAEHya6UX0D05ynIdwtHDNGzMAE+wnUuI4YMGSyOtTUJPNi0NxcAiurML6qyMTE\nU7s81fJyLjx0/DgnBuXlQGgoJwQzZgDe3p3EoCMqGyHdPoSkjnASoVnUt9QjtSiVDxNl3c1CkFMQ\nwt3C8d3c7+Dv4A9DfcMHH4gg+gGFogm1tWmoqjp+L4l8GUJhIESimfD0/BlC4cQ+N64bUKqqgOTk\ntjBRcTEQHMyJwbPPAuPHA3oD63l36zno6+vD1NSUX29sbISJCTfYSCAQoLa2dmAMJM9BLTRKG3G6\n5DQfJvrz1p+Y5DCJDxNNdpwMIwPV95QnCIBrS1FTcwo1NSdQXZ2M+voL95LIMyASPQILi2nQ19ei\nwZC1tUBKSluYKDcXmDq1LWcwcSLQx5fyAWvZrS5IHAaGFnkL0krS+CTy2dKzGDdsHMJHhGOG2wxM\ncZ4CU0PTBx+IIPoBmawGNTUnUV2djJqaE5BIrsDc3A+WliGwsgqFhcUUjRmJ3CMaGoCTJ9vCRJmZ\nQEBAW84gIADo52EBJA5Er5ApZDhfdp73DM6UnMFom9GYMWIGwt3CMd1lOoRGWvTPR2g1UmkFampS\nUF3NeQaNjdkQCgNhZRUKS8sQWFgEaZdn0NQEnDnTFia6eBHw82vLGQQFAcaqLZclcSB6hEwhw6Vb\nl5BckIzjBceRWpQKNys3XgxCXENgZaxl5XyE1tLSchvV1SdQU5OM6uoTaGoqgKXlVN4zEAr9B2Qq\nzH6juRk4e5ZLIiclAenpgI9PW5ho2jTAbGAHdpI4EF0iaZEgrTQNKYUpSC1ORVpJGlytXPny0lC3\nUNia2qrbTGKQ0NxcgurqZF4QWlpuw9JyOu8ZaF0CuboaOHUKSE3lcgcXLwJjxnDlpeHhXDLZwkKt\nJmq0OCQkJGDDhg2Qy+VYt24dNm7cqLR9+/bt+Pzzz8EYg1AoxHfffYdx48YpG0ji0CPKJeVILUpF\nanEqUotSkXknExPsJ2C6y3RMd5mOqc5TYW1irW4ziUEAN+jsBp88rq5Ohlxed88r4DwDrpW1FpU6\nl5S0CUFqKnD9OhAYCEyfzgnB5MmAULPCsBorDnK5HJ6enkhMTISjoyMCAgIQFRUFLy8vfp/Tp0/D\n29sblpaWSEhIwKZNm3DmzBllA0kcOsEYw/Wq67wYpBSm4Fb9LUx1norpLtMR7BIMfwd/amVNDAit\nbaxbhaCm5gQYk93zCkJhZRUKU1MtmpJVoQCuXWsTgpQUQCJpE4Lp07n8gaFml26rfCa43pKeng4P\nDw+4ubkBAJYvX47Y2FglcZgyZQr/efLkySgpKVGVOVqNXCHH5duXkVqUipSiFKQWpUIgECDYJRjT\nXabjpYCX4GvnS4POiAGBMQUkkkzeM6ipOQGBwAhWVqEQicLh5rYJJiYe2iMGLS3A+fNtQnDyJGBl\nxQlBaCjw5puAp+cDB53pGioTh9LSUjg7tzVTc3JyQlpaWrf7//zzz5gzZ06X2zZt2sR/DgsLQ1hY\nWH+ZqZE0ShuRXprOC8HpktNwFDpiust0zB89H5898hncrNy055+P0GoUimbU119sN84gBYaG1rC0\nDIGNzTyMHPkFjI1d1W1mz6mtBU6fbvMMzp0DRo/mPILVq4HvvwcctG/OELFYDLFY3G/HU1lYKSYm\nBgkJCfjxxx8BANu2bUNaWhq++eabTvsmJSXhxRdfxMmTJyESKc+xOhjCShUNFThZfJILExWl4s/b\nf2Ks3Vg+RDTVeSqGmg1Vt5nEIKGl5RZqak6htvY0ampOob7+EkxNPWFhMeVeEjlEu3oTlZVxItDq\nGeTmAv7+bWGioCDA0lLdVvY7GhtWcnR0RHFxMb9eXFwMJyenTvtdvnwZf/3rX5GQkNBJGHQRxhiK\naop4ryClKAXFNcUIcgpCsEswPprxESY7TaYBZ8SAwJgM9fVXUFvbJgYyWfU9IZiKESM+hIVFAPT1\ntWTeeMaA7Gzl5HF1NVdKGhwMfPstMGlSvw8400VU5jnIZDJ4enri2LFjcHBwQGBgYKeEdFFREWbM\nmIFt27YhKCioawO13HNQMAUy7mQo5QukcinvFUx3mY7x9uNhoKdFZXyE1iKVVqK29gxqa0+hpuYU\n6urOwcjIGZaWU2BhMRWWllNhYjIaAm3poCuVcmWkrUKQmsqNJ2hNHAcHcyWmA9yXSBPQ2GolAIiP\nj+dLWdeuXYs33ngDW7ZsAQCsX78e69atw969e+Hi4gIAMDQ0RHp6urKBWiYOlY2VSC9Nx5mSM0gr\nTcOZkjMYajoUwa7BmO7MlZV6WGtRso7QWhhToKHh2j0hOI3a2lNobi6FhUUgLCymwsJiCiwsgmBo\nqEUe+61b3MjjtDTu67lzgLs7JwStizNNHAVouDj0B5osDlK5FJdvX+ZFIK00DTfrbsLfwR+TnSYj\nyDEIQU5BGGY+TN2mEoMAmawOdXXp9/IFp1BbmwYDAxEsLafyXoGZma/2jC9obOS8gvZiUFfHjSlo\nXYKCgEEQju4NJA4DBGMMJbUlSh7BpVuX4GblhiCnIEx2nIwgpyB4D/WmklJC5XADza7zeYLa2lNo\nbMyDubkfny+wsJiCIUO05MWEMSAvr00E0tK45nReXpwAtApBF5PaEF1D4qAiJC0SnCs7pyQGcibn\nRWCy42QEOAbAwki9Q+SJwYFc3oj6+vNKVUR6eoa8R2BhMQXm5n7Q09OSRGtVFdd/qFUM0tMBU9M2\nEZg8mWtbbUqFGb2FxKEfUDAFrt29hrSSNJwpPYO0kjTkVuZi3LBxSmJAYwuIgaC1/URd3dl7yePT\nqK+/AjMzn3ZewVQYG2tJbF0mA65cUQ4PlZZyVUPtxUALxxZoMiQOvaBcUo600jTeIzhbehY2pjZK\n4aHxw8bTZDbEgNDcXIa6urP3xOAs6urOQV/fBEJhACwsJsPCYiqEQn/o62vJW3RJiXJ46MIFwNVV\nOU/g49PnyWyI+0Pi8ABa5C24dOuSUniooqECAY4BvBhMdpxMg8yIAUEqrUBd3TklIWCsGUJhQLvF\nH0ZGw9Vtas+QSLjWE+3FoLlZOU8QEKCTg8w0HRKHdjDGUFhTqCQEl29fxijrUXz10GSnyRhjOwZ6\n2lLHTWgtMlkd6usvKAmBVHoXQuFEJSEwNtaScKVCAeTkKIeHcnIAX19lMRgxgpLGGgCJQzt+ufgL\n3jz+plJ4aJLDJJgP0ZLRnYTWolA0ob7+TyUhaGoqgLn5OCUhMDX11I4BZq1CcOEC5xmcP8+VldrY\nKOcJJkxQ+YxmRO8gcWiHXCGHnkBPO97CCK2FMRkkksx2QnAWDQ3XYGrqqSQEZma+0NPT7LbOAAC5\nnGs50SoCFy4Aly4BtrZcxdCkSdzi5wcMpfCrtkDiQBAqhDEFGhtzlYRAIrkMIyNnJSEwN5+gHXMc\nyxf3xhYAABEOSURBVGRAVlabCJw/D/z5J2Bv3yYCEydyizVNDqXNkDgQRD/BGENzc5GSENTVnYeh\noTUvBBYWATA3nwgDAy0Y3yKVcgPJ2gvBlSuAk1ObCLR6BFY0f7iuQeJAEL2AMRkaGrJRX38R9fWX\n+EUgMFTyCIRCfwwZogWhlOZmICNDOUeQmQm4uSkLwYQJap/bmBgYSBwI4gHIZHWQSC4riYBEchVG\nRk4wN5+gtAwZYq/5OaumJs4DaJ8jyMoCRo5UzhGMHw+YUzHGYIXEgSDuwRhDS0uZkgjU119Cc3MZ\nzMx8OwjBWO2Yo6ChAbh8WVkIcnK4HkPtcwTjx1OrCUIJEgdiUMKFhXLaiQAXHgIAc3O/diLgB1PT\nURAINHw0LmNcS4nLl5WX69e5+Qjah4bGjaPyUeKBkDgQOo9cXo/6+o5hoUwYGTl2ERYarvlhIYmE\nywd0FAJDQ+7B337x9gaMqI0L8fCQOBA6AxcWutlFWKgUZmY+SiJgZjYWBgZCdZt8fxQKoKCgswiU\nlHDeQHsRGDsWGKYl7bUJrYDEgdBKFIoWNDbmtvMIWsNCig5hoQn3RhVreFiopoZLErcXgStXuIlo\nOnoDo0ZxXgJBqBASB0KjYUyOxsZ8SCQZkEgyIZFkoKEhE42N+TAycoG5+VglMRgyxEGzw0IyGZCb\n29kbqKjgegx19AZoljJCTZA4EBoBYwo0NRUoCQD3NQdDhtjDzMwXZmY+/FdT0zHQ09PwpGp5eWcR\nyMri5h3o6A24uw/KSewJzYXEgRhQuFHEJfzDv00MsmBgIFISAVNTH5iZeWl+yWhdHXDtGvfgbx8a\namzsLAI+PoBQw3MdBAESB0JFMMYgld7mH/7tv+rrm3YQAG4xMNDwnv0VFZwAZGUBV6+2fb17F/D0\n5OYr9vHhxgyMGwc4O1PraUJrIXEg+oxUWtFJACSSDAC45wkoh4QMDW3UbPF9YAy4dUv54d/6tbGR\nKw318lL+6uoK6Our23KC6FdIHIgeI5NVQyK52ikvIJc3Kj3820RgmOYmhxUKoKioaxEwNOxaBBwc\nyBMgBg0kDoQSjMnQ2HgdDQ3ZaGzMRkND26JQNMDUdEynvICRkZPmioBMBuTndxaAa9e4SqCOAuDl\nRXMOEARIHAYtUuldpQd/qxA0Nd3AkCHDYWo6BqamnjAx8YSpKbdodJloUxPXM6ijCOTnc2/8HUVg\nzBial5gg7gOJgw7DDRTL7+QBNDZmQ6GQ8gLQupiYeMLExENzJ51hDLh5kxsn0Lpcu8aJQHExVw7a\nUQQ8PQETDb0egtBgSBy0HK4q6E4nD6Ch4RqamopgbOwMU9MxSh6Aqamn5uYDGOOqf9oLQE4O9zUv\nj3vQjx7NjRIeNYp7+Ht7Ax4eNGqYIPoREgctQaFoQmNjnpIH0NBwDY2N2QD0lLyAViEwMRkJPT0N\nbbpWXd21AOTmctvbC8CoUW3rFAoiiAGBxEGDUCia0dR0414oqHXJRUPDNbS0lMHY2K2TB2Bi4glD\nQ1vN9AIkku4FoKGhewGwsaGqIIJQMxotDgkJCdiwYQPkcjnWrVuHjRs3Km2/du0ann76aVy8eBEf\nffQRXn311c4Gapg4yGTVSg//pqa2zy0tt2Fs7Axj45EwMRkJExMPmJh4wNTUE8bG7tDT08CwSVMT\nl/TtSgSqqrjZxboSAHt7EgCC0GA0Vhzkcjk8PT2RmJgIR0dHBAQEICoqCl5eXvw+5eXlKCwsxL59\n+yASiTRCHFrbRnd88LeuKxRN7R7+3NK6bmzsopndQyUSrnX0jRtc3L+9ANy6xc0z3JUAODlRvyCC\n0FL6+uxU2ZMsPT0dHh4ecHNzAwAsX74csbGxSuIwdOhQDB06FHFxcaoyo0sUCimamwvbPfjz2onB\ndejrmys9+K2tI+55ASNhaGineSEgmYybI+D6dU4AbtxQ/lxTwwnAiBGcJ+DpCcyfzwmAqytgoIGC\nRhCEWlHZU6G0tBTOzs78upOTE9LS0np1rE2bNvGfw8LCEBYW9sCfkcvruw3/NDeXwsjIQemt39Jy\n2r11dxgYWPTKTpXBGNchtKsH//Xr3PSSw4ZxD393d+7rnDnc1xEjuBAQeQAEodOIxWKIxeJ+O57K\nxKE/367bi8P9qKw8jIKC99HYmAe5vBbGxu68B2Bm5gsbm4X3BMANenpD+s2+fkEi6frB3/rZyEj5\n4e/vDyxbxn12caGpJAlikNPxxfm9997r0/FUJg6Ojo4oLi7m14uLi+Hk5KSq0wEATE294O7+CUxM\nRt6bS1iD3pZlMm6gV1cP/+vXubbRrW/6rSIQGtq2TiWgBEEMICoTB39/f+Tm5qKgoAAODg7YuXMn\noqKiuty3vxLOxsYuMDZ26ZdjPTQNDdzDv6io7WtREVBYyD38y8q48E77h//cuW2fhw2j0A9BEBqD\nSktZ4+Pj+VLWtWvX4o033sCWLVsAAOvXr8etW7cQEBCA2tpa6OnpQSgU4urVqzA3b5scRiNKWeVy\nrqqn9YHfUQCKi7k3f2dnbnFxUV7c3bmvQzQslEUQhM6isaWs/YXKxYExrpqn4wO//XpZGTewq/Vh\n31EAnJ25TqD05k8QhIZA4vAgmpu5ap6u3vhb1xnr/LBvv+7oSAlfgiC0ChKH9pw/D2zfriwAFRVc\ny+eu3vZbP1ta0mhfgiB0Co0dBKcWBAJg+HBg8uQ2ARg+nKaAJAiCeEh0y3MgCIIgAPT92UkZVIIg\nCKITJA4E8f/buduQpvo/DODXwnwjlWSmmZaSNt2yzVwKFT1pBKUrsyeDrLAHipAioqCXkWgRJNSL\niuiZjOKGUtJKKjJLylQiTbKaaI+UuWpLm67v/ULcv3XuzMr9V6frA+fFdn6D78XGLnbOziEiBZYD\nEREpsByIiEiB5UBERAosByIiUmA5EBGRAsuBiIgUWA5ERKTAciAiIgWWAxERKbAciIhIgeVAREQK\nLAciIlJgORARkQLLgYiIFFgORESkwHIgIiIFlgMRESmwHIiISIHlQERECiwHIiJSYDkQEZECy8HL\nrl275u0RPErN+dScDWC+v51Hy6GkpATR0dGIiopCXl7ef67Jzs5GVFQUDAYDqqurPTnOb0ntH1A1\n51NzNoD5/nYeKwen04n169ejpKQEdXV1OHXqFB48eOC25sKFC3j06BEaGhpw4MABrF271lPjEBHR\nD/BYOdy+fRuRkZEIDw9H//79sXjxYpw7d85tzfnz57Fs2TIAQGJiIqxWK169euWpkYiIqLfEQ86c\nOSMrV650PT5+/LisX7/ebU1KSoqUl5e7HiclJUllZaXbGgDcuHHjxu0ntl/hAw/RaDS9Wtf1/f/t\n1329n4iIPM9jh5WGDx+O5uZm1+Pm5maEhob2uObp06cYPny4p0YiIqJe8lg5mEwmNDQ0oLGxEQ6H\nA6dPn4bZbHZbYzabcezYMQBARUUF/P39ERQU5KmRiIiolzx2WMnHxwd79+7FzJkz4XQ6kZWVhZiY\nGOzfvx8AsGbNGsyaNQsXLlxAZGQk/Pz8cPjwYU+NQ0REP+KXzlj0gaamJpk6darodDrR6/WSn58v\nIiItLS2SnJwsUVFRMmPGDGltbXW9JicnRyIjI0Wr1crFixe9Nfp3tbW1SUJCghgMBomJiZGtW7eK\niDqyfamzs1OMRqOkpKSIiLryjRw5UmJjY8VoNMr48eNFRF35WltbJT09XaKjoyUmJkYqKipUk6++\nvl6MRqNrGzhwoOTn56smX05Ojuh0OhkzZoxkZGRIe3t7n2bzejm8ePFCqqurRUTkw4cPMnr0aKmr\nq5PNmzdLXl6eiIjk5ubKli1bRESktrZWDAaDOBwOsVgsMmrUKHE6nV6b/3vsdruIiHR0dEhiYqKU\nlZWpJlu33bt3y5IlSyQ1NVVERFX5wsPDpaWlxe05NeXLzMyUQ4cOiUjXZ9RqtaoqXzen0ynBwcHS\n1NSkinwWi0UiIiKkvb1dREQWLlwoR44c6dNsXi+Hr82ZM0cuX74sWq1WXr58KSJdBaLVakWkq/1y\nc3Nd62fOnCm3bt3yyqw/wm63i8lkkvv376sqW3NzsyQlJcmVK1dcvxzUlC88PFzevHnj9pxa8lmt\nVomIiFA8r5Z8X7p48aJMmjRJRNSRr6WlRUaPHi1v376Vjo4OSUlJkUuXLvVptt/q3kqNjY2orq5G\nYmIiXr165To5HRQU5Lo47vnz527/egoNDcWzZ8+8Mm9vfP78GUajEUFBQZg2bRr0er1qsgHAxo0b\nsWvXLvTr97+PkpryaTQaJCcnw2Qy4eDBgwDUk89isSAwMBArVqzAuHHjsGrVKtjtdtXk+1JBQQEy\nMjIAqOP9Gzx4MDZt2oQRI0YgJCQE/v7+mDFjRp9m+23KwWazIT09Hfn5+RgwYIDbPo1G0+N1E729\npsIb+vXrh5qaGjx9+hTXr1/H1atX3fb/ydmKioowdOhQxMXFffN6lD85HwCUl5ejuroaxcXF2Ldv\nH8rKytz2/8n5Ojs7UVVVhXXr1qGqqgp+fn7Izc11W/Mn5+vmcDhQWFiIBQsWKPb9qfkeP36MPXv2\noLGxEc+fP4fNZsOJEyfc1vxqtt+iHDo6OpCeno6lS5di7ty5ALpa7+XLlwCAFy9eYOjQoQD+3Gsj\nBg0ahNmzZ+Pu3buqyXbz5k2cP38eERERyMjIwJUrV7B06VLV5AOAYcOGAQACAwORlpaG27dvqyZf\naGgoQkNDMX78eADA/PnzUVVVheDgYFXk61ZcXIz4+HgEBgYCUMd3S2VlJSZMmICAgAD4+Phg3rx5\nuHXrVp++d14vBxFBVlYWdDodNmzY4HrebDbj6NGjAICjR4+6SsNsNqOgoAAOhwMWiwUNDQ1ISEjw\nyuzf8+bNG1itVgBAW1sbLl++jLi4OFVkA4CcnBw0NzfDYrGgoKAA06dPx/Hjx1WT7+PHj/jw4QMA\nwG6349KlS4iNjVVNvuDgYISFheHhw4cAgNLSUuj1eqSmpqoiX7dTp065DikB6vhuiY6ORkVFBdra\n2iAiKC0thU6n69v3zkPnS3qtrKxMNBqNGAwG11/OiouLpaWlRZKSkv7zL1k7duyQUaNGiVarlZKS\nEi9O37N79+5JXFycGAwGiY2NlZ07d4qIqCLb165du+b6t5Ja8j158kQMBoMYDAbR6/WSk5MjIurJ\nJyJSU1MjJpNJxo4dK2lpaWK1WlWVz2azSUBAgLx//971nFry5eXluf7KmpmZKQ6Ho0+zaUR48yIi\nInLn9cNKRET0+2E5EBGRAsuBiIgUWA5Ev8jpdHp7BKI+57G7shKpxfbt23Hy5EkEBgYiLCwM8fHx\nKCoqgtFoxI0bN7BkyRJMnjwZmzZtgs1mw5AhQ3DkyBEEBwd7e3Sin8ZyIOrBnTt38M8//+DevXtw\nOBwYN24c4uPjAXRdvHnnzh10dnZi8uTJKCwsREBAAE6fPo1t27bh0KFDXp6e6OexHIh6UF5ejrlz\n58LX1xe+vr5ITU117Vu0aBEAoL6+HrW1tUhOTgbQdZgpJCTEK/MS9RWWA1EPNBrNN+8b5efnB6Dr\nKn+9Xo+bN2/+P0cj8iiekCbqwcSJE1FYWIhPnz7BZrOhqKjIta+7NLRaLV6/fo2KigoAXYeb6urq\nvDIvUV/hLweiHphMJpjNZowdOxZBQUGIjY3FoEGD3O546evri7NnzyI7Oxvv3r1DZ2cnNm7cCJ1O\n5+XpiX4eb59B9B12ux1+fn74+PEjpkyZgoMHD8JoNHp7LCKP4i8Hou9YvXo16urq0N7ejuXLl7MY\n6K/AXw5ERKTAE9JERKTAciAiIgWWAxERKbAciIhIgeVAREQKLAciIlL4FybzSFlPT742AAAAAElF\nTkSuQmCC\n",
"text": [
"<matplotlib.figure.Figure at 0xbee63cc>"
]
}
],
"prompt_number": 13
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"# plot Prob(admit=1) isolating GPA and Prestige\n",
"isolate_and_plot('gpa')"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAYYAAAEVCAYAAADkckIIAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzsnXdcVGf2/9+oIL3bKIKKIogFC3ZEjSX22NYS15qY9ttk\nk92U3RRTNslm3d1vNm4Sk41RUbErxh6M2AW7YsFKsaD03mbm+f1xYWAEbMwww/i8X695wcy93Hvu\nzPB87jnnec6xEEIIJBKJRCIpo4GxDZBIJBKJaSGFQSKRSCQ6SGGQSCQSiQ5SGCQSiUSigxQGiUQi\nkegghUEikUgkOkhhMGEaNGjA9evXa9yemppKQEAAxcXFBjn+gwgKCmL//v1P9LePyt27dwkMDKSk\npKTGfV5++WU+++yzWp1n1qxZfPDBB0/89w4ODiQkJNTKBnNiwYIFzJgxw6g2JCUl4eDggJyN/2RI\nYdAzvr6+2Nra4uDgQPPmzZk9ezb5+fkGOdeXX37J7Nmzady4sUGO/yDi4uIIDQ0FajcQrF27lj59\n+mBnZ8fAgQN1tjVr1oyBAwfyww8/1Pj33333He+///4TnbscCwsLLCwsHmnfsLAwfvrpJ53XcnNz\n8fX1rZUN5sSjvpeGpGXLluTm5mptqe5zk9SMFAY9Y2FhwdatW8nNzeXkyZMcP3682jtalUpVq/MU\nFxezfPlynn/++Vodx9i4ubnx5ptv8u6771a7ffr06SxevNjgdjzqnaUpDHpPExqNRi/HkZ/b4yGF\nwYB4eHgwfPhwzp8/Dyihm2+//Za2bdvi7+8PwI8//kjbtm1xc3Nj7Nix3LlzR+cY27Zto02bNjRp\n0oS3335bO4DFxMTg7OyMh4eHdt+ff/6ZwMBAHB0dadOmTZU77X/84x94eHjg5eXFkiVLdLbNmjWL\nV155hREjRuDg4ED//v1JSUnh9ddfx8XFhYCAAE6fPq3d39fXlz179rBz506++OIL1qxZg4ODA8HB\nwY/1Hg0ePJiJEyfSokWLareHhIRw/fp1kpOTq91eOQyUlpbGqFGjcHFxwc3NjdDQUO37dfHiRcLC\nwnBxcSEoKIhffvml2uNlZmYyatQomjZtiqurK6NHj+bWrVsA/PWvf+XAgQO89tprODg48Ic//AHQ\nDcnNmjWLV199lVGjRuHo6EivXr10wnW7d+/G398fZ2dnXn31VQYMGFDjnWxhYSEzZ87E1dWVwMBA\nvvrqK7y9vbXbfX19+fLLL+nQoQOurq7MmTNHG1bMysqq8Tqq48svv8TPzw9HR0c6dOjA5s2btduW\nLl1Kv379+POf/4yrqyutW7dm586d2u03btxgwIABODo6MnToUNLS0mo8T3R0NF5eXnzxxRc0adKE\nVq1asWrVKu32WbNm8fLLLzNixAjs7e2Jjo7m9u3bTJgwgaZNm9K6dWu++eYb7f6xsbF0794dJycn\nmjdvzltvvQVAQkICDRo0QK1W1/i5PeyzWLJkCYGBgbi6ujJ8+HCSkpJqvC6zQ0j0iq+vr4iKihJC\nCJGUlCQ6dOggPvzwQyGEEBYWFmLo0KEiMzNTFBUViT179gh3d3dx6tQpUVxcLP7f//t/IjQ0VHss\nCwsLMWjQIJGZmSmSkpJEu3btxP/+9z8hhBCLFi0SI0eO1Dn3tm3bxPXr14UQQuzbt0/Y2tqKkydP\nCiGE2LFjh2jWrJk4f/68yM/PF1OnThUWFhbi2rVrQgghZs6cKdzd3cXJkydFUVGRGDRokPDx8RHh\n4eFCo9GI999/XwwcOFDnOvfs2SOEEGLBggVixowZOra8/PLLwtnZudpH586dq7xvP/74owgLC6v2\nPe3UqZPYsmVLtdtmzZolPvjgAyGEEO+++6546aWXhEqlEiqVShw8eFAIIURJSYlo06aN+OKLL0Rp\naan47bffhIODg4iPj9ce4/333xdCCJGeni42btwoCgsLRW5urpg0aZIYN26c9nxhYWHip59+0rHh\n/vfRzc1NHDt2TKhUKjF9+nQxZcoUIYQQqampwtHRUWzatEmo1Wrx9ddfC0tLyyrHK+edd94RYWFh\nIisrS9y8eVN07NhReHt7a7f7+PiIjh07ips3b4qMjAzRt2/fR76O+1m3bp24c+eOEEKINWvWCDs7\nO5GSkiKEEOLnn38WlpaW4n//+5/QaDTiu+++Ex4eHtq/7dWrl3jrrbdESUmJ2L9/v3BwcKjyfShn\n7969olGjRtr99+3bJ+zs7LSfxcyZM4WTk5M4fPiwEEKIgoIC0bVrV/Hpp5+K0tJScf36ddG6dWux\na9cu7blXrFghhBAiPz9fHD16VAghxI0bN4SFhYVQq9XVfm4P+yw2b94s/Pz8xKVLl4RarRafffaZ\n6NOnT43vn7khhUHP+Pj4CHt7e+Hs7Cx8fHzEq6++KoqKioQQygCyd+9e7b5z5swR77zzjvZ5Xl6e\nsLS0FImJidr9y/8BhBDi22+/FYMHDxZCCPHZZ59pB5yaGDdunPj666+FEELMnj1bvPfee9ptly9f\n1hnQZs2aJV588UXt9m+++UYEBgZqn589e1Y4Oztrn1cWho8++kg8//zzj/Du1MyDhKFv374iPDy8\n2m2VheHDDz8UY8eOFVevXtXZZ//+/aJ58+Y6r02dOlUsWLBAe4zyAfV+Tp06JVxcXLTPw8LCtOJc\nzv3v4wsvvKDdtn37dtG+fXshhBDLli2rMrh4e3vXKAytW7cWu3fv1j7/3//+J7y8vLTPfX19xeLF\ni3XO1aZNm0e6jofRpUsXERkZKYRQhMHPz0+7LT8/X1hYWIi7d++KxMRE0ahRI1FQUKDdPm3atBq/\nD+XCUHn/yZMni08//VQIoQjDzJkztduOHj0qWrZsqXOMzz//XMyePVsIIURoaKj46KOPRGpqqs4+\n1QlD5c/tYZ/F8OHDdT4XtVotbG1tRVJSUrXXZW7IUJKesbCwIDIykszMTBISEli0aJFOcrhyKODO\nnTv4+Phon9vZ2eHm5qbj8lfev2XLlty+fRsAFxcXcnNzdc69Y8cOevXqhZubGy4uLmzfvp309HTt\nue4/1v00bdpU+7u1tbXOcxsbG/Ly8h79jdAjubm5ODs717hdlIWL/vznP+Pn58fQoUNp06YNf//7\n3wG4ffu2zrUD+Pj4aN/LyhQUFDB//nx8fX1xcnJiwIABZGdn6+QgHhavbtasmfb3yu/b7du38fLy\n0tn3/ueVud/u6vat6fvxKNdRmeXLlxMcHIyLiwsuLi7ExcVpvzsAzZs31/5ua2sLQF5eHrdv38bF\nxQUbGxvt9srf6eqobv/yEKqFhYXOdSYmJmrPUf744osvuHfvHgA//fQTly9fJiAggJCQELZt21bj\neSt/bg/7LBITE7Vh1PLQJPDAcJw5IYWhjqn85fTw8NCZ5pifn096ejqenp7a1yrHNZOSkrTbOnXq\nxOXLl7XbiouLmTBhAm+//Tb37t0jMzOTESNGaAeCFi1aVDmWIa6pnJdeegkHB4dqHx07dnykY4CS\npL969SqdO3d+qB329vYsXLiQa9eusWXLFv71r3/x22+/4enpSXJyss6gmJiYqPM+l5//n//8J5cv\nXyY2Npbs7Gz27duHUDzrB9r5KHh4eHDz5k3tcyGEzvP7adGihU5upbo8S03fj4ddR2USExN58cUX\n+e9//0tGRgaZmZkEBQU9UkK+RYsWZGZmUlBQoHO8B71P1e1fOVdW+W9btmxJq1atyMzM1D5ycnLY\nunUrAH5+fqxatYrU1FTeeecdJk6cSGFhYZVz3m/Pwz6Lli1b8sMPP+icNz8/n169ej30PTEHpDAY\nkalTp/Lzzz9z5swZiouL+ctf/kKvXr107uYXLlxIVlYWycnJ/Oc//+F3v/sdoCRls7KytHeIJSUl\nlJSU4O7uToMGDdixYwe7d+/WHmfy5MksXbqUixcvUlBQwMcff6xjy6MMAjXRvHlzEhISdI7x/fff\nk5ubW+3j3Llz2v00Gg1FRUWUlpai0WgoLi6mtLRUuz02NhZfX98qd/zV2b1161auXr2KEAJHR0ca\nNmxIw4YN6dmzJ7a2tnz11VeUlpYSHR3N1q1bmTJlivYY5cfJy8vDxsYGJycnMjIyqrxPzZo149q1\nazW+Fw96H0eMGMG5c+eIjIxEpVLx3//+l5SUlBr3nzx5Ml988QVZWVncunWLRYsW6QxwQgi+/fZb\nbt26RUZGBn/729+034+HXUdl8vPzsbCwwN3dHY1Gw88//0xcXFyN+1fGx8eH7t2789FHH1FaWsrB\ngwe1g/aDKN//wIEDbNu2jUmTJmmvqTIhISE4ODjw1VdfUVhYiFqtJi4ujuPHjwOwYsUKUlNTAXBy\ncsLCwoIGDaoOa/d/biNHjnzgZ/HSSy/x+eefc+HCBQCys7NZt27dI70n5oAUhjrk/ruWwYMH8+mn\nnzJhwgQ8PDy4ceMGq1ev1tln7NixdOvWjeDgYEaNGsWcOXMAsLKyYtasWaxYsQJQFln95z//YfLk\nybi6uhIREcHYsWO1xxk+fDhvvPEGgwYNol27dgwePFjHnvvn8lc3t7+mu8Dyf2o3Nze6d+/+WO/J\n8uXLsbW15ZVXXuHAgQPY2Ngwf/587faVK1fy8ssv1/j3le28evUqQ4YMwcHBgT59+mhnmlhaWvLL\nL7+wY8cOmjRpwmuvvUZ4eDjt2rWrcow33niDwsJC3N3d6dOnD88++6zOdb/++uusX78eV1dX3njj\njQfaU/k1AHd3d9atW8fbb7+Nu7s7Fy9epHv37jWuQ/nwww/x8vKiVatWDB06lEmTJmFlZaVz3GnT\npmlDZ23bttWu6XjYdVQmMDCQt956i969e9O8eXPi4uLo16/fI10TwKpVq4iJicHV1ZVPPvmEmTNn\nVnuecpo3b46LiwseHh7MmDGDxYsXV/tZgDLja+vWrZw+fZrWrVvTpEkTXnzxRXJycgDYtWsXQUFB\nODg48Mc//pHVq1dr388HfW5ubm4P/CzGjRvHO++8w5QpU3BycqJjx47s2rXrgddlVtRpRkOiV1JT\nU0X79u21yW1z4+7duyIgIEAUFxcb2xSDoFarhYeHh4iOjn6k/b/99lsxYMAA7fPKEwDqC3v37tVJ\noJsKj/tZmDsG9Rh27txJ+/btadu2rTYReD/R0dEEBwcTFBREWFiYIc0xO8rvdIyx8rkuaNq0KRcu\nXNC5S67v7N69m6ysLIqLi/n8888Baoxbp6SkcOjQITQaDfHx8fzrX/9i/PjxdWmuWfM4n8XTRiND\nHVitVvPaa68RFRWFp6cnPXr0YMyYMQQEBGj3ycrK4tVXX2XXrl14eXk9cGGMRGIOHDlyhGnTplFS\nUqJdSFaTsJeUlPDSSy9x48YNnJ2dmTp16gPDavUFU1mF/DifxdOGhRCGqTJ15MgRPv74Y+0KyS+/\n/BJAp/TBt99+S0pKCp988okhTJBIJBLJE2Awj+HWrVtV5mDHxMTo7HPlyhVKS0sZOHAgubm5vP76\n61WKsZnK3YVEIpHUN570vt9gOYZHGdBLS0s5efIk27dvZ9euXXz66adcuXKlyn6ibDqhqTw++ugj\no9tQH2wyVbukTdKmp8Gu2mAwj6F8UVE5ycnJVVYaent74+7ujo2NDTY2NoSGhnLmzBnatm1rKLMk\nEolE8hAM5jF0796dK1eukJCQQElJCWvWrGHMmDE6+4wdO5aDBw+iVqspKCggJiaGwMBAQ5kkkUgk\nTwW1dBgM5zE0atSIRYsWMWzYMNRqNXPnziUgIEBbW3/+/Pm0b9+e4cOH06lTJxo0aMALL7xQL4TB\nFKfVmqJNYJp2SZseDWnTo2MqdiUlwbJlsHRp7Y5jsFlJ+sLCwqLW8TKJRCIxVwoKYNMm+PlnOHUK\npkyBWbMgJOTJx06DeQyGxtXVlczMTGObUStcXFzIyMgwthkSiaSeIQQcPaqIwfr10LMnvPACjB0L\n1ta1P369FYbMzMx670nIqbgSieRxuHULwsOVUJEQMHs2nDsHlQoF64V6KwwSiUTyNFBUBFu2KN5B\nTAxMnAhLlkDv3mCoe0spDBKJRGJiCAEnTihisGYNdOmi5A02bICyPkkGRQqDRCKRmAh378KKFYog\nFBYqYnDiBDykKZ7ekcIgkUgkRqSkBLZtU8Rg/3547jn473+hf3+opudQnWB2jXoaNmxIcHAwHTt2\nZPLkydW2+XtcykvyltO3b99aH1MikTzdnD4Nb7wBXl7wf/8H48fDzZuKQAwYYDxRgHq8jqGm1x0c\nHMjNzQXg+eefp1u3bvzxj3/UblepVDRq9HiOUuVj6hO5RkMiebpIS4NVq5TBPyMDZs5UHm3a6P9c\ntRlfzM5jqEz//v25evUq+/bto3///owdO5agoCA0Gg1//vOfCQkJoXPnzvzwww8A3Llzh9DQUK3H\ncfDgQd59910KCwsJDg7WVn61t7cHlH7Fr7zyCgEBAQwdOpSRI0eyYcMGAE6cOEFYWBjdu3dn+PDh\nD+ztK5FIzBeVCrZuhQkTwM8PYmNh4UK4cQM++cQwolBrhIlTk4k1vW5vby+EEKK0tFSMGTNGfP/9\n9yI6OlrY2dmJhIQEIYQQixcvFp999pkQQoiioiLRvXt3cePGDfHPf/5T/O1vfxNCKK3+cnNzdY55\n/znWrVsnRowYIYQQIiUlRbi4uIgNGzaIkpIS0bt3b5GWliaEEGL16tVizpw5j3wNEomk/hMXJ8Sf\n/iRE8+ZC9O4txA8/CJGVVXfnr834YnbJ5/K7e4DQ0FDmzJnDoUOHCAkJwacstb97927OnTvH+vXr\nAcjJyeHq1av06NGDOXPmUFpayrhx4+jcufMDz3Xw4EEmT54MQLNmzRg4cCAA8fHxnD9/nmeeeQZQ\nutl5eHgY5HolEonpkJkJq1croaJbt+D3v4foaPD3N7Zlj4fZCYONjQ2nTp2q8rqdnZ3O80WLFjFk\nyJAq+x04cICtW7cya9Ys3nzzzSqNgyrzoBhehw4dOHz48GNaL5FI6htqNfz2m7LobMcOGDZMCREN\nGQINGxrbuifDrHMMNTFs2DC+/fZbVCoVAJcvX6agoICkpCSaNGnCvHnzmDt3rlZgLC0ttftWpm/f\nvmzYsAEhBHfv3iU6OhoAf39/UlNTOXr0KKA0JLpw4ULdXJxEIqkTbtyAjz6C1q3h3XehXz+4fl1Z\nkDZ8eP0VBTBDj6G6+kMWFhY6r8+bN4+EhAS6du2KEIKmTZuyadMmoqOj+cc//oGlpSUODg4sX74c\ngBdffJFOnTrRrVs3wsPDtceaMGECe/bsITAwEG9vb7p27YqTkxOWlpasX7+eP/zhD2RnZ6NSqfjj\nH/9YL0qKSySSmikshI0bFe/g7FmYNg0iI5WVyeaE2U1XrWvy8/Oxs7MjPT2dnj17cvjwYZo2bfpI\nf2sq1yCRSGpGCDh+XBGDtWshJATmzIExY6BxY2NbVzO1GV/MzmOoa0aNGkVWVhYlJSV8+OGHjywK\nEonEtElNhZUrFUHIz1fE4PRp8PY2tmWGR3oMRsQcrkEiMSdUKti9WxGDqCjFK5gzB0JDjbsS+UmQ\nHoNEIpHUgitXlCmmy5YpHsGcOfDTT+DkZGzLjIMUBolE8lSSnw/r1ineQXw8zJiheAsdOhjbMuMj\nQ0lGxByuQSKpT5S3xFyyRGmJ2a8fzJ0LI0aAlZWxrdMvMpQkkUgkDyAlRWmJuWQJaDRKqOj8eZAF\nCapHCoNEIjFLSkth+3ZFDPbvV8pa/+9/0KeP4VpimgtSGCQSiVlx8aIiBuHh0Lat4h2sXAllRZEl\nj0A9m4BVP1i0aBHdu3fH2tqa2bNnG9scicTsycmp8AYGD1bKUezbBwcOwOzZUhQeF+kxGABPT08+\n+OADdu3apZcOchKJpCpCKAP/kiWweTMMGgR/+YtSp+gxe3FJ7kO+fQbgueeeA+D48ePcvHnTyNZI\nJObFrVuwfLkiCFZWyqyir74CWXRAf5itMOgruVSb2aRyKqpEoh9KSpQuaEuWwOHDMGkSrFih1C2S\niWT9Y7bCYApjcnWVXiUSyaNz4YKyAnnFCmjfXvEO1qyB+9qrSPSMQZPPO3fupH379rRt25a///3v\nVbZHR0fj5OREcHAwwcHBfPbZZ4Y0p86RHoNE8vjk5MCPP0Lv3kqzm8aN4eBBJZn8+99LUagLDOYx\nqNVqXnvtNaKiovD09KRHjx6MGTOGgIAAnf0GDBjAli1bDGWGUZEeg0TyaAihDP4//VSRSH7/faUb\nmkwk1z0Ge8tjY2Px8/PD19cXgClTphAZGVlFGMzxrlqtVlNaWopKpUKtVlNcXEyjRo1oWJ9bOkkk\nBuDOHaVw3ZIligDMnQt//zs0a2Zsy55uDCYMt27dwrtS4XIvLy9iYmJ09rGwsODw4cN07twZT09P\nFi5cWG2XswULFmh/DwsLIywszFBm64VPP/2UTz75RPt8xYoVLFiwgA8//NCIVkkkpkFpKWzbpngH\nBw/CxInKLKOePWUiuTZER0dr2wvXFoMV0duwYQM7d+7kxx9/BJTBMSYmhm+++Ua7T25uLg0bNsTW\n1pYdO3bw+uuvc/nyZV0DZRE9icQsuH9F8ty5yuwimTMwDLUZXwyWfPb09CQ5OVn7PDk5GS8vL519\nHBwcsLW1BeDZZ5+ltLSUjIwMQ5kkkUjqmNxcxTPo00fJGzRsqNQtOnAAZs2SomCqGEwYunfvzpUr\nV0hISKCkpIQ1a9YwZswYnX3u3r2rVbTY2FiEELi6uhrKJIlEUgcIAYcOKTWKWraEX36B996D5GT4\n8kto187YFkoehsFyDI0aNWLRokUMGzYMtVrN3LlzCQgIYPHixQDMnz+f9evX891339GoUSNsbW1Z\nvXq1ocyRSCQGJiWlYkUyKKGizz+H5s2Na5fk8ZGNeoyIOVyD5Onm/tLWEyYonkLv3jKRbGxkox6J\nRFKnXLqk9EhevhzatJGlrc0NKQwSieSRyMuDtWsV7+DqVWUV8t69SqkKiXkhQ0lGxByuQWLeCKEU\nrVuyBDZuhP79K3okW1oa2zrJg5ChJIlEolfu3KlIJFtYKGJw8aJMJD8tSGGQSCSAUtp62zZFDMpX\nJC9dCr16yUTy04Zs7WkASkpKmDt3Lr6+vjg6OhIcHMzOnTuNbZZEUi3nz8Nbb4G3N/zf/ymCcPNm\nRYVTKQpPH9JjMAAqlYqWLVuyf/9+WrZsybZt25g8eTLnzp3Dx8fH2OZJJGRnw+rVindw86ayCvng\nQaVUhUQik891ROfOnVmwYIG27SfUv2uQ1G80GqWnwZIlymrkIUOUaaZDhyqlKiTmhUw+V4PFx/rx\nf8VHtR+47969y+XLl+nQoYMeLJJIHo+kJKW09c8/K+sM5s6Ff/8b3N2NbZnEVJEeg4EpLS3l2Wef\npW3btnz33Xc62+rLNUjqH0VFEBmpeAfHj8OUKYp30LWrzBk8LdRmfJHCYEA0Gg3Tpk0jLy+PyMjI\nKo166sM1SOoXp04pYhARAcHBihiMGwc2Nsa2TFLXyFCSCSKEYO7cuaSmprJ9+3bZvU1iMNLTYdUq\nRRAyM2H2bMVLKGueKJE8NlIYDMTLL7/MpUuXiIqKonHjxsY2R2JmqNUQFaWIwa5dMHIkLFwIAwdC\nAzkJXVJLZCjJACQmJtKqVSusra11PIUffviBqVOnap+b8jVITJPr15Uk8tKlyirkOXOU/IGLi7Et\nk5gaMsdQTzGHa5AYnoIC2LBB8Q7On4fp05VwUadOxrZMYsrIHINEYmYIAbGxihisW6esQH7tNRg9\nGqysjG2dxNyRwiCRmBC3b0N4uBIqUquVFcnnzoGnp7EtkzxNSGGQSIxMUZGyEvnnn+HoUaUL2k8/\nyTpFEuMhhUEiMQJCwIkTihisWQNduijewfr1YGtrbOskTztSGCSSOiQlBVasUEJFhYWKGJw4AbK2\nosSUkMIgkRiYkhLYulXxDg4ehOeeg2+/VbqhyVCRRN8IITh552StjiGFQSIxEKdOKZ5BRAR06KB4\nBxERSiE7iUTfpOSlsOLsCpadWUZeSV6tjiXXMRgRc7gGiS6pqbBypSIIWVmKGPz+99C6tbEtk5gj\nRaoifon/hWVnlnEo+RDPtX+OmZ1n0t+nPw0bNJQL3Ooj5nANEigthe3bFTHYuxfGjFEWoA0YIMtT\nSPSPEIJjt4+x9PRS1p5fS+fmnZnZeSbjA8Zjb1XhjsoFbibI888/z549e8jPz8fd3Z25c+fy17/+\n1dhmSfTI2bOKGKxcCe3aKWKwbBk4OhrbMok5civnljZUVKIuYVaXWZx48QQ+zvqfuSA9BgNx/vx5\n2rRpg7W1NfHx8QwYMIClS5cyfPhw7T6mfg2SqqSlKXmCpUuVsNHMmcrDz8/YlknMkcLSQiLjI1l2\nZhlHbx5lYuBEZnaeSV/vvlg8ZOaC9BhMkPu7tTVq1IimTZsayRpJbVCpYOdORQyiopRKpn//u1LJ\nVFZTl+gbIQRHbx5l2ZllrLuwjm4tujGryyw2TN6ArWXdLHIxW2GIjtbPPMCwsCe/o3/llVdYtmwZ\nxcXFLFq0iK5du+rFJkndcP68IgYrVkCrVkoi+aefwMnJ2JZJzJGbOTcJPxPO0jNLAZjVeRan55/G\n28m77o0RBmTHjh3C399f+Pn5iS+//LLG/WJjY0XDhg3Fhg0bqmyryUQDm643NBqN2Lt3r3BzcxMx\nMTE62+rLNTxNZGQI8d//CtGjhxAeHkK8+64Qly4Z2yqJuZJfki9WnFkhhiwfIly+dBHzf5kvjiQf\nERqNptbHrs34YrAcg1qtxt/fn6ioKDw9PenRowcREREEBARU2W/IkCHY2toye/ZsJkyYoLO9vuYY\n7ufll1/G2tqaf//739rX6ts1mCulpUqzm+XLYfduGD5c8Q6GDJGhIon+EUJwKPkQy84sY8OFDfT0\n6smszrMY4z8GG0v99WA1yRxDbGwsfn5++Jb1F5wyZQqRkZFVhOGbb75h4sSJHDt2zFCmmASlpaW4\nubkZ2wxJGUIoC9CWL1eSyX5+ynqDxYtl0xuJYUjMSiT8bDjLzizDsoEls7rMIu6VODwcPIxtWhUM\nJgy3bt2JgVQtAAAgAElEQVTC27siNubl5UVMTEyVfSIjI/ntt984duxYjVn2BQsWaH8PCwsjLCzM\nECbrjdTUVPbs2cPo0aOxtrYmKiqKdevWERUVZWzTnnpu3lSmly5frtQqmjEDDh2Ss4okhiG/JJ+N\nFzey9MxSzqSc4XdBv2PV+FV09+j+0FlFj0t0dDTR0dF6OZbBhOFRLvqNN97gyy+/1Lo8Nbk9lYWh\nPmBhYcH333/Pyy+/jBCCdu3aER4eTo8ePYxt2lNJXh5s2qSIwYkTMHGi4hn07StrFUn0j0ZoOJB4\ngGVnlrHp0ib6evfl5e4vM7rdaBo3Mlz/9/tvmj/++OMnPpbBhMHT05Pk5GTt8+TkZLy8vHT2OXHi\nBFOmTAEgLS2NHTt2YGlpyZgxYwxlVp3g7u6uN+WWPBlqNURHK2IQGQn9+sELL8CWLWCjvzCuRKLl\nRuYNlp9ZzrIzy7CzsmNW51l8Pvhzmts3N7Zpj43Bks8qlQp/f3/27NmDh4cHISEh1Safy5k9ezaj\nR49m/PjxugaaSfK5OszhGkyNCxcUMVi5Epo2VUJFU6dCs2bGtkxijmQWZrL2/FpWnFvBpbRLTA2a\nyqwuswhuHqz3UNHjYpLJ50aNGrFo0SKGDRuGWq1m7ty5BAQEsHjxYgDmz59vqFNLnjLu3YPVqxVB\nuHMHnn8eduyAoCBjWyYxR0rUJWy/sp3ws+FEXY9iWJthvN3nbYb7DceyoaWxzdMLsiSGETGHazAW\n5e0ww8Nh/34YPVqZVTRokJxiKtE/QgiO3DxC+Nlw1p1fR1DTIJ7v9DwTAyfibO1sbPOqxSQ9BolE\n3wgBhw8rnsH69RAcrISKVq4EBwdjWycxR65mXGXF2RWsOLsCy4aWzOg0w2CF60wJKQwSk+f6dcUz\nCA8HS0ulaN3p0+BthEoBEvMnvSCdNefXEH42nOuZ15kaNJXVE1fTrUU3o+cN6goZSjIi5nANhiIr\nC9auVcQgPh6mTFFCRd26ySmmEv1TpCpi6+WtrDi7guiEaJ5t+ywzOs1gSOsh9TZvUJvxRQqDETGH\na9AnlUtT7NoFQ4cqoaLhw8HKytjWScwNjdBwKOkQ4WfD2XBxA12ad2FGpxmMDxiPY+P631RD5hgk\n9RZZmkJS18SnxRN+NpyV51ZiZ2nHjE4zjFfF1ER5qDCcP3+e/fv3k5CQgIWFBb6+vvTv379KvwGJ\n5HFITFSEYMUKKCiQpSkkhuVe/j3WxCl5g+ScZKZ1nMam322ic7POT03e4HGoMZQUHh7ON998g5ub\nGyEhIXh4eCCE4M6dO8TGxpKWlsbrr7/O888/b1gD63ko6cqVK3Ts2JFJkyYRHh6us62+XIO+SEuD\ndetg1Sq4eFEpTTF9urIqWf5vSvRNYWkhW+K3EH42nINJBxntP5rnOz7P4NaDadTA/IMlBgklZWZm\nsmfPHhxqmAeYk5PD0qVLn+ikTxOvvvoqISEhT+1dSX6+UoZi5Uo4cABGjIC334Zhw2TeQKJ/NELD\nvoR9rDi3gk0XN9HdozszOs1g9cTV2FvZG9u8eoNMPhuQ1atXs2nTJgIDA7l69epT4zGUlsKvvyqe\nwdat0Lu34hmMHSvXG0gMw4XUC0re4OxKXG1cmdFpBlM7TjXJktZ1RZ0nnz/55BM+/PDDJzphnaGv\nO/QnfGNzcnL46KOP2Lt3Lz/88IN+bDFhyhefrVqlhIv8/GDaNPjXv5SaRRKJvknJSyHiXAQrzq0g\nJS+F6R2ns23aNjo262hs0+o9TyQMP/74o+kLg5HvxD/44APmzZuHh4eHWYeRzp9XwkQREUrV0unT\n4ehRaN3a2JZJzJGsoiw2XtxIRFwEx28fZ6z/WL565ivCfMNo2EDWQtEXNQpDTbkFgMLCQoMYYy6c\nPn2aPXv2cOrUKQCzCxclJSlF61auhIwMpXrppk3QubNMIkv0T0FpAb/E/0JEXAR7E/byTOtnmN9t\nPlumbNFrK0xJBTUKg4uLC7GxsTRvXrWWuLesRfBA9u3bR0JCAi1btgQgLy8PtVrNxYsXOX78uJGt\nezIyMipmFMXFwYQJ8J//QP/+0KCBsa2TmBsl6hJ2X9tNRFwE2y5vo6dXT6YFTWPZuGU4WTsZ2zyz\np0ZhmDFjBklJSdUKw9SpUw1qVH3nxRdf1L5HQggWLlxIQkIC33//vZEtezwKCpQZRatWwb59ygrk\nN99UfjY2XCMqyVOKWqNmf+J+IuIi2HhxI+3d2zM1aCr/HvZvmtrJRFVdImcl1QEff/wx165dY/ny\n5Tqvm+I1qFQQFaWEiX75BXr1UpLI48aBY/2vEiAxMYQQHL99nIi4CNacX0NTu6ZMDZrK7zr8zuwr\nmBoag9dKysnJwdHRkezsbJyc6taNMwdhqAlTuQYhlITxqlVK4bpWrZQk8uTJsvOZxDBcSL1ARFwE\nq+NWAzCt4zSmdJhCQJPqOzxKHh+DC0PXrl05efIkwcHB2oRqXSGFwXBcvKh4BqtWKYvNpk9XvIM2\nbYxmksSMSchKYHXcaiLiIkgvSOd3Qb9jatDUp6qcdV1i8HUM9X0AllRw86YytXTVKqUl5tSpsGED\ndOkiZxRJ9M/dvLusPb+WiLgIrmRcYULABP4z/D/09+lPAws5a8FUqbcFQ1xcXOr9XYZLHZUPTUlR\nBv+1a5UZRePHKwvPQkNlG0yJ/skqymLTxU1ExEUQeyuW0f6jeT/0/Xrd2+Bpo94KQ0ZGhrFNMGnu\n3YONGxUxOHkSRo2CP/1J6XEgZxRJ9E1BaQFbL28lIi6C3278xqBWg5jXdR6bp2zG1tLW2OZJHpN6\nKwySqqSlKQvN1q6FY8eUgnV/+INSsM5GrgOS6JlSdal2rcHWy1sJ8QxhatBUfh77M87WzsY2T1IL\npDDUczIyYPNmRQyOHFHWGLz0EkRGgq28UZPoGY3QaNcabLiwgXZu7ZgaNJV/Dv0nzezlFDZz4ZFm\nJV2+fJl27doRHx+Pv79/Xdilxdgzd0yRrCxl4F+7Fg4ehCFDlKmlI0eCnZ2xrZOYGxqhIeZmDOsu\nrGPt+bW42boxNWgqU4Km4Ovsa2zzJDVg0OmqX3/9Na+//vpDXzMUUhgUcnKUVchr1yqrkAcNUsRg\n1ChZylqifzRCw5HkI6y7sI4NFzfgYOXApA6TmBw4mQ5NZffG+oBBhaG6tQtdunTh9OnTT3TCx+Vp\nFobcXKWfwZo1sHcvDBigiMGYMXIVskT/aISGQ0mHtGLgYu3CpMBJTOowicAmgcY2T/KYGGQdQ0RE\nBKtWreLGjRuMHj1a+3pubi5ubm5PdDLJw8nPh23bFDGIilLaXk6eDEuXgrPM50n0jFqj5mDSQdZd\nWMfGixtxt3VnUuAkomZEyVXITzE1CkOfPn1o0aIFqamp/OlPf9Iqj4ODA507d64zA58GCgpgxw5F\nDHbtUjqeTZ4MP/4Irq7Gtk5ibpQXqysXg+b2zZkUOIm9M/fi7163OUSJaVJvi+jVd4qKFDFYu1b5\n2aOHIgbPPQfu7sa2TmJuqDQq9iXsY/3F9Wy8uBFPB08mBU5iYuBE2rq1NbZ5EgNgkBxD3759OXTo\nEPb29lVWGFtYWJCTk/PQg+/cuZM33ngDtVrNvHnzeOedd3S2R0ZG8uGHH9KgQQMaNGjAP/7xDwYN\nGlTlXOYiDMXFikewdq0SLgoOht/9ThED2f5Som9UGhXRCdGsu7COTRc30dKpJZMCJzEhcAJ+rn7G\nNk9iYAxeRO9JUKvV+Pv7ExUVhaenJz169CAiIoKAgIq4ZX5+PnZl8yvPnTvHc889x9WrV3UNrOfC\nUFICv/6qiMEvv0DHjooYjB8P1bS6kEhqRam6lL0Je1l3YR2bL22mlXMrJgZOZGLgRFq7yH6rTxMG\nL6KXmZlJcnIyKpVK+1rXrl0f+DexsbH4+fnh6+sLwJQpU4iMjNQRBrtKk+7z8vJwN5MYSmGhIgYb\nNypiEBCgiMEXX4CHh7Gtk5gbpepS9tzYw7oL64i8FImfqx8TAydy7IVjcp2B5Il4qDB88MEHLF26\nlNatW9OgUg/HvXv3PvDvbt26pdMC1MvLi5iYmCr7bd68mffee487d+6we/fuao+1YMEC7e9hYWGE\nhYU9zOw6JztbCQ9t2gS7d0O3bkqI6LPPwMvL2NZJzI0SdQlR16NYd2EdW+K30M6tHZMCJ/HRgI9o\n6dTS2OZJjEB0dDTR0dF6OdZDQ0nt2rUjLi4OKyurxzrwhg0b2LlzJz/++CMAK1asICYmhm+++aba\n/Q8cOMC8efOIj4/XNdCEQ0l37yqLzjZuhEOHlHUG48fD6NEygSzRP8WqYn69/ivrLqzjl/hfCGgS\noOQMAibg7ST7sEt0MWgoqUOHDmRmZtLsMVt5eXp6kpycrH2enJyM1wNunfv3749KpSI9Pd2k10kk\nJChewaZNcPasUpto9mwlhyBXIEv0TZGqiN3XdrP+wnq2Xt5Kh6YdmBQ4ic8HfY6no6exzZOYKQ8V\nhr/85S8EBwcTFBRE47J6zRYWFmzZsuWBf9e9e3euXLlCQkICHh4erFmzhoiICJ19rl27RuvWrbGw\nsODkyZMAJicKQiidzjZuVMQgKQnGjoV33oHBg8Ha2tgWSsyN3OJcdl7dyeb4zWy/sp1OzToxKXAS\nXz7zJR4OMkklMTwPFYbf//73vPvuuwQFBWlzDI/SIKdRo0YsWrSIYcOGoVarmTt3LgEBASxevBiA\n+fPns2HDBpYvX46lpSX29vasXr26lpejH4RQylZv2qQIQkFBRXObvn2hkaxJK9EzKXkpbInfwuZL\nmzmYdJC+Lfsyzn8cC4cspIVDC2ObJ3nKeGiOoUePHhw7dqyu7KlCXeUYVCo4cEARgs2blSql48cr\nj27dZNtLif6JT4tn86XNRMZHcjHtIs/6PctY/7E82/ZZHBvLYliS2mHQdQxvvvkmjRs3ZsyYMdpQ\nEjx8uqq+MKQwFBUp00o3bVKmlfr4KELw3HPKFFOJRJ9ohIZjt46xOX4zmy9tJqc4h7H+YxnXfhxh\nvmFYNXy8CR4SSXUIoSE39xhOTr0MJwxhYWHVho4eNl1VX+hbGHJyYPt2RQx27YLOnRUxGDdOEQaJ\nRJ8Uq4rZm7CXzZc2syV+Cy42LoxrP46x/mPp7tGdBhYNHn4QieQhqNWFZGZGkZ6+hfT0X7C0dCck\n5LzprXzWF/oQhtTUimmlBw5A//6KVzBmjCxFIdE/2UXZbL+yncj4SHZe3UlQ0yDG+o9lbPuxtHNr\nZ2zzJGZCSck90tO3kZ4eSWbmXhwcuuLuPhY3t9HY2LQxTCjpn//8p7JDDcH1N99884lO+Lg86cUl\nJVVMKz11Sul7PH680gdZ9jKQ6JtbObeU5HH8Zo4kHyHUJ5Rx7ccxut1o2fJSoheEEBQWxpOWtoW0\ntEgKCs7j4jIUd/cxuLqOwNJStxSzQdYx5ObmYmFhQXx8PMeOHWPMmDEIIdi6dSshISFPdDJDk5gI\nK1cqYnDjhrLQ7K234JlnwMbG2NZJzAkhBBfTLrL5kpIvuJpxlZHtRvJi1xfZMHkD9lb2xjZRYgYI\noSY7+zDp6YoYaDSFuLmNwdf3Q5ydw2jQoPHDD/IEPDSU1L9/f7Zv345D2eqt3NxcRowYwYEDBwxi\nUBUDH0P1fv1VCRk99xyEhspppRL9otaoOXrzKJHxkWy+tJkiVZE2XxDqE4plQ0tjmygxA9TqPDIy\ndpflC7bRuLEXbm5jcHcfi7198CMtFwADr3y+d+8elpYVX3hLS0vu3bv3RCczNEOGKA+JRF8UqYrY\nc30Pm+OV5HEzu2aMaz+O1RNXE9z80f9JJZIHUVx8m/T0X0hL20J29gEcHXvh7j4GX99PsLau+9pX\nj7TALSQkhPHjxyOEYPPmzcycObMubJNIjEJmYSbbrmxj86XNRF2PonPzzozzH8d7/d6TpaslekEI\nQX5+HOnpkaSlbaGw8Cqurs/SvPnvCQxcRaNGTka175FmJZ04cYIDBw5gYWFBaGgowcHBdWEbYNpF\n9CTmw9WMq2y/sp0t8VuIvRXLoFaDGOs/llHtRtHEromxzZOYARpNKdnZB0hLiyQ9XSkpVB4icnLq\nT4MG+g1F1kmjnrt371JUVKR1nVu2rBv3RgqDxBAUq4o5kHSAbVe2sf3KdnKKcxjRdgSj2o5iaJuh\n2FnZPfwgEslDUKmyycjYSVpaJBkZO7Gx8cPdfQxubmOxswsyaCjSoMKwZcsW3nrrLW7fvk3Tpk1J\nTEwkICCA8+fPP9EJH9tAKQwSPXEr5xbbr2xn+9Xt/HbjNwKbBDKy7UhGtB1Bl+Zd5GIziV4oKkos\nyxdEkpMTg5NT/zIxGEXjxnVXEdegwtCpUyd+++03hgwZwqlTp9i7dy/h4eEsWbLkiU742AZKYZA8\nIeWziLZf3c62y9tIzklmWJthjGg7guF+w3G3lU0zJLVHCDU5ObFkZOwgPX0LxcU3cXMbhZvbGFxd\nh9KwoXGmLht0VpKlpSXu7u5oNBrUajUDBw7k9ddff6KTSSSGJq0gjV1Xd7HtyjZ2XduFt6M3I9qO\n4L8j/ktPr540aiDnMEtqT0nJXTIydpGRsYOMjN00buyJq+tw/Py+wcmpNxYW9ft79lDrXVxcyM3N\npX///kyfPp2mTZtiby8X70hMAyEEp1NOa3MF51PPM9B3ICPbjuSrIV/h5Sj7qkpqj+IVxJQJwQ4K\nC6/i7DwYN7dnadPmHzRubF7fs4eGkvLz87G2tkaj0bBy5UpycnKYPn16nTXUkaEkyf3kFufy6/Vf\nlXzBle3YW9kzst1IRviNINQnlMaNDLMaVPJ0oXgFO8vE4NeyhWbP4ur6LI6OffQ+i0jfGCTHIIR4\naMb8UfapLVIYJEII4tPj2X5lO9uubCP2Viy9vXprE8dt3doa20SJGSCESusVpKfvoKjoOi4ug3F1\nfRZX1+F1mjh+YjQapTjc7t1Y/OUv+s8xhIWFMWrUKMaOHUu7droVIePj49m8eTPbtm1j//79T3Ri\nieRBFJYWsi9xnzZEVKIuYUTbEfwh5A8Mbj1Y1iKS6IXi4jtkZu4iPX0HmZm/Ym3dElfXZ/Hz+zeO\njr1N3isAICUFdu9W+gj8+iu4uipVQ2tBjR5DcXExK1euJCIigri4OBwcHBBCkJeXR1BQENOnT2fa\ntGlYWRm2uYj0GJ4eErMStdNJ9yXso3PzzozwG8HIdiPp2LSjLD8hqTVCqMjOPqLNFRQVJeDi8kwl\nr6Ae9NQuLoZDhxQh2LVLqR46aJAiBkOHgq8vUAcL3NRqNWlpaQC4u7vTsGHDJzrZkyCFwXwpVZdy\nOPmwdjrp3fy7DPcbzsi2IxnaZiiuNq4PP4hE8hCKi+9ocwWKV+BbJgTP4ujYy/S9AiHgypUKIdi/\nX2kxOWyY8ujZs9qKoQYRhsLCQr7//nuuXr1Kp06dmDt3Lo2MUK5UCoP5UJ4r+PXar/x6/Vf2Je7D\nz9VPmyvo4dGDhg3q7qZDYp5U9QoS7/MKWhjbxIeTkwN79lSIQUlJhRA88ww8wuQfgwjD5MmTsbKy\nol+/fuzYsQNfX1++/vrrJzpJbZDCUL9JzU8l6noUv17/lajrUQgEQ1oPYUjrIQxuPZimdrKFnqT2\nFBffruQVRGFt3QpX12dxc1O8ApNfV6DRwIkTFUJw+jT07l0hBh06wGOGUg0iDB07duTcuXMAqFQq\nevTowalTp57oJLVBCkP9orC0kINJB/n1uuIVXM+8zgCfAYoYtBmCv5u/zBVIao1GU0xOzlGtGBQV\nJePqOqTMKxiGlVVzY5v4cG7frkgaR0VBkyYVQhAaCra2tTq8QVY+Vw4bGSOEJKkfaISGMylntEJw\n9OZROjbtyJA2Q/jm2W/o6dlTNrCR1BohNOTlnSYzcw+ZmVHk5BzG1jYAV9ehtG37LY6OIabvFRQV\nwcGDFV7BzZsweLAiBF99Bd7exrZQS40eQ8OGDbGtpFiFhYXYlPXHtLCwICcnp24MlB6DyZGcnawV\ngj3X9+Bi46IND4X5huFkbdxa8pL6jxCCoqJrZGZGlYnBb1hZNcXFZTDOzoNxdg7D0tLF2GY+GCEg\nPr5CCA4eVEJC5V5Bjx4GbTNZJ2W3jYUUBuOTU5xDdEK0IgbXfiW9MJ3BrQZrw0Mtneq+w5TE/Cgp\nSSEz8zetGAihwsVlMC4uz+DiMrh+LDDLytJNGms0ukljl7oTMykMEr1Sqi4l9las1is4k3KGXl69\ntEIgS1RL9IFKlUNW1j6ysvaQmbmH4uKbODuHacXAxqYe5KMKC+HwYdi7F377Dc6dg759K8QgIOCx\nk8b6QgqDpFYIIbicflkrBPsS9uHr7MuQNkp4qH/L/thY2hjbTEk9pzxhXJ4nyM8/h4NDiNYjcHDo\navp5gpISiIlRRGDvXjh+HDp1goEDlUffvmBjGv8rUhgkj01qfip7buzRhoc0QqMVgsGtBtPMvpmx\nTZTUc5SE8RkyM6PIytpDdvYhbG0DtB6Bo2MfGjY0jUG0RlQqZfAv9wiOHoX27SuEoF8/cHAwtpXV\nIoVB8lDySvI4nHxYEYNrv3It85qcRirRKxUJ4z1lYrAXS0t3XFyeqT8JY7Uazpyp8AgOHlRKTAwc\nqJSdCA0FZ2djW/lImKww7Ny5kzfeeAO1Ws28efN45513dLavXLmSr776CiEEDg4OfPfdd3Tq1EnX\nQCkMT0ROcQ6Hkg6xL3Ef0QnRxN2Lo2uLrgxsNZAhrYfIaaQSvVBScrdMCPaQlbUHjaa0zCMYXJYw\nNvE+BRoNnD9f4RHs3w/Nm1d4BGFh4F4/O/2ZpDCo1Wr8/f2JiorC09OTHj16EBERQUBAgHafI0eO\nEBgYiJOTEzt37mTBggUcPXpU10ApDI9EVlEWB5MOsi9xH/sS9nEh9QLdPboT5hvGAJ8B9PLqJfME\nklqjUuWSnb1PO3OocsLY2XkwtrbtTdvzFAIuX67wCKKjwdGxwiMIC4MW9aBkxiNg0NaeT0psbCx+\nfn74llX6mzJlCpGRkTrC0Lt3b+3vPXv25ObNm9Uea8GCBdrfw8LCCAsLM4TJ9YqMwgwOJB5QhCBx\nH5fTLxPiGUKYTxgLhy4kxDME60bWxjZTUs9RqbLIzj5EVtZ+srP3k58fV5YwHoy//0+mnzAWAm7c\nqPAI9u5V1g4MGgQjR8LChdDSPKZbR0dHEx0drZdjGcxjWL9+Pbt27eLHH38EYMWKFcTExPDNN99U\nu//ChQu5fPkyP/zwg66B0mMAIL0gnf2J+4lOjGZfwj6uZ16nl1cvBvgMIMw3jB6ePbBqaNgS6BLz\np6TkLtnZB7RCUFh4DUfHnjg59cfJqT+Ojr1NP2GcnKwIQPmjuLjCIxg4EFq3NtoU0rrEJD2Gx3En\n9+7dy5IlSzh06JChzKl33Mu/pwhBQjT7EveRlJ1EH+8+DPAZwPejvqdbi24yRyCpNUVFiVoRyMra\nT2npPZyc+uHk1J927b7H3r4rDRqY+A3H3bu6HkFmZkWO4J13wN//qRACfWIwYfD09CQ5OVn7PDk5\nGS+vqomos2fP8sILL7Bz505c6nBVoKmRkpfCvoR92mTx7dzb9GvZjwE+A/hpzE90bdGVRg1M2GWX\nmDxCCAoKLpGdfUArBEKU4OQUipNTfzw9X8POLggLCxMvfZ6UpDSqOXhQyRHcvq3MFho0CF57DYKC\noIFcgFkbDBZKUqlU+Pv7s2fPHjw8PAgJCamSfE5KSmLQoEGsWLGCXr16VW+gmYaSbuXc0uYH9iXs\n427+XUJ9QhngM4ABPgPo0ryL7E0gqRVCqMnLO1MmAgfIzj5Aw4Z2ODmF4uzcHyenUGxs2pp2slit\nhrg4RQTKxaCoSFk/0LcvDBgAwcFQh83D6gsmOSsJYMeOHdrpqnPnzuW9995j8eLFAMyfP5958+ax\nadMmWpYlfywtLYmNjdU10EyEITk7WRsW2pe4j4zCDEJ9QgnzCWOA7wA6Nu0ohUBSKzSaYnJzj5eF\nhg6Qk3MYKysPnJ1DtTkCa2sTT7QWFEBsbIUQHDkCzZopQlD+8POToaFHwGSFQR/UR2Eo71R2OPkw\nB5MOEp0QTV5JniIEZdNHOzTtIOsNSWqFWp1PTs4RrRDk5h7Dxsa/khD0w8rKxBsh3btX4QkcOqTU\nGurUSfEG+vWDPn2gqYlfg4kihcHI5JXkcezWMQ4nH+bwzcMcvXkUx8aO9PHuQ1/vvoT5hhHgHmDa\nLrvE5CktzSQ7+yDZ2YoQ5OWdw8EhWJsjcHLqQ6NGJlzyvLx38cGDFY9795TBv1wIevSodYMaiYIU\nhjpECEFCVgKHkw9z5OYRDicfJj49ns7NOtPHuw99vPvQ26s3LRzMY5GMxHgUF9/WCkFW1n6Kim7g\n6NhLmyNwcOhp2lNHS0rg1KkKETh0SCkwV54f6NdP6U8g8wMGQQqDASlSFXHyzkkdIQDo692X3l69\n6ePdh64tutK4UWOj2Sip/2g0ReTmniIn5wg5OUfJyTmKWp2Pk1MfnJ0H4OTUv2zqqAlPUc7KUnIC\n5aGh48eVfEB5bqBvX5PqUmbuSGHQI7dzb3Mk+QiHbx7mcPJhzt49S3v39oo34NWH3t698XHykWEh\nyROjFJtL0ApATs5R8vPPYWvbHkfH3jg69sLRsRc2Nn6m/T1LStL1Bq5fh+7dK4SgVy9wMuHQlpkj\nheEJKVWXcvbuWR1vILckVxsO6uPdhx4ePbCzsjPI+SVPB2p1Hrm5x3WEACx0RMDBoRsNG5rw96y0\nVJk2evhwhRAUF1eEhPr1U6aNWpqwR/OUIYXhEUkrSOPozaNKkjj5MCfunMDHyUdHCNq5tTPtuzSJ\nSSOEhsLCK2UCoISFCgquYG/fSSsCjo69adzY23S/Z+X1hWJjKx6nTys1hXr3ltNG6wlSGKpBIzRc\nSL2g4w2k5KUQ4hmiDQv19OqJs3X9qK0uMU1UqixycmIr5QZiaNTIUccbsLfvQoMGJpyDSkuDY8d0\nheDqiXYAABJRSURBVMDKCnr2hJAQ5dGtmwwL1TOkMJQRnxbP6rjVHLl5hKM3j9LEromON9ChSQe5\niEzyxAihJj//vE5IqLg4GQeHbpVCQj1p3NiEZ6QVFiozhSqLQGqqkhsoF4GQEPD0NLalkloihaGM\n/Yn72X5lO729etPbuzdN7eTCGMmTU1Jyj5ycGK0I5OYew8qqhVYEnJx6l9UWMtEaVmo1XLpUIQAx\nMcrzwMAKAejZUykyJ2sLmR1SGCSSWqJWF5Kff47c3Fhyco6SnX0ElSodB4eQSmGhECwt3YxtavUI\nAbdu6YrAiRNKOYnKItClC1jLPh1PA1IYJJLHQK0uIC/vDHl5J8nNPUFu7gkKC69ga9sOB4ceWiFQ\nupGZ6J10drayTiAmpkIMVKoKAQgJUcJDbiYqZBKDI4VBIqkBtTqfvLzT5Oae0ApBYeE1bG0DcHDo\nhoNDV+ztu2Fv35EGDUz0TrqkBM6erfAEYmOVZjTBwboJYh8fOUtIokUKg0SC0o84L+80eXknyjyB\nkxQV3cDOrgMODt2wt1eEwM4uyHRnCRUXK83pT51SHsePK4Xl/Px0k8MdOigtKiWSGpDCIHnqUKly\nyjyAk2XewAmKipKxt++IvX3XMm+gG7a2HUy3jEROjrI+oFwETp9WGtW3aaN4A126KNNEu3UDe3tj\nWyupZ0hhkJg1KlVWJQFQfpaU3MbOrpNWAOztu2JrG2C6InDnToUAlItASgp07FghAsHBSvcxGxMu\njCepN0hhkJgNpaUZOknhvLyTlJTcxd6+y32egL9pThPVaODataoioFLpCkBwMLRrJyuLSgyGFAZJ\nvaSk5F6lxLAiBKWl6djbB+skhm1t25lmH+KSEt18wOnTcOYMuLrqCkBwMHh5ycSwpE6RwiAxaTSa\nIvLzL5Cff5a8vHPk558lP/8cGk0J9vadKyWGu5VVFDXBKaI5OcqgX1kE4uOhdWtdEejSRREGicTI\nSGGQmARCaCgqSiQ//5yOCBQVJWJj44e9fSfs7DpiZ9cJe/uOWFl5mmYhuZSUqqGg27eV+H9lL6Bj\nR5kPkJgsUhgkdY5KlaVz95+Xd5b8/DgaNXLCzq5jmQgoQmBr60+DBlbGNrkq+flw4YISDoqLUx5n\nzihTRisLQHk+QE4PldQjpDBIDIZGU0phYXzZwH9O+1OlysTOLkh7918uApaWLsY2uSrFxUrYp3zw\nj4tTxODOHaVOUIcOijcQFKQ0ovf2lvkASb1HCoOk1gghKCm5rTP45+WdpbDwCtbWLSuFgBQBsLb2\nNb1cgEoFV6/qegBxcZCQAK1aVQz+5ULQpo30AiRmixQGyWOhVueRn39eKwDl+QALi4Y6ISB7+07Y\n2gaaXsN5jQYSE3Xv/uPilMVhLVpUCED5o107aGyiK50lEgMhhUFSLSpVNgUFlygouEhBwSXy8y9S\nUHCe4uLb2NoGaO/+y39aWTUztsm6CKEkfe/3AC5eBGdn3bv/oCAICAA7E26PKZHUIVIYnmKUENCd\nssG/sgBcRK3OwcbGHzu7AGxtA7C1bY+tbSC2tm1Nb3FYWlpVDyAuTukhXPnuv0MH5eEsO+9JJA9C\nCsNTgBAqCguvVeMBXKJBA2tsbdvfJwABNG7sZVp5AI0Gbt5UEsHlj4sXFQEoKqrqAXToAE1lsyWJ\n5EmQwmBGqNX5ZYO/IgDlg39R0TWsrDywtQ24TwDam17zmNxc3cG//HHlitI32N+/4tG+vbIewNNT\nzgSSSPSIFIZ6hhCC0tLUasI/lygtvYeNTdtKAqDc/dvYtDOtJLBarcz2qU4AsrOhbVtdAfD3V5LA\njo7GtlwieSowWWHYuXMnb7zxBmq1mnnz5vHOO+/obL906RKzZ8/m1KlT/O1vf+Ott96qamA9FgYh\n1BQVJVYb/wdRdtevKwDKNFATqguUkVH94H/9uhLmuX/w9/dX6gLJHsISiVExSWFQq9X4+/sTFRWF\np6cnPXr0ICIigoCAAO0+qampJCYmsnnzZlxcXOqlMAihoqgokcLCq/c9rlBUlIClZZNqBcDSsqnp\nlIMoKVEG+uoEoLi4+sG/bVuwtTW25RKJpAZqM3YabGpKbGwsfn5++Pr6AjBlyhQiIyN1hKFJkyY0\nadKEbdu2GcoMvaDRlFBUlFDD4J9E48YtsLHx0z6cnQdgY+OHtXVrGjY0kcFTCKUG0JUrVQf/pCRl\ntW95uCckBGbMUJ43by5j/xLJU4bBhOHWrVt4e3trn3t5eRETE/NEx1qwYIH297CwMMLCwmppXVU0\nmiIKC29UGfgLC69SXHyLxo29dAZ/F5ch2Nq2xdq6lem0iSwuVuL+164pHkDlnzduKHP8K8f++/VT\nfrZpA1YmWMtIIpE8MtHR0URHR+vlWAYTBn2GSSoLQ21QqwsoKrquHfgLCq5ofy8pScHa2qds4G+L\nrW173NxGld35+5hGETghlJj//YN++e/37il3/q1bK4N969YQGqr8bN0aHByMfQUSicRA3H/T/PHH\nHz/xsQwmDJ6eniQnJ2ufJycn4+XlZajTaVGr8ygsvFblrr+w8CqlpelYW7fS3vXb23emSZMJZYN/\nS9NY9KVSQXJy9Xf9168r+5QP+m3aQM+eMG2a8tzbW9b+kUgktcZgo0j37t25cuUKCQkJeHh4sGbN\nGiIiIqrdV1/J5ZSUZVy+/DI2Nm3+f3v3GhJVGsYB/D9abZbd7GKbuZurRmrNTKZJu3QRoZslgQUZ\nmd0sumLLEm0QRUW0u/RBdr8UlNCV6AJpWUF0My3sohUZYQulY9ZKWaOW6+g8+2GY2U6NejTPObPL\n/weHYTqv8Fef3kff856jZ/Lv0yceQ4akISAgEl99FeIbO37q6rxP+n/+6WoKwcHKyT811fUaHg4M\nGMA1fyLSlKbbVc+fP+/Zrrps2TL8/PPP2Lt3LwBg5cqVePnyJeLj42G32+Hn54c+ffqgrKwMgYGB\n/wbswJV1p9MBk8nf+Lt9HQ7XHb4VFa61/U+XfOrrlcs97kn/u++AESP4wDci+mI+uV21q/jkdlW7\n3fV0z4oK769//eXazfPtt67DPfG7J3/u9CEijbExdCWn07Wts7VJv6LCte/fPel/883nryEhXOsn\nIkOxMXREY6NrHf/5c+8Tv83menJnWxN/UBB/4icin8bG4ObeztnWT/u1ta5HNrQ26YeG8g+8E9F/\nHhuD2/79wE8/eZ/w3a9Dh/I5PkT0v8fG4OZ0ctInIsKXNYb/1yzKpkBE9MU4kxIRkQIbAxERKbAx\nEBGRAhsDEREpsDEQEZECGwMRESmwMRARkQIbAxERKbAxEBGRAhsDEREpsDEQEZECGwMRESmwMRAR\nkQIbAxERKbAxEBGRAhsDEREpsDEQEZECGwMRESmwMRARkQIbAxERKbAxEBGRAhtDJ1y9etXoCJ/x\nxUyAb+ZiJnWYST1fzdVZmjaGCxcuYNSoUYiMjMQvv/zidcz69esRGRkJi8WCkpISLeN0GV8sAl/M\nBPhmLmZSh5nU89VcnaVZY2hpacHatWtx4cIFlJWV4dixY3j8+LFiTH5+Pp4+fYry8nLs27cPq1at\n0ioOERGppFljKC4uRkREBEaMGIHu3btj/vz5OHPmjGJMbm4uMjIyAAAJCQl4+/YtXr16pVUkIiJS\nQzRy4sQJWb58uef9oUOHZO3atYoxs2bNksLCQs/7pKQkuXPnjmIMAB48ePDg0Ymjs7pBIyaTSdU4\n19zf+sd9ep6IiLSl2VJSSEgIKisrPe8rKysxfPjwNsfYbDaEhIRoFYmIiFTQrDHExcWhvLwcz549\nQ1NTE44fP46UlBTFmJSUFBw8eBAAcOvWLfTv3x/BwcFaRSIiIhU0W0rq1q0b/vjjD0ybNg0tLS1Y\ntmwZoqKisHfvXgDAypUrMXPmTOTn5yMiIgK9e/dGTk6OVnGIiEitTl+d6EIVFRUyZcoUiY6OlpiY\nGMnOzv5szOHDh8VsNsuYMWPk+++/l/v37xueya24uFj8/f3l1KlTmmbqSK4rV66I1WqVmJgYmTx5\nsuGZampqZNq0aWKxWCQmJkZycnI0zfThwwcZP368WCwWiYqKkk2bNnkdt27dOomIiBCz2Sz37t0z\nPJPeda726ySib52rzaVnnavJpHeduzU3N4vVapVZs2Z5Pd/ROveJxlBdXS0lJSUiIlJXVycjR46U\nsrIyxZiioiJ5+/atiIicP39eEhISDM8k4vqGJCYmSnJyspw8eVLTTGpz1dbWSnR0tFRWVoqIq1iN\nzrR161bPf6SamhoJCgoSh8Ohaa6GhgYREXE4HJKQkCAFBQWK8+fOnZMZM2aIiMitW7c0ryk1mfSu\nczWZRPSvczW59K5zNZmMqHMRkT179siCBQtk9uzZn53rTJ37xCMxhg4dCqvVCgAIDAxEVFQUXrx4\noRgzYcIE9OvXD4DrngebzWZ4JgD4/fffMXfuXAwePFjTPB3JdfToUaSmpnou9g8aNMjwTF9//TXs\ndjsAwG63Y+DAgejWTbOVTABAr169AABNTU1oaWlBUFCQ4rwR99G0l0nvOleTCdC/ztXk0rvO1WQy\nos5tNhvy8/OxfPlyr7s4O1PnPtEYPvbs2TOUlJQgISGh1TH79+/HzJkzDc9UVVWFM2fOeO7YVrtF\nV+tc5eXlePPmDRITExEXF4dDhw4ZnikzMxOPHj3CsGHDYLFYkJ2drXkWp9MJq9WK4OBgJCYmIjo6\nWnG+qqoKoaGhnvfDhw/XfCJuL9PH9KpzNV8nI+q8vVxG1Hl7mYyo8w0bNuC3336Dn5/36bwzde5T\njaG+vh5z585FdnY2AgMDvY65cuUKDhw40Oqzl/TMlJWVhd27d8NkMkFcy3K6ZGovl8PhwL1795Cf\nn4+LFy9ix44dKC8vNzTTrl27YLVa8eLFC5SWlmLNmjWoq6vTNI+fnx9KS0ths9lw/fp1r8+z+fR7\npvWkpyYToG+dt5fJqDpvL5cRdd5eJr3r/OzZsxgyZAjGjh3b5velo3XuM43B4XAgNTUVCxcuxJw5\nc7yOefDgATIzM5Gbm4sBAwYYnunu3buYP38+wsLCcOrUKaxevRq5ubmG5woNDcXUqVMREBCAgQMH\nYtKkSbh//76hmYqKijBv3jwAQHh4OMLCwvDkyRNNM7n169cPycnJuHPnjuLfjbyPprVMgP513l4m\no+q8vVxG1Hl7mfSu86KiIuTm5iIsLAxpaWm4fPkyFi1apBjTqTrvsqsfX8DpdEp6erpkZWW1Oub5\n8+cSHh4uN2/e9JlMH1u8eLEuuzXU5Hr8+LEkJSVJc3OzNDQ0yOjRo+XRo0eGZtqwYYNs27ZNRERe\nvnwpISEh8vr1a80y1dTUSG1trYiIvH//XiZOnCiXLl1SjPn4otzNmzc1v9CrJpPeda4m08f0qnM1\nufSuczWZ9K7zj129etXrrqTO1Lm2V0VUKiwsxOHDh2E2mzF27FgArl/JKioqALjuedi+fTtqa2s9\n65zdu3dHcXGxoZmMoCbXqFGjMH36dJjNZvj5+SEzM7PNtWw9Mm3evBlLliyBxWKB0+nEr7/+6vUi\nZ1eprq5GRkYGnE4nnE4n0tPTkZSUZOh9NGoy6V3najIZQU0uvetcTSa96/xT7iWiL61zkwgfRkRE\nRP/ymWsMRETkG9gYiIhIgY2BiIgU2BiIiEjBJ3YlEfmqHTt24MiRIxg8eDBCQ0Mxbtw45OXlwWq1\n4tq1a2hubsaBAwcQHx+P4uJiZGVlobGxEQEBAcjJycHIkSON/hSIOoyNgagVt2/fxunTp/HgwQM0\nNTUhNjYW48aNg8lkwocPH1BSUoKCggIsXboUDx8+RFRUFAoKCuDv749Lly5h8+bNOHnypNGfBlGH\nsTEQtaKwsBBz5sxBjx490KNHD8yePdtzLi0tDQAwceJE2O122O12vHv3DosWLcLTp09hMpngcDiM\nik70RXiNgagV7mcDqbVlyxYkJSXh4cOHyMvLQ2Njo4bpiLTDxkDUih9++AF5eXn4+++/UV9fj7Nn\nz3rOHT9+HABw48YN9O/fH3379oXdbsewYcMAgH+NkP7TuJRE1Iq4uDikpKTAbDYjODgYY8aM8fyt\nhJ49eyI2NtZz8RkANm7ciIyMDOzcuRPJycm6P4adqKvwkRhEbWhoaEDv3r3x/v17TJ48Gfv27cOP\nP/6IPXv2IDY21uh4RJrgbwxEbVixYgXKysrQ2NiIxYsXex4SSPR/xt8YiIhIgRefiYhIgY2BiIgU\n2BiIiEiBjYGIiBTYGIiISIGNgYiIFP4B0G8Uzg3sdTwAAAAASUVORK5CYII=\n",
"text": [
"<matplotlib.figure.Figure at 0xbdc7dcc>"
]
}
],
"prompt_number": 14
},
{
"cell_type": "raw",
"metadata": {},
"source": [
" "
]
}
],
"metadata": {}
}
]
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment