Skip to content

Instantly share code, notes, and snippets.

@903124
Created August 17, 2020 16:34
Show Gist options
  • Save 903124/5961252d89b6a695c28e9007e948fef0 to your computer and use it in GitHub Desktop.
Save 903124/5961252d89b6a695c28e9007e948fef0 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {
"ExecuteTime": {
"end_time": "2020-08-17T16:33:59.612164Z",
"start_time": "2020-08-17T16:33:35.849252Z"
}
},
"outputs": [],
"source": [
"import pandas as pd\n",
"from scipy import stats\n",
"import numpy as np\n",
"import matplotlib.pyplot as plt\n",
"import seaborn as sns"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {
"ExecuteTime": {
"end_time": "2020-08-17T16:33:59.988413Z",
"start_time": "2020-08-17T16:33:59.614166Z"
}
},
"outputs": [],
"source": [
"rapm_data = pd.read_csv('EH_rapm_sk_stats_ev_regular.csv')\n",
"gar_data = pd.read_csv('EH_gar_sk_stats_regular.csv')\n",
"xgar_data = pd.read_csv('EH_xgar_sk_stats_regular.csv')"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {
"ExecuteTime": {
"end_time": "2020-08-17T16:34:02.692298Z",
"start_time": "2020-08-17T16:34:02.647329Z"
}
},
"outputs": [],
"source": [
"metric_df = pd.merge(rapm_data,gar_data,on=['Season','Player'])\n",
"metric_df = pd.merge(metric_df,xgar_data,on=['Season','Player'])"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {
"ExecuteTime": {
"end_time": "2020-08-17T16:34:04.089248Z",
"start_time": "2020-08-17T16:34:03.455474Z"
}
},
"outputs": [],
"source": [
"\n",
"metrics = ['G±/60','xG±/60','C±/60','GAR/60','xGAR/60']\n",
"metric_matrix = np.zeros((len(metrics),len(metrics)))\n",
"year_range = range(7,19)\n",
"for year in year_range:\n",
" df_year_y_1 = metric_df[metric_df.Season == \"{0:0=2d}\".format(year) + '-' + \"{0:0=2d}\".format(year+1)]\n",
" df_year_y_2 = metric_df[metric_df.Season == \"{0:0=2d}\".format(year+1) + '-' + \"{0:0=2d}\".format(year+2)]\n",
" \n",
" merge_df = pd.merge(df_year_y_1,df_year_y_2,on='Player')\n",
" for i,metric_y_1 in enumerate(metrics):\n",
" for j,metric_y_2 in enumerate(metrics):\n",
" slope, intercept, r_value, p_value, std_err = stats.linregress(merge_df[metric_y_1 + '_x'], merge_df[metric_y_2 + '_y'])\n",
" metric_matrix[i][j] += r_value**2 / len(year_range)"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {
"ExecuteTime": {
"end_time": "2020-08-17T16:34:04.700744Z",
"start_time": "2020-08-17T16:34:04.183313Z"
}
},
"outputs": [
{
"data": {
"text/plain": [
"Text(0.5, 1.0, 'Average R squared of metrics 2007-2019\\n (min TOI 50 mins)')"
]
},
"execution_count": 5,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXIAAAElCAYAAAD5r2lGAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy86wFpkAAAACXBIWXMAAAsTAAALEwEAmpwYAABo3ElEQVR4nO2dZ5gUxdaA3zMzG2BhWTaTQQkKCiLJBAhGQC+Y8ZoT6jVfrzlhwnTN+plzuJgFBREUFSQoiIAEkQwLbCTDxpnz/aje3ZmNA2xvwHp5+qG761TVOds9p6tPVVeJqmKxWCyWhounrhWwWCwWy75hHbnFYrE0cKwjt1gslgaOdeQWi8XSwLGO3GKxWBo41pFbLBZLA8c6ckudIiJrROT4vcz7kIhki0h6Teu1N4jIeSIyua71sPz92K8cuYj8KCJbRCSqrnXZV0RktIgUishOEdkqIjNF5Mi61qu+ICJtgJuBrqqa6nJd7UVERcRXlZyqfqCqJ9ZgvUeIyBQR2SwiWSLyiYi0CEoXEXlMRHKc7XERkTJ6/yAiu0Xkz+AHpojc6dxbxVuuiAREJLESXYaJyM/OvZguIq+JSNOg9CgReVNEtjvp/y6T/1URWebUcXGZtCgReVpENjq/3/8TkYga+BP+bdhvHLmItAf6Awr8w4Xyq/wRu8RHqtoESAR+AD6pAx3Cppb/Ru2AHFXNrMU6K8Ul25sDrwLtMfbuAN4KSh8FjAB6AN2BU4Arg9L/B/wOJAB3AZ+KSBKAqo5R1SbFG/AY8KOqZleiSzPgIaAlcDDQGngiKH000MnRcxBwq4icHJS+APgXMK+Csm8HegOHAJ2Bw4G7K9HDUhGqul9swL3ADOAp4GvnXBSwFTgkSC4JyAWSneNTgPmO3Eyge5DsGuA2YCGQD/gwN91KzI9qCXBakLwXeBLIBlYD12IeLD4nvRnwBrAJ2ID5YXgrsWc08H7QcVenrKRK5DsCPwHbnPo/Cko7AfjTSXvBkbu8knral9H5EmCpY+8q4Mog2WOBNOdvlA68h2kcFP+NcoCPgfigPBcAa520u5y/8fGV2NQMeBfIcvLc7ZR/vHMNA8BO4O0K8hbrdiuQ6fzNRwBDgb+AzcCdQfKV6g2sc/4mO53tSOBizP32tFPWQ865n4PK7AZMcdIziusD+gJzge3O+afCvMcPB3YEHc8ERgUdXwbMdvY7Y+7ZpkHp04GrKihXHLsv2oPf2+nAH0HHG4ATg44fBMZWkO9n4OIy5+YCZwUd/xNYX9c+pSFt+02LHLgQ+MDZThKRFFXNBz4Hzg2SOxv4SVUzReRw4E1MKyYBeAUYXyY0cy4wDIhT1SLMDd8f42TuB94Pet29AhgCHIb50Y0oo+M7QBHG6fYETgQur84wEYl07MsBtlQi9iAwGdOKaw087+RNBD7DOMFER/+jq6sziEzMwy4W49Sfdv5uxaQC8ZiW2CjgeozdAzGtty3Ai44uXYGXMM68JeZv3rqKup/H/J0PcMq7ELhEVb/D/J03qmlRXlxJ/lQgGmiFedC/BpwP9MJcw3tF5ABHtlK9gQHO/3FOfbOc436Yh1sy8HBwxU7Y4TtgklNeR+B7J/lZ4FlVjQUOxDw0wmEAsDjouBumpVvMAudccdoqVd1RSXow/YEUzH0SLiW6iEhzjI2V6VId4mzBx61FpNke6PP3pq6fJDWxAccAhUCic/wncJOzfzzmhi6WnQFc6Oy/BDxYpqxlwEBnfw1waTV1zweGO/tTCW2xHo/TusX8UPKBRkHp5wI/VFLuaKAA86bgxzjxY6vQ413Ma3jrMucvxGmlOceCaamG1SKvoJ4vgRuc/WMdHaOD0pcCxwUdt3CujQ/jTMcGpcU4+cu1yDFvN/mYGHjxuSsxr//FdadV8fc4FtNq9zrHTR27+gXJ/AaMCEPvcn8TTOt7XZk6L8ZpkTvX9vdKdJuGaQQk7sE93h3Tsu8fdM4PHBR03MnRUzAPy9llyniYit9e3qjofBW6nIB50HV2jts49UaXkVlTQd6KWuQPYX6XSZiH7y9OeS3C1envvu0vLfKLgMlaGt/70DkHxrk2EpF+ItIO01r+wklrB9zsdOBsFZGtmJuyZVDZ64MrEpELRWR+kPwhmJYuTr71leRtB0QAm4LyvoJpzVXGx6oah3kILMK0JCvjVswP+FcRWSwil1akk5pfzvoK8leIiAwRkdlOh9tWTGgiuEMsS1Xzgo7bAV8E2bgU43BSKtBlF+YBVRGJQCQmpFLMWkzrOlxyVNXv7Oc6/2cEpecCTcLQuzKq+ju2wbz9VMRlmNDHnyIyR0ROqaIcRKQj8A3mATo9KGkn5k2pmFhgp3ONy6YVpwe30BGRRsBZmLfF4nP9gzpBF5eRPwLz+zpTVf8K0qO4/ErrqoKHMbH8+Zhw0ZeYh2i96P9oCNRFB16N4tyIZwPeoGFoUUCciPRQ1QUi8jGmhZSBiZ8X32DrgYdV9eFyBZdSMj2k8yB4DTgOmKWqfhGZT+lr4SZCQwVtgvbXY1qYiWpCNGGjqtkiciUwR0Q+VNVNFcikY0I7iMgxwHciMs3RqUQPZ1RDsF67gMZBx6lBslGY1+0LgXGqWigiXxL6Glx2+sz1mLeYGWV1FJFNmI6y4uPGmPBKRWRjfsztMH0RAG0xsVg3qErvdpXkKWt72fLOrShBVZcD54qIBxNr/lREEpwHW0V1f4d5c3yvTPJiTEfnr85xD0pDL4uBA0SkadD93gPjhIM5HdPS/zFIv+mUPuCCdekJjMf8nb4Pkt/iXNsemD6BsrpUiarmYvqTrnXqGQX8FvQQtlTD/tAiH4FpOXXFtLYPwziL6RgHBObmPQc4j9Ab+TXgKqe1LiIS4wyzakrFxGB+vFkAInIJpkVezMfADSLSSkTiMJ2AADjOdzLwpIjEiohHRA4UkYHhGKmqfwLfYlre5RCRs0Sk+CGyxdHTD0wAuonI6c7IiusJctaYVtAAEWnrxCTvCEqLxDwUs4AiERmCietXxcvAw8XOT0SSRGS4k/YpcIqIHOPE/R+gknvQ+RF/7JTV1Cnv38D71dS/t1SldxamY/WAyjJXwNdAqojc6Ayvayoi/ZyyzxeRJFUNYEJnYK5VCCLSCvNG+aKqvlxBHe8C/3but5aY4ZhvAzit5fnAfSISLSKnYcIzZePgFwHvOq34ShGRQzDx/utU9atKdLlbRJqLyEGYRsXbQfkjRSQa0wiIcHTyFNspIi2d3+ARwD3AfVXpYylDXcd29nXD3FxPVnD+bMxIiuLRFyswLY/IMnInA3MwP6hNmCF+TZ20NZSJ32JeAzdjWoxPEToCxIcZxZCDGbVyE6ZVKU56M0xcPg0zguR3YGQldo0mKHbtnOuHaUEnVyD/OKa1uhPzSj+qjI1/UcGoFSf9Rcf+FZgfYEk8GLgG8yazFTMqZSzwkJN2LGXi1BjH/G9MX8MOR5cxQekXYUaBhDNqpTnGcWdhWrj3Ap7K6i6TNyTduTYKtA869zNwfph6P+DosRU4gjIjVByZkHOYh/z3mAdrOnC7c/59TNhgJ6bVOqISG+4jdLTMTkzopDhdnOu+2dkex7nXnPT2mJZ2rmNX2Xu5FU7nexi/s7coHSVUvC0OSo/CDBwoHonz7zL5f3RsCd6OddIGOPfBbkfP8+rarzS0rdjBWFzAacG+rKqVvZrXCSLyI+Yh8Xpd62KxWPad/SG0Um8QkUYiMlREfM5r8X2UdqxaLBaLK1hHXrMIZljZFkzYZCkmHGCxWCyuYUMrFovF0sCxLXKLxWJp4FhH/jdBRP4nIiP2Mu/LInJPDavU4BGRb0TkouolqyzjHyIytqZ0svw9sY78b4CIdMd8oDFub/Kr6lWq+uBe1Ls46AtBv4jkBR3f6ci0FpEPxEzDuktEfi37paOYKWQ7VlJH8HS/xdsBQemVTuW6r6jqEFV9p3rJKssYDxziXCOLZa+wjvzvwZXAB1rLHSKq2k1Lp0mdDlyrpVOnjhGReMxY7gLMBEuJmHH4H4rImXtQ1UdB5TZR1VVBaZVO5VqP+B9mwjGLZa+wjvzvwRDMR0AAiMjFIjJDzGT+W0VklYgc5ZxfLyKZwSEDEXlbRB5y9o8VkTQRudmR2+R84bo33IT5sOQyVU1X1VxV/R/mo6snRUSqzl41IlI8t/V9TtmfAX8AZ1QiP1rM4g3vi8gOEflDRDqLyB2OretF5MQg+R9F5HJn/2IxCy/8V8ziCKud7wgISl/llLtaRM4LqvpHzAybFsteYR35fo6IxAAdMF/MBdMPM896AmbagrFAH8x0q+cDL4hIufk2HFIxX6m2wkwA9aKYqUz3lBOAz9R8qh7Mx5h5VTqHWc6pYib1WiwiVwed35OpXEvKwnzB2hzTkv8W8ztphfm685Uq8vbD/J0TMV9ZvuF8dh4DPAcMUdWmwFGYz+eLWQq0F5Gyk1xZLGFhHfn+T5zzf9mZ6Far6ltq5jT5CDOR1gOqmq+qkzHhjgrj0phpBx5Q1UJVnYhpVXfZC90SMdMilGVTUHp1fIyZWycJM73AvSJSPFlVE8y0BMFsw0xpWxnTVfVbNRObfeKU+6iqFmIedu3FzKNTEWtV9TXnb/oOZirc4tkTA5hYeCNV3aSqwRNKFV+bysq1WKrEOvL9n63O/2WdV9npXFHVyqZ4LUuOhs7guLsK2arIxji7srQISq8SVV2iqhtV1a+qMzGLNhTH18OayrUMZf8G2Vp+KtzKbC1ZBFpVdxfLqpnV8BzgKsw0xhOciaWKKb42W6vQy2KpFOvI93McJ7KS8MMUtcl3wBnFs+AFcTZmkqy/ymepluKFFSBoKteg9LCnV61JnFb+CZiH1J+YmTeLORizCMP22tbLsn9gHfnfg4mYJczqG09jWshviEiqM7XpuZjRJbeEM8pGRIY7U6eKiPTFTNM7DvZoKldXEZEUZ7x4DGZO+p2ETls7ELNwhMWyV1hH/vfgVeC8fR0FUtOoag5mmb5ozOIROZipZC9Q1Y/CLGYkZvrdHZg5sR8rM7Z7JGaF9i3Ao5iVbbJqxoKw8WDmCt+ImW52IGZF+WLOpepOVIulSuxcK38TRORDzNJxX9a1LpZSRORUzIPr7LrWxdJwsY7cYrFYGjg2tGKxWCwNHOvILRaLpYFjHbnFYrE0cHx1rUBlTIrvst8F7wNlP0TfT2jatF4NhqkR1qftd7cfAHHN6loDdxi6ddk+34Ry9RFhX3R9aXa9uulti9xisVgaOPW2RW6xWCy1iXjqVSN7j7CO3GKxWACPr+EGKKwjt1gsFmyL3GKxWBo89WwGiz3COnKLxWLBtsgtFoulwWMducVisTRwrCO3WCyWBo4dtWKxWCwNHNsit1gslgaOdeQWi8XSwLHDDy0Wi6WBY1vkFovF0sCxnZ0Wi8XSwLEtcovFYmngWEdusVgsDZyG7MgbblCoEhKP60//XybRf+5kOtxwRYUyBz9yF/3nTubo6eOJ7d615Pwhz49h0LKZHD3jq3J52l5xPv1/mcTRM7+m8+hbXNO/MhKP68+AOZMYOG8yB9xYsV1dH7uLgfMmc8yM8cT2KLXr0BfGcNzymfSfGWpXp9uvZfCSaRwz/UuOmf4lSScMcNWGsjQ/9hh6/fgNvad/S+t/VWzTAfffRe/p33L45HHEHGJsimyRyqEfvUOvqRM4/LuvaHnpBSXyvrhmHPLBG/SeNolDPngDX7PYWrElmBYn9eeUPydx6vLJdL2tYrt6PXsXpy6fzJAF42nes/Ra9XtjDKdnzGToH6HXqs2ZJzN00dec619KfK9DXNW/Mty4BwHajTqfAXMm0X/W13S5v/Z/W8WIR8Le6hv7lyP3eOj6+L3MPftyfj5yGC3OOIWYLgeGiCQeP4DGB7Zneu8TWXTTPXR9cnRJ2oYPP+e3sy4vV2z8Mf1IHnIcP/c/lRlHncKaF95w25JQPB66/fde5px5OdP6DaPlmafQpIxdSScMoPEB7fnp8BNZdMM9HBJkV9qHnzPnzPJ2Aaz+v7f5uf8Ifu4/gqwp09y0IhSPhwMfupfFF17Bb4NPIWn4MBp3CrWp+aABNOrQjrn9T2L5bffSccx9AKjfz6oHH+O3wcNYMHwkLS46ryRv639dwdYZs5k74GS2zphd6QPCLcTjofeL9/LDkMuZ0HUY7c49hdiDQ+1qOWQATTu156tOJ/LrqHvo89LokrRVb3/ODyeXv1bbFv3F9NOvI3PaHLdNqBiX7sH4/v1IGXocPx99KtOPPIXVz9fybysI68jrCXG9urN79Vpy16ahhYWkfz6BlCHHhcikDD2OjWO/BGDb3AVExMYSlZIEwJZZcyncsq1cuW0uPZfVz76KFhQCUJC92V1DyhDXqzu7V5XatemzCaQMLW/XBseurXMX4GsWZNfMiu2qS5oe1p28NevIW2dsyho/kfgTQ21KOPE4Mj8bB8CO3xfgi40lIjmJwswsdi1aAoB/1y5yV6wkMjWlJE/Gp18CkPHplyScdHztGQUk9O3OzhVr2bU6jUBhIWvHTqD18FC7Wg0/jtXvGh1zfllAZFws0anmWmVNn0vB5vLXavufq9jx12rX9a8Mt+7Bdpeey8qnXyVQR7+tYDxeT9hbfcM1jUSkmYicIyL/FpGbnP04t+oDiGqRQu6G9JLjvI0ZRLVIqUYmvZxMWWIObE/zI3tzxJSP6fvVe8T2PLRmFa+G6BYp5AXpnFuBXWVl8jamE12NXQDtRp3HMTPGc+gLY2o1DBGVmkL+xk0lxwWb0olKDdU3MgyZqNatiOl2MDt+X2DyJCZQmJkFQGFmFhEJ8W6ZUCGNWqWwa33pddidlkHjVqE6N26Vwu4QmfRyMvUNt+7BmI7tiT+qN0d99zH9JrxHs1r+bQVjW+RlEJELgXnAsUBjIAYYBPzmpFWWb5SIzBWRuRPzt+5NxeXPqe65TNlifV4imsUy+4SzWXbf4xz25jN7rtu+UOEXZ9XbpdXYtfaN//HjYSfw8zHDyU/P5OCHb98HJfeQCk0K1bdis0tlPI0bc/Arz7Fq9CP4d+6qWf32lnCuw15cqzrHpXtQvF4i4mKZefzZ/HnP4/R8+5m913EfaciO3K1RK3cBvVR1a/BJEWkO/AK8W1EmVX0VeBVgUnyXPb6z8zem06hVaslxdMsU8tMzK5TZWiKTWk6mLHkbM8j4egoA2+b9AYEAEQnNKczZsqcq7hV5G9OJDrKrUcsU8jdlVikTjl0FWTkl++vf/YTeY1+uIY2rJ39TBlEtW5QcR7ZIJT8jM2wZ8fno+upzZH35FTmTppTIFGTnlIRfIpKTKMyp3Vf13LR0YtqUXofGrVPI3Rhq1+60dBqHyKSWk6lvuHUP5m3MIP2r0t+WBgJEJjSnoJZ+W8F46l/EJGzcUl0o97gGIEDFbbEaYdu8P2h8QHsatW2NRESQevowMidNDZHJ/GYqLUeOAKBZ7x4Ubt9BfkZWleVmTviO+AFHAND4wPZIZEStOXEwdsUc2J5G7YxdLc4YRsY3oXZlfDOVVo5dcb17UBSGXcXxS4CUU45nx9LlNa57ZexY8AfR7dsR1aYVEhFB0j+GsnlKqE05U6aSfMZwAJr27IF/x46SsEmnJx5i9/KVbHjt7ZA8m6dMJeXMEQCknDmCnMnfu25LMDlz/qBpp/bEtG+NJyKCdiOHsWF8qF0bxk+lw4VGx4R+PSjctoO89KqvVV3j1j2YMeE7EpzfVsyB7fFERNSJEwfwioS9VYeInCwiy0RkhYiUe9UVkYNEZJaI5IvIfypI94rI7yLydTi6u9UifxiYJyKTgfXOubbACcCDLtWJ+v0sufUBen/6OuL1kvbBZ+z8cwVtLh4JwPq3x5I15ScSTxjIgN+m4M/N5Y9r7yzJ3+O1J2l+dF8iE5pz7KKfWP7o82x4/1PSPviMQ58fw9EzviJQUMgf/6rFEIRj1+JbHqDvZ6+D10va+8autpcYu9a9NZasyT+RfMJABv4+hcDuXBZeU2rXYa8/Sfwxxq5Bi41dae99ykEP3ELsIQehQO66DSy68d7aM8rvZ+U9D3LI+28gXg8ZH33G7r9WkHr+OQCkv/8RW6b+RPzgAfT+eTKB3Dz+utnYFNvncFLOHMGupcvoOekLANY89jRbfpjG+hdf4+CXniZ15Bnkb9jE0qtvrD2bMNdq7rUPMOhbcw+uevMzti1ZQccrzbVa8cpYNk78iZZDB3Lqiin4d+cy+5LSa3XUh0+ScmxfohKbM2L9Tyy873lWvfkprUccT+/n7yEqKZ6BE15h6/ylFY5ucdMuN+7B9e9/RvcXxtB/5lcECgtZWMu/rWAia6gTU0S8wIsYf5cGzBGR8aq6JEhsM3A9MKKSYm4AlgJhdVyJW7E5J4xyEtAK0wpPA75V1bAet3sTWqnvBAJ1rYE7NG1a/2KG+8r6tP3u9gMgrllda+AOQ7cu2+eb8KC3zwz7ov958aeV1iciRwKjVfUk5/gOAFV9pALZ0cBOVf1v0LnWwDuYBvG/VfWU6vRx7ctOx2GPdRSLBTq5VZfFYrHsK+GETIoRkVHAqKBTrzp9fGAar+uD0tKAfnugyjPArUDTcDO4NWrlfRFJdPZPAhYDjwHzReQsN+q0WCyWfWFPYuSq+qqq9g7aXg0qKowhPhUjIqcAmar6257o7laLvIeqZjv79wH9VXWN49y/Bz5xqV6LxWLZK7w1N6wwDWgTdNwa2Bhm3qOBf4jIUCAaiBWR91X1/KoyuTVqxeOEU8CMVFkH4Dh3O1GXxWKpd3gl/K0a5gCdRKSDiEQCI4Hx4eigqneoamtVbe/km1qdEwf3nOr9wA8i8iIwA/hERMYBg4FJLtVpsVgse01NjVpR1SIRuRb4FvACb6rqYhG5ykl/WURSgbmYUSkBEbkR6Kqq2/emTlccuap+LCLzgCuAzk49RwL/U9Vv3ajTYrFY9oUaDK2gqhOBiWXOvRy0n44JuVRVxo/Aj+HU54ojF5Fzgcmqepsb5VssFktNE0bIpN7iVmilHSacEoHp3PwG+FXr/YQSFovl70pNtshrG1c6O1X1UVUdDAwFFgCXYr70/FBELhSR+j3Vm8Vi+dtRk5/o1zaujiBR1R3AF86GiHQFhmAmzTrJzbotFotlT6iPDjpcXB8KKCLNiz/LV9UlIrJJVZ90u16LxWLZEyIbcJC8NiZuLDv9XO1OR2exWCxh4PVI2Ft9ozY+zilrdf37K1gslr89NrRSAUErATUv3lfVdwlzzgGLxWKpTerhUpxh42aLvIPzf5SzX+zAG+5jz2Kx7LfYFnkFqOr9ACIyvHi/OMmtOi0Wi2VvqalP9OuCehsjT0houE/Hyuj9/oi6VsEV5KAeda1CjbPxwqfrWgVXiOsSX9cq1FsasB+vFUd+QTXHFovFUuc05NCKWwtLDBeRawBUdZGI/CIiq0RkFXCQG3VaLBbLvlCD09jWOm61yG/FzKVbTBTQB4gB3gI+dalei8Vi2Ss8DbhF7pYjj1TV4DXrflbVHCBHRGJcqtNisVj2mvrY0g4Xtxx58+ADVb026DDJpTotFotlr4lowJ2dbqn+i4hcUfakiFwJ/OpSnRaLxbLXeDwS9lbfcKtFfhPwpYj8E5jnnOuFiZWPcKlOi8Vi2WtsaKUMqpoJHCUig4FuzukJqjrVjfosFotlX6mHDe2wcXs+8qmAdd4Wi6XeY1vkFovF0sBpyMMPG3A/rcVisdQcEZ7wt+oQkZNFZJmIrBCR2ytIP0hEZolIvoj8J+h8GxH5QUSWishiEbkhHN1ti9xisViouU/0RcQLvAicAKQBc0RkvKouCRLbDFxP+cEfRcDNqjpPRJoCv4nIlDJ5y2Fb5BaLxYLp7Ax3q4a+wApVXaWqBcBYYHiwgKpmquocoLDM+U2qOs/Z3wEsBVpVq3vYVlosFst+zJ7MtSIio0RkbtA2KqioVkDwl+1phOGMyyIi7YGewC/VydrQisVisQCePWjWquqrwKuVJFfUZt+jdRhEpAnwGXCjqm6vTt46covFYgEiam4geRrQJui4NbAx3MwiEoFx4h+o6ufh5LGO3GKxWKjRceRzgE4i0gHYgJkJ9p/hZBQRAd4AlqrqU+FWuN858tj+x9D27rsQr4esjz8l/dXXysm0vecumg0cQCA3j9W33cHuJUuITE2lwxOPEZGUCIEAWR99TMY77wHQ+rZbiBs0CC0sJH/dOlbffif+HTtq1a7pCzN4+P2FBALKmQPbMerULiHpqsrD7y9k2oIMoqO8PHJFL7q1jwNg+64C7n7zd5anbUeAhy8/nJ6dEnj20yV8//smPCLEx0bxyBWHk9K8Ua3ZNG3Wch5+chKBQICzhh/OqIv6l7fpyW/4aeZyoqMjePTeEXQ7qCWr1mZz052flMit37iF60cN4uJzj+TGOz9h9dpsAHbszKNpk2jGfXB1rdkEENXzCJpdcTN4POyeMo6dn71bTib2ipuJ7nUUmp/H1mcfoHDVMoiIJHHMK0hEJHi95M38nh3/M/dv039eSXS/ARBQ/Ns2s/W5Bwhszq5VuzydexE5/GoQD0W/TqLox4/LyUT842q8B/WBwnzyP34S3bCiNFE8RF//HLo9h/y37jPywy7He3A/8BcRyNlIwcdPQd6u2jIphJpqkKtqkYhcC3wLeIE3VXWxiFzlpL8sIqnAXCAWCIjIjUBXoDtm8Z0/RGS+U+Sdqjqxqjr3L0fu8dBu9L38dfGlFKRn0PWzT9g6dSp5K1aWiDQbOICodu344/iTiDmsB+0euI+lZ56D+v2sf+Qxdi9Zgicmhm5ffMa2GTPJW7GS7TNmkvbfp8Dvp/UtN9PiqlGkPfFkrZnlDygPvLuAN289mpT4Rpx13w8MPrwFHVvFlshMW5jB2oxdfPvECSxYuYX7357Px6OPBeDh9xfS/9AUnruuHwVFAfLyiwC4bFgnbjizKwDvTl7J/335J/df0rN2bPIHeODxibz1wgWkJMdy5kWvMbh/FzoekFxq08zlrFm/mcmfXc+CRWmMfmwCn7x1BQe0Syxxzn5/gAHDnuSEYw8G4JkxZ5Xkf/SZb2nSJKpW7CnB46HZlbeSc9+1+HMySfrvO+T9Op2i9atLRKJ6HYWvRRsyrzqDiM6H0Ozq28i+5VIoLCDnnn+hebng9ZL46Gvk/TaLwr8WsfOL99nx4SsAxJxyNk3PuZxtLz1ae3aJh8jTriH/tTvRbdlEX/cc/iWz0cx1paYf1AdPYkvyHr8UT9uDiDztWvJfuLEk3XfMCAKZ65HoxiXn/H/No/CbNyEQIGLIpUQMOscc1wE1uUKQ43gnljn3ctB+OibkUpaf2YsF6verUSsx3buTv3Yd+evT0MJCNk+YSPPjjguRiTv+OHK+HAfArvkL8DaNJSIpicKsLHYvMUM1A7t2kbtyJZEpKQBs/3kG+P0A7Jy/gMjU1Fq0Chau3Ezb5BjaJMcQ6fMw9IjWfD9vU4jM9/M2MfzoNogIh3WMZ/vuQjK35rEzt5C5y3I4c2A7ACJ9HmJjIgFo0iiiJH9ufhFSi1+2LVy8gXat42nTKp7ICB/DTjyE76ctC7Vp2jJGDO1hbDq0Ddt35JGZHfomNGvOKtq0jqdVi7iQ86rKN98t5pQTD3XblBAiOnWjKD0Nf8ZGKCoid/pkovsOCJGJ7juA3B/Mb7zwr0V4YpriaZ5g9M7LBUC8PvD6KO4j09zSVqpENQKt3TXMPW26oNmb0M3p4C+iaMFPeLsdGSLj7XokRfO+ByCw7k+kURNoatYIlWaJeA/qQ9Gvk0LyBJbPg0CgNE9cYi1YUzE1OPyw1tmvWuSRqSkUbCp1cAXp6cT0CF0YODIlVKYwPZ2IlBQKs7JKZVq1onHXg9m5YEG5OpLOPIPNE6p8y6lxMrbk0SKhNOSRGt+IBSu3hMpszqVFfKhMxuZcfF4TNrnjtXksW7eNbh3iuPP87jSOMpf+6U8WM27Gepo28vHOHaGhDTfJyNpOakrpG0VKciwLF6eFymSGyqQmx5KRuZ3kxKYl5yZMWcQpJx5Srvy5v68lIT6G9m0TXNC+crwJSfizM0qO/TmZRHbuVkYmOVQmOxNvQjKBLTng8ZD05Lt4W7Rm18RPKfxrcYlc0/OvpvGgoQR27STn7toNF0mzBHRb6W9Et2XjaRMa3vM0S8C/NUhmaxaeZgkEdmwm4tQrKZj4BhLVmMrw9TmRogXTal75MGnIc624tWaniEg/ETldRE5z9uvmz1S25VLhwKBSGU/jxnR84TnWP/wIgZ2hsboWV1+JFhWRM/4rFxTdM8L5Y4pAkV9ZsmYr5x7XgS8eGkyjKB+vffVXicxNZ3Xjx2dO5pSj2vD+d6vcU7gMFTUopYxVFbU5g2+jgsIipk5bxsnHdSsn9/XkRZxyUu22xg3lr0w5W6u6BwMBsm46n4zLTiGyc1d8bQ8oEdnx/ktkXHYquT9NImbYWRUU4ibhjKiryHbFc3BfdOfW0Hh5GXyDR6IBP/7f626OvQiPhL3VN2rckYvIicByYDQwFBgG3A8sd9KqylsyyP6LbVv3uO6C9AwiW7QoOY5MTaUwM7NKmYggGfH56PjCc+SM/4otk6eE5Es4bQRxgwax6uZb9livfSWleTSbcnJLjtM355LcPDpUJr4RmzaXlWlEanwjUuIb0eNA84p7Up+WLFm7tVwdpxzZhilzNrhjQAWkJseSnlE6PDYjczvJSU2rlEkvIzNt5gq6HdSCxIQmIfmKivxM+XEpQ48v7+Ddxp+TiTcxpeTYm5BMYHNWqEx2GZnEZPxlZHTXTvL/mEfU4aHhC4Dcad8SfeTgGta8anRbNtKsdHEvaZaIbt8cIhPYlo3EBcnEJaHbN+Nt1w1v1yOIvv0dos67Hc+BPYgceWuJnLfX8XgP7kfB/x5335Aq8IiEvdU33GiRPwscr6pDVPVyZzsZM+/As1VlVNVXVbW3qvY+rVncHle8648/iGrfjsjWrZCICOKHDWXL96FP+K3fTyVhhPlaNuawHvh37CgJq7Qf8xC5K1eS8dbbIXli+x9Di1GXs/yqqwnk5e2xXvvKoQc0Z23GTtKydlFQFGDi7DQG92wRIjO4ZwvGzViPqjJ/xWaaNo4gOS6apLhoWsQ3YtUmE1uetTiLA1saZ7gmfWdJ/qnzNtGhZagjddWmri1Zsz6H9Ru2UFBYxITJixjcP/RVfXD/Lnw5cYGx6Y/1NG0SFRpWmfwHwyqIgc+cs4oD2iWSmtLMdTvKUrh8Cb4WbfAmtwSfj0b9TyTv1+khMnm/TqfRoKEARHQ+hMCunQS25OCJjUNinIdSZBRRPfpSlLYWAG+L0mHJ0X0HULRhTa3YU0wgbRmS2BJpngJeH74eA/EvmR0i418yG9/hpk/K0/YgE9ffsZnCSW+RN+YC8h69iPwPHiWwcgEFY43T9nTuRcSxZ5H/9mgozK9Vm8rSkB25GzFyH2ZAfFk2ABEVnK85/H7W3f8gXd58A7wesj/9jLwVK0g69xwAsv73Edt+/IlmAwdw6PeTzfDD2+8EoEmvw0k8bQS7/1xGt/FfAJD25NNs+2ka7e67B09kJF3eNr3pO+cvYO29o101JRif18M9F/bgssdnEFA4Y0A7OrWOZexUMxJi5OAODOyRwrQF6Zx4yxSiI72Mufzwkvx3X9CdW16aS6E/QJukGMZcYdKe/HgxazbtQDxCy4TG3H/xYbVnk8/LvbcM5fLr38MfUM44tSedDkzmf5/NAeDcM/ow8OhO/DRzOSec/hyNoiMYc0/pdBW5eQXM/GUVD9xxarmyJ05exLAK4ua1QsDPtlefIGH0c2b44fdfUbR+FY1PPh2A3ZM+J/+3GUT3Porklz83ww+ffxAAT/NEmt94n/nEUDzkzviO/Lk/AxB74TX4WrUDDeDPTGdrbY5YAQgEKBj3f0Rd/jB4PBTNmYxmrMV3hHkgFc2eSODPXwkc1Ifo296EgnwKPql+GHTkiGvAF0H0FWMA8K/7k8LPn3fVlMqojw46XERruPdbRO4AzsZMFFM830Bb4BzgY1V9JJxy5nQ6qHa75WuB3u+PqGsVXEEO6lG9UANj44VP17UKrhDXJb6uVXCFxo9P2mcv/PnKG8P2Oacf+Ey98vo13iJX1UdE5EvMbF9HYnpA0oDzqpuK0WKxWOqKhtwid2vNzqWY6RcBEJEEVc1xoy6LxWKpCXx7MmtWPcONUSuPikiis99LRFYBs0VkrYgMrOn6LBaLpSbw7MG/+oYbGg1T1eJJIP4LnKOqnTCjVmrvu3aLxWLZA+yolVAiRMSnqkVAI2cVDFT1LxGp5YkvLBaLJTzqo4MOFzcc+YvARBF5FJgkIs8AnwPHAfNdqM9isVj2GY/Uv5BJuLgxauV5EVkEXAV0duroDHwJPFTT9VksFktN0JA7O2vckYvIkcCPqvpDTZdtsVgsbtGQQytuPIIuAn4TkbEicrEzgbrFYrHUazziCXurb7gRWrkKQEQOAoYAb4tIM+AHYBIwQ1X9NV2vxWKx7AuePV/Pod7g2nzkqvon8CfwtIg0AgYBZwFPAb3dqtdisVj2hoYcWqmVhSVUNRez7FHtrshgsVgsYVIfQybh4qrmIjKrqmOLxWKpL/g8nrC3+obbGkVXc2yxWCz1gprs7BSRk0VkmYisEJHbK0g/SERmiUi+iPxnT/JWhCuhFREpXm02pnhfVadR8epdFovFUufUVIxcRLyYDyNPwMz8OkdExpeZ/XUzcD0wYi/ylsOtGPklzv8Jzr4C0whvqUmLxWKpdWpw1EpfYIWqrgIQkbGYab1LnLGqZgKZIjJsT/NWhFvT2F7iKDGveL84yY36LBaLZV/Zkxa5iIwCRgWdelVVX3X2W1G6qA6YlnW/MIveq7xuj1op+5cJ+y8Vn+LuqnB1gfQ4uq5VcAW56eG6VqHG2RRR/zq0aoLoo1rVtQr1Fp/HG7as47RfrSS5Ij8XbiN2r/K67chvq+bYYrFY6gVSc8MP04A2QcetgY1u5nW12aGqk6s6tlgslvpCDS4sMQfoJCIdRCQSGAmMD1ONvcrrxqRZxwAHqOq7zvGnQPGKrw+p6tSartNisVj2lZpqkatqkYhcC3wLeIE3VXWxiFzlpL/szEE1F4gFAiJyI9BVVbdXlLe6Ot0IrdwPXBd03AW4GIgB7gSsI7dYLPWOmvyyU1XLfcmuqi8H7adjwiZh5a0ON0IrsWXGPC5X1d+cceRNXajPYrFY9hnBE/ZW33CjRR4XfKCqpwcdprhQn8VisewzXqmVqadcwY1Hy58VDHJHRE4BlrlQn8Visewz++V85CLiUdXAXpR5EzBBRM4E5jnnegFHAafsRXkWi8XiOjU4/LDWqUrzec6ybXuEqq4AugO/A+2dbZpzruWeq2ixWCzuU4PDD2udqjS6EnhWRF4TkeZ7Uqiq5mM+X80E/gOMBZ4AHtlbRS0Wi8VNRDxhb/WNSkMrqvqLiPQDrgLmisg3QCAo/fpqyu4HPAbMxIxW+QDYP79Rt1gsDZ76GPsOl+q6aeOBPkAW8BtBjjwMCoFcoBFmHvLVexlzt1gsFtfxSsOd36mqzs6rgFswIZHLVHVPZy6cA4zDPAgSgFdE5ExVPXNvlbVYLBa3qI/jw8OlqhZ5f+BIZ97cveEyVZ3r7KcDw0Xkgr0sy2KxWFxlvwytqOp5+1JwkBMPPvfevpRpsVgsblEfOzHDpeF+ymSxWCw1SH0cVhgu1pFbLBYL4PU0XHdYpeZi3jUWquohtaTPPtOo39Ek3nAb4vGy/evP2fr+G+VkEm64nZgj+xPIyyNzzN0U/LUUb3IKKXePwRufCBpg+/hP2fbJBwBEduxC0i33IJFRqN9P9pMPkb90Ua3aNW3Gnzz82DgCgQBnndaPUZcNDklXVR5+bBw//byU6OhIHn3wHLodbCZXe/u9aXzy+S+IQOdOLXjkgXOIiorgm8kLeOGlyaxcncknH1zPod3aVFS1a5zU9QiePfsmvOLh9RnjeWxyaOStS0o73rrwbg5v04W7xr/Mk999WJK2+qEv2JG3C38gQFHAT59HzYqCPVp34uV/3ka0L5KigJ9//e8J5qytcrnDGifysCOIveQm8HjI/X48u74sH1Fseum/iep5JFqQz7YXHqRo9TKIiCT+gZeQiEjwesmfNZWdH78ekq/xP/5J7IXXk3HJSeiObbVlEgDTl21hzPhVBFQ5s08KVwwKvV9UlTHjVzFt2RaiIzyMObsz3Vo1AeC4R+cQE+XFK4LXI3x6/WEl+d6fsZEPZm7C6xEGHtycW4Z2qE2zSthfOztR1YCILBCRtqq6rraU2ms8HpL+fRcbbxpFUWY6rV8fy66ff6BwzaoSkcZH9CeyTTvWjRxGVLfuJP3nbjaMOg/8frJf+C8Ffy1FGjWm9ZsfsXvOLArXrCLhX/9my1svs3v2zzQ+oj8J//o3G6+7tNbM8vsDPDDmC956ZRQpKc0485/PMvjYrnQ8MLVEZtrPf7JmXRaTv7qdBX+sY/RDn/HJBzeQkbGNdz+czsQvbiU6OoIbbnmXCZPmc/rwPnTumMrzT1/EfQ9+Wmu2FOMRDy+O/A8nPHc9aVsymXP7W4xfOJ2l6WtKZDbv3s71Hz/FiB4DKyxj0NPXkLMr1Jk9ftq13D/hDSYtnsWQbkfy+OnXMujpf7lpSigeD7GX/4ctD1yPf3MmCY++Rd7c6fjT1pSIRPY8El+LNmRfdxYRnboRO+pWNt9xGRQWsOX+a9G8XPB6iX/oVSJ+n0XhcjMdtSchmajuffFnbao9exz8AeXBL1fyxuWHkNIskrNfmM+grgl0TGlcIjNt2RbWZucx6ZZeLFi3gwe+WMFH1x5Wkv7OqENpHhM6xO+XlVv5fkkO427qSaTPQ87OgtoyqRwNubMzHM1bAItF5HsRGV+8ua3Y3hB18KEUpq2jaGMaFBWx87tviDlmUIhM4/6D2DHJqJ+/eCGeJk3xJiTiz8mm4K+lAGjubgrXrMaX6EzWqIo0jgHA06QJRdlZtWcUsHDROtq1SaBN6wQiI3wMO/kwvv8xdK75739YzIhTeyMiHNa9Hdt35JGZtR0wD4K8/EKKivzk5RaSnBQLwIEHpHBA++RataWYvu27siIrjdXZGyn0FzF27hSG9xgQIpO1Ywtz1y6l0F8UdrmKEhttrlWzRk3YuK12r1VEx67409PwZ26EoiLyZkwhuk+oXdF9BpD7o5luunD5YjyNm+CJSzD65+UaIa8P8Ya2s2IvvpEd771QJ0uYL1y/g7YJ0bRJiCbS52FojySmLskJkZm6eDPDeyWbe7BdLNtz/WRur9oxj52VzhXHtiHSZ1xRQpNI12yojv19Gtv796ZgETkIGI5ZFVox686NV9Wle1NeOPiSkinKTC85LsrKILpr91CZxDIymRn4EpPx52SXyqS2JLLzQeQtWQhA9nOP0eKpV0i85j/gETZcVbujKDMyt5GaGldynJIcx8I/1paXSSmVSU1pRkbmNg7t1oZLLzqWQSc9RFR0BEcf2ZljjupSS5pXTqu4JNZvKR3ZmrYlk34duoWdX1WZfP1zKMor07/gtZ/HAXDjJ8/w7XXP8N/Tr8PjEY56YlQ1JdUsnvgk/NmldvlzMonoFGqXJyEJf06QzOZMPAlJBLbmgMdDwmNv401tze5vPytpjUf17o9/cxZFa1fUjiFlyNxWQGpcVMlxSrMoFq7bESKTsT2f1Galjji1WSSZ2/NJjo1EgMteX4QInNOvBWf3M2+Ta7Jz+W31Np79dg2RPg+3DuvAoW3qZtmC/bpFrqo/VbRVlUdEbsPMryLAr5iPgwT4n4jcXkW+USIyV0Tmjk3fvGeWmAIqMqBamWAJadSI1IefJufZx9DduwCIHXEOOc89ztozTiD7+SdIuuOBPddtH6joUywpY4dW0EwTEbZt3833Pyzi+4l3Mn3KveTmFjDu69/cUjVsyuoPFdtZGUf/dxS9HrmIIS/cxDUDz6R/x8MAuHrA6dz06bO0vWs4N33yLG9ccFcNaRwmFd6D5YQqkHGEAgFybrmQrCv/QUTHrvjaHACRUcSccTE7P6ps0Xb3qejSlDW1QhnH1g//1Z3Pb+jJq5d248NZG5mzyoTEigLK9twixl7Tg1uGdeCmD/5kz789rBka8lwr1WokIkeIyBwR2SkiBSLiF5Ht1WS7DOijqo+q6vvO9ijQ10mrEFV9VVV7q2rvkanxlYlVSlFmBr7k0rixLymFouzQ75mKssrIJKeUtqC8PlIfepodkyewa9r3JTJNh/yDXT99B8Cuqd8SfXDt9v2mpjQjPX1ryXFG5laSk2NDZZLjSM8olUnP2EZyUiwzZy+ndasE4uObEBHh5cTjDuX3BWtqR/EqSNuSSZvmpWGd1s2T9ygMsmmbeYPK2rGFL+b/RN/2XQG46IihfP77DwB8Mu97+rbrWoNaV08gJxNvYqld3oRkAluyysskBMnEJxPYnB0io7t3UrB4HpE9j8CX2hpvcgsS//s+Sf/3BZ6EJBIffwdP3J7/RvaWlGaRpG/NLznO2GZa2sGkxkaRvq00lJK+rYAkRyY51rTmE5pEcny3BP5Yb1rzqc0iOeGQBESE7m2a4hFhy67wQ2k1iQQCYW/1jXAeLS8A5wLLMfOmXO6cq4oAFU9Z24I9m69lj8j/cxERbdrha9EKfD6aHD+EXTN+DJHZ/fMPND35HwBEdetOYOfOkrBK8h33U7B2Fds+ejckjz87i+ievQFo1KsfhWm12+97aLc2rFmXzfq0HAoKi5gwaT6DB4a+rg8+titffjUXVWX+wrU0bRJNclIsLVPjWLBwLbm5Bagqs35ZzoEd6n6hpjlrl9IpuQ3tE1oQ4fUxsvcJjF84Pay8jSOjaRLVuGT/xIP7smij6dDeuDWbgZ0OB2Bwl94sz1rvjgGVULhiKd4WbfAmtwCfj+ijTyB/Tqhd+XOn0+jYoQBEdOpGYPdOAltzkNg4pLEZ5UFkFFHd++DfsJaidSvJumwoWf86jax/nUYgJ4vsWy8isHUv3lr3kkNbN2VtTi5pm/MoKAowcUEWgw4OfZAM6hrPuN8yzT24djtNo70kx0ayu8DPrnzjnHcX+Jnx11Y6pZp+jOO6JTB7pWmdr87KpdAfoHlMHQ0D1ED4Wz0jrL+Yqq4QEa+q+oG3RGRmNVluBL4XkeVA8S+pLdARuHZvla0Wv5/sp8bQ4qmXzfDDCV9QuHolscPPAmD7uE/YPWs6jY8cQNuPJhLIyyNrzN0ARHfvSdOT/0H+ir9o/dYnAGx+5Tl2z55O1uOjSbzhdvB60YJ8Mh/fq26Dvcbn83LvHadx+dWv4Q8oZ4zoQ6eOqfzvY3MZzj37KAb2P5iffv6TE055lEbREYx54BwAenRvx0kndOe0kU/j83o4+KBWnHPmEQBM+f4PHnz0SzZv2cmV177BwV1a8sbLtRNT9gf8XDv2v3x73bN4PR7enPk1Szat5sr+pwHwyvQvSImNZ+7tbxMbHUNAA9w4eCRdHxhJYpM4vrjyMfO38Xj5cM5kvl0yG4ArPniEZ8++CZ/HS15hAaM+qOWZkwN+tr/+X5rf/awZfjj1a4rSVtPoRGNX7uQvyJ83k8jDjyLxhU/R/Dy2/d9DAHibJ9Ls2nvA4wUR8mZ+T/5vM2pX/0rweYW7hx/I5W8sIhCA0/uk0Ck1hrGzzQiakUe0YOBBzZm2bAsnPf4b0ZEexpzVCYCcHYVc954ZAlrkh1N6JtG/i5kZ+/TeKdz96XJOfWoeEV7hkbM7Vxh2qxXqoYMOF6kuHiUi04Djgdcxc6ZsAi5W1R7V5PNgQimtMEHBNGCO8zColpXHHFo3gTIXOfC7MXWtgivITQ/XtQo1zqbaHexSaySfH36HckPCM+KNfff+hd+E73MihlRZn4icDDwLeIHXndBycLo46UOB3RifOs9JuwkT+VDgD+ASVc2rqr5wQisXOHLXAruANsAZ1WVS1YCqzlbVz4D5GIde98MlLBaLpSICgfC3KhARL/AiMAToCpwrImU7a4YAnZxtFPCSk7cVcD3Q2/kQ0wuMrE71cEatrMW0qFuo6v2q+m9nObeqDPlBRBKd/QuAiY7iH4nIddXVabFYLLVOzcXI+wIrVHWVqhZgRvANLyMzHHhXDbOBOBFp4aT5gEYi4gMaY4ZuV0k4o1ZOxbSoJznHh4XxQVCSqhZ3w1+PmQ73csyqQVdUV6fFYrHUOoGisLfgodLOFty51IrSvkEwYeVWZWqrUEZVNwD/BdZhwtjbVHVydaqHE1oZjXnCbAVQ1fmYBZWrotB5RQDYiQnJAORjXhUsFoulfrEHoZXgodLOFjzIv6L4edn4e4UyzvrIw4EOmJF/MSJyfnWqhzNqpUhVt+1hT/JNwGQR+QxYDEwVkUmYxSre2pOCLBaLpVaouVEraZi+xGJaUz48UpnM8ZhlMbMARORz4Cjg/aoqrLRFLiITRaQDsEhE/gl4RaSTiDyPWVC5UlT1R6fyTZi1O3/DtMavU9X/VpXXYrFY6oSai5HPATqJSAcRicR0VpYNR48HLhTDEZgQyiZMSOUIEWnsjGw5Dqh2WpOqWuRvA98C7wGHYBzxh865B6sqVETuACap6kvVKWCxWCz1ghpqkatqkYhci/GVXuBNVV3srIOMqr6MGQAyFFiBGX54iZP2i4h8CswDioDfgWrnZqhqqbePRWQCcC9wMsahF8d5rgGeqqLc1cANItIDWAB8A0xW1S3VKWSxWCx1gfoLw5atLtCsqhMxzjr43MtB+4rxoxXlvQ+4L2xlqD5GXojpqIwCmhDmBJqqOhYz5AYR6Yl5EHzujK/8DtNa/3VPFLVYLBZXacBfdlbqyJ0vk57CxHIOV9Xde1OBqv6OeT14RERigRMwXy1ZR26xWOoP+6MjB+4CzlLVxVXIVIuIdCsuQ1W3i8ifzteeFovFUn9owI680lErqtp/X524Q9kFC8svYGixWCx1zf4+++E+UrZfoI6mNrNYLJYqqIfzjIeLa45cRO51dlOK91W1dpfWsVgslnAJ1M2CFjWBmy3y4kUlC4P2LRaLpX5SD0Mm4eKaI1fVdwBE5Ibi/eIkt+q0WCyWvcaGVqpkr2Lkq/8qqF6ogeHvf1tdq+AKGR2a1bUKNU7etirn8W+w5Iyt9mvvBknSiBooxLbIQxGRjkCKqs7AzBVQfH4AZhJ1i8ViqV80YEcezjS2e8MzwA4AVQ1eIXY35pN/i8ViqV/U0ApBdYFboZX2qrqw7ElVnSsi7V2q02KxWPaeorCWE66XuOXIo6tIa+RSnRaLxbL31MOWdri4FVqZIyLllnQTkcswc5NbLBZL/SKg4W/1DLda5DcCX4jIeZQ67t5AJHCaS3VaLBbL3tOAW+SuOHJVzQCOEpFBmEUpACao6lQ36rNYLJZ9xjryilHVH4Af3KzDYrFYagTb2WmxWCwNnHoY+w4X68gtFosFbGjFYrFYGjwN2JG7NfzQYrFYGhSqGvZWHSJysogsE5EVInJ7BekiIs856QtF5PCgtDgR+VRE/hSRpSJyZHX12Ra5xWKxQI21yJ1F5l/ErE+chvmuZryqLgkSGwJ0crZ+wEvO/wDPYhaoP1NEIoHG1dVpHbnFYrFATY5a6QusUNVVACIyFhgOBDvy4cC7apr3s51WeAtgFzAAuBhAVQuAaqeCtaEVi8VigT2aNEtERonI3KAteFbXVsD6oOM05xxhyBwAZAFvicjvIvK6iMRUp7p15BaLxQJ75MhV9VVV7R20vRpUUkVrLpQNrFcm4wMOB15S1Z6YFnq5GHtZrCO3WCwWqMm5VtKANkHHrYGNYcqkAWmq+otz/lOMY6+S/S5GnjCoP50fvgvxetjw/iesff61cjKdH76LxOMH4s/NY8l1t7PjDxO66vrMGBJPOJaC7BxmDzy1RP6QV58mpmMHAHyxTSnavoNfBo+oFXuKaXzkMSTffAd4vGwb9ylb3nm9nEzSzXcSc/QANC+X9PvvJH/ZUnwpqaSOfgRvQiKosu2Lj9k69v2SPHFnn0fc2f9E/X52/fwT2c8/WWs2RfY4giaX3AgeL3nfj2f3uPfKyTS55CYiex4F+Xls/78HKVr9F0RE0vz+l8AXgXi95M/+gV2fmL9HzFmXEX3ccALbtwCw638vU/D7rFqzCSC691E0/9ct4PGw65sv2f7RW+Vkmv/rVqL7Ho3m55HzxH0UrvgTIiJJeeoNJCISvF5yp3/HtndfBqDRgONpdsFVRLTtQMZ1F1Dw15JyZbpNRI9+NLnoRsTjJXfqV+SOL3+9Yi66iaieR6L5eex46SGK1pjrFXff/yEREeDxkv/LD+z+9A0AGp95GdGD/4EWX6+xr1Awv3avVwk1N/xwDtBJRDoAG4CRwD/LyIwHrnXi5/2Abaq6CUBE1otIF1VdhlmYp9qLvX85co+HLo/dy+9nXULexgz6Tv6U7G+nsuuvlSUiCccNoPEB7ZnZ70Rie/XgoMdHM2fI2QBsHPs56994n24vPBZS7KJRN5Xsd7r/Noq276wde4rxeEi+9W42XHs5hRkZtHvnI3ZN+4GC1aV2xRw1gMi27Vhz+slEH9Kd5NvvY/0lI9GiIrKeeZz8ZUuRxo1p9+6n7P5lFgWrV9KoV19iBg5m7bkj0MJCvM3ja88m8dD0spvZ8tANBHIyaf7Im+TPnY5/w5oSkcieR+JNbcPm68/C16kbTS+/lS13XQ6FBWy9/1o0Pxe8Xpo/8Ar582dRtHwxALsnjCX3qw9rz5ZgPB6aX3c7mbddjT87g9QXPmD3rJ8oWreqRCS67zH4WrVl08XDiTz4UOKvv5OM6y+EwgIybxmF5uWC10fK02+SO2cGBUv/oHDNSrLvv5n4G++uG7vEQ9NL/8PWh53rNeYNCn4rc70OOxJfi9ZsvvFsfB270eTyW9h69xXmej14HTjXK+7+lymYP5uiFeZ65U4cS+7X/6sbu4Kpoc5OVS0SkWuBbwEv8KaqLhaRq5z0l4GJwFBgBWbBnUuCirgO+MAZsbKqTFqF7FeOvNnh3cldvZbctWkAZHwxgaSTjwtx5ElDjmPTx18CsP23BfiaxRKZnERBZhZbZ88luk3ZPolQUv4xhN9Ov8g1GyoiutuhFK5fR+EGY9f2Kd8QM3BwqCMfOJjtE8YBkLdoId6mTfEmJOLPycafkw2A7t5NwZpV+JKSKVi9krgzRrLlndfRwkIA/Fs2U1v4OnalKD2NQKZ548yf+R1RfQawO8gxRPUeQN60bwAoWr4YiWmCJy6BwNYc48QBvD6zhTG2tzaI7HIIRRvX40/fAMDuH7+l8VHHsj3IkTc6ciC7vvsagIKlf+Bp0hRPfCKBzdnGiQPi8yG+UruK1q2uZUtC8XXsij/oeuXN/I7I3v3JDXbkvfuTN20SAEUrFiONS68XZa9XfVyDvQY/CFLViRhnHXzu5aB9Ba6pJO98zGyxYbNfxcijUlPI25Becpy3KYOoFinlZTaWyuRvTC8nUxlxR/SmICuH3NVra0bhMPElpVCUUapzUUY6EUnJZWSSKQyWyczAlxxql69FS6K6HEzeYrN4U0S79jQ6rBdt3hpL61feIarrIdQW3vgkAjmZJceBnEw88UkhMp74JALZGUEyWaUy4qH54++Q+PpECv74laIVpW+fjU86k/gn3qPp1XchMU3dNaQM3sRk/FmlOhdlZ+BNDLXLl5iMP7P0WvmzM/AlOtfT4yH15bG0+uR78ubNpuDPRbWid3V44pPw5wRdi81ZeCu4XmVlQq7Xo2+T+OoECv+YE3K9Gp10Js0fe5cmV95Z69crhAa81JsrjlxEfCJypYhMcr5aWiAi34jIVSIS4UadTsXlz5VpqUkYMpWRevoppH/x9d5otm9UoHM5lauxSxo1puVjz5L11CMEdu0y57xePE1jWX/JSLKf/S8txzxVk1pXTTjXocJ+fUdGA2y59SJyrhpOxIFd8bY5AIDdkz8n57oz2XzrhQS2ZNPkwutrVu/qCGu8QkXX0xEKBEi/aiQbzj2JyC6HENH+wBpXscYofxNWLqMBttx+MTn/GoHvwIPxtjbXK3fK52y+/iy23H4Rga05xJx/nbs6V0UDXljCrRb5e8BhwGhMHGgYcD/QA3i/skzBYzMn5G7d40rzN6UT3Sq15Di6RQr56ZkhMnmb0oluWSoT1TK1nEyFunm9JA07gYwvJ1YrW9MUZabjSynV2ZeSSlF2ZhmZDCKCZZJTKMpyZLw+Wj72DNsnfc3OH74LKXfnD1MAyFvyB6oBvHHNXbSkFH9OJp6E0rcKT0IygS3ZITKBnCw8iSlBMknlZHT3TgqWzCPysCPM8bYtZjV0VXK/H0fEgQe7aEV5/FmZeJNKdfYlpuDPyQqRKcrKwJtceq28Fcjorp3kLZhLdO+j3FU4TAKbs/AmBF2L+CT8Za/X5sxyMhVdr8IlvxN5mPmIMfh65U0dR0THri5aUQ22RV6Ow1X1alWdrappzjZbVa8GelaWKXhs5rBGcXtc6fbf/6DRAe2JbtsaiYgg5bRhZH0bupZF1qSptDh7BACxvXpQtH0HBZlZFZQWSvyAo9i9fBX5mzKqla1p8pYsIqJtO3wtW4EvgtgThrBrWug077umTSV22HAAog/pTmDnjpLYeOo9D1KwZhVbP3wnJM/OH6fSuI/5QUW0bYdERODfuqUWLIKilUvxtWiDJ6kFeH1EHXU8+XOnh8jkz51O9IAhAPg6dUN37yKwNQdpGoc0bmKEIqKIPLQP/g0m3OWJSyjJH9X3WIrWr6I2KVi2mIhWbfGmtgSfj8bHnkTurB9DZHJn/UTM8acAEHnwoQR27SSwORtPs+ZIjLFLIqOIPrwfhevX1Kr+lVG0cine1NYl1yv6qOMp+O3nEJmC334mesDJAPg6Vna9Iok8tDf+jRVcrz4Da/16BaN+DXurb7jV2blFRM4CPlPVAICIeICzANc8hfr9LLv9AXp+9Dri9bLxw8/YtWwFrS4aCcCGd8aS891PJB4/kKN+nUJgdy6Lb7izJP8hLz9J86P7EhHfnGPm/8Sqx59n44efApBy2lDSv5jglupV4/eT9fjDtH7uNfB62D7+CwpWraDZ6ecAsO3zj9g1YxoxRw+g/ReT0Lw80h+4C4DoHocTO2w4+cuX0faDzwHIefEZds2cxrbxn5N670O0GzsOLSwkffSdlapQ4wT87HjzSeLuegbxeMj94Wv8aauJPsGsBJg35QsKfp9J5OFHkfDcJ2hBPtv/7yEAPM0TiL3mXsTjARHyZk2lYN4MAJqcfw2+9p1BFX/WJna8+lilKrhl1+YXHiP5kf8zww+/HUfh2lU0OeVMAHZ+/Sl5v/5Mo37H0OKd8Wh+Hpv/OxoAb3wiCbc+AB4PiIfd06aQ94t5uDU6ehDNr7kNb7PmJD30HAUrl5F1R4V9Za7ZtfOtp2h259OIx0te8fU6fgQAed99aa7XYUcS/+wnZvjhyw8D5no1vfoec708HvJnfU/BvJkAxJx3Db52nUqu187XH689m8pSWP9a2uEi4czktceFirQHHgMGYxy3AHHAVOB2Va22C/675C7177G3j7Rt561rFVwhrkOzulahxsnbklfXKrhCo4RGda2CKySNnVlR78QeUfT+BWH7HN/57+1zfTWJW2t2rgHOARCRBMwDI7vKTBaLxVKX1MOQSbi44shFpC2Qqap5wGbgYme+3SXAa6pa5Ea9FovFstf4G25oxa3OzolBZT+KGbXyC9AHeLWyTBaLxVJXaEDD3uobbnV2elR1t7N/PNDH6fR8X0QWuFSnxWKx7D0NOLTiVot8vYgMdvbX4Mzy5cTLLRaLpd6hhYGwt/qGWy3yK4B3RGQ0sA2YLyK/A82Bf7tUp8Visew99fBDn3Bxa9TKOmCQiBwMdAbexlm7rnhcucVisdQrGnBoxa1RK98Ck4BvVHWcG3VYLBZLTVIfOzHDxa0Y+UWYD4FGi8g8EXlJRIaLSBOX6rNYLJZ9w6/hb/UMt0Ir6ZhwytvOp/n9gCHArSKSC0xW1Tr8FtdisVhCqY+dmOHi+sISTkx8lrPdKyKtgYFu12uxWCx7RAP+IMg1Ry4irYAWwEJVLRCRZOBG4GJVbelWvRaLxbI32Bh5GUTkRmA+8DwwW0QuApYCjYBebtRpsVgs+4SNkZdjFNBFVTc7866sAAao6myX6rNYLJZ9w7bIy5GnqpuhZEz5X9aJWyyW+kxNLiwhIieLyDIRWSEit1eQLiLynJO+0JlUMDjdKyK/i0hYa0u61SJvLSLPBR0nBx+rai0vpGixWCzVUOivkWJExAu8CJyA8yGkiIxX1SVBYkOATs7WD3jJ+b+YGzDh6Nhw6nTLkd9S5vg3l+qxWCyWGqEGOzv7AitUdRWAiIwFhmOm8S5mOPCumpV9ZotInIi0UNVNzsi+YcDDhDmliVvjyN+p6LyIRAOnhlPG/riaTutBbetaBVeIPr1HXatQ4xRMXlK9UAMk6qx+1Qv9Xam5TsxWwPqg4zRCW9uVybQCNgHPALcCTcOt0K0YeQlOrGeIiLwLrMVZOchisVjqE3syH7mIjBKRuUHbqKCiKloGruxTokIZETkFsyjPHkUx3BxHPgD4J+YV4VfgaKBD0DzlFovFUm8IpxOzRFb1VSpfJCcNZ+puh9bAxjBlzgT+ISJDgWggVkTeV9Xzq9LHrXHkaZiVgWYAXVX1DCDXOnGLxVJfqcEVguYAnUSkg4hEAiOB8WVkxgMXOqNXjgC2qeomVb1DVVuransn39TqnDi41yL/DBiBCaP4RWQc5V8tLBaLpd7gr6G5VlS1SESuBb4FvMCbqrpYRK5y0l/GLIc5FPONzW7gkn2p063OzhucrzsHAecCTwDNROQcYIKq7nSjXovFYtlbtAYXllDViRhnHXzu5aB9Ba6ppowfgR/Dqc+1GLmj6FRgqohEACdjnPqLQKJb9VosFsvesCcx8vqGWwtLDAdaq+qLzqmfgWRn/yY36rRYLJZ9wU6aVZ5bCQ3uRwG9MdPXXuxSnRaLxbLX1OQn+rWNW6GVSFUNHuz+s6rmADkiEuNSnRaLxbLX+IvsfORlaR58oKrXBh0muVSnxWKx7DX1saUdLm6FVn4RkSvKnhSRKzEfB1ksFku9ogbHkdc6brXIbwK+FJF/AvOcc70wsfIRLtVpsVgse019dNDh4tY48kzgKBEZDHRzTk9Q1alu1GexWCz7SkMOrbi6+LLjuK3ztlgs9Z6a/CCotnHVkVssFktDIVBDn+jXBdaRWywWCxCwMXKLxWJp2NgYeT2i8ZHHkHzzHeDxsm3cp2x55/VyMkk330nM0QPQvFzS77+T/GVL8aWkkjr6EbwJiaDKti8+ZuvY9wFoMeZJItp1AMDbpCn+nTtYd97ptWqXp3MvIodfDeKh6NdJFP34cTmZiH9cjfegPlCYT/7HT6IbVpQmiofo659Dt+eQ/9Z9Rn7Y5XgP7gf+IgI5Gyn4+CnI21VbJjF9YQZjPviDQEA5c2A7rjilc0i6qjLmgz+YtiCD6EgvY644nG7t4wDYvquAe96cz/IN2xHgocsPp2fHeJ4Yu4gf5qcT4fXQJjmGMZf3JDYmstZsAvAceDi+k64Ajwf/71Pwz/i0nIzvpFF4OvWCwnwKxz2Lpq8sTRQPkZc/he7YTOHYB4z8wHPx9jwJ3b0NgKKp7xJYUbsrKE6ft4GHX59rrtcJHRl1xiEh6arKw6/PYdpvG4mO8vLI9UfR7cAEALbvLODuF2exfN1WRODha4+i50FJ3PTENFZv2G5kdhUQGxPJl8+cUqt2lehvW+T1BI+H5FvvZsO1l1OYkUG7dz5i17QfKFhd+iOJOWoAkW3bseb0k4k+pDvJt9/H+ktGokVFZD3zOPnLliKNG9Pu3U/Z/cssClavZNOdN5fkT7zxVgI7d9SuXeIh8rRryH/tTnRbNtHXPYd/yWw0c12JiOegPngSW5L3+KV42h5E5GnXkv/CjSXpvmNGEMhcj0Q3Ljnn/2sehd+8CYEAEUMuJWLQOea4FvAHlAffXcAbtx5NSnwjzh79I4N6ptKxVelas9MWZrA2fSeTHj+eBSu38MA7C/jovoEAjPngD445NJlnr+tLQVGAvPwiAI7qlsxNZ3XF5/Xw348W8+rXy/nPOd0q1MEVxINvyFUUvn8Puj2HyMufIrDsFzS79ENnT8deSEJLCl64EmnVhYhhV1Pwxn9K0r39TkWz0yCqcUjRRb+Mwz/ri1ozJRi/P8ADr/zKm/cfT0pCY8665RsG921NxzZxJTLTftvI2k07+Pal4Sz4K5v7X/6Fj58YCsDDb8yh/+Etee62gRQU+snLNwsdP33LgJL8j745l6a1/NANpiG3yF1f6q02ie52KIXr11G4IQ2KCtk+5RtiBg4OkYkZOJjtE8YBkLdoId6mTfEmJOLPySZ/2VIAdPduCtaswpeUXK6OpsefxI5vJ5Y77yaeNl3Q7E3o5nTwF1G04Ce83Y4MkfF2PZKied8DEFj3J9KoCTSNB0CaJeI9qA9Fv04KyRNYPg+cnvrAuj+RuNqblHLhqi20TWlCm+QYIn0ehvZrzdR56SEyU+elM/zotogIh3WMZ/vuQjK35rEzt5C5y3I4c2A7ACJ9npJW99GHJuPzmtu6x4HNydiSW2s2AUirTuiWTejWDAgU4V88DU+X0OUaPV2OwL/ADObSDcsgKgaaOB9DN03A06kP/t8n16re1bFweQ5tWzSlTWpTIiO8DD2mHd//sj5E5vtf1zP82APM9eqSxPZdhWRu3s3O3QXMXZzBmcd3BCAywktsk1CHrapMmrGWYf3b15ZJ5bAfBJVBRE7CfPjTCrOgxEZgnKpOqirfvuJLSqEoo9QZFGWk0+iQ7mVkkikMlsnMwJecgj8nu1SmRUuiuhxM3uKFIXkb9eyFPyeHwvVrXbKgYqRZArotq+RYt2XjadMlRMbTLAH/1iCZrVl4miUQ2LGZiFOvpGDiG0iZFl4wvj4nUrRgWs0rXwmZW3JJjW9UcpwSH83ClVtCZDK25JKaUCqTGh9N5pZcvB4hvmkkd74+j2XrttO1fRx3nn8ojaNCb+fPp69lSN/W7hpSBmmagG4rvZd0ew6eVp3Ly2wPktmRY87t3ELESVdQ9N1bSGQjyuLrMwxv90EENq6gaMobtRoGy9i8mxaJpdMkpSbEsGB5djUyjcnYnIvPI8Q3i+aO52aybM0Wuh2YwJ2X96ZxdESJ7NwlmSTERdO+ZSx1RUMetVLjLXIReQa4AfgJeByzqMRPwPUi8mw1eUsWNP0oa0tVopUVUO6UllvytII1T4OEpFFjWj72LFlPPUJgV+gPpemJw9gxuXZb445WFZyrfi1XVcVzcF9059bQeHkZfINHogE//t9rb8h/uetC+UtTUbtHEPwBZcnabYwc3IHPHxxE4ygvr339V4jcy+OX4fV4OPWo2nXkYV2rikRQPJ36oLu2oZtWlkstmvsN+c+PouCVG2DnFnwnXFYTyoZPRdcrTJmigLJk5WbOHdKZL54+hUbRPl77bHGI3ITpaxjWv0NNabtXBAIa9lbfcKNFPlRVO5c9KSIfAX9hnHyFBC9o+lefrnv81yrKTMeXklpy7EtJpSg7s4xMBhEpqeQVyySnUJTlyHh9tHzsGbZP+pqdP3wXWrjXS5NBx7PuwrP2VK19RrdlI81K5xqTZono9s0hMoFt2UhckExcErp9M75D++PtegTeg/oiEREQ1ZjIkbdSMPZxALy9jsd7cD/yX729doxxSIlvRPrm0rBHxuY8kuNCW6GpzRuRnlMqk745j6Tm0QimBd/jQBM6OrFPS16bsLxE7suf1/Hj/HTeuu1opKIHt4vojmykWWmISmIT0B2h10q35yCxiSV+T5oaGe/BR+Pt0hdvp17gi4SoxkSM+DeFXz4Fu7aW5PfP+5aIc++tBWtKSUlozKbs0oZNes4ukuMbVSOzm+T4RogIKQmN6dHZ3J8nHdmW1z4vdeRF/gBTZq3jsyeHumxF1TTg74FciZHniUjfCs73gRL/6Qp5SxYR0bYdvpatwBdB7AlD2DXthxCZXdOmEjtsOADRh3QnsHNHSVgl9Z4HKViziq0fvlOu7MZ9j6Rg7WqKMjPcNKFCAmnLkMSWSPMU8Prw9RiIf8nsEBn/ktn4Dj8OAE/bg9DcXbBjM4WT3iJvzAXkPXoR+R88SmDlghIn7unci4hjzyL/7dFQmF+rNh3aIY61GTtJy9pFQVGAib+kMahnaojMoJ6pjJuxDlVl/orNNG3kIzkumqS4aFrEN2b1JtPpPHtJFh1bNgXMSJjXJyzn/248gkZRtd+XrxuWI/EtkbgU8PjwdhtA4K/QeeICf/2Ct4fpu5FWXSB/N+zcQtHUd8l/5hLyn7ucws8eJ7B6oXHiUBpDBzwHHYlm1m5479BOCazdtIO0jB0UFPqZ+PNaBvdtEyIzuG9rxv24ylyvZVk0jYkgOb4xSc0b0SIxhlUbzIibWQvTObBNs5J8sxZsokPrWFIT63aG60Ag/K2+4cadfjHwkog0BdKcc22A7bi9qITfT9bjD9P6udfA62H7+C8oWLWCZqefA8C2zz9i14xpxBw9gPZfTELz8kh/4C4AonscTuyw4eQvX0bbDz4HIOfFZ9g108SNm544pNY7OUsIBCgY939EXf4weDwUzZmMZqzFd4RpwRTNnkjgz18JHNSH6NvehIJ8Cj55qtpiI0dcA74Ioq8YA4B/3Z8Ufv68q6YU4/N6uPuC7lz+xEwCAeX0Ae3o1DqWsVNXAzBycAcG9khh2sIMTrplCtFRPsZc3rMk/13nH8otL/9GYVGANsmNefjywwF46L2FFBQFuOyJGQD0ODCe0RcfVis2AaABir55mYjz7gfx4J//HZq1Dm+vkwHw/zaJwPK5eDr2JvLaV83ww/FVRhwBiDj+EiSlA6Do1kwKJ7xYbZ6axOf1cM8Vfbns/u8J+JUzju9Ip7ZxjJ1kQlojT+7MwF6tmPbbBk686ktzva4/qiT/3Vf04ZanfjbXK6VJSNqE6Ws4pY7DKlA/HXS4iFYUrKyJgkVSMZ2dAqSpano1WULYm9BKfaf1oLZ1rYIrRJ/eo65VqHEKJi+paxVcIeqsftULNUDk4Lv3OYb2W5eDwvY5vZb9Wbsxu2pwo7MzUkREVdNV9TegKXCeiJxc03VZLBZLTVGToRUROVlElonIChEp1wElhuec9IUicrhzvo2I/CAiS0VksYhU2qcYjBsx8jlAnKPULcDDQCPgZhF5xIX6LBaLZZ+pKUcuIl7gRWAI0BU4V0S6lhEbAnRytlHAS875IuBmVT0YOAK4poK85XDDkXtVtXjs4DnAcar6EEbxYS7UZ7FYLPtMDbbI+wIrVHWVqhYAY4HhZWSGA++qYTYQJyItVHWTqs4DUNUdwFJMiLpK3HDk20WkeBKGbCDa2fe5VJ/FYrHsM3viyIO/eXG2UUFFtQKCP3tNo7wzrlZGRNoDPYFfqtPdjVErVwEfiMgCIBOYKyI/Ad2BMS7UZ7FYLPvMnoxaCf7mpQL27gu+IBkRaQJ8Btyoqtur06fGHbmqFgfuTwQ6AwswT5t/q+rWmq7PYrFYaoKiohorKg0z5LqY1phpSsKSEZEIjBP/QFU/D6fCGnfkIvIq8A3wnap+U9PlWywWixvU4FDsOUAnEekAbABGAv8sIzMeuFZExgL9gG2quknMp8hvAEtVtfqPQRzcCK28CZwM/FtECoDJwCRVXeBCXRaLxVIj1NQHQapaJCLXAt8CXuBNVV0sIlc56S8DE4GhwApgN3CJk/1o4ALgDxGZ75y7U1Wr/BrRjdDKbGA2MFpEEjAhlptFpDswD+PUy6+KYLFYLHVITX7Z6TjeiWXOvRy0r8A1FeT7mUqmVasKVyejUNUc4H/Ohoj0wTh2i8ViqVc05E/03ZqP3As0V9Vs5zgSuAjT4XmwG3VaLBbLvtCQHbkbn+iPBDYDC0XkJxEZBKzCxIPOq+n6LBaLpSYoKgp/q2+40SK/G+ilqiucYYizgJGqWjeLDVosFksYNOQWuRuOvEBVVwCo6jwRWW2duMViqe/Uw4V/wsYNR54sIv8OOm4SfLwnYyMtFoultrAt8lBew0xdW9FxA37mWSyW/RnryINQ1fsrS3OGH1osFku9oz52YoaLaysElVRg5tIdCZyL+Qy1t6sV7gUiMsqZBGe/Yn+0a3+0CfZPu/ZHm+orrjhyEWmHcdznYiZKbwf0VtU1NV5ZDSAic+vjA2Zf2R/t2h9tgv3Trv3RpvqKG+PIZ2I+TY0AzlTVXsCO+urELRaLpaHjxkIPWZjOzRQgyTlnOzktFovFJWrckavqcOBQzARZ94vIaqC5iPSt6bpqkP01jrc/2rU/2gT7p137o031ktro7EzBrN05Emijqm2qyWKxWCyWPcB1Rx5SmUg7VV1baxVaLBbL3wA3OjuPEZELg44/FZGpIjIVOLCm67NYLJa/O250dt4PzA067gLcAowGbnWhPsCEcETkQxFZJSK/icgsETltD/LfISLnOftni8gSEVksIh8GyVwkIsud7SI37KhCv04i8rWIrHTs+0FEBuxhGa+IyNHO/nUissyx8fEgmTtEZIWTdlJN2xFUT6qIjHXsWSIiE0Wkc5h5zxWRu5z9Y0VkvmPHT0EyJzs2rBCR2120o8r7TkSeFZENIuIJOnexiGQ5ev8pIjeVKbOFiEx29tuKyGQRWer8ndo75zuIyC/OvfiRM1W0G/ZVe9+JyDgRmVXm3GjH7vmO3ueWST9SRF5z9rs7f7fFIvKHiEQ753s5xytE5DkR2eMFF/42qGqNbsCcMsefB+3PqOn6nHIFM8viVUHn2gHXlZEbDVxcSRk/YEbZdAJ+x8ynDpDs/B+PmY43Hmju7Dd3w54KdIsG/gL+EXTukIpsAd4Gjq2knPmYpacGAd8BUWVs7IpZLDsK6ACsBLy1dL0OA/oHHV8MjK4k/ztALyAOWAK0LWOH19H9ACDSsalrbd93mIbSOsyKWceWse0FZz8ByMb0HxWnXwLc7Oz/CJzg7DcBGjv7H2NmFQV4Gbi6Lu475xqsB5YCHYLOjwb+4+x3ArYDEUHp9wNnYL4uXwj0CPp7eJ39X4Ejnb/zN8CQmrZxf9ncaJHHBR+o6ulBhyku1AcwGDPrYvBSSmtV9flwMotILBCpqlnAFcCLqrrFKSfTETsJmKKqm520KZi1SWsUEekjIgtFJFpEYkRkMXAlMEtVxxfLqeoiVX17D8o9GPhLVf3A1cCjqprvlFVs43BgrKrmq+pqzHqCbow2GgQUlrle81V1enUZnVbZYZhRUf/ENBTWOWUU29EXWKGqq1S1ABiLsa2mqe6+GwQsAl7CfBxXDjWraK0AWgSdPhn4RsxX0T5VneLI7lTV3c7fYDDwqSP/DjBiXwzZh/vuDOArzN94ZCU2LsesS9k86PRxmMbEicBCddb0VdUcVfWLSAsgVlVnqfHq7+6rjfszbjjyP0VkWNmTInIKsMyF+gC6YX7Y5RCRQ53Xu/nAVcADxcdi1hQFOB743tnvDHQWkRkiMltEip11K0zLo5g051yNoqpzMCtsPwQ8DryPaeVVaB+AiJwUZOM/gNed41+CxIYAk5z9zkB/59X8JymdA6dWbMS06n4re1JEEoLseAC4KuhaHeqI9QQWOD/uzpihrT86r/3FfTO1ZUel953DuZhlDr8AThGRiLICItIW0/Jd6Bx7gS6qugRj31YR+VxEfheRJ5z0BGCrqhbPDrLP9u3NfVfGxv9RycNKzLoEy4sftCKSiHmQb3NsVBH5VkTmiUhx+LWVY1cxbl3D/QI3Zj+8CZggImdSehP0Ao4CTnGhvnKIyIvAMZjWUh9MCw4RGQ2sqaAlezLwlrPvw7wKHgu0BqaLyCFUvCCqW0N+HgDmAHnA9cATwYki8oWj41+qerqqfotZsRsReRt4W1V/LFPmSZSu1O3DtI6OAPoAH4vIAdSujeUrMq3Tw8DEkYH2qjq6jNjJmNdsMHb0wrTuGgGzRGQ2dWRH8H2HWQ19KHCTqu5wHqonAhMc8XPErJ7VBbhCVfOc8/2A4gewD+iPeXitAz7ChGVKWshB1IR9e3TfiRla3BH4WVVVRIpE5BBVXeRkuUlErsCEuILfXk8EJjv7PszfrA+m1f69iPyGCcW4YeN+iRsfBK0AumPizO2dbZpzrmVN1+ewGDg8SIdrMD/upEpzhNIXE48D8+Qfp6qFTnhhGebmTQOCx8C3Bjbuo96VEY+JhzbFtNbK2nca5gcdH05hItIYiFPVYn3TMCEJVdVfgQCQSO3ZuBjjgPeGYCeQBkxS1V1q1oedBvSgdu2o7L47GWgG/CEiazDOKrjF+pGqdsM46idFJNU5H/zmlAb87oSIioAvnfqygTgRKW6I1ZR9e3rfnYNpEKx2bGxPaHjlaVXt4si9W9yJWYGNP6lqtqruxkzvcbhzvnVQWW7+3ho+bgXfMbHB2zCto8bA85h4mxt1CaYVc3XQubaY1nd1ebth4sLFxycD7zj7iZhX9ATMzbsac+M2d/bjXbJnPCb+exfwAqa1uYLQTqcBwI9hljcMExMvPr4KeMDZ7+zYKM7fIrizcxXudXb+gmmJFp/rAwysJl8zTOuv+PhgTEjM59xjizBhG5+jewdKOzu71eZ9hxNqCDofA2Q6el6M09nppD0LPOLsz8TEhsF02i4Akpzjt4BrnP1PCO3s/Fdt33eYjt4jg9I6YPomIKiz0zkeh4m5i2NT8TcszTFv7o2d6/YdMMxJm4N5ayzu7Bzqxu9tf9jcK9jcuC84F3sRcAfgcbG+FpgOl9WY1vUPmJbAoZjRGhVtCcB/CO2FF+ApzGiIP4p/LE7apc6NvQK4xCU7LsQZ6eP8kH/BdGwdhGmtrHL+ppOB4x25kyqx7xcn/QVCR01EYmKgi5wf0eCgtLswIz6W4eIoAczb2cdOXYsxIYd+VVyrQ4EzKTOSBTO0dYljy41B54diRlysBO6q5fvuIswC5LFlZD937smLCXXkLYF0zFvE1DJ5TsDEz//AjEiKdM4f4NS3AuPUo2rzvsO0vjfgOOSgcuY513E0oY68l3NP9cWE/oLznO/cA4uAx4PO93bOrXTuYdkXG/fnzbUvO51xrQ9jbsQmwN2qOtaVyvYBEZkCXKiqm+paF7cQkXlAP1UtrGtd9gUReR14XVVn17UubiAi5wOtVfXRutbFLUTkbkyrvd75goaMm458AeZ16kFMy/cVTE/1ma5UaLFYLH9T3HTkvVV1bplzF6jqe65UaLFYLH9TanXSLIvFYrHUPG58EGSxWCyWWsQ6covFYmngWEducQUx/CwiQ4LOnS0ik6rKtwflfygiVwcd93PmCqn2a2URGeB8Dl7kfIFssTRobIzc4hrO1AafYD4x92LGg5+sqiv3oiyvmgm/io9TMOOa+wE5mHHPN6nqz0Ey7TFjlo8tU1Z7IBbzDcF4Vf0Ui6UB48ZcKxYLYGbKE5GvMF/4xmA+QrrLmQDLh/m4Z5zjWN9zZACuVdWZInIscB+wCTMHS9egsjNE5L+YCZ7mYGbQK3Hi1ei1BkBEAvtoosVSL7CO3OI292O+9isAvsZ8uXipiMQBv4rId5hP109Q1TwR6YT5vL23k78vcIiaeW/K8jLmK8pjg+Qtlr8d1pFbXEVVd4nIR8BO4GzgVBH5j5McjZmbZCPwgogcBvgx878U82slThxVDYjIK0BvNTMnAiWz9BXPs9LWmRYX4FlVfat8SRZLw8Y6ckttEHA2Ac5Q1ZB56Z3phTMwMxd6MNOoFrMrzLJLUDNLX6Uxcotlf8OOWrHUJt8C1zkr3CAiPZ3zzYBNqhoALsB0jFosljCxjtxSmzwIRAALRWSRcwzwf8BFzqIQnam+Fb5POMuapQFnAa84y5pZLA0WO/zQYrFYGji2RW6xWCwNHOvILRaLpYFjHbnFYrE0cKwjt1gslgaOdeQWi8XSwLGO3GKxWBo41pFbLBZLA+f/ATc7pbGWnG2fAAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 432x288 with 2 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"sns.heatmap(metric_matrix, annot=True, fmt=\"1.3f\", xticklabels=metrics, yticklabels=metrics,cmap='RdYlGn')\n",
"plt.xlabel('Year Y+1')\n",
"plt.ylabel('Year Y')\n",
"plt.title('Average R squared of metrics 2007-2019\\n (min TOI 50 mins)')"
]
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.7.3"
}
},
"nbformat": 4,
"nbformat_minor": 4
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment