{ "cells": [ { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "from scipy import stats" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "# matplotlib for plotting\n", "import matplotlib.pyplot as plt\n", "%matplotlib inline\n", "from IPython.core.pylabtools import figsize\n", "import matplotlib\n", "\n", "import json\n", "s = json.load(open('style/bmh_matplotlibrc.json'))\n", "matplotlib.rcParams.update(s)\n", "matplotlib.rcParams['figure.figsize'] = (10, 3)\n", "matplotlib.rcParams['font.size'] = 14" ] }, { "cell_type": "code", "execution_count": 17, "metadata": {}, "outputs": [], "source": [ "figsize(12, 6)\n", "x = np.linspace(-5, 5, 10000)\n", "y = stats.norm.pdf(x)\n", "n_tests = 1000\n", "\n", "# Random values\n", "x_random = np.random.randn(1000)\n", "\n", "# Non significant random values at p = 0.05\n", "x_random_nonsig = [x for x in x_random if abs(x) < stats.norm.ppf(0.975)]\n", "y_constant_nonsig = [0.0 for _ in range(len(x_random_nonsig))]\n", "\n", "# Significant random values at p = 0.05\n", "x_random_sig = [x for x in x_random if abs(x) > stats.norm.ppf(0.975)]\n", "y_constant_sig = [0.00 for _ in range(len(x_random_sig))]" ] }, { "cell_type": "code", "execution_count": 18, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAsAAAAGQCAYAAACgSdaFAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvNQv5yAAAIABJREFUeJzs3Xd4FFXbwOHfpEEKAUIJgSSEBAgQaaEXKQb4gtKRJl0URLoCKrygL0XwBVERKSqCSEcBKUrVUEMNICUgLUgXQk+AQDLfH7MbdpPdJLspm/Lc17XXJjNnzjxzZnb27JkzZxRVVRFCCCGEECKvsLN1AEIIIYQQQmQlqQALIYQQQog8RSrAQgghhBAiT5EKsBBCCCGEyFOkAiyEEEIIIfIUqQALIYQQQog8RSrAQgghhBAiT5EKsBBCCCGEyFOkAqyjKEqUoihqktcTRVH+URRlpaIojW0dY1roY7d1HOmlKIqdoihvKIryq6IoVxVFeaooyh1FUQ4pijJRUZTiKSybK8oguzD4bPjZOpbsIj1lkh2Oz+yyT1OKQ1GUtoqi7FEU5YHBObladok9J8qpZZdT404r3XfdLkVR7iuK8kj3PTdIURSr62iW5qkoykITdSDD12nrtzB7crB1ANnQZuCG7u/CQHWgE9BJUZT3VFX9wmaR5RGKongDa4EaQAJwANgJFADqAf8BhiuK0k9V1ZU2CzQXUBQlCigNlFFVNcq20eQOUqbppyhKdeBn3b9/ANd1f9+xTURCZA5FUb4B3gWeANuBZ0AIMAsIURSlk6qq8VmY5x7gnInp101My9GkApzcVFVVw/T/KIriCHwFDASmKoqySlXVK7YKLrdTFMUD2AX4AWHAm6qqXjSY7wi8D0wGliuKEq+q6i82CDUvCQEcgau2DiQbyellkl3iNxdHO7Tvp09VVR1rOENRlOwSuxDpoihKR7SK6g2gkaqqZ3XTPYE/gfbAYLQ6SFbl+b2qqgut2Z6cRrpApEJV1WdoFa6HgBPQwrYR5XrfoFV+DwItDSu/oO0PVVWnAu8BCvCDoihFszzKPERV1fOqqp7WfRYEOb9Mskv8KcTho3s/a8EyQuQ0H+neP9BXVAFUVb2J1ugG8KGFXSEyI8/cSVVVeakqQBSgAk3MzD+km/+RiXl1gGm6NDeBOOAa2iW8umbyU7XiVwG6AOHAI7SK9nagYQqxVgbWoF0OjAEigLeS5mtiudLAbOAC8BS4i/aL8I0U1mUYZx/dNsag/bqcDxTTzcsP/Bf4G+2yyz9orbSOFuyDACBet85aqaRVgL90aT9JIeb+wBEgFogGVgMvmckzEPgRuKTbhw91x8UaoKOZZVyB0WgV9gfAY+Ak8Anglkp59gP265ZTgbq693/NlRtgryt7FQiy9hjU7Us1hZefic+Gn4l8LDqmsPK4t2bfmCi3u2iXAgskmdfGYLtbJpnnrlvmLmCXUpmktUytLYNUts+i8klln1YDfuXF+eUw2pUYo9gzaL8axYH2uTFXfgtTi93gMzlSF8M9tM/kBWAV8GqStFl97k5zbNacX9JwnKS039O0LqACVp6nrN2m1PZ5CtsbjdaNrpQu/2O6fRUL7EhpX2X2C/DWbdNTwNlMmiu6NPUzO09goW56H1uVSZbvA1sHkF1epF4BPqub/6aJeduA57oP13rdyfO4Lv1zoJOJZfQn9Qlolb4dwAog0uAArmdiuca6D68KnAaWoXUViAdmYP4Lqg7al7iKdsJdDmxFO+GrwCJASSHOz3QxbQF+QesPpOq22Q2t39AdtC/c39C+OFXgWwv2wTDdMifSmH6kLn2EmZi/0JVLmK6czuimxyQ98aH9qNBXRCN12/gz2hdVLLDJxPq90U7c+i+DzWgVh2sGZVPYTGxf62LbCSxF+wIuCOzTzW9rZptf080/mJ5jEGiIdsJ7pEvzs+5//auoic+GX3qPKaw47q3ZN2bKbo0uj9ZJpn9lENfnSea11k1fbeZ84WdpmVpTBqlslzXHrrl9+gpahUSf1zK0HzTP0SqKKVaALd2mpHGgdX1YiNYHUQV2G5TfWynFrptXGu28qKJVSH9HOy736PZLmA3P3ZbGZvH5JQ3Hirn9btG6sP48ZdU2pbTPU9hWf90yt4ATuv35B9rn46rBvqpuSRlm1IsX55aIFNLoz1mDMjtPXlSAf0SrS3wLTAT+D4Mf/7npZfMAssuLFCrAQJDuwxMHeJuYHwp4mpjeWrdMNOCSZJ7+JBoN1DCYbqc78FRga5JlnHnx6+1TDCoXaBVjfaVTTbJcfrQWWX2l0N5g3ktoLR8qMMDENujjvAFUNJhemBcn8+No/XYLGsyvhtZylgCUTuM+WKTL74c0pm+sSx8POJiIOQatD5R+ugJM0c37B8hvMO8H3XRTLfxuJK+UKcBeXlRmXQzmOQM/YdBqZSK2e0BtE+t6BxMVLoP5K82cvCw+BpMc935p+Gz4GUyz6pjCuuPeon2TwnYM1uXzZZLpJ3TH903gWJJ5X5opb7PlllqZWlMGqWyXxeVjZp+68KIi8l+Mzy/10Sptyc4v6dkmc2VFCq1RKSxjh3Y1TEW7iTbpj88CQEh6PzdWHsMWxYaV55c0HCum9rvF68KK81R6tsncPk9lWzsb7Ku/gPIG81zRfriowIpU8llokI8lrxRjBYbq0q1JIY3+x/n0NG6z1Xmmsp0ngcqWHGs54WXzALLLCxMVYLRKXku0lsMEYKAV+S7R5ftakun6A2uwiWWK6+Y9weASE9BTN/0cBhUOg/mf6/NNMl2/XBQmLlmh9QtSgbMm5unj7G9i3nBeVEArmpj/q25+rzSW1e+69FPSmL6CQXzFTcSc7KSBdmnuvG5+d4PpG3XTqqVx3S116cMx8esY7QR7E+1HQGETsY0xk28htBa4OAxaYQ2OxydorRYe6T0Gkxz3fmn4bPgZTLPqmLLyuLdo36SwHRV1+ZwwmOapm7YUrbUzIcmxpG8NDExruaVWptaUQSrbZXH5mNmnvXTTzpg5pj/Tx55R22SurLCuAtxON/0iZi7/Wni8ZOS526LYsPL8YuV+t3hdWHGeSs82pfaZMrOt/9MtcxvTP3JCdPMvpZLPWxhfyUnrq2gq+Y7RrX9xCmkm69LMS+M2W50n2vf5EKCSbl94obXkH9WlvwmUSu/nKju9pBN0cn8ajNN5B+1yfmm0voFzzC2kKEpRRVH6KIoyXVGU73Vj6i1Eaw0DKG9m0Q1JJ6iq+i/apeV8QBGDWY1178tV00OY/GRmHfrllqimbxxZgHaAl1UUpZSZPDaZmKYfKuWSqqqRJubrO+CXNJNneimpzF+cdIKu3Jbp/m1iMOuA7n2uoijNFUXJl0rer+ref1FVNcHEemLQujU4ALVMLL/aVKaqqt5DayFyBN5IMrsr2jGxTlXVZMNBpeMYtEZ6jylLjntL941JumP0GhCkKEoJ3eQQ3fs23UvRT9PdNf0ScFVV1TPWrDMVlpRBSjKkfHixT1eYOqbRfiSkJqO2yRqhuvclqqo+TutCWXTutjS29J5fLGHxuqw8T2XlNoE2jCZoIzvdNDH/b927W0qZqKr6vaqqfax43U4lPv33l2rBNqXG6jxVVf1SVdWvVVU9papqjKqq11VV3QjURuvyUpwXN9jlCjIMWnL6cYAVoATQCO1y7yJFURqoqppsfDxFUQag9ZlxSSFfdzPT/zEz/QHaL+n8BtO8de8XzSwTZWa6vgJicjlVVZ8oinJNl64UpocXMjX026MU5hnOz29mflL6E4ZnGtPrH4aRgOnxQVMrJ2+DadOAl9EqP1uAp4qiHEW7TLZYVdXjSfLw1y+nKMq0VOIsZmLapRTSL0T7EukNzDSY3ttgvpF0HoPWSO8xZclxb+m+Scl2tNbrELQWPsMKsP7Loxnaj6RXDJbJDJaUQUoyqnz0+9TcsZnSMauXUdtkjdK69zQP2J+F525LY0vv+cUS1q5rIZadp7JsmxRFUYBg3b/JGkJ09BXff9OzrnR4mCQOU/TzHqaQJlPzVFU1TlGUKWhXdF9NLX1OIhXg5JKOA+yFVimuDCxRFKWuqrteoJtfE5iD1kd4FNqNFFeAWFVVVUVRPkX71WSytdJMS4u1zP3qS8uvwhRbU1OJM6O24TDQA200hLSorXs/pqrqcyvWl1geqqrGAs0URamD1lrTAO2hG3WA0YqifKyq6gSDZe117zsw/8NDL1nFIZVWoK1ox1CwoiiVVVU9rihKoC6WGyRpjU/vMWildB1Tlhz3VuyblJiqAJ9TVfUfAEVRzqNVgMG4cpzhMuqzn8HlA+b3aarxZvD5LFNls3N3Uuk6v2TRuiw6T6VjPdYoi9ZN4x9VVW+YSVNP934kpYwURXkL7QZXS41MpRU4SvdeOoU0+uEAo1JIk9l5wosfbuauEOdIUgFOhaqq1xVF6YzWib420B3jX5Svo50gZ6qqOt1EFmUzMBx9K5qfmfllzEzXt9D6m5qpKEp+tP4+huuwhQ1o/ZgrKopSS1XVg+YS6n7h99L9u95MMj+0u4pNTQftcrgRVVX3ow1NhqIoTmiX974DPlEUZYXBpfDLuvdVqqp+Yy5Oa6iqmqAoyk9oX7590Mah7qObvdhE95esPAb1svyYsmDfpERfmQ1RFCUA7YtibpL5AxRFKceLCvAfGRF/ZsuA8tF/Hsx9efplRJyZSF9pCkxj+qz83FgaW6adXzJqXVacp7Jym2rq3h+kkEb//ZHag5Qa8qJV2xKf8OKqpin6ineQoijOZhpFaiVJm5rMyBNedOd5lGKqHEb6AKeBqqqn0cY6Be3LxPCHg4fu/TJJKIpSDGiegaHs0L13VRTF3sT87qks1y1J7Hq90b4IzqmqarMKsK57if7xp9/oKlHmDEXro/cQ7eEZpiQrD125ddH9G5ZKPHGq9kScfWjlU8Vg9u+6904p5ZEOC3Xv3XWVmR5JphtKzzEYp3u39MewTY+pVPZNSstdRbvJy5cXg8IbdnHQ/90frcJ32or4rS3TDGNl+ezUvXc2M0h+twwKL7Ns1r33SOXcoZeV525LY8vs80tGrWuh7j0t56ms3CZ9/19fU+cnRVFaofV5P4k2LJhZuv68ihWvqFTyvYw2MogTJspEUZTGaN30bqDdOJiqzMhTp7Pu3WyjVI6kZoM78bLDi9THAS7Gi7E2+xpMf0837TDGg4UXQGuZ1N8x/EmS/EzeTW0iHj+DaYbDFE3EeJiihrwYf1RNkpfhkFXTMB7QvxIvBiw3OwyamRib6OaHmZn/ialtT2U/FEVrLVHRWt78ksx3BD5AG3kiAehsLmZdeTQ0mK4Ak3TzrmBwNzbaoyMDTeTlj9a/WMVg2DK0H4/6h6PMxcSoDLplkw6fleJ+T5J2jy69ftiag2bSWXUM6tL8qZvXOoU4TB2LVh1TVh73Fu2bNJTrN7plHuuOI8M71Yvojiv9WLhfpzXOtJapNWWQyvZYXD5mytnVYL+Nw/j8UocX579ksVu7TSlMX4jlo0DYobVqqWitegWTzE861FhWnrstjc2q80sajhVzsVm9LtJ+nrJ6PVZ8JsIM9t+YJPNC0IahfIqZh51k1QvtKoSKNq5+WYPpxXkxXvIwE8tNQeuWkGzEJGvyRBu2tBVJRpdC+xH/Hi8eUPV/tiyvDC9/WweQXV6kUgHWpRmvS3Me3bizaDc76CsC/6Ld3b8GbYzI62hPS8uQk6huetKB6peifeGm9iCMurx4aME5tJt8NutOAiqpPAjDTIxNyOAKsG45X158WTxHO8EuRftSiuZF5babmeX1Jz79gzD+1C2vH7c4FmicZBn9UC/n0Tr7L0FrDdSXzzIT6/HmxdPoHqCNhbwMrW+c/qEbNyzZ70nSvm2wLSrmv4CsOgZ1yw41iP9n4Hvdq0gajkWLj6nUtt/UuqzZN6mUa3uDMj1sYn6Ewfx2lnw+01Km1pRBKttjzbFrbp82RxvCSgVOoX1u/kD7HOrPL3GWnCdSWZ+56QuxsAKsm1eGFw8teoA2is8ytAdqGD1sgqw/d6c5NtXK80sajhVzsVm9LtJ4nkrPelLa5ybSKsB9tHO/fmjQCLSHjugr4E8w0Xhiixfa1WUV7Xt9ve44vK+btgbTQ57qPx8LMyJPXgzTF43WMrwKrQ+3/oEh8cBoW5dVhpe9rQPILi/SVgF240ULST+D6cWBebx4HOxltL53XpipBGLlSVQ3ryraF91dtMrcUeCd1PJF69s3B+3O/adov4J3oHUVSFb5TUN+TciECrBuWXu0y2nr0Vq943TbexitFbdECsvqT8QKWuvYMV056Z9Ul2xAb7Rfv3PRKt63dOXzD9pd9Z0w8yQctJbQQbpyvKOL8zraiXYayR83aUkF2J0XT/1Lcexfa45B3XJ2wH/QfkzpKz5JK6ApHYsWHVPWHPfW7psU1lGIFy0an5mYrx8/NB4oZMXnM8UytaYMUtkei8snlfiDgXW8OL8cAQag3TyjAtcyYr+mMn0hVlSAdfMLoI2Hehiti1QsL55U+H9J0mb1uTvNsenSW3R+ScOxklJsVq0LC85T1q4ntX2eJK1+fPiTaN8BH6D98HiK9l3yI1DBknLL7Bdaf/09aD8K9I8fH4T57x3952NhRuSJ9uPsS7QHlVxFO2891pXbDxg88CU3vRTdxgshhBBmKYrSE61Vf4Oqqq1tHY8QpiiKor9R/UdVVfvYOByRjclNcEIIIQBQFKW4oiilTUyvi9Y6B6ZvcBIiu9CPAJG7btgSGU6GQRNCCKFXBdiqKMoJtG4tcWg3JlXXzf9JVdXUho0SwpakAizSRLpACCGEAEBRFG+0PqqN0frBFkDrQ3gUreV3sSpfGiKb0g3fdx/tUcwFVFV9auOQRDYmFWAhhBBCCJGnSB9gIYQQQgiRp2R6H+CwsDA1X758mb0aofP48WOcnZ1tHYbIRLKP8wbZz3mD7OfcT/Zx1oqNjb0dEhJSLLV0mV4BzpcvHxUqVMjs1QidiIgIKe9cTvZx3iD7OW+Q/Zz7yT7OWhEREZfSkk66QAghhBBCiDxFKsBCCCGEECJPkQqwEEIIIYTIU6QCLIQQQggh8hSpAAshhBBCiDxFHoUshBAiW3j69ClPnjxBURRbh5JlChcuzIMHD2wdhshEso8zjqqq5M+fn4wYXlcqwEIIIWwuJiYGAHd39zxVAba3t8fV1dXWYYhMJPs446iqSmxsLM+fP093mUoXCCGEEDYXHx+Pq6trnqr8CiEsoygKrq6uxMfHpzsvqQALIYQQQog8RSrAQgghhBAiT5EKsBBCCCGEyFOkAiyEENnAvXtP2Lo1im+/PcbkyeH8+OMl5sw5wrZtUURHP7Z1eCIHW7p0KT4+Phme76BBg+jatWvi/wkJCYwYMYKAgAA8PDzYvXt3sjS5UVhYGB4eHty/f9/WoQgLyCgQQuRA9kePGv0fX62ajSIR6REb+4yffz7DsmWR7N9/3USKywDY2SnUqlWCbt0q0qlTBZyd5dSdXQwaNIhly5YxZswYRo4cmTh99+7dtGnThrNnz1KkSBFtYmys8cIuLhkSw6JFi5g/fz4XLlzAzs4OHx8fWrZsydixYwFo3749zZs3z5B1GZoyZQqqqib+v3XrVpYuXcq6devw8/OjcOHCVKlSxShNVhg0aBDR0dEsX77cbJoPPviA7du3c/DgwWQ3Xt67d4+KFSsydepUevfundnhChuRs6gQOZBUeHO2Z8/iWbDgBDNmHOTff7VKkZOTHcHBJahUqQjFirlw7do1HB0LERkZzaFDN9i//zr7919n4sS9jB5dh759K+PgIBfxkrp79y4XLlzg3r17FCpUCH9/fwoXLpyp68yfPz8zZ86kT58+FC1a1HzCDKrwGlq8eDEfffQRkydPplGjRsTFxXH69GkOHjyYmMbZ2RlnZ+cMX7e7u7vR/xcuXMDT05M6deokTnNycsrw9WaEXr168d1337F3714aNGhgNG/VqlU4ODjQoUMHG0UnsoKcPYUQIgsdP36LZs1W8OGHO/j331iqVy/ON98059y5/vz22+tMn96UDz6oQ+/epZk+vSkbN77OuXP9mTOnOdWrFyc6+gkffLCDxo2XceLELVtvTrZy9+5dDh06xNOnT/Hw8ODp06ccOnSIu3fvZup6GzZsiK+vL9OnT08x3d69e2nWrBleXl4EBgYyZswY4uLiEue3bt2akSNHMnHiRMqWLUv58uUZN24cCQkJZvPctGkTrVq1ok+fPvj7+1OhQgXatWvH5MmTE9OY6gLxxRdfEBgYiI+PDwMHDuSzzz6jatWqifP1XRfmzp1LUFAQZcqUYdCgQcQatGIbdm8YNGgQY8eO5cqVK3h4eCTmlbQLhKqqzJo1i5o1a1KiRAmCgoKYMGFC4vz//ve/1K5dm5IlS1K1alU+/vhjnjx5kjh/6tSp1K9fn19++YXg4GB8fX3p0aMH0dHRifOXLVvGli1b8PDwSOyKkVRQUBDVq1dnyZIlyeYtXryYtm3bUqBAAQBmzpxJgwYN8Pb2JigoiBEjRqT4YItFixbh5+dnNM1UN4l9+/bx6quvUrJkSYKCghg1ahQPHz5MnL97926aNWuGt7c3fn5+NGvWjDNnzphdr7CMVICFECILqKrK/Pl/ERKyguPHb+Pr687ixa+xbVsXunWriJub+ZayAgWc6NKlItu2dWHRotfw83MnMjKaZs1W8O23x7L8EnN2deHCBVxdXXFxcUFRFFxcXHB1deXChQuZul47Ozs+/vhjFixYwMWLF02muXbtGp07d6ZKlSqEhYUxc+ZMVq9ezdSpU43SrVq1Cnt7ezZt2sRnn33G3LlzWbNmjdl1Fy9enIiICKKiotIc7y+//ML//vc/xo4dy59//kn58uWZPXt2snTh4eFERkayevVq5s+fz8aNG5k3b57JPKdMmcKoUaMoWbIkkZGRbN++3WS6iRMnMn36dIYPH87evXtZsGABpUqVSpzv4uLCzJkz2bdvH9OnT2f16tV8/vnnRnlcvnyZtWvX8tNPP/Hzzz/z119/MWnSJAAGDx5Mu3btaNy4MZGRkURGRlK7dm2TsfTo0YN169YZVWaPHTvG8ePH6dGjR+I0e3t7pk6dyp49e5g3bx779+/no48+MlO6aXP8+HFef/11Wrduza5du1i4cCFHjhxh+PDhADx79ozu3bvTsGFDdu3axZYtW+jfvz92dlJtyyhSkkIIkcmePYtnxIg/GDUqjOfPE+jXrzK7d7/Bq68GWPTgB0VRaNUqgF27utOnz0vExSXw4Yc7GDlSyzevu3fvXrJL/c7Ozty7dy/T1928eXPq1KmTWBFL6ocffsDT05Pp06cTGBjI//3f/zF+/Hh+/PFHo1ZVfctw2bJlad++PQ0bNmTHjh1m1zt69GgKFy5McHAwNWvWZMCAASxfvpxnz56ZXWbevHl069aNXr16UbZsWUaMGEGNGjWSpStQoACff/45gYGBvPLKK7Rt25adO3eazNPd3R03Nzfs7e3x9PQ02RXk0aNHzJkzh48//pgePXrg7+9P7dq16devX2KaUaNGUbduXXx9fWnevDnvvfceq1evNsrn+fPnfPPNNwQFBVG7dm169+6dGJebmxvOzs7ky5cPT09PPD09zXbD6NixI4BR/osXL6Z8+fLUrVs3cdqgQYN4+eWXKV26NA0bNuTjjz9m7dq16frhOXPmTDp16sTAgQMJCAigVq1aTJs2jTVr1nD37l3u3bvHw4cPCQ0NpUyZMpQvX57OnTtTrlw5q9cpjEkFWAghMtHTp8/p2/d3Fi06Sf789syd24Jp05qm2OKbGldXR2bMeIXvvgslXz57Fiw4Ts+eG3j69HkGRp7zFCpUiMePjUfMePz4MYUKFcqS9X/yySf8+uuvHDlyJNm8v//+m1q1ahm14NWtW5e4uDijVuOgoCCj5by8vLh9+7bZdZYoUYItW7awe/du3nnnHVRV5b333qNZs2ZGFWtDZ8+eJTg42GiaqQpwYGAgDg4vbhUqUaIEt25Z3+3mzJkzPH36lEaNGplN8+uvv9KyZUsqVKiAj49PYrcKQz4+Pkb9j62Ny93dnTZt2iR2g3jy5Am//PKLUesvaN0X2rVrR1BQEL6+vrz55ps8fvw4xf2SmmPHjrFs2TJ8fHwSX61btwbg4sWLFCtWjM6dO9OhQwe6du3K7NmzuXr1qtXrE8lJBViIHMhp4UKjl8ienjx5To8eG/nttwsULJiPdes60rlzhQzLv2PH8qxd2wEPj/xs3hxF796/5elKsL+/PzExMcTGxqKqKrGxscTExODv758l6w8ODqZ169b897//TTZPVVWUp09Rbt9OfOkZXgVwdHQ0Wk5RlBT7AOtVqlSJt956i2+//ZbVq1dz/Phx1q5dazZ9Wq48WBuLOam1mB48eJC33nqLV155haVLlxIWFsaYMWOStWYbVsr1cVnbGtuzZ08OHz5MZGQkGzZsICYmhi5duiTOj4qKolu3blSqVImFCxfyxx9/8OWXXwIY9d82ZGdnlyyepNuQkJBAnz592LFjR+Jr586dHDp0iEqVKgEwd+5cNm/eTN26ddm4cSO1a9cmLCzMqu0UyckoEELkQK7vvWf0f1yfPrYJRJgVH59A//6b2b79EkWK5GfNmva89FKxDF9PnTperF3bnnbt1rBlSxT9+m3ixx9fxd4+77VvFC5cmJo1a3LhwgXu3LlDoUKFqFixYqaPAmFo3Lhx1KtXL1kf2MDAQNauWAGXLiW2Au+LjMTJySnZDVPpFRgYCEBMTIzJ+eXKlSMiIoLu3bsnTouIiMjQGMzFlS9fPnbu3ElAQECy+fv378fLy4tRo0YlTkva+psWjo6OxMfHpyltvXr1KFeuHEuWLOHEiROEhoZSrNiLz2lERATx8fFMmjQpcb/99ttvKeZZtGhRHj16RExMDK6urgCcOHHCKE2VKlU4ffp0qj/OKleuTOXKlRk+fDjt27dnxYoVNGnSJE3bJlKW986QQgiRyVRV5aOPdrJhw3nc3Z349dcOmVL51XvppWKsXdueQoXy8dtvFxg7dlemrSu7K1y4MDVq1CAkJIQaNWpkaeUXtFbo3r17J7vN8q4UAAAgAElEQVRZ7M033+RGdDTvffUVZy5dYvO+fUyYMIHevXvjko7h0d5//32mTZvGvn37uHz5MgcPHmTgwIG4uLjQtGlTk8sMGDCAZcuWsXjxYs6fP8/MmTM5fPiwRf3RrVGgQAEGDBjAhAkTWLJkCRcvXuTw4cP88MMPAJQtW5br16+zatUqoqKi+OGHH/jll18sXo+vry+RkZGcPXuW6OjoFPtDA3Tv3p0lS5awa9euZN0fAgICeP78OfPmzePSpUusWrWKb7/9NsX8atasibOzMxMnTuTChQts2LCBhUmu1I0YMYIDBw4watQojh8/zoULF9i0aRPv6Ro3Lly4wIQJEzhw4ABXrlxhx44dnD59OvHHjUg/qQALIUQGmzv3KN9//xdOTnYsXtyKSpVSGBs2g7z0UjEWL26Fk5Md3357jHnzjqa+kMgUo0aNSnaZvmTJkvw8ZQp/nTtHwwEDGDRtGh06dOCDDz5I17qaNGnC4cOH6devH7Vq1aJXr16ANtJD2bJlTS7TsWNHRo4cyYQJE2jSpAmRkZH07duXfPnypSuWtBg/fjzDhg1j+vTp1K1bl969e3Pt2jUAQkNDGTJkCGPGjOHll18mLCzMqtEWevXqRfny5QkJCaFcuXLs378/xfRdu3YlNjYWLy8vQkJCjOZVrVqVSZMmMWvWLOrVq8fSpUtNdnExVLRoUebMmcO2bdto2LAhK1asSLYdlStXZv369Zw/f57XXnuNxo0bM2nSJDw9PQFtNIyzZ8/Sp08fatasyZAhQ+jWrRuDBw+2uDyEaUpmD58THh6uVqiQcX3eRMoiIiKS3dwgcpeIiAhCmjUzmnb3zh0bRSOSCg+/Sps2q4mPV/n++1A6dChvVT7WfpZXrTrNgAFbsLNT+PXX9jRo4G3V+rPagwcPkj1YIbexT3JzXHz16kaXyW2pZ8+ePH/+nGXLltk6lFwnu+zj3CSl80VERMThkJCQmqnlIS3AQgiRQW7ciOHNN38nPl5lyJBgqyu/6dGpUwWGD69BQoLKW29tSnzSnBB6sbGxzJo1K7GbwIwZM/jtt9+M+gQLkdtJBVgIITJAQoJK//6buHkzloYNSzFuXH2bxTJmTD0aNCjFzZux9O+/ifh4GSNYvKAoCtu3b6dVq1Y0adKENWvWMHfuXFq1amXr0ITIMjIKhBBCZIDZs4+we/dVihd34fvvW+LgYLv2BQcHO777LpQmTZaxc+cVvvnmCEOHJh/nVeRNzs7OKT5dToi8QFqAhRAinU6dimbSpL0AzJwZQvHi1t/Vn1FKlHBl1iytr/inn4YTGRlt44iEECL7kAqwEEKkQ1xcPAMHbiYuLoFevYJo0aKMrUNK1KyZHz17BhEXl8CgQVt59ixtY6MKIURuJxVgIYRIh2++ieD48dv4+bkzadLLtg4nmYkTG+LjU4CjR/9l5szMf9iBEELkBFIBFkIIK0VF3WfatAMAfPFFCG5uTjaOKDl393zMnKl1hfj88wNERd23cURCCGF7UgEWQggrqKrK6NFhPHkSz+uvB9K4sY+tQzKrcWMfOncO5MmTeD74YAeZPf67EEJkd1IBFkIIK6xff55t2y7h7u7ExIkNbR1Oqv7734a4uzuxdWsUGzdesHU4QghhU1IBFkIIC8XGPmPMmJ0AjB9fH0/P7P+UJ09PV/7zn3oAfPjhDmJintk4oryjdevWjB49OsPz9fDw4Ndff038/++//6ZFixZ4eXlRtWpVk2lyowEDBtCjRw9bhyFyGBkHWIgcKGbGDFuHkKfNmXOUa9ceUaVKMXr3fsnW4aRZ376VWbz4FH/9dYs5c44wcmRtW4eU492+fZupU6eydetWbt68ScGCBalYsSLDhg2jadOmACxatAgHBwcSnj7N0HVHRkZSqFChxP8nT56Mi4sL+/fvx8XFxWSarODh4cGCBQto27atyflxcXEEBQUxYMAARo4cmWz+/PnzGT9+PJGRkbn+8djCdqQCLEQOFNenj61DyLNu3ozhq68OAdoIC/b2OedCmr29HRMnNqRt2zV89dVhevYMyhGt15ZwXLUK54kTsbt6lYRSpXg8bhzPOnXKtPX17t2bx48fM3PmTPz9/bl16xZ79uzhzp07iWkKFy4MgFqgQIau29PT0+j/ixcv0rJlS3x9fc2myQ6cnJzo1KkTS5cu5f3330dRFKP5S5YsoU2bNlL5FZkq55y5hRAiG5g6dT+PHj0jNLQML7+cfW98M+fll30IDS1DTMwzPvtsv63DyVCOq1bhOnw49leuoKgq9leu4Dp8OI6rVmXK+u7fv094eDjjx4+ncePG+Pj4EBwczJAhQ+jYsWNiuqRdIP7991/eeOMNSpYsSd26dVmyZAn169dn6tSpiWk8PDxYuHAhffr0wdvbm+rVq7Ny5Uqj9Rt2b/Dw8ODEiRNMmzYNDw+PxLySdoG4fv06/fv3JyAggFKlStGoUSN27doFaBXo7t27U6FCBby9vWnSpAmbN282WmfVqlWZPn06I0aMwNfXl6CgIGbOnGk0H6Bv3754eHgk/p9Uz549iYqKYvfu3UbTT5w4wdGjR+nZsyegtbD369ePoKAgSpUqRf369Vm+fHlKu4WWLVvy0UcfGU1L2k0iISGBL774gurVq1OyZEkaNGjAL7/8kjhfVVWmTp1KlSpVKFGiBBUrVmTQoEEprlfkLFIBFkKINDp1KpqffjqJvb3CJ580sHU4Vvv44wbY2yv89NNJTp/OPU+Ic544EeXxY6NpyuPHOE+cmCnrc3V1xc3NjU2bNvHkyZM0L/fuu+9y+fJl1q5dy/z581m1ahVXrlxJlm7atGm8+uqr7Ny5k/bt2zNkyBAuX75sMs/IyEjKlSvHoEGDiIyMZPDgwcnSxMTE0Lp1a/755x9++ukndu/ezahRo4zmN2vWjNWrV7Nz505at25Nr169+Pvvv43ymTNnDpUqVSIsLIxhw4bxySefcOCANhzg9u3bAfjyyy+JjIxM/D+pihUrUqNGDZYsWWI0ffHixQQEBFC/fn0Anjx5QrVq1Vi+fDl79uzh7bffZujQockqzpaaMGECy5cvZ/r06YSHhzNs2DCGDh3Ktm3bAFizZg1z587l888/59ChQyxZsoTg4OB0rVNkL1IBFkKINJo8eS8JCSp9+1amfHkPW4djtcBAD3r3fon4eJUJE/baOpwMY3f1qkXT08vBwYFZs2axcuVKypQpQ4sWLRg3bhyHDh0yu8zZs2f5448/mDFjBrVr1yYoKIhZs2YRGxubLG2XLl3o3Lkz/v7+jBkzBgcHB8LDw03m6+npib29Pa6urnh6euLm5pYszc8//8y///6b2OJcpkwZWrduzcsvaw9weemll+jbty+VKlXC39+f999/nypVqrBu3TqjfJo2bcrbb7+Nv78//fv3x9/fn507tZtCixYtCkDBggXx9PRM/N+Unj17sn79eh48eADA06dPWbVqlVFLrbe3N0OGDKFy5cr4+fnRt29fWrZsyerVq83mm5qHDx8yb948vv76a0JCQihdujSdO3eme/fu/PDDDwBcuXKFEiVK0LRpU7y9vQkODqZfv35Wr1NkP1IBFkKINIiIuMnvv1/ExcWBUaNy/s1jo0fXwcXFgU2bLhIRcdPW4WSIhFKlLJqeEdq0acOpU6dYunQpzZo148CBA7Ro0YIZZm5UPXv2LHZ2dlSvXj1xmre3NyVKlEiWNigoKPFvBwcHihQpwq1bt6yO9fjx41SqVIkiRYqYnB8TE8PHH39M3bp1KVOmDD4+Phw9ejRZ67RhXAAlSpTg9u3bFsfToUMH7O3tE7sebNy4kYcPH9K1a9fENM+fP2fatGk0bNiQgIAAfHx8+P333022mKdVZGQkT58+pWPHjvj4+CS+Fi1axMWLFwFo3749jx49olq1agwdOpRff/2VuLg4q9cpsh+pAAuRA9kfPWr0Eplv6tR9ALz9dlWKFXOxcTTpV7y4C2+/rfXP1G9bTvd43DhUZ2ejaaqzM4/HjcvU9ebPn5+mTZsyevRoNm/eTI8ePfjss8+SV5hiY1H1XSVMtPgm5eBgfJ+6oigkJCRYHWdqD0AZP348v/76K2PGjGHDhg3s2LGD4OBgnj0zHjIvo+Jyc3Ojbdu2id0gFi9eTIsWLYxu3Pvqq6+YN28eQ4cOZe3atezYsYPQ0NBkMRmys7NLtq2G6fWxLl++nB07diS+9u7dyypdf3EfHx8OHjzI9OnTcXV1ZezYsYSEhJhsqRc5k1SAhciB3F95xeglMteBA9fZtu0Sbm6ODB6ce/oBDh4cjJubI9u2XeLAgeu2DifdnnXqRMyXXxLv7Y2qKMR7exPz5ZeZOgqEKYGBgTx//jxZv2D7M2eo6OBAQkICxzdsSJx+9epVbty4kelxValShVOnThEdbbrf9759++jSpQtt2rQhKCiIkiVLEhUVZfF6HB0diY+PT1Panj17EhERwebNm9m5c2fizW+GMbVs2ZLOnTtTuXJlypQpw7lz51LMs0iRIty8+eKqhqqqnDx5MvH/SpUq4eTkxJUrV/D39zd6eXt7J6ZzdnYmNDSUKVOmsGXLFk6ePJli9xaRs0gFWAghUjFlitZCOmBANYoUcU4ldc5RpIgzAwZUA3JPK/CzTp148Ndf3IuO5sFff2Vq5ffOnTu0bduWlStXcvLkSS5dusTatWv5+uuvadSokclhvMr5+BBSqxbDv/ySgwcPcvLkSQYPHoyLi0uy4cAy2uuvv07RokXp0aMH4eHhXLp0id9//z1xFIiAgAA2btzIsWPHOHXqFAMGDLDo5j49X19fdu7cyc2bN7l3716KaWvXrk1gYCADBw7E09OTZs2aGc0PCAggLCyM/fv3c+bMGd5//32uptKnu1GjRmzdupXNmzdz9uxZPvroI6MfGO7u7gwcOJCxY8eydOlSLl68yF9//cX8+fNZtGgRoLVGL168mFOnTnHp0iWWLVuGo6Mj/v7+FpeHyJ6kAiyEECnYu/cqO3ZcpkABJwYNqp76AjnMu+9Wp0ABJ8LCLhMenjk3i+VWrq6u1KxZk3nz5tGqVSvq16/PpEmT6NixY+LNVKbMGT2aUkWL0qZNG9588006depE0aJFyZ8/f6bHu379ery8vOjWrVvi0Gv6ivfkyZMpVqwYr732Gp06daJmzZrUq1fP4vVMnDiR3bt3U7lyZRo3bpxq+h49enDv3j26deuGvb290bwPPviAatWq0alTJ1q3bk3BggVp3759ivn16tWLLl268O677/Lqq69SuHBhQkNDjdKMHz+e999/n6+++oq6devSsWNHNm7cSOnSpQHtJr4ff/yRV199lYYNG/L777+zePFioxZikbMpqfUJSq/w8HC1QoUKmboO8UJERIQM1ZLLRUREEJKkleSuwaD7ImN16LCGsLDLjB5dmw8/rJtl683Kz/LUqfv43/8O0LixD2vWpFy5yCwPHjzI9Q8+sD9yxOj/+OrViYmJ4cmTJ1SqVInvvvuONm3a2Cg6kVliYmJwdc1dD5yxtZTOFxEREYdDQkJqppaHtAALIYQZR4/+S1jYZdzcHHnnnWq2DifTvPNONdzcHNmx4zJHj/5r63ByvR1HjrBxzx6ioqKIiIigX79+FClShJCQEFuHJkSeIRVgIYQwQ//I4z59KlOoUOZenralQoXy06dPZeDFNovM8+z5cyYtWEDDhg3p27cv+fPnZ8OGDdJKKEQWckg9iRBC5D3nz99j3bpzODnZMXBg7m391Rs4sBrffnuUdevOcf78PQICCtk6pFyrWa1aNKtVK7ELhFR8hch60gIshBAmfP31YVQVOneugJdX8qdq5TZeXm506VIRVdW2XQghcjOpAAshRBLXrz9i+fJIFAWGDKlh63CyzJAhwSgKLF8eyfXrj2wdjhBCZBqpAAshRBJz5x4lLi6BVq0CKFeusK3DyTJlyxamdeuyxMUlMHeuPGFQCJF7SQVYCCEMPHwYx8KFJwAYNizVkXRynaFDtRbvRYtO8uhRXCqphRAiZ5IKsBBCGFi+PJKHD+OoW7ckwcGetg4nywUHe1K7thf37z9lxYrTtg5HCCEyhVSAhRBCJyFB5dtvjwEwYEBVG0djO/oxj7/99hgJCZn7sCQhhLAFqQALIYTO9u2XOH/+Ht7eBXjttQBbh2MzrVoFUKqUG2fP3mX79ku2DkcIITKcjAMsRA70tFcvW4eQK+lv/HrrrSo4OOTd9gEHBzvefrsqn3yyh3nzjtK8uZ+tQ8oV1CJFbB2CEEIn757hhcjBYr/80ugl0u/MmTv8+ec/ODs70LNnkK3DsbmePYNwdnbgjz/+4fTpaFuHk6MNGjSIrl27kuDra/TKqHz1EhISGDFiBAEBAXh4eLB79+5kaXKjjNrGzCirhIQEPvzwQ6N9IrIHaQEWQghI7PvbpUsFChfOvY89TqvChfPTpUsFFi48wbffHmPGjFdsHVK2dfv2baZOncrWrVu5efMmBQsWpGLFigwbNoymTZsyZcoUVDXj+1InzXfr1q0sXbqUdevW4efnR+HChalSpUqmrDurDRo0iOjoaJYvX55sXkaVb9J8WrduTcWKFfnf//5ndZ5bt25l5cqVRvskp5o/fz5ff/01N2/epEKFCnz66afUq1cvXctMnTo1WfkWL16c06cz/wZcaQEWQuR59+49YcWKSAD698+7N78lpS+LFStOc+/eExtHk3317t2biIgIZs6cycGDB1m2bBkhISHcuXMHAHd3dwoWLJjh602a74ULF/D09KROnTp4enri5OSUaevOTjJqGzOjrC5cuEDx4sWN9ok14uJsOyTh6tWr+eijjxgxYgRhYWHUrl2bzp07c+XKlXQvU65cOSIjIxNfWdVKLhVgIUSet3RpJLGxz2nc2IcKFaSfpl6FCkVo3NiHx4+f22xItMIeHla9CjRtajbPAk2bml3OUvfv3yc8PJzx48fTuHFjfHx8CA4OZsiQIXTs2BFIfmk9JiaGgQMH4uPjQ/Xq1fniiy/o2rUrgwYNSkzTunVrRo4cycSJEylbtizly5dn3LhxJCQkJKYxzHfQoEGMHTuWK1eu4OHhQdWqVU2uW1VVZs2aRc2aNSlRogRBQUFMmDABgG3btvHqq69SpkwZ/P396dixI2fOnDHa3rTEldI69PNnzpxJcHAwJUuWpEGDBqxcudLisjdVDvoY33//ff7zn//g7+9PuXLlmDt3Lk+fPmXUqFH4+flRuXJlVqxYYTafQYMGsWfPHr7//ns8PDzw8PDgn3/+Ye/evTRv3hwfHx9Kly5Ns2bNOHXqlNm4xo4dy9WrV432ydOnT/noo48IDAzEy8uL5s2bs2/fvmTl/P777zNu3DjKlStHy5YtTa4jOjoaDw8PZs+eTUhICF5eXtSqVYs//vjD6vI0Zfbs2XTr1o3evXsTGBjIZ599hqenJz/88EO6l7G3t8fT0zPxVbRo0QyN3RypAAsh8jRVVfnxR+3BF2+9VcXG0WQ/fftWBmDBguO54lJ6RnN1dcXNzY1Nmzbx5EnaWsnHjRvHnj17WLRoEcuXL+fEiROEh4cnS7dq1Srs7e3ZtGkTn332GXPnzmXNmjUm85wyZQqjRo2iZMmSREZGsn37dpPpJk6cyPTp0xk+fDh79+5lwYIFlCpVCoDY2Fjeeecdtm3bxvr163F3d+eNN95I1vqYWlwprQNg8uTJLF68mP/973+Eh4czYsQI3nvvPbZs2ZKm8kuLVatW4ebmxtatWxk2bBhjxoyhR48eBAQE8Mcff9C1a1eGDRvG9evXTS4/ZcoUatWqxRtvvJHYMunp6Un37t2pW7cuO3fuZOvWrbzzzjvY29ubzWPUqFF4eXkZ7ZOPP/6YtWvX8vXXXxMWFkalSpXo1KkTN27cSLYNqqqyceNGZs+ebXIdf/31FwDff/8948ePZ9euXVSqVIn+/fvz+PFjo7QzZszAx8cnxZep4zAuLo5jx47RNMmPyqZNm3LgwAGTcVmyzKVLlwgKCqJatWr069ePqKgok3lmNOkDLITI0/bsucrZs3fx8nLl//6vjK3DyXZatixDiRKu/P33XcLDr1G/fqnUF8pDHBwcmDVrFsOHD+fHH3+kcuXK1KlTh7Zt21KzZvInCT569IglS5Ywe/ZsmjZtSkxMDDNnzuSll15KljYwMJAxY8YAULZsWRYtWsSOHTsSW5YNubu74+bmltiaZsqjR4+YM2cOn376KT169ADA39+f2rVrA9CmTRuj9LNmzaJ06dJERERQt27dNMWV2jpiYmKYPXs2v/zyS2Jf0NKlS3P48GG+//57WrRokUJpp12FChX48MMPAa0l9quvvsLR0ZF33nkHgNGjRzNz5kwOHDhA27Ztky3v7u6Ok5MTLi4uieV59+5d7t+/T2hoKGXKaOeK8uXLm43B1D6JiYlhwYIFfPXVV4nbOmPGDHbt2sX8+fMZO3Zs4vK+vr5MmjQpxe08fvw49vb2rFixgnLlygHwySefUKNGDc6ePUuVKi9+1Pft25d27dqlmJ+Xl1eyadHR0cTHx1O8eHGj6cWKFWPHjh0m80nrMjVq1GDWrFmUL1+eW7du8fnnnxMaGsrevXvxsOKKjCWkAixEDpT0Uu1dXV9DYbkFC44D0KNHUJ4e+swcR0d7evQIYvr0A/zww3GpAJvQpk0bWrRoQXh4OAcPHmT79u188803/Oc//+G9995LTGd/5Aj/nD/Ps2fPqJX/xY2Wrq6uVKxYMVm+QUHGo5F4eXlx+/Ztq+M8c+YMT58+pVGjRibnX7x4kU8//ZTDhw9z+/ZtVFUlISEhWZ/NlOJKbR1nzpzhyZMndOrUCUVREqc/e/YM3wwYGcNUjIqiULRoUSpVqpQ4zdHRkUKFCllUnoULF6Zbt268/vrrNGrUiEaNGtG2bVu8vb3TnEdUVBTPnj2jTp06idPs7e2pVatWsu4m+i4TKTlx4gShoaGJlV/Qts1c/Om5Cc9wf4F29SzpNEuXad68udH8mjVrEhwczLJly4y6BGUGqQALIfKsW7di2bDhPHZ2igx9loJevYKYMeMg69ef49atWIoVc8mydWfGj7uHf/6Z4Xnmz5+fpk2b0rRpU0aPHs3QoUP57LPPGDx4sFE6fTeS1CoOkLwioyiKUV9bS6XWhaVbt254eXkxY8YMvLy8cHBwoF69esm6QKQUV2rr0KdbunRpsoqjuYqbNUzF6ODgkGyapeX5zTffMHDgQLZv386mTZuYPHkyP/30EyEhIWla3pL97+rqmmqa48eP0759e6NpR44cIX/+/JQtW9Zo+owZM/jiiy9SzG/lypXJRnYoUqQI9vb23Lx502j67du3KVasmMl8rFkGwM3NjQoVKnDhwoUU48wI0twhhMizli49xbNnCbRo4Ye3dwFbh5NteXsXoEULP549S2DpUtM3/AhjgYGBPH/+PFm/YP9SpXB0cOCwwTBPsbGxREZGZklM+fLlY+fOncnm3blzh7///pv33nuPJk2aEBgYyKNHj3j+/HmGrcNw/uXLl/H39zd6+fj4WLVdmcXR0ZH4+Phk01966SWGDRvG+vXradCggcmh2cwpU6YMTk5ORje9xcfHc/DgQQIDAy2K78mTJ5w7dy5ZJX7u3Lm0b98eFxfjH6p9+/Zlx44dKb6qVauWbD1OTk5UrVqVsLAwo+n6kR1MsWYZ/TadPXvWbDeejCQtwEKIPCkh4cXNb336JO9/KYz17VuZTZsusnDhCYYMqYGdXeotWHnBnTt36Nu3L927dycoKAg3NzeOHDnC119/TaNGjXB3dzdK7+bsTI/QUD757jsK16iBu7s7s2fPTlfLbloVKFCAAQMGMGHCBJycnKhfvz537tzh2LFj9OnThyJFirBo0SJKlSrF9evX+fjjj5O1mqZnHW+++SYFChRg8ODBjB8/HlVVqV+/PjExMRw6dAhFUejTp4/ZvB8+fMjx48eNpmXmEG++vr5ERETwzz//4OrqyoMHD1i0aBGhoaGULFmSqKgoTp06Rd++fdOcp6urK3379mXChAkUKVIEX19f5syZw61bt+jXr59F8elHn1i9ejWNGzemaNGifP7551y8eNHk6Azp6QLx7rvvMnDgQIKDg6lTpw4LFizgxo0bRtv+3Xff8f3337N///40LzNu3DhCQ0Px9vbm1q1bTJ8+nZiYGLp162ZVnJaQCrAQIk8KC/uHqKgHeHsXICSktK3DyfZeecUXX193Ll16wB9/XKJZMz9bh5QtuLq6UrNmTebNm8eFCxeIi4vDy8uLjh07MnLkSJPLTHrnHWKfPKF79+64uLjw7rvvcuvWLfLnz/wHsIwfP55ChQoxffp0rl27RrFixejatSt2dnbMnz+fDz/8kAYNGlCmTBkmTpyYYoXU0nXojRkzhmLFijFr1ixGjhxJgQIFeOmllxg6dGiK+YaHh9O4cWOjaa1bt8bNzc3iGNNi8ODBvPvuu9SrV4/Hjx+zZcsWzp8/z5tvvkl0dDTFihXj9ddfZ9iwYRbl+8knnyTmf//+fSpXrsyqVasoUaKERfkcP34cf39/Ro8ezVtvvcW9e/do2rQp27Zty/AW1A4dOnD37l0+//xzbt68ScWKFVmxYoVRq310dDRnz561aJlr167x9ttvEx0dTdGiRalRowZbtmzJkqsBSmYPaxMeHq5WqFAhU9chXoiIiCA4ONjWYYhMFBERQUizZkbT5CY4y/XqtZENG84zdmw93n+/lq3DSSY7fpa/+OIgEyeG89pr/vz0U6sMzfvBgwfJWktzG/sjR4z+j69enZiYGBwcHKhatSqDBw9O1mdY5HwxMTFp6s9rqVGjRnH79m0WLFiQ4XlndymdLyIiIg6HhIQkH4IlCekDLITIc65ff8Tvv1/AwcGO7t0rpb6AAOCNNyphb6+weXMU//4ba+twcqxjZ8+yaim8YlMAACAASURBVPt2Lly4wIkTJ3j33Xd59OhRspuZhEjJ8ePHk43IIdJOKsBCiDxnxYrTxMerhIZqY9yKtPH0dKVFCz+eP0+w2ZPhcotZP/9MkyZN6Nq1K7du3WLDhg1GD4sQIiWqqnLq1CmpAKeD9AEWQuQpqqomjmQgrb+W6969Er//fpElS04yeHD1NA3nJIxVLVeOHXPmJHaByIzL4yJ3UxSFf/75x9Zh5GjSAiyEyFP277/OuXP38PR0kZvfrNC8uR/Fi7vw9993OXToRuoLCCFENiQVYCFEnqJv/e3SpaI8+c0Kjo72dOmi3di8eLGMCSyEyJnk7C+EyDNiYp6xdq02TM8bbyR/9KxIG33XkTVr/iYm5pmNoxFCCMtJBVgIkWesX3+OR4+eUatWCcqX97B1ODlW+fIe1K7txaNHz1i37pytwxFCCItJBVgIkWcsWaJdsn/jDbn5Lb30rcCLF5/MsDwze1x6IUTOp6pqhpwrpAIsRA70vGpVo5dIXVTUffbsuYqzswPt25ezdTg5Xrt25XB1dSQ8/Brnzt1Nd34uLi7cv38/V1eCVRcXo5cQwjKqqnL//v0MGTlFhkETIgd6+Oeftg4hx1m6NBKANm3K4u6ez8bR5HwFCjjRtm05li49xdKlpxg/vkG68nNwcMDNzY0HDx7k3qHVvLyM/3/wgOjoaIoUKWKbeESWkH2ccVRVxc3NDQeH9FdfpQIshMj14uMTWLZM3/1Bbn7LKD16VGLp0lOsWHGasWPrYW+fvouKDg4OFCxYMIOiyxnOnTtHmTJlbB2GyESyj7Mn6QIhhMj1du68wtWrjyhd2p0GDbxtHU6uUaeOF35+7ly/HsOuXVdsHY4QQqSZVICFELnesmVa94du3SpiZ5dLL6/bgKIodO6sjQm8cqU8GlkIkXNIBVgIkas9fBjHxo3nAejaVbo/ZDR9BXj9+vMyJrAQIseQCrAQIlfbsOE8jx8/p169kvj6uts6nFzH378QtWqVICbmWeIPDSGEyO7kJjghcqACTZsa/S+jQpinvzTfqVMFG0eSe3XtWpGDB2+wYsXpxBZhkZx8boXIPqQCLEQO5HDsmK1DyBGuX3/Ezp2XcXKyo127srYOJ9dq164cH320gx07LnP9+iO8vNxsHVK2JJ9bIbIP6QIhhMi1fv75b1QVWrQoQ6FC+W0dTq5VuHB+WrQoQ0KCys8//23rcIQQIlVSARZC5FqrVmndH7p0kcvymU1fxjIahBAiJ5AKsBAiVzp16jYnTtymUKF8NGtW2tbh5HrNm/tRuHB+Tp68zYkTt2wdjhBCpEgqwEKIXGnlyjOA1j81Xz653SGzOTnZ0759OQBWrJBWYCFE9iYVYCFErpOQoCZ2f5BRCbKOvhvEzz+f4fnzBBtHI4QQ5kkFWAiR6+zefYXr12Pw9XWnTh0vW4eTZ9SsWYKAgELcvBnLjh2XbR2OEEKYJRVgIUSuo+/+0KlTIIoijz7OKoqi0KlTIKC1AgshRHYlFWAhRK4SG/uMdevOAtL9wRY6dtQqwBs3ak/gE0KI7EgqwEKIXGXTpos8evSM4GBPypUrbOtw8pyAgEJUr16cR4+esXnzRVuHI4QQJkkFWAiRq8jNb7anbwVevVoeiiGEyJ6kAiyEyDXu3HnM9u3/YG+vJA7JJbJe+/blUBTYsuUi9+8/tXU4QgiRjFSAhRC5xvr153n+PIHGjX0oVszF1uHkWV5ebjRs6E1cXAIbNpy3dThCCJGMVICFELnGmjXaJff27cvbOBLRsaO2D2Q0CCFEdiQVYCFyoLt37hi9BNy8GcPu3VdxdLSjVasAW4eT57VpUxZHRzt27brCzZsxtg4nW5DPrRDZh1SAhRC5wrp150hIUGnWrDQFC+azdTh5XqFC+WnWrDQJCSpr1561dThCCGFEKsBCiFxBP+KAdH/IPjp00EaD+OUXGQ1CCJG9SAVYCJHjXbnykP37r+Ps7EBoaBlbhyN0QkPL4OrqyKFDN4iKum/rcIQQIpFUgIUQOZ7+Envz5n64uTnZOBqh5+rqSMuW/oCMCSyEyF6kAiyEyPH0FeAOHaT7Q3bz+usyGoQQIvtxsHUAQgjLuQwfbvR/7Jdf2igS27t48R4RETdxc3OkeXM/W4cjkmjSxJfChfNz+vQdTp26TaVKRW0dks3I51aI7EMqwELkQPkWLTL6Py9/ka5Zo7X+tmzpj7OznNKyGycne9q2LcvChSf4+eczjB+fdyvA8rkVIvtIUxcIRVHaKIoi3yxCiGxHXwGW7g/Z1+uvvxgNQlVVG0cjhBBp7wM8EbiuKMosRVHqZGZAQgiRVmfO3OHkydsULJiPpk19bR2OMKNu3ZKULOnG5csPOXDghq3DEUKItFWAVVWtCjQDHgO/KIpyRlGU/yiK4peJsQkhRIr0jz5u1SoAJyd7G0cjzLGzU2jXrhzwYp8JIYQtpXkUCFVVj6mqOgrwAQYBnYDziqLsVBSlu6IoMqKEECLLqKqa2P2hfftyNo5GpEa/j9atO0d8fIKNoxFC5HUWVVoVRQkAxgNzgPy6v78DBgM/Z3h0QghhxsmTtzl79i5FiuSnUSMfW4cjUhEc7Enp0u7cuBHDvn3XbB2OECKPS+tNcIMURdkH7Ac8gZ6qqgaqqjpZVdWfgBCgRSbGKYQQRvQPVmjTphwODnIBKrtTFMNuEGdtHI0QIq9L67dGS+BzoKSqqu+qqrrPcKaqqrFAh4wOTgghTDHs/tChg3R/yCkMu0E8fy7dIIQQtpPWCnCYqqqrVFWNM5yoKMp7+r9VVd2SoZEJIYQZERE3uXTpAV5ertStW9LW4Yg0qly5GAEBhbh9+zG7d1+xdThCiP9v777jpKruPo5/zpSd7RUEdlmkC1hQVEQUQRDFgigGFbAbjUk0MUZTNDwxxhIVCY/Rx8TELjZUFBUVC/aCAWwU6dWls5StU87zx+wOu8DCLOzs3Zn5vl+v+5q5d9p3d+be+c25556bxKItgP+ngeV/aqogIiLRqtv9we1W94d4oW4QItJS7PWbwxgz2BgzGPAYY06una+Zfgpsb56YIiJhoZDllVfU/SFenXtu+IQlr7++BL8/6HAaEUlW+zq72yM1lz7g0TrLLbAWuC4WoUREGvLllz9SUlJGcXEWxxzT1uk40ki9ehXQo0c+CxZs5oMPVjF0aEenI4lIEtprC7C1tpO1thMwqfZ6zdTZWtvfWju1mXKKiADUGfu3O8YYh9PI/qhtBa5tyRcRaW7RngnuklgHERHZl0AgxKuv6uQX8a72vXvjjSVUVQUcTiMiyajBLhDGmPnW2p4111cR7vawG2tthxhlE5EGbHv/facjOOLTT9ewYUMFXbrkcsQRrZ2OI/upa9c8Dj+8Fd99t5H331/J6ad3djpSs0jW9VakJdpbH+Cr6ly/KNZBRCR6wSOPdDqCI2pHf1D3h/h37rnd+e67jUyZsihpCuBkXW9FWqIGC2Br7Sd1rn/YPHFERPasujrIa68tBtT9IRGcc043brvtM956aykVFQHS0vZ1TLaISNOJ9lTINxhjjqy53s8Ys9IYs9QYc3xs44mIhH344SpKS6vo2bOAnj0LnI4jB6hjxxz69GnDjh1+3nlnudNxRCTJRDuC/G+AZTXX7wImAHcAE2MRSkRkVzu7P6j1N1HsPCnGQoeTiEiyibYAzrHWbjXGZAG9gX9Yax8BDoldNBGRsMrKAG+8sQTYOYSWxL/aAnj69OXs2FHtcBoRSSbRdrpaZYzpDxwKfGStDRpjsgGdxkfEASmPP15vvvqyyxzJ0VzefXcFO3b46d27NV265DodR5pI+/ZZ9O3bjpkzS3j77eWcd15i/7hJtvVWpCWLtgC+CXgRqAbOq1l2FjAzFqFEZO8ybrih3nyif5HW7iJX62/iOffcbsycWcIrryxM+AI42dZbkZYs2hNhTLPWFlprO1prZ9UsngycHbtoIiJQVubn7bfDhyCo/2/iOfvsrhgTbuXftq3K6TgikiSi7QOMMSbHGNPXGDPYGDMYGFAziYjEzFtvLaO8PMCxx7aluDjb6TjSxNq1y6R//yKqqoK8+eZSp+OISJKIdhi0y4AfgdeAR+pM/4lZMhER4JVXwt0fRo5M7N3jyay2ZX/KlEUOJxGRZBFtC/AdwE+stW2stZ3qTMlx+h4RccS2bVW8885yjAnvKpfENHx4V1wuw4wZKyktrXQ6jogkgWgLYA8wPZZBRER2NW3aUqqrQ5xwQhHt2mU6HUdipHXrdAYMaI/fH+L115c4HUdEkkC0BfDdwJ+MMVH3GRYROVA7T36h7g+JTt0gRKQ5NeZMcH8CttecBjkyxTCbiCSxzZsr+OCDVbjdhuHDuzgdR2Js+PCueDwuPvpoFRs3ljsdR0QSXLTjAF8U0xQiIrt47bUlBAIhBg/uQKtW6U7HkRjLy0tl0KBi3n13Ba+/voTLLjvc6UgiksCiHQf4w4amWAcUkeT0yivhXeHq/pA8at/r2q4vIiKxEu0waD5jzB3GmKXGmK01y041xlwb23gikozWrSvj449X4/W6OOssdX9IFmee2ZmUFBeffrqGtWvLnI4jIgks2j7AfwcOA8YCtmbZXODnsQglIslt6tTFhEKWIUMOJifH53QcaSbZ2T6GDu2IteHPgIhIrERbAJ8LjLHWfg6EAKy1a4CiWAUTkeRVOxKATn6RfM45R90gRCT2oi2Aq9nlgDljTGtgU5MnEpGktnr1dr744kfS0jwMG9bJ6TjSzIYN60R6uoeZM0tYtWqb03FEJEFFOwrEZOAJY8xvAIwx7YCJwHOxCiYiDSubMMHpCDHz6qvh1t+hQzuSmZnicBppbhkZXk47rRNTpixiypRF/OpXRzsdqckk8norEm+iLYBvBu4BvgPSgUXAv4G/xCiXiOxF9WWXOR0hZtT9QUaO7F5TAC9MqAI4kddbkXgT7TBo1dba6621mUAbIMta+xtrbXVs44lIMlm2rJTZs9eRmell6NCOTscRhwwZcjBZWSl8880GFi/e4nQcEUlA0Q6D1ssY8zNjzB+BkUDP2MYSkWRU2/p7xhmdSUuLdgeVJJrUVE9k+DudGllEYmGvBbAJe5Rw14ebgbOBW4BvjTGPGWNMM2QUkSRRe+S/uj/Iued2A2DKFI0GISJNb18twFcDg4B+1tqDrbXHW2s7AMcDA4CfxTifiCSJ+fM3MW/eJnJzfQwa1MHpOOKwgQOLyc9PZcGCzcybpwGHRKRp7asAvhj4lbX2q7oLa+avr7ldRJqZ++uv602JoHZX9/DhXUlJcTucRpzm9boZPrwrkDitwIm43orEq311susFfNjAbR8CTzVtHBGJRvbgwfXmt2ze7FCSpmGtjRQ5tbu+RUaO7MYTT3zPyy8v5Oab+xHvve4Sbb0ViWf7agF2W2u37+mGmuXRnkhDRKRB3323gSVLSmndOo0TT2zvdBxpIfr3L6JNm3SWLdvK11+vdzqOiCSQfbUAe40xJwMN/ezWYdoicsBefjnc/WHEiG54PPpdLWFut4sRI7rx8MPf8PLLCznqqDZORxKRBLGvAnY98Og+bhcR2W/WWo3+IA0aObI7Dz/8DVOmLOIvfzkRlyu+u0GISMuw1wLYWtuxmXKISJL66qu1rF69ncLCTPr2bed0HGlhjj22Le3bZ7F69XZmziyhX79CpyOJSALQvkYRcVRt6++553ZT657sxhijMYFFpMmpABYRxwSDIV59Ndz/V90fpCG1n41XXllEIBByOI2IJAIVwCLimM8+W8O6deV06pTDkUce5HQcaaGOOKI1XbrksmFDBZ9+usbpOCKSAFQAi4hjakd/GDmye9yP8SqxE+4GEW4Fru0yIyJyIKIqgI0xE4wxR8Y6jIgkD78/yNSpiwGd/EL2beTI8GfktdcWU10ddDiNiMS7aFuAvcDbxpjvjTG/N8ZopHoROSAffLCKLVsqOeSQfHr2LHA6jrRwPXoU0KtXAaWlVXzwwUqn44hInIuqALbWXgcUAn8AjgTmG2PeNcZcYozJjGVAEUlMr7yi7g/SOLXdIKZMWeRwEhGJd1H3AbbWBq21r1trRwP9gNbA48BaY8x/jDFFMcooIgmmsjLA668vAdT9QaJX+1l5440lVFQEHE4jIvEs6gLYGJNtjLnSGDMD+Aj4EhgA9AR2AG/GJqKIJJr33lvB9u3V9O7dmq5d85yOI3Gic+dcjjrqIHbs8PPuu8udjiMicWxfp0IGwBjzInAa4cL3n8Ar1tqqOrffAGyNSUIR2U3VJZc4HeGAvPRS7ckvNPavNM6553Znzpz1vPzyQoYP7+p0nEaJ9/VWJJFEVQADXwDXWmvX7ulGa23IGNOm6WKJyN6UT5zodIT9tm1bFW+9tRTQyS+k8c45pxv/8z+fMH36crZvryYrK8XpSFGL5/VWJNFE2wViwJ6KX2PMy7XXrbXlTZZKRBLWtGlLqawM0r9/Ie3bZzkdR+JM+/ZZ9OtXSEVFgDffXOp0HBGJU9EWwCc3sHxQE+UQkSQxefIPAPzkJ4c4nETi1U9+Et5zUPtZEhFprL12gTDG3FZzNaXO9VqdgRUxSSUiCWndujI+/HAVXq+LESM0+oPsnxEjuvGHP3zEBx+sZMOGclq3Tnc6kojEmX21ABfXTK4614uB9sAqYFRM04lIQpkyZRGhkOWUUw4mLy/V6TgSpwoK0hgypAPBoNWYwCKyX/baAmytvRzAGPOZtfbfzRNJRBLViy+q+4M0jVGjevD228t54YUFXH11b6fjiEicabAANsZ0tNYur5l9zxjTeU/3s9bqKASRZpaXn19vfsvmzQ4lid6SJaXMnr2OzEwvw4btcXMiErVhwzqRmell9ux1LFlSSpcuuU5H2qd4XG9FEtXeukB8V+f6YmBRzWXdSfueRCQqkycvAGD48K6kpUU7AqPInqWneznrrC7Azj0LIiLRarAAttZm1bnusta6ay7rTu7miSki8cxaq+4P0uRGjeoBhH9cWWsdTiMi8STqUyGLiOyvOXPWs3TpVtq0Seekk9o7HUcSxEkntadNm3SWLt3K7NnrnI4jInFkb32APwb2+ZPaWntSkyYSkYRTO17rued2x+3W725pGm63i5Eju/PQQ18zefIPHH10W6cjiUic2FtHvP80WwoRSViBQIgpUxYCMGqUuj9I0xo1qgcPPfQ1U6Ys5PbbB+Dx6AeWiOxbgwWwtfaJ5gwiIonpo49WsX59OV265HLkkQc5HUcSTO/erenWLY9Fi7bwwQcrOeWUjk5HEpE40OBPZWPMxXWuX9HQ1DwxRSRe1T34zRjjcBpJNMaYyJ4FnRpZRKK1ty4Qo4Gnaq5f3MB9LPBokyYSkYRRXu7n9deXAOr+ILHzk58cwp13fsG0aUspK/OTkeF1OpKItHB76wJxRp3rJzdPHBFJJG+9tYwdO/z06dOGzp1b/okKJD517JjDsce25auv1vLmm0s11J6I7FPURwsYY3KNMWONMTfVXOrbTET26vnnwye/OP/8Hg4nkURX+xl74YUFDicRkXgQVQFsjBkMLAd+BRwLXAcsN8YMiV00EYln69aV8f77K/B6XZx3Xnen40iCO+ecbng8LmbMWMmGDeVOxxGRFi7aFuAHgKuttcdZa8+31vYDrgIejF00EYlnkyf/QDBoOfXUjhQUpDkdRxJcQUEaQ4Z0IBi0vPTSQqfjiEgLF20BXAi8tMuyKYBGHReR3Vhree65+QBccEFPh9NIsqj9rNV+9kREGrK3USDqehL4JXB/nWU/r1kuIs0s0Lu30xH26rvvNjBv3iby81M59dSOTseRJDFsWCdyc318++0G5s7dyKGHtnI6Uj0tfb0VSSbRngrZBfzcGPM7YA1QBLQBvoh5QhHZzfYZM5yOsFfPPhs+EOm887qTkuJ2OI0ki9RUDyNHdufRR7/j2Wfnc/vtA5yOVE9LX29FkkljToX871gGEZHE4PcHeeml8AkJLrxQ3R+keY0e3ZNHH/2OyZMX8Oc/98fr1Q8wEdmdToUsIk3qvfdWsHFjBYcckq9TH0uz69OnDd2757Fw4Rbef38lp53WyelIItICNWYc4DbGmOHGmMt1KmQRacizz4YPQLrwwh469bE0O2MMo0eH9zw888w8h9OISEsV7TjA5wBLgNuAfxEeB/hfNHyKZBFJQlu2VPL228twuQyjRunkF+KMUaN64HIZ3nprGZs3VzgdR0RaoGhbgG8HLrfWHgWU1VxeDcyKWTIRiTsvv7yQ6uoQAwcWU1iY6XQcSVKFhZkMGlSM3x/i5Zc1JrCI7C7aYdA6WGsn77LsCWAtcGPTRhKRfck6+eR68y3l6PLa7g+1u6BFnDJ6dC/ef38lzz47n5/+tGUMP9ZS11uRZBRtAbzeGNPGWruO8CmQjwc2Ajq8VsQBnm++cTrCbhYu3Mzs2evIzPRyxhmdnY4jSe6MMzqTnZ3CnDnrmT9/Ez17FjgdqUWutyLJKtouEP8GTqy5/ndgBvAN8H+xCCUi8ae29XfEiG6kp3sdTiPJLi3Nw7nndgd0ZjgR2V1UBbC19m5r7Us1158EugNHW2vHxTKciMQHvz8YKYAvuuhQh9OIhNV2xXnhhQUEAiGH04hIS9KYYdDcxpgTjDGjgGJARxaICADvvLOc9evL6dYtj7592zodRwSAY49tS5cuuaxbV86MGSudjiMiLUi0w6AdASwCJgM31VwuMsa0jCMLRMRRTz01Fwi3/mrsX2kpjDGMGRNuBa79jIqIQPQtwI8CDwJF1tq+QBHwQM1yEUliP/64g3feWYHH4+LCCzX2r7Qso0f3wu0Ojwm8fn2503FEpIWItgDuDky01lqAmsv/BbrFKpiIxIfnnptPKGQZNqwTrVunOx1HpJ62bTM47bROBAIhHQwnIhHRFsDTgLN3WTYceKNp44hIPAmFLJMmhU83e/HFOvhNWqbaz+ZTT82lph1HRJJcg+MAG2OeAmq3FG7gOWPMLGAV4YPgjgZejXlCEWmxPvtsDcuWbaWwMJPBgzs4HUdkj4YMOZh27TJYsqSUTz9dw4kntnc6kog4bG8twIuBJTXT98CdwNvAvJrLOwEdVSCSxGoPLBozphdud9SDyog0K4/HxZgxvQAdDCciYQ22AFtr/9KcQUQkvpSWVvLaa4sBGDtWpz6Wlu3iiw9lwoSvmDp1MX/720Dy8lKdjiQiDmrMOMAnG2MeNca8XXM5OJbBRKRle/HFH6isDDJwYDEHH5zjdByRverQIZtBgzpQVRXkhRcWOB1HRBwW7TjAPwWeB9YCLwMlwDPGmKtimE1EWihrLU8/HT74TWd+k3hRezDck0/qYDiRZNdgF4hd/A4Yaq39pnaBMeZ54CXg37EIJiIt1+zZ6/j22w3k5aVy5pmdnY4jEpUzzuhMq1ZpzJ+/iVmz1nHMMTproUiyirYLRAHhg9/q+gHIb9o4IhKNLZs315ua26OPfgfAmDE9SU2N9ne0iLNSUtxceGG4v/qTT37f7K/v9HorIjtFWwB/AkwwxqQDGGMygHuBz2IVTERaptLSSqZMWQjAZZcd7nAakca56KLwaBBTpixi27Yqh9OIiFOiLYCvAQ4Hthpj1gGlQG/gZ7EKJiIt07PPzqeyMsigQcV06ZLrdByRRunePZ8TTiiirMyvg+FEktg+C2BjjAHSgFOAToTPANfJWjvQWvtjjPOJSAtireXxx8O7ji+/XK2/Ep+uuCL82X3kke90MJxIktpnAWzDW4fvgJC1drW1dqa1dnXso4lIS/Pxx6tZtGgL7dplcPrpOvhN4tOZZ3ahTZt0fvhhM59+usbpOCLigGi7QMwBuscyiIi0fLUHv1188aF4PDrzm8SnlBQ3l1xyGACPPPKtw2lExAnRHr79AfCWMeZxYBUQ2WdkrX206WOJyN6kX399vfnyiRNj/ppr15YxbdpS3G4TKR5E4tWllx7GhAlf8cYbSykp2UG7dpkxf00n1lsR2bNoC+ATgGXAwF2WW0AFsEgz8z35ZL355vgiffrpuQQCIc46qwuFhbEvFkRiqbAwkzPO6Mxrry3hySfn8vvfHxfz13RivRWRPYtqH6a19uQGJp0OWSQJBIMhnnhCB79JYrnyyiOA8JjAfn/Q4TQi0pz2WgAbY9KNMXcaY6YaY241xviaK5iItBxvvbWMNWt20KlTDgMHFjsdR6RJDBjQnm7d8igpCXfvEZHksa8W4AcID3u2APgJMD7miUSkxfnXv74G4Kc/PQKXyzicRqRpGGMircC1B3iKSHLYVwF8OnCqtfZ3NdfPin0kEWlJ5s7dyCefrCEz08vYsb2cjiPSpC68sAcZGV4+/ng1CxZscjqOiDSTfRXAGdbaEgBr7SogJ/aRRKQl+ec/w62/Y8b0IjtbvaAksWRn+xg16hAA/v1vDYkmkiz2VQB7jDEnG2MGG2MG7zpfs0xEEtTGjeW8+OIPGANXXdXb6TgiMXH11eHP9nPPzWfz5gqH04hIc9jXMGjrqT/M2aZd5i2g00GJJKgnn5xLVVWQU0/tSJcuuU7HEYmJHj0KGDy4A++/v5Innvie3/zmWKcjiUiM7bUF2Frb0VrbaS+Til+RBOX3ByNnyfrZz450OI1IbP3850cB8J//fEt1tYZEE0l0OpepiOzR1KlLKCkp45BD8hk0SEOfSWIbPLgDhxyST0lJGVOnLnY6jojEmApgEdmj2qHPfvaz3hijoc8ksRljuOaa8J6Ohx6ag7XW4UQiEksqgEVkN//971r++9+15OT4GDWqh9NxRJrF+ef3ID8/lTlz1vPllyVOxxGRGFIBLCK7+cc/ZgFw2WWHkZHhdTiNSPNIS/NETvX9f/83x+E0IhJLKoBFpJ4lS0p5/fUlpKS4dPCbJJ0rrzwCr9fFtGlLWbFiq9NxRCRG9jUMmoi0QNvefz9mz/3gg7OxNrw7uG3bVnHUAgAAIABJREFUjJi9jkhL1LZtBiNHduf55xfw0ENf87e/DWyy547leisijaMCWCQOBY+MTcvs+vXlPPvsfACuvbZPTF5DpKW79to+PP/8Ap56ai433dSXgoK0JnneWK23ItJ46gIhIhH//vc3VFUFOf30TnTvnu90HBFHHHpoK4YO7UhFRYCHH/7G6TgiEgMqgEUEgB07qiMnvrjuuqMdTiPirN/85hgg/KNw+/Zqh9OISFNTASwiADz99DxKS6vo27cd/foVOh1HxFH9+hVy3HHtKC2t4sknv3c6jog0MRXAIoLfH4wM+3Tdder7KwJw/fXhVuD/+785VFUFHE4jIk1JB8GJxKGUxx+vN1992WUH9HwvvPADq1dvp1u3PE4/vfMBPZdIohg6tCM9exYwf/4mXnjhBy6++NADer6mXm9FZP+pABaJQxk33FBv/kC+SAOBEH//+1cA3HDDMbhcOu2xCIDLZfj1r4/mmmum849/zGLMmJ643fu/47Qp11sROTDqAiGS5KZMWcjSpVvp1CmH8847xOk4Ii3KyJHd6dAhm8WLS3nttSVOxxGRJqICWCSJBYMhxo8Pt/5ef/0xeDzaJIjU5fG4Iv3i7713JqGQdTiRiDQFfduJJLGpUxezaNEWiouzuOCCHk7HEWmRLrqoF4WFmcyfv4nXXlvsdBwRaQIqgEWSVChkue++na2/KSluhxOJtEw+n4cbbgiPCHHPPWoFFkkEKoBFktS0aUuZN28T7dplMGZMT6fjiLRoY8fubAWeOlWtwCLxTgWwSBIKhSz33vslAL/+9TH4fBoQRmRv6rcCf6lWYJE4pwJYJAm9+uoivvtuI+3aZRzw2KYiyWLs2F4UFWWyYMFmtQKLxDkVwCJJJhAIcdddXwBw0019SUtT669INMKtwMcCagUWiXcqgEWSzPPPL2Dx4lI6dsxm7NheTscRiStjxvSMtAK/9NIPTscRkf2kAlgkiVRVBbj77nDf3z/+sR9er0Z+EGkMn8/D739/HAB33vkF1dVBhxOJyP5QASySRJ54Yi6rV2+nZ88CRo7s7nQckbh04YU9OeSQfFas2Mbjj3/vdBwR2Q8qgEWSRFmZn/vumwnALbf0w+3W6i+yPzweF+PGHQ/A+PEz2b692uFEItJY+gYUSRIPPTSHDRsq6NOnDaef3tnpOCJx7fTTO9O3bzs2bqzgwQdnOx1HRBpJh3+LxKGyCRMadf+1a8v43/+dBcCtt56AMSYWsUSShjGGW289gTPOeJEHH5zDFVccwUEHpe/1MY1db0UkdlQAi8Sh6ssua9T977rrC8rK/JxxRmdOPLF9bEKJJJl+/QoZNqwTb721jPHjZ3LPPYP2ev/GrrciEjvqAiGS4ObO3cikSfPweFzceusJTscRSSjjxvXH5TI89th3LFiwyek4IhIlFcAiCcxay7hxHxMKWa644nC6ds1zOpJIQunZs4BLLz2MYNByyy0fY61OjiESD1QAiySwd99dwQcfrCInx8fvftfX6TgiCemPfzyOnBwfM2asZPr05U7HEZEoqAAWSVB+f5Bx4z4G4Le/PZb8/DSHE4kkplat0iMnx7jllo90cgyROKACWCQOub/+ut60J//61zcsXLiFTp1yuOqqI5o5oUhyufLKw+nWLY+lS7fyr399s8f7RLPeikjz0CgQInEoe/DgevNbNm+uN//jjzu4557wKY/vvnsgPp9WdZFY8nrd3HHHAM4/fyrjx8/kggt67DYs2r7WWxFpPmoBFklA48Z9zI4dfs46qwunnNLR6TgiSeGUUzpy6qkd2b69mj//+ROn44jIXqgAFkkwH364iilTFpGW5uGOOwY4HUckqdx550n4fG6ef34BH320yuk4ItIAFcAiCaS6OsjvfvcBED7wrbg429lAIkmmc+dcbrwxPOLKjTfOoLIy4HAiEdkTFcAiCeSBB2azaNEWunTJ5Ze/PMrpOCJJ6brr+tC9ex6LF5cyceJ/nY4jInugAlgkQSxcuDly4Nv48SfrwDcRh6SkuJkwIXzA28SJ/2XRoi0OJxKRXakAFkkAoZDlV796j+rqEBdd1IuBA4udjiSS1Pr3L2LMmF5UV4e48cYZOkOcSAujAlgkATzyyLfMnFlC27YZ/PWvOvBNpCX4y19OoKAglY8/Xs1jj33ndBwRqUMFsEicW0EOt932GQD33juInByfw4lEBKCgII177z0ZgD//+VOWkudwIhGppQJYJI6FMFzF2ZSV+RkxoitnntnF6UgiUsc553TjnHO6UVbm5wpGEMI4HUlEUAEsEtceoC/v0IX8/FTuvnuQ03FEZA/uvXcQrVun8SEdeZBjnY4jIqgAFolbc2nN7zkFgIkTh+x22lURaRkKCtK4777wqBC/ZyiLyHc4kYioABaJQ9W4uYiRVOLlCmZz1lnq+iDSkp11VhfG8i0VeBnDeVTjdjqSSFJTASwSh/6Hk/madnRmMxN5y+k4IhKFB5jGwZTyX4q4hcFOxxFJaiqAReLMrFml3MMJuIzlkWEVpFxyodORRCQKaZeczxOnl+E2lvGcwLvvLnc6kkjSUgEsEkd+/HEHd931AxbDjTcdR+9n7qZ84kSnY4lIFMonTuSISXdz85/6A/DLX77D5s3VDqcSSU4qgEXihN8f5Kc/fYvSUj8DBxZz0019nY4kIvvh178+hoEDi9mwoYK77vqBQCDkdCSRpKMCWCRO3HHHF3zxxY8UFKTw8MOn4XZr9RWJRy6X4aGHTqV16zTmzNnK7bd/7nQkkaSjb1CRODBt2hLuv38WbrfhT386hNatNeSZSDxr2zaDRx89HZcL7r9/Fq+8ssjpSCJJRQWwSAs3b95GrrlmOgDjxvXn8MNzHE4kIk3hhBPac801nQC47rp3mTdvk8OJRJKHCmCRFmzjxnLGjHmdHTv8nHded667ro/TkUSkCZ17biGjRh1CWZmfSy55ndLSSqcjiSQFFcAiLVR1dZBLL53GypXb6NOnDffffwrGGADy8vPrTSLS8u1pvTXG8Pe/D+bww1uxdOlWLrnkDaqqAg4nFUl8KoBFWiBrLb/97Qw+//xH2rXL4KmnziQtzeN0LBGJgfR0L5MmDadt2ww++WQN11//PtZap2OJJDQVwCIt0F13fcmkSfNITXXz1FNn0a5dptORRCSG2rfP4tlnh5OR4eX55xdw991fOh1JJKGpABZpYR599FvGj5+Jy2V49NHT6dOnjdORRKQZ9O59EP/5zzBcLsM998zk6afnOh1JJGGpABZpQaZOXcxNN30AwIQJJzNsWGdnA4lIszrttE7cffdAAK6//n0NjyYSIyqARVqId99dztVXv4W1cPPN/bjkksOcjiQiDrjyyiP4/e+PIxSyXH3127zzznKnI4kkHBXAIi3AjBkrufjiN6iuDnH11b357W+PdTqSiDjod7/ryy9/eRSBQIhLL32DTz5Z7XQkkYSiAljEYR99tIqxY1+jqirIFVcczl13nRQZ7kxEkpMxhttuO5HLLjuMysogF144lY8+WuV0LJGEoQJYxEEffLCS0aNfo7IyyCWXHMo99wxS8SsiQLgIHj/+ZEaP7kl5eYALL5yq7hAiTUQFsIhDXn11ERdcMJWKigAXXdSLCRMG43Kp+BWRnVwuwz/+cQqXX344lZVBLrrodV5/fYnTsUTingpgEQc8/vh3XHHFm/j9Ia655kgmThyi4ldE9sjlMowfP4if//xI/P4Ql18+TUOkiRwgFcAizchayz33fMkNN8zAWvjTn47njjsGqPgVkb0yxnD77QO48ca+BIOWX/3qPe6443OdMU5kP6kAFmkmFRUBrr76bf72ty8xBu6772RuuOFY9fkVkagYY7j55n5MmDAYt9tw331fcc0106mqCjgdTSTuqAAWaQZr15YxfPhLvPTSQjIzvUyadBaXX36407FEJA5ddtlhPPNM+LTJkyf/wIgRUygp2eF0LJG44nE6gEii+/LLEq64YholJWUUF2fx7LPD6dWr1QE9Z6B37yZKJyLNpSnX26FDO/LGG+cxevRrzJxZwuDBz/HYY2fQr19hk72GSCJTASwSI6GQ5YEHZvPXv35GMGg57rh2PPnkmbRunX7Az719xowmSCgizamp19sjjjiIGTNGc+WVb/Lpp2s4++yXue22E/nZz3qra5XIPqgLhEgMbN5cwdixr3HrrZ8SDFquvbYPU6eObJLiV0Sk1kEHpTNlyrn84hfhs8bdfPNHjB79GuvXlzsdTaRFUwEs0sSmTVtC//6TePvt5eTm+njmmbO47bYT8XrdTkcTkQTk8bi4/fYBPPbY6eTk+Jg+fTknnjiJt99e5nQ0kRZLBbBIEyktreTnP5/ORRe9wfr15fTvX8iHH45m2LDOTkcTkSQwYkQ3PvlkDCed1J6NGysYPfo1rr32XTZvrnA6mkiLowJY5ABZa5k8eQHHH/80zz+/gLQ0D3feeRJTp55HcXG20/FEJIkUFWXx8svnctttJ5KS4uKZZ+bRr9/TvPDCAo0ZLFKHCmCRAzBv3iaGD3+Jn/1sOuvWldO3bzs+/HA011xzpE5uISKOcLkM117bh48/HssJJxSxcWMF11wznZEjX2H+/E1OxxNpETQKhMh+WLu2jPvum8njj39PMGgpKEjl1ltPZPTons1S+GadfHK9eY0KIdLyNfd6261bHlOnjuSZZ+YzbtzHfPjhKgYMeIaLLz6UP/zhONq0yYjp64u0ZCqARRph69Yq7r9/Fv/619eUlwdwuQxXXnk4t9xyPLm5qc2Ww/PNN832WiLSNJxYb40xjB3bi9NO68g998zksce+44knvuell37gF784imuuObJZt10iLYW6QIhEYePGcu688wuOPPJx/v73/1JeHuDMMzvz8cdjuPfek/UFIiItWqtW6dxzzyA+/XQsw4Z1YscOP/fcM5PevR/njjs+14FyknTUAiyyF8uXb+XBB2czadI8KiuDAJx4YhHjxvXn2GPbOZxORKRxunfP55lnhvP552u4556ZfPjhKu677yv++c+vGTu2Fz/96RF07ZrndEyRmFMBLLKLQCDE9OnLeOKJubz77nJqD5w+7bSO/PrXx+hUoyIS944/vogpU85l5swSxo+fybvvruDhh7/h4Ye/YciQg7nqqiMYMuRg3G7tKJbEpAJYpMbChZt58cWFTJo0l5KSMgBSUlyMHNmda689ml69ChxOKCLStPr2bccLL4zg++838O9/f8vkyQt4770VvPfeCgoLMzn//EO44IKeHHJIvtNRRZqUCmBJaitWbGXKlEW8/PJCvv9+Y2R51665XHLJYYwe3ZOCgjQHE4qIxN5hh7Xmf/93CH/+c3+efnoejz32HStWbGPixFlMnDiLo446iHPO6c6ZZ3amc+dcp+OKHDAVwJJUAoEQX31VwjvvLGf69OXMm7dzTMzs7BTOOqsrF1xwCCee2B5jNI6viCSX/Pw0fvWro7nuuj58+WUJzz03nylTFjFnznrmzFnPn//8CT165HPmmV047bROHHnkQXg86iYh8UcFsCS0YDDE3Lkb+fzzH/nsszV89NFqtm6tityekeFl2LBOjBzZncGDO+DzaZUQETHG0K9fIf36FXLXXQOZPn0Zb7yxlOnTl7FgwWYWLNjMffd9RVZWCieeWMRJJxUzYEAxPXrk6yRAEhf0bS8Jw1rL2rVlfPvtBr75Zj3//e9avvyyhO3bq+vdr2vXXIYO7cjQoR05/vhCFb0iInuRluZhxIhujBjRjerqIJ99toZp05by/vsrWLp0K2++uYw331wGQE6Ojz592nD00W05+ug2HH10G1q1Snf4LxDZnb75JS5t2VLJokVbWLw4PH3//Ua+/XYD69eX73bfgw/O5vjjwy0ZAwa0p1Mn9V8TEdkfKSluBg3qwKBBHQBYvXo7H320io8+WsXHH6+mpKSMGTNWMmPGyshj2rXLoEePAnr23Dl165ZHVlaKU3+GiApgaZm2b69mzZrtrF69nTVrdtRM21m2bCuLF29h48Y9D9qenZ1C794HcfjhrTnqqIPo16+QoqKsZk4vIpIc2rfPYsyYXowZ0wuANWu2M2vWOmbNWsusWWv5+uv1lJSURQrjulq1SuPgg7Pp1CmXgw/OpmPHHDp0yKJNmwzats0gKytFx2JIzKgAlkbxTp5M2l//imvNGkJFRVSMG4d/1KgG7x8KWcrK/GzfXl1v2ry5go0bd06bNu28XLeujG3bqht8Tgj33e3SJZeuXfPo2jWXnj0LOPLIg+jQIbvBDaZ38mTSb7gBU1a2c6HPB1U7+wTjclF12WVUjB+/8zF/+ANmyxYAbEYGpKSE510uCIUiD7V5eZT/7W/4R43a+X9avRrcbggGCbVvj//UU/FOn15vee2lzcjAVFTUf86MDEx5OZHBiIHBDfxPli1bxmeffcayZcsoLS2loqICn89HcXExXbp0Ydu2bSxYsICVK1eyZcsWrLX4fD6stXi9XoqKijjooINYvHgxy5YtY9u2bQSDQdxuN7m5uaSlpbFhwwY2bNiA3+8nLS2NgoICvF4vGzZsoLS0lMrKSgDcbjcejwe32421lmAwSCAQIBgM1vnX+/B4PASDQUKhEKFQiGAwiK3zt9aVkpJCeno61lqqq6uprKyM3Dc9PZ2cnBx27NjB9u3bG/gPNS+32w1Q728+0Oez1hKq8/k4ED6fD4Cqup//A+DxeAgEAo16jNfrJTMzk0AgQGVlJcFgEJfLFfnMQPh9d7vdBAIB/H4/xhhSU1Mjy30+H8FgkIqKCqy1pKen4/P5CAQCWGsj9w2FQpHnLS8vj2Stfazf78flcuHz+cjPzycvL4+2bdtSXFyMtZatW7dSVVVFampq5H+XkpKCy+VizZo1rF+/Ho/HQ5s2bfD7/fz444/4/X7S08O7/8vLy0kFioBLgRFAbn5+g+vzvtiUFMjMxGzZQqioCP+pp5IyZUr9bRVEtnd1t0/13oM623SbmwvV1fW3kQDGhLdBdbZlDW37i4qyKCrK4uyzu9b8f0OsXLmN+fM3MX/+ZubP38S8eRtZtmxrZPs/a9a6Pf6N6eke2rYNF8Nt2mTQunU6OTk+8vJSycsLX+bmptZc+sjKSsHnc+930dzY7zeJb6ahL5um8vnnn9sePXrE9DV2FQyGCARCkZph56WtN7+nZbv+Pxp67L7nd3++aB+7r2x7e825c+dy6KGHEgpZAoFQ5DIY3H0+GNx5GQjsvAyF6s6HqKwMUlUVwD/rW4JvvUNVACrxhCe3jx2HH0VFQVuqqgJUVATYvr2aHTv8NZfVu/1d0UhNdVNUlEX79lkUFmbSvn0WRUWZHHxwDl275lJYmNmojZx38mQyfvELTBTFiAWqrriCwHHHkXHttRi/P+rXsV4vVRdfjO/ZZ8PF7B6eOxbtGbNnzeKFF17A5/OxZMkS5s+fj9vtplu3bqxbtw6v14vb7Wbbtm38+OOPGGMoKysjEAjg8XgoLCzEWsuKFSvIz89n27ZtlJWVUVVVFSks6xZytZ+92sJVxAl1i26Xy4UxJvI59Xq9uFwuAoEAxphIEVxdXY0xZp8/JrKyssjJycFaS0ZGBmlpabhcLrZu3YrP54tMq1evJhgMkpOTQyAQYN26dRhjaN26NVu3bqWsppjMycmBdesYDdwNNPXO/2i2LdbrpeyBByJFnXfyZDKuv36P26p9vl5aGmUTJ+53gRgKWUpKdrBixTaWL9/K8uXhyzVrtrN2bRnr1pVRXt64H1QAbrchPd1LZqaXjAwv6enhy9rr6ekeUlLc+HxuvN7aSxdpC+eT9fqr+AKV+AiQQpAUr5vQFZdgTzwBt9vgdrtwuUzN9V3nXZHlLperzn0Mc+fO4/DDD8UYgzHU++6qO197ffdle75PeL7+9T3dJxqN+T6N9q5er8uRE6nMnj171pAhQ47Z1/0SsgB+7rn5/OIX7zTra0rDMjO9ZGWlkJmZQlZWeMrLS6VVqzQKCtIil61b116mk5+f2qS7vrKPOAL36tVR39+63YTatWvUY+o+NppCuyk98I9/UFZWxpYtW5gzZ07ki72srIwOHTqwePFirLWkpKRQWlqK3+9nx44d+P1+srOzyczMZNOmTRhjqK6uxlpLIBAgEAhEnqt2eW3LbnV1dZO1boo0BbfbHflMejyeeq3JKSkp+P1+QqFQpEV5b3w+HxkZGXg8HtLT00lNTY08f2VlJQcddBDr16+nvLyctLQ03G43VVVVlJaWYowhLS0tct9QKERqair5paV8RLgV2CnB9u3Z9u23QOO3i3t7rqZmrWX79mrWri2LFMSbNlWyZUslW7dWsWVLZWQqLQ3P79hRTXW1fpC3FA89NJQLLujZ7K8bbQGckF0gXC5DSkr4V0fdX0R7mm9oWTSP3X1+7/fbV46meM3q6ip8Ph/GhH95ejyuyKXLFe28C48n/CvW4zGkpnpITfWQc/8E0vCTSqDe5CNI4PlJpKa68fk8kSI3OzuFjAxviziVpmvNmsY9IBhs/GPqPLa5WKB082bW3ncfRUVFrFmzhoqKinBrE7B27VpSU1Mj3QVqv/zLy8trooZ3OVdVVUUet23bNjweT6TYLS8vj3zxG2MizxHrH88ijVX3M1m3dbf2s1+7LJrPbnhPWLj12O/34/V6I599v9+PtTbSbaO2oK6uro48trq6OlKAA+EuFl270m7x4qb8kxut7nZtv7dxTfT4vTHGkJ3tIzvbR/fu0Z+Fzu8PUl4eoKzMH5nKy8OXO3ZUU1ERoLo6RHV1MDJVVQVx3TuBatxU4a659FCFGz9uKk8/o2ZvaXjPaChk680Hg7ZmWSiyvO58VVUVXm8K1tbfc7vrHmJr7R7vU3cvb/15u9fniVZjtuWNed6WPhxeQhbA55/fg/PPb95W55Zi9uzZ9OnTJybPnf3yz/fYWhBs355tQzvG5DWbSqioqHEtHQfQAhzp29scavqZtm3blm3btpGenk5aWlqkL25WVhaVlZWkpqZGWoArKyvxer1UVVXhdrsJhUKRx/n9/sh9a1uBPR5PzUuFi16Xq/6PS5GWou5n0uVyRQpQY0xkPhQKRfXZDTcAhPd21HYhqv0h6PV6I32Ra/uuezweUlJSKC8vx+VykZKSEvnRWNvPPhAI8KPbTXsH95yEiorqXT+QFuC6z9VSeL1ucnLc5OT4GvW47GcXNPz9NumhA8oUy+9l2X/ON81J3KgYNw6bVv+0wDYtjYpx4xxKFL2KceOwNcXivlig6tJLw4/xehv1OtbrperSS3f7P9V97qZSmxOgf//+bNmyBZ/PR+vWrdm8eTNbt26lU6dObNy4kdatW9OuXTu8Xm/kQKJwP7Xwblufz0e7du2oqqoiLy8vcnBadXU1GRkZkcK3tgiu+xwiTqn9cQbs9lkM94EM77GoLX5rC9e6P+QaUnvQpdfrJSUlBZ/PR0pKCtXV1aSlpVFdXU2bNm0iPya9Xi+pqamR183Kyoq8PoS7VOzYsYPxeXns/RDf/RPNtsV6vfW213vapkf9enGy7Y9WPH+/yf5x33rrrTF9gdWrV9/aqlWrmL6G7FRSUkK7du1i8tyhQw8l2KED7q+/xmzfTqh9e8rvuisujpINHXoowU6d8MyYAXUPavP56rfWulxUXX45FePHhx/TsSOeTz+FmhZVm5GBzcwMz7tc9fYH2bw8yidMoOo3v9n5f9q2LdxKay2h9u2p/slPMBs21Ftee2kzMsJZ6j5nRgbsqa+i2x3JCZCXl0dRURGbN2/G7/eTm5tLdnY2Pp+Pww47jJNOOoni4mLcbnfk6Pnc3FwKCwspKCggMzOTbt26MWDAALxeb6SVKzMzk7S0NDp06EDHjh0jR9e7XC6ys7MpLi6msLAQl8tFMBiMtHjVtob5fL7IwUhQf1ebz+eL9KusPYBpb7viUlJSyMrKwufzRV6vVnp6Oq1atYr0YW4Jav+upuoq4na79/k/aoy6o3A0BY/H0+gDIr1eLzk5OXi93sjftevfmZKSQkpKSqSYrO1fW3tgWm5uLunp6ZHPXW5uLvn5+aSmppKamkpBQQGtWrUiMzOTrKwssrOzI0VtWloaKSkp9boHpaWl0aZNG9q1a0enTp3o06cPxcXFpKenk56eTnFxMe3bt6d169a0bduWbt26kZ6eHvmxeNhhh1FUVBQZMaKwsJDWrVvj8XjIycmhols3vD160L2kBPcBfFZtSgo2JwcqKyPbFtfy5fW3VSkpke1d7fap7vZ61226zcsLNxTseuBvbYt5nW1ZvGz7oxXL77dYfi/L7kpKSko6d+788L7ul5AHwSUz7WpJfHqPk4Pe5+Sg9znx6T1uXtEeBKf9lyIiIiKSVFQAi4iIiEhSUQEsIiIiIklFBbCIiIiIJBUVwCIiIiKSVFQAi4iIiEhSUQEsIiIiIklFBbCIiIiIJBUVwCIiIiKSVFQAi4iIiEhSifmpkN97770NwIqYvohEbN68uVV+fv5Gp3NI7Og9Tg56n5OD3ufEp/e42R08ZMiQ1vu6U8wLYGlexpj/Wmv3eQ5siV96j5OD3ufkoPc58ek9bpnUBUJEREREkooKYBERERFJKiqAE8/DTgeQmNN7nBz0PicHvc+JT+9xC6Q+wCIiIiKSVNQCLCIiIiJJRQVwgjLG3GiMscaYVk5nkaZnjLnXGLPAGPOtMWaKMSbX6UzSNIwxw4wxPxhjFhtj/uB0Hml6xphiY8wMY8x8Y8xcY8yvnc4ksWOMcRtj5hhjXnc6i+ykAjgBGWOKgaHASqezSMy8AxxmrT0CWAj80eE80gSMMW7gQeB0oBcw2hjTy9lUEgMB4LfW2p5AP+CXep8T2q+B+U6HkPpUACemvwO/A9TBO0FZa6dbawM1s18A7Z3MI02mL7DYWrvUWlsNPAeMcDiTNDFrbYm1dnbN9e2Ei6MiZ1NJLBiHBbaoAAAFd0lEQVRj2gNnAv9xOovUpwI4wRhjzgbWWGu/cTqLNJsrgDedDiFNoghYVWd+NSqMEpoxpiNwFPCls0kkRiYSbpAKOR1E6vM4HUAazxjzLtB2DzfdAtwMnNq8iSQW9vY+W2tfrbnPLYR3p05qzmwSM2YPy7QnJ0EZYzKBl4DrrbXbnM4jTcsYcxaw3lo7yxgzyOk8Up8K4DhkrT1lT8uNMYcDnYBvjDEQ3i0+2xjT11q7thkjShNo6H2uZYy5FDgLGGI1nmGiWA0U15lvD/zoUBaJIWOMl3DxO8la+7LTeSQmTgDONsacAaQC2caYp621FzmcS9A4wAnNGLMcOMZau9HpLNK0jDHDgAnAQGvtBqfzSNMwxngIH9Q4BFgDfAWMsdbOdTSYNCkTbqF4Athsrb3e6TwSezUtwDdaa89yOouEqQ+wSHx6AMgC3jHGfG2M+afTgeTA1RzYeC3wNuEDo15Q8ZuQTgAuBgbXrL9f17QSikgzUQuwiIiIiCQVtQCLiIiISFJRASwiIiIiSUUFsIiIiIgkFRXAIiIiIpJUVACLiIiISFJRASwiIiIiSUUFsIiIiIgkFRXAIiIJzBjjczqDiEhLowJYRCQKxphBxhjbwPT4Xh7X1hjzhDHmR2NMlTGmxBjzujHmoDr3yTDG/M0Ys6TmPmuNMS8ZY9rUuU9fY8y7xpgdxpgyY8x7xpi+u7zW48aY1caY440xnxljKoB76tx+lTHmG2NMpTFmozHmEWNMftP+p0REWj6P0wFEROLEbOD4XZadDNxJ+LTFDXkKOBi4CVgFtAGGAOkAxpgU4B3gSOAu4AsgBzgNyAPWGWOOAD4E5gGXARb4A/ChMaaftfabOq+XAzwHjAduBipqXudvwG+B+2uyFAG3A4cZY/pba4ON+m+IiMQxnQpZRGQ/GGO6AV8C7wOjbAMbU2PMDuBma+39Ddx+BfAIMMJaO7WB+7wInAJ0tNaW1izLBpYDH1hrR9Ysexy4FDjHWvtqncd3BJYAf7HW3lZn+QnAJ8C51tpXov3bRUTinbpAiIg0kjEmD3gdWAxcbK21xhi3McZTZzI1d/8KuMkY82tjzOF1ltc6FVjbUPFb4yTg9driF8Bauw2YCgzc5b6Bmmx1DSW8vZ9UNyPhAn5bzfOLiCQNFcAiIo1QUzi+CKQCZ1trK2puWgL460yX1iy/gHCh+jvgW2CNMeZ/jDG1298CYM0+XjYfKNnD8rWEu0nUtX4P3Rlq+xsv3iWjH8iuySAikjTUB1hEpHEeBI4FTrDWrq2zfDhQd8SFZQDW2vXAL4FfGmMOIVwY/wXYADwEbAQO28drbgba7mF525rb6tpTV4xNNZenAlv2cruISFJQASwiEiVjzG+AKwm3/H5X97Zd5/fEWvsDcLMx5hp2Fr3TgQuNMcOtta818NAPgTONMVnW2u01WbIIF90fRBH9HSAEdLDWvhPF/UVEEpoKYBGRKBhj+hMeWeFJYLMxpl+dmzdYa5fs4TE5wLvAJGAB4S4HIwh3W5hec7engauAZ40xdxHul5tFeBSIidbaBcBfgbOA94wxdxNu5f094ZEkIge1NcRau6TmcQ/UtEJ/CFQCxYT7B//HWjujEf8OEZG4pgJYRCQ63QkfN3FZzVTXE3tYBuEiczbhAvdgwq2wPwBja0dpsNb6jTGnAn8Grq653AR8Sk33Bmvtt8aYQcAdNa9lCA+XNnCXIdAaZK292Rgzn5ruGISL6FXAe8CiaJ5DRCRRaBg0EREREUkqGgVCRERERJKKCmARERERSSoqgEVEREQkqagAFhEREZGkogJYRERERJKKCmARERERSSoqgEVEREQkqagAFhEREZGkogJYRERERJLK/wOqf2v/MJXBzAAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "\n", "plt.plot(x_random_nonsig, y_constant_nonsig, 'ko', label = 'Non Significant Values', alpha = 0.2)\n", "plt.plot(x_random_sig, y_constant_sig, 'ro', label = 'Significant Values')\n", "plt.plot(x, y, '-', color = 'navy'); plt.xlabel('z-score', size =16); plt.ylabel('Probability Density');\n", "plt.vlines([stats.norm.ppf(0.975), stats.norm.ppf(0.025)],linestyles='--',\n", " color = 'red', linewidth=4,\n", " ymin = 0, ymax = y.max(), label='Significance Limits for $p = 0.05$')\n", "plt.yticks([])\n", "plt.legend(prop={'size': 14});\n", "plt.title('Random Observations with significance level $p=0.05$', size = 22);" ] }, { "cell_type": "code", "execution_count": 20, "metadata": {}, "outputs": [], "source": [ "# Non significant random values at p = 0.05 with correction\n", "x_random_nonsig = [x for x in x_random if abs(x) < stats.norm.ppf(1 - (0.025/n_tests))]\n", "y_constant_nonsig = [0.0 for _ in range(len(x_random_nonsig))]\n", "\n", "# Significant random values at p = 0.05 with correction\n", "x_random_sig = [x for x in x_random if abs(x) > stats.norm.ppf(1 - (0.025/n_tests))]\n", "y_constant_sig = [0.0 for _ in range(len(x_random_sig))]" ] }, { "cell_type": "code", "execution_count": 21, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAsAAAAGOCAYAAACZlbVuAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvNQv5yAAAIABJREFUeJzsnXlcVUX7wL8DirK4gBqigCgqKu77lhvoT0s0c88NyzRzy1xKe11yeV9Ns3LJpSw1QYvc99RCRHEljRLNVBT3LTdQQTi/P865l3sv97LDVZnv53M/cOaZ5TkzZ+Y8Z85zZoSiKEgkEolEIpFIJPkFG2srIJFIJBKJRCKR5CXSAJZIJBKJRCKR5CukASyRSCQSiUQiyVdIA1gikUgkEolEkq+QBrBEIpFIJBKJJF8hDWCJRCKRSCQSSb5CGsB5gBAiRgihmPyeCCEuCSF+EkK0tLaOGUEI0UrTPdTaumQXIUQpIcSnQogjQog7QoinQoirQojNQoi+QgiLfUPXhnmp78uKEGKFVp+B1tbleUIIEajVy4ospM3xOhVCTDUzhunGsX+EEMuEEBVzqrzsIoQYJIQ4LoSIM9C1uLX1ymuyO1YJIZppbRsthLgvhEgQQlwXQuwRQowRQrySk/q+LGSn/0ryjgLWViCfsQu4rv3vDNQBugPdhRAfKoryhdU0y0cIIXoAywEn4B5wEHgAeAKvAQHAaCFEZ0VRLltN0RccIUQr4Ddgn6IorayrzcuBEMILuABcVBTFywoqnAPCDY5LAvWBd4E+Qoi2iqIctIJeeoQQHYFvgCfAbuCuJkqwmlIvGEIIJ+B7oJsWdAEIBeKA0kATwA+YKoTwVxTlsDX0tBa6hwpFUYS1dZFkHWkA5y2zFEUJ1R0IIQoCXwFDgVlCiBBpcOUuQohuwFrtcBLwmaIoCQZyT+A71MF9vxCinqIod1PnJMkhJgCzgGvWVuQ5YwNwCLhvbUVMCFcUJdAwQAhRGFgJ9AAWoT7YW5Pu2t+RiqJ8Y1VNrE/VzCYQQtihTtY0BU4DgxVF2W8SpzDQH/gUcMsBPV82ntf+KzFAukBYEUVREoExwEPADmhnXY1eboQQpVBnfgUwVlGUGYbGL4CiKJeA14HDgBewIK/1zE8oinJNUZTTiqLIG4UBiqLc1+rluX8wUBTlCfCJdlhbCFHMmvoAHtrfs1bV4jlAu4ZOZzLZZFTjNwZoZmr8avk+URRlGVAbiM62oi8ZL1L/zc9IA9jKKIryGPhbO3Q1lQshygkhJgghfhNCxGq+qne147fM5WnoqyuEKCiE+EQIcVrz17sphFitzXSaRQjxhhDigOY/968QYndG/JSFEE2FEOs0HzGdr9jPQojGFuKHanq20nzNdmrl3RdC7BJC1DaI218IcVQI8Ug7/9VCiNLp6WTCcKAoEAVYdDdRFOWpFheglxDCO41zHiyE+F0IEa/5Eq8XQlS3ENdHCLFSCHFRq5+HQvUP3yCE6GohTSMhxFohxGUtzS2h+ik3txBf7/MnhHhHCHFYCPFA5wMphHis5VMyjXM6psV/3SCsmhBimhDioFB9pXW6bBdCtDeTRyiq+wNAS2HsNxpqEM+iv6pQ6addJ/9q1+85IcQiIYSHaXwz599TCBGhXTMPhRB706i3TLeNmTzGaOXPNiOL1GSpXhULIT7TZKMNwlL5EGr/X9AOy5nUaYwFnSoKIYKFEDeEOnacFkJ8JNLwcc8i1w3+L2hGj4JCiOEG1+NjofqVzhJCuJiJ76U7L+06eF8IcULrZ/8KITaZ9jPdtQS01oJ+M6ifqSZxqwohlgshLmjX1b9C9WvtZO7kRMp3HF5CHR9/09IoQhunhPF41kIIsU0IcVsIkSyEeCMv68IgbaZ8gIUQRYGR2uGH6b39UhTlhqIoZ0zyyG6/TTVuaeHptoFB+b2EEL9o9f9UqN/bfCNUFyJL5+4hhJgnhDgl1HvfA61dvtbVr9B84U31Nq1nkY4PsMjevbKeUO8Bd7Rr56QQ4h1L5yVJA0VR5C+Xf6hP0grQyoL8rCZ/24zsP5rsH1R/trWoPnjPtPCvzKRppckOAHtQZ5i3or6Wua7JLgHFzaQdr8l16YNRDcYk4EstPNRMuqFaHAU4oqU7oh0nAe+aSROqyT/Xzuewdn7RWvh9oBIwF9V/bzewzuAc/gTsMtEOv2vpxmYw/p9a/JEm4br6+UI7t1BgDXBGC48DmpukqYHqZ6xo57cO+BmIAOKBnWbKHwMka79jQIgWPymNOtXptkCLE6a1xTGgmFa/CjDKwjn7avJrgK1B+Lda+ClgB/AjcNSgvA9N8vkY2KnJrgMrDH4fG8RbocUJNEkvgCBNlgD8oul+Xgu7AzRI4/ynaee/T9NVd009BZpkt20s1F0tLY/jJuElSOkbzzDpd8BxTVbTICxQC1thEDZI00sBHpnU6Vwzdfolah86p9Xdb0Ci7vrI5Bg21VQfE3kbTX7TjKywVraub2wFftKuMQV1fKxgksbLQLYS1Z/3F9Q+EEPK+FDBpH5WkDI+7DSonzcM4vXSrgPdGPKzdp3owqaZOQddmQu0v4dR+9V+XbuRMp4t0to7CnVc2A28npd1YdofMtHOnbU0/2LQ/zORPrv91uy4lYk2KIjafxXUvhum1ZOu/98F6pspv51WhwpwBVivXReRmj5TtXhvkNK/FIz7oGFfDcRCfyF798pZWr3qrq0DBrqMyWx75fef1RXIDz/SMIBRDY5n2kXtbkbeAPA1E14J1YhVgEYmslYGneIo8IqBrBgpN9xPTNLV0XRJBAJMZGMN8gw1kdXS0iQB3U1kvbTwBKC6iUzXqZOBbgbhNtqgoGgd/TpQzUDuQoqx2S+DbVDQYNBpkcE032vxV5qE6+ohzjAv1MH/f6Q8YBQ2kH2nhU8wU44TqY2y9qQMxqbt2wx1sE4AKlvQ7R7Q0ExZunwjLZzzZ5p8jkl4S8DLTPxGBrq4m8h012GoubK0OCswbwC/T4rx7GsQbgvMJ8UgKGTh/O8A9UyuqWWabLdJmky1TRrnIoCb2nXmYhDeTcv/D+1vF5NrOQm4AQiD8EDM3EAxMIYyUKcKquFqYyBrQcoDlEdGzktLN9WCPiWATqgz0wowOo1rKhooaxBuT4qxEmHhPBVU48nbQFYI2KbJvjFTXiiWx9uaqIbuQ6CDicyXlDG1tYksRgtPRDNm0yhXQfWbTat/5VVdKGTOAJ6mpdmb0TQ53G/NjluZaINZWpx9pB6Phmuyf4ACBuGepDwA/8dQZiCvZ07fNOohEPP9Jbv3SgWTiTKgLykPQQ5Zabf8+rO6AvnhhxkDGHUViA6ohlwyMDQL+b6LeWOlFSmGZQ0z6Xpo8l9Nwpeb67QG8mOYN4B16YIspPvR3ABt0KmDzaSpbdDhU91MgNGa7LsM1pWrQX4+GUyjG0y3m4Tr8plrJo0t6oybAvQxCNfdpGpnsOzDWvwOFuS6B5LPLeg20UI6G+CyFqemGd2varJUD11p6DpTSzPMwnUYmkbaFZg3gHV1aG42xI4UQ6WPiUx3/sPNpHtFkz0BCma1bdKpC921bvhAt1gLC9D+LjKQddXC1pjkE2iuL5I5A/gIBka1gXy7Ju+fifOaalC35n5XgZ5m0tmjGpsK0NaMvCTqbLaC6m9qep4K0NFMuoaa7LwZWSiWDWBd+5gdb0l5WFlnEh6jhS9Lo4505f5iQW6NukjTUDMTX3etrsloGpP02e23ZsetjLQB6sNkvFbHr1iIs1XLI8AgTPdmc20mzjPNesVy/83uvfJnC+lOafIMTe7In/qTPsB5y28GvkJ3UW9E5VCNnMWWEgkhCgshOgshZgghlmq+bitIWaKmsoWklxRFiTITrvsoooxJeEvt72oL+VkK16VbYUH+nfa3lQX5TjNh/6Qj133gYnoOOUl6S9ykqg9FUZJQX02B8fke0f4uEUK0FUIUslio6p/bAHVW4hcL0fZpf5tYkK83F6goSjLwg3YYaCJuh/pF9zFFUf4yo1cRzbdullDXBtVdh620KJauw0whhHAHKqA+wP1gKlfUDxeDtMNWpnKNrWbS3UR9tVsIdeZSR4bbJgPs0f76GYT5AbGKomwBYgF/ExnA3myUaYntinZ3NMFS/88I51Bfw+t+P6O+UXID5gkhXjOJXw91Fv2qoii7TTNTFOU2sEU7bGWmvGeY7/+ZPgeh+j3r3oD8bCFalvpVBuM8N3WRG+RQv81O/bZGfcjYp/V1c5hrX903DN9moOzskt17ZapxTeO5uAZeNOQyaHmLbh1ggbqWYgtUn7BVQohmiqL8Y5pACNEE1UfMPY18i1oIv2Qh/IH2t7BJuK6MCxbSxVgIL5tOunMm8UxJtfSboiiPhBAW5agzJZD6HCxxF3VgtkGdDT6TdnRAnTEEuGVBnl49GbbZHOBVVIPnF+CpEOIE6oC82uRBpTzqNVIUeGZQD+YoZSH8YhppVqD66PYRQoxXFOWZFj7AQG6EEKIz6uCc6iMdAyxdh5lFd51cU9QVBsyR3jWV1rXvjPF1k5m2SQ+dIesP6oc1qO5K3xvIA4UQHoqixJK7BnBm+39GSLUMGujHqd3AFiFEE0VRdA8V6Y0NkHZbXjO4PvUoivJA6xeZeVgpQco1ejMX+lV6cZ6nurCEbqzLygYXOdFvs1O/FbS/r2fgwz/D9i2n/c3sahlZIbv3ytzo0/kWaQDnLabrALuhGsU1gCAhRGPDGRshhAPqh2uuqK9OFqPOjD5UFCVZCNFOS29pJE/OlbOwjKVBJ72Z1DT11GYts4WiKIlCiChUH6xGqB9HpEdD7e/xrBZrUH484C+EaIQ649AMdRaiETBeCDFFUZRpWnRb7e99YGM6Zdw2W7C6uoh5pRTljBAiQiu/A6rRUgz1A5gEUmawAf3MzhrU2ZVZqP7ZMUCcdh0OBpaSfjtnFF0+ad3E0iwrM9dMJtsmvbzOCyEuABWFutJKG0201+BvoFbebtRZ8/OKoqRlFGWVPOv/iqJECCGWobomjSNlLd7stmVOnoOuXyVh+W1WeljsVxmI8zzVhSV0Y11dIYSt9kYro+REv81O/era9wzqGrxpYbgaS3rGcm6QK/dKSeaQBrAVURTlmlB3JfsD1djqg/HA3ALV+D2uKMogM1nk9NajV1Cfor1IeRI1xCuNdN5aWnPpyhvEsyabUQ3g/kKIuRZeDwMghKgPVEMdcLZbiOYFnLQQDqpfpBGKumPSYa0MO+At1F2rpgohflTUJYViteiJ5mbbcogVqEZeIOpr116oswc/K6mXPuqI9pGOoigTzOSV09ehbsa/jBCikKIuS2dKjl9TGWybjLAXdUUCf1KW5NIZwDoXCX9Sbma5MftrDXR933DzBV1blscyeTU+3EY1nuxRfcQfpRM/p3me6sISv6G+XSuO+nHjhkyktUq/NUA3bkZlcty8BPhov9zeiOpFuVfmC6QPsJVR1EXKv9YOpwohDB9KdK+bYzGP2XWAs4HOP6qPBbmlcF26/hbkA7W/oVnQKSdZiPqBRHXUmSqzaD6gC7XDEHOuKRqp6kMIYQv01A5D01JGUZQERVFWoM5WCNQv1FEU5Qrq6hclhbqdcG6wFtUY6CjUtUctuj+QxnWo1ZWldXJ1m4xk6kFbUXdDPI86PvU1U2ZBUq790MzknQkdzLZNBtEZtH6oM8B/KYpyXcv3OuoHK36k+AJnxgDOUp3mEbr1sg0Ny+PacVkhhJ9pAiFECdSPAyGXxwfNfUD3ANItrbi5xHNTF5ZQFOUBKZv/fC7MrEtsiBDiFSGEj5bW2v12D+oKC/5CWzs4g+zS/pqbZLJEIoDJ/TojvCj3ynyBNICfD2aiGmbeQD+DcJ1PUhshRBVdoBDCRggxGfVVbU6yCHVWqp/pxyxCXaS/voV081E/0OgthOhikq476qoTiVo8q6F9GKFbOWOuUDcIsTOMo7223or6+vsSMCyNLN8XBhsrCNUR71PUGdErqMsa6WTv624UJuVVQF1+CYx92yZpf1drri6m6eyEEJ0038tMo93oNqB+mT0FdTb4Oml/ZNNVCKHfrEWruwWk+N6ZopvFqJiFG8U87e90k2vfFnUpqXKo9WXpY6YMk4W2SY9fUa+xN1E/StljIt+D+mZHt+LAr5nI+xaqEewqhHDORLpcRbsOh2iHm3Xh2ivtJdrhV5rbly5NYVS3LifgkKIoB/JA1WmoY9FX2gedRq+ctbHVT5jZ3CW7PId1YYlPUd+ElAfChZnNY7Tx523UtdUNZ/zzrN+aoijKDdR7WHFgs2H5Bno4CyEGGY5jms6PUDc9mqDpapjGQwhRzyQr3diW2a2mX4h7Zb7B2stQ5Icf6WyEocWZrMU5h/EahZtJWbppB+rM3TnUTjIbM8tMkc7yU6SxlBIwQZMlo264EYTqopEEfGUpX9T1H5M1+SEtnW4pryTML2UWmla9kMZSM+mdYzrt0ZuU5Ybuohq8wRhvMPI74JmWXqRshPGblv60Fh4PtDRJc8KgfTdp9bOXlMX3Uy07BHxooM8Z7Vr4WavXe1r4exmtMzP5tzU4FwWT5fQM4hVAXRBeQf3YYjPqh5lXtHrUXRcrzKTVpTuF+mX4t8A4A/kKTR5okk6Qshb0U1TDfA0pyyzdJY0F9TPQF72y0zYZqNuTBvVquqZ2JwPZCQvpA9Oo0/WaLEbT9VvU7wvSrFMD+VRNPjUT56NL8w/Gi///hPGGKL9isP61ltZw84dH2vXzIylL7l0kjc0f0tDJbFuT/rjSC7WP6upwu6bPAdQHDMWwPi1dN5kt10p1keHxwCRdEYPrTNc3NqL2yb2kLOd2H4N1e8mlfpuJNihIylJiz0jZQGgLaj/X9ekqJuk6GJzTZVI2wzmOwUYYBvHnaXFvot6TvwW+zWD/zY175QrS6PPyZ+F6sbYC+eFHxgxgJ1J2MHrHINwO+Ah1x6LHqAP0ZtQZylbksAGsyd8EDqJu9HBfG/DaZCDfZqiD5g1UA/2GNpCY3UggA53a4qCYni4ZaBNXYLo2QP6LOqt2DdUY7k8auyCRclMQ2mB2EvWGehd1VtXc2ssdUWd/ftfa8CnqDPMvqB8M2Vgoq5Y2uP6jtf9D1K2zN6G+snMxp1sG68CGlHU5FdJY+xf1hjgb1RB/otXVGlS/uUAsD/ZeqDek66QY86EG8hVYGLi1+u2H+sHiPa3OLqC6DJndxCG988e8AZyltkmnbnU3yESgiImsKCk7sn1uIX1adVpCuyZiDfKJMZBbrFNNPpWsG8Cmv2eoRsAe4G0s9BtUw2QE6pJzD7Vr6LR2TZWwcN3kigGsxfFGfXA7hTrOxaEaabtQtwIuk951k5VyrVAXWTKADdK/ql1rZzRdE1D78m7gAwv65ni/zWgbGMTthGqwX9V0voN6D/0edTe3gmbSlEd1fTurtcl97fpYiMFGTFpce9QdTM9r+RvpTxr9V5Pn9L1yBdIAzvRPaJUnkUgkEolEIpHkC6QPsEQikUgkEokkXyENYIlEIpFIJBJJvkIawBKJRCKRSCSSfIU0gCUSiUQikUgk+QppAEskEolEIpFI8hW5vptQaGioUqhQodwuRqLx+PFj7O3tra2GJBeRbZw/kO2cP5Dt/PIj2zhviY+Pv+3n51cqvXi5bgAXKlSIKlVSbcgiySUiIyNlfb/kyDbOH8h2zh/Idn75kW2ct0RGRmZo507pAiGRSCQSiUQiyVdIA1gikUgkEolEkq+QBrBEIpFIJBKJJF8hDWCJRCKRSCQSSb5CGsASiUQikUgkknxFrq8CIZFIJBLJi8DTp0958uQJQog8K9PZ2ZkHDx7kWXmSvEe2cc6hKAqFCxcmJ5bXlQawRCKRSPI9cXFxABQtWjRPDWBbW1scHR3zrDxJ3iPbOOdQFIX4+HiePXuW7TqVLhASiUQiyfckJSXh6OiYp8avRCLJHEIIHB0dSUpKynZe0gCWSCQSiUQikeQrpAEskUgkEolEIslXSANYIpFIJBKJRJKvkAawRCKRPAfcu/eE3btjWLbsJDNnRrBy5UUWL/6dPXtiuHPnsbXVk0iyTHBwMB4eHjme77Bhw+jVq5f+ODk5mdGjR+Pt7Y2Liwvh4eGp4ryMhIaG4uLiwv37962tyguFXAXiBeXE5RNGx7Xda1tJE4lEklXi4xP5+eczrFkTzeHD18zEiAXAxkbQoEFpeveuSvfuVbC3l0O3RGXYsGGsWbOGiRMnMnbsWH14eHg4nTp14uzZs5QoUSJXdVi1ahXLly/n/Pnz2NjY4OHhQYcOHfjkk08A6NKlC23bts3xcv/3v/+hKIr+ePfu3QQHB7N582a8vLxwdnamZs2aRnHygmHDhnHnzh3Wrl1rFB6fEK///z8T/0NYaBhHjx5N9eHlvXv3qFq1KrNmzWLAgAF5onN+RI6iLyjS4JVIXlwSE5P4/vs/mTfvKDdvqjdFOzsb6tYtTbVqJShVyoGrV69SsGBxoqPvcOzYdQ4fvsbhw9eYPv0g48c3YuDAGhQoIF/iPW/8+++/nD9/nnv37lG8eHEqVKiAs7NzrpZZuHBh5s+fT2BgICVLlszVskxZvXo1EyZMYObMmbRo0YKEhAROnz7N0aNH9XHs7e2xt7fP8bKLFi1qdHz+/HlcXV1p1KiRPszOzi7Hy80qDnYO+v/fGfgOK75bwcGDB2nWrJlRvJCQEAoUKMCbb76Z1yrmK+ToKZFIJHlIVNQt/P1/5OOP93HzZjx16rzCokVt+eefwWzf3o25c1vz0UeNGDCgHHPntmbbtm78889gFi9uS506r3DnzhM++mgfLVuu4c8/b1n7dCQG/Pvvvxw7doynT5/i4uLC06dPOXbsGP/++2+ultu8eXM8PT2ZO3dumvEOHjyIv78/bm5u+Pj4MHHiRBISEvTygIAAxo4dy/Tp06lYsSKVK1dm0qRJJCcnW8xz586ddOzYkcDAQCpUqECVKlV44403mDlzpj6OOReIL774Ah8fHzw8PBg6dCizZ8+mVq1aernOdWHJkiX4+vpSvnx5hg0bRnx8fKo4uv8/+eQTLl++jIuLiz4vUxcIRVFYuHAh9evXp3Tp0vj6+jJt2jS9/NNPP6Vhw4aUKVOGWrVqMWXKFJ48eaKXz5o1i6ZNm7Ju3Trq1q2Lp6cnffv25c6dO3r5mjVr+OWXX3BxccHFxYWIiIhU9ebr60udOnUICgpKJVu9ejWdO3emSJEiAMyfP59mzZrh7u6Or68vo0ePTnNjjVWrVuHl5WUUZs5N4tChQ7z22muUKVMGX19fxo0bx8OHD/Xy8PBw/P39cXd3x8vLC39/f86cOWOx3BcNaQBLJBJJHqAoCsuX/4Gf349ERd3G07Moq1e/zp49PenduypOTpZnqooUsaNnz6rs2dOTVatex8urKNHRd/D3/5Fly07m+SteiXnOnz+Po6MjDg4OCCFwcHDA0dGR8+fP52q5NjY2TJkyhe+//54LFy6YjXP16lV69OhBzZo1CQ0NZf78+axfv97I+AN19tHW1padO3cye/ZslixZwoYNGyyW/corrxAZGUlMTEyG9V23bh2fffYZn3zyCb/99huVK1fm66+/ThUvIiKC6Oho1q9fz/Lly9m2bRtLly41m+f//vc/xo0bR5kyZYiOjmbv3r1m402fPp25c+fywQcfcPDgQb7//nvKli2rlzs4ODB//nwOHTrE3LlzWb9+PZ9//rlRHrGxsWzcuJEffviBn3/+mT/++IMZM2YAMHz4cN544w1atmxJdHQ00dHR1KtXz6wuffv2ZfPmzUbG7MmTJ4mKiqJv3776MFtbW2bNmsWBAwdYunQphw8fZsKECRZqN2NERUXRrVs3AgIC2L9/PytWrOD333/ngw8+ACAxMZE+ffrQvHlz9u/fzy+//MLgwYOxsXl5zMaX50wkEonkOSUxMYnRo39l3LhQnj1L5p13ahAe/havveadqY0XhBB07OjN/v19CAysTkJCMh9/vI+xY9V8Jdbl3r17qV7129vbc+/evVwvu23btjRq1EhviJny3Xff4erqyty5c/Hx8eH//u//mDx5Mt9++63RrKpuZrhixYp06dKF5s2bs2/fPovljh8/HmdnZ+rWrUv9+vUZMmQIa9euJTEx0WKapUuX0rt3b/r370/FihUZPXq0WSOxSJEifP755/j4+NCmTRs6d+5MWFiY2TyLFi2Kk5MTtra2uLq6mnUFefToEYsXL2bKlCn07duXChUq0LBhQ9555x19nHHjxtG4cWM8PT1p27YtH374IevXrzfK59mzZyxatAhfX18aNmzIgAED9Ho5OTlhb29PoUKFcHV1xdXV1aIbRteuXQGM8l+9ejWVK1emcePG+rBhw4bx6quvUq5cOZo3b86UKVPYuHFjth5858+fT/fu3Rk6dCje3t40aNCAOXPmsGHDBv7991/u3bvHw4cPad++PeXLl6dy5cr06NGDSpUqZbnM5w1pAEskEkku8vTpMwYO3MGqVX9RuLAtS5a0Y86c1mnO+KaHo2NB5s1rwzfftKdQIVu+/z6Kfv228vTpsxzUXJJZihcvzuPHxit2PH78mOLFi+dJ+VOnTmXTpk38/vvvqWR///03DRo0MJrBa9y4MQkJCUazxr6+vkbp3NzcuH37tsUyS5cuzS+//EJ4eDjvvfceiqLw4Ycf4u/vb2RYG3L27Fnq1q1rFGbOAPbx8aFAgZRPlUqXLs2tW1l3+zlz5gxPnz6lRYsWFuNs2rSJDh06UKVKFTw8PPRuFYZ4eHgY+R9nVa+iRYvSqVMnvRvEkydPWLdundHsL6juC2+88Qa+vr54enry9ttv8/jx4zTbJT1OnjzJmjVr8PDw0P8CAgIAuHDhAqVKlaJHjx68+eab9OrVi6+//porV65kubznEWkAv6CsOLTC6CeRSJ4/njx5Rt++29i+/TzFihVi8+au9OhRJcfy79q1Mhs3vomLS2F27YphwIDt0gi2IhUqVCAuLo74+HgURSE+Pp64uDgqVKiQJ+XXrVuXgIAAPv3001QyRVEsvm0wDC9YsGAqWVo+wDqqVavGoEGDWLZsGevXrycqKoqNGzdajJ+RNx9Z1cUS6c2YHj16lEHPIFS0AAAgAElEQVSDBtGmTRuCg4MJDQ1l4sSJqWazDY1ynV4ZmY29/ei20Q+gX79+HD9+nOjoaLZu3UpcXBw9e/bUp4mJiaF3795Uq1aNFStW8Ouvv/Lll18CGPlvG2JjY5NKH9NzSE5OJjAwkH379ul/YWFhHDt2jGrVqgGwZMkSdu3aRePGjdm2bRsNGzYkNDQ03fN8UZCrQLygfLj+Q6PjwMaB1lFEIpGYJSkpmcGDd7F370VKlCjMhg1dqF69VI6X06iRGxs3duGNNzbwyy8xvPPOTlaufA1bWzm/kdc4OztTv359zp8/z927dylevDhVq1bN9VUgDJk0aRJNmjRJ5QPr4+PDxo0bSU5O1s8CHzp0CDs7u1QfTGUXHx8fAOLi4szKK1WqRGRkJH369NGHRUZG5qgOlvQqVKgQYWFheHt7p5IfPnwYNzc3xo0bpw8znf3NCAULFiQpKSlVeOy9WKPjkk4ladKkCZUqVSIoKIg///yT9u3bU6pUyjgRGRlJUlISM2bM0Lfb9u3b0yy/ZMmSPHr0iLi4OBwdHQH4888/jeLUrFmT06dPp/twVqNGDWrUqMEHH3xAly5d+PHHH2nVqlWaaV4U5AgpkUgkOYyiKEyYEMbWrecoWtSOTZvezBXjV0f16qXYuLELxYsXYvv283zyyf5cK0uSNs7OztSrVw8/Pz/q1auXp8YvqLPQAwYMSPWx2Ntvv83169cZO3YsZ86c4ZdffmHatGkMGjQIBwcHC7mlz5gxY5gzZw6HDh0iNjaWo0ePMnToUBwcHGjdurXZNEOGDGHNmjWsXr2ac+fOMX/+fI4fP54pf/isUKRIEYYMGcK0adMICgriwoULHD9+nO+++w6AihUrcu3aNUJCQoiJieG7775j3bp1mS7H09OT6Ohozp49y507d9L0hwbo06cPQUFB7N+/P5X7g7e3N8+ePWPp0qVcvHiRkJAQli1blmZ+9evXx97enunTp3P+/Hk2btzIihUrjOKMHj2aI0eOMG7cOKKiojh//jw7d+7kww/VybXz588zbdo0jhw5wuXLl9m3bx+nT5/WP9y8DEgDWCKRSHKYJUtO8O23f2BnZ8Pq1R2pVi3312atXr0Uq1d3xM7OhmXLTrJ06Yn0E0leSsaNG5fqNX2ZMmX46aef+OOPP2jZsiUjRozgzTffZNKkSdkqq1WrVhw/fpx33nmHBg0a0L9/f0Bd6aFixYpm03Tt2pWxY8cybdo0WrVqRXR0NAMHDqRQoULZ0iUjTJ48mVGjRjF37lwaN27MgAEDuHr1KgDt27dnxIgRTJw4kVdffZXQ0NAsrbbQv39/KleujJ+fH5UqVeLYsWNpxu/Vqxfx8fG4ubnh5+dnJKtVqxYzZsxg4cKFNGnShODgYLMuLoaULFmSxYsXs2fPHpo3b05QUFCq86hRowZbtmzh3LlzvP7667Rs2ZIZM2bg6uoKqKthnD17lsDAQOrXr8+IESPo3bs3w4cPz3R9PK+I3F4+JyIiQqlSJed83iQqLuNdjI7vfnYXUF+XmH5cIHm5kG38fBMRcYVOndaTlKTw7bftefPNylnKJ6vtHBJymiFDfsHGRrBpUxeaNXPPUvn5jQcPHqTaWCEvMHxNnZ/p168fz549Y82aNdZWJcfRtfHvl40/TqzjXsdKGr34pNVfIyMjj/v5+dVPLw85AyyRSCQ5xPXrcbz99g6SkhRGjKibZeM3O3TvXoUPPqhHcrLCoEE79TvNSSTPC/Hx8SxcuFDvJjBv3jy2b99u5BMskeQ20gCWSCSSHCA5WWHw4J3cuBFP8+ZlmTSpqdV0mTixCc2aleXGjXgGD95JUpJcI1jy/CCEYO/evXTs2JFWrVqxYcMGlixZQseOHa2tmiQfIVeBkEgkkhzg669/Jzz8Cq+84sC333agQAHrzS8UKGDDN9+0p1WrNYSFXWbRot8ZOdL8blQSSV5jb2+f5u5yEkleIGeAJRKJJJucOnWHGTMOAjB/vh+vvJL1r+pzitKlHVm40B+A//43gujoO1bWSCKRSJ4fpAEskUgk2SAhIYmhQ3eRkJBM//6+tGtX3toq6fH396JfP18SEpIZNmw3iYmp1yaVSCSS/Ig0gCUSiSQbLFoUSVTUbby8ijJjxqvWVicV06c3x8OjCCdO3GT+/NzfbEAikUheBKQBLJFIJFkkJuY+c+YcAeCLL/xwcrKzskapKVq0EPPnq64Qn39+hJiY+1bWSCKRSKyPNIAlEokkCyiKwvjxoTx5kkS3bj60bOlhbZUs0rKlBz16+PDkSRIffbSP3F7/XSKRSJ53pAEskUgkWWDLlnPs2XORokXtmD69ubXVSZdPP21O0aJ27N4dw7Zt562tjkQikVgVaQBLJBJJJomPT2TixDAAJk9uiqvr87+Tl6urI//5TxMAPv54H3FxiVbWSJJXBAQEMH78+BzP18XFhU2bNumP//77b9q1a4ebmxu1atUyG+dlZMiQIfTt29faakgyiVwH+AVl3pvzrK2CRJJvWbz4BFevPqJmzVIMGFDd2upkmIEDa7B69Sn++OMWixf/ztixDa2tkiSb3L59m1mzZrF7925u3LhBsWLFqFq1KqNGjaJ169YArFq1igIFcv52Hx0dTfHixfXHM2fOxMHBgcOHD+Pg4GA2Tl7g4uLC999/T+fOnc3KExIS8PX1ZciQIYwdOzaVfPny5UyePJno6Ogc3R7bo/jz6yaVH5EG8AtKYONAa6sgkeRLbtyI46uvjgHqCgu2ti/OizRbWxumT29O584b+Oqr4/Tr5/tCzF6/SIScCWH6welceXiFskXKMqnpJLr7dM+18gYMGMDjx4+ZP38+FSpU4NatWxw4cIC7d+/q4zg7O+dK2a6urkbHFy5coEOHDnh6elqM8zxgZ2dH9+7dCQ4OZsyYMQghjORBQUF06tQpR41fgJJOJXM0P0n2eHFGbolEInkOmDXrMI8eJdK+fXleffXFm9F59VUP2rcvT1xcIrNnH7a2Oi8VIWdC+GDvB1x+eBkFhcsPL/PB3g8IOROSK+Xdv3+fiIgIJk+eTMuWLfHw8KBu3bqMGDGCrl276uOZukDcvHmTt956izJlylCzZk2CgoJo2rQps2bN0sdxcXFhxYoVBAYG4u7uTp06dfjpp5+Myjd0b3BxceHPP/9kzpw5uLi46PMydYG4du0agwcPxtvbm7Jly9KiRQv2798PqAZ0nz59qFKlCu7u7rRq1Ypdu3YZlVmrVi3mzp3L6NGj8fT0xNfXl/nz5xvJAQYOHIiLi4v+2JR+/foRExNDeHi4Ufiff/7JiRMn6NevH6DOsL/zzjv4+vpStmxZmjZtytq1a9NqFjp06MCECROMwkzdJJKTk/niiy+oU6cOZcqUoVmzZqxbt04vVxSFWbNmUbNmTUqXLk3VqlUZNmxYmuVKMoc0gCUSiSSDnDp1hx9++AtbW8HUqc2srU6WmTKlGba2gh9++IvTp+UOcTnF9IPTefzssVHY42ePmX5weq6U5+joiJOTEzt37uTJkycZTvf+++8TGxvLxo0bCQoKIiQkhMuXL6eKN2fOHF577TXCwsLo0qULI0aMIDY21mye0dHRVKpUiWHDhhEdHc3w4cNTxYmLiyMgIIBLly7xww8/EB4ezrhx44zk/v7+rF+/nrCwMAICAujfvz9///23UT6LFy+mWrVqhIaGMmrUKKZOncqRI+pyhHv37gXgyy+/JDo6Wn9sStWqValXrx5BQUFG4atXr8bb25umTZsC8OTJE2rXrs3atWs5cOAA7777LiNHjkxlOGeWadOmsXbtWubOnUtERASjRo1i5MiR7NmzB4ANGzawZMkSPv/8c44dO0ZQUBB169bNVpkSY6QBLJFIJBlk5syDJCcrDBxYg8qVXaytTpbx8XFhwIDqJCUpTJt20NrqvDRceXglU+HZpUCBAixcuJCffvqJ8uXL065dOyZNmsSxY8cspjl79iy//vor8+bNo2HDhtSoUYOFCxcSHx+fKm7Pnj3p0aMHFSpUYOLEiRQoUICIiAiz+bq6umJra4ujoyOurq44OTmlivPzzz9z8+ZN/Yxz+fLlCQgI4NVX1Q1kqlevzsCBA6lWrRoVKlRgzJgx1KxZk82bNxvl07p1a959910qVKjA4MGDqVChAmFh6kepJUuqbgbFihXD1dVVf2yOfv36sWXLFh48eADA06dPCQkJMZqpdXd3Z8SIEdSoUQMvLy8GDhxIhw4dWL9+vcV80+Phw4csXbqUBQsW4OfnR7ly5ejRowd9+vThu+++A+Dy5cuULl2a1q1b4+7uTt26dXnnnXeyXKYkNdIAlkgkkgwQGXmDHTsu4OBQgHHjXvyPx8aPb4SDQwF27rxAZOQNa6vzUlC2SNlMhecEnTp14tSpUwQHB+Pv78+RI0do164d8+aZ/1D67Nmz2NjYUKdOHX2Yu7s7pUuXThXX19dX/3+BAgUoUaIEt27dyrKuUVFRVKtWjRIlSpiVx8XFMWXKFBo3bkz58uXx8PDgxIkTqWanDfUCKF26NLdv3860Pm+++Sa2trZ614Nt27bx8OFDevXqpY/z7Nkz5syZQ/PmzfH29sbDw4MdO3aYnTHPKNHR0Tx9+pSuXbvi4eGh/61atYoLFy4A0KVLFx49ekTt2rUZOXIkmzZtIiEhIctlSlIjDeAXlBOXTxj9JBJJ7jJr1iEA3n23FqVKOVhZm+zzyisOvPuu6h+pOzdJ9pjUdBL2BeyNwuwL2DOp6aRcLbdw4cK0bt2a8ePHs2vXLvr27cvs2bPNGkyZ2QTFdOUIIQTJyclZ1jO9sidPnsymTZuYOHEiW7duZd++fdStW5fEROMl+3JKLycnJzp37qx3g1i9ejXt2rUz+nDvq6++YunSpYwcOZKNGzeyb98+2rdvn0onQ2xsbFKda2JiIknJScQnxPM4QXWTWbt2Lfv27dP/Dh48SEiI6i/u4eHB0aNHmTt3Lo6OjnzyySf4+fmZnamXZA25CsQLSpv5bYyO735210JMiUSSXY4cucaePRdxcirI8OEvjx/e8OF1Wb78D/bsuciRI9do2NDN2iq90OhWe8jLVSDM4ePjw7Nnz3jy5Al2dsbbc1euXJnk5GROnDhB/fr1Abhy5QrXr1/Pdb1q1qxJSEgId+7cMTsLfOjQIXr27EmnTp0A1f82JiaGihUrZqqcggULkpSUlKG4/fr1o3379uzatYuwsLBUPsGHDh2iQ4cO9OjRA1CN+H/++SfN1S1KlCjBjRspb1UUReGvv/6ilHspztw8g00JGwoWLMjly5dp1szytwT29va0b9+e9u3bM2LECKpXr86xY8do0aJFhs5NkjZyBlgikUjS4X//U2dIhwypTYkS9unEfnEoUcKeIUNqA3IWOKfo7tOdPwb+wZ2Rd/hj4B+5avzevXuXzp0789NPP/HXX39x8eJFNm7cyIIFC2jRooXZZbwqVapEmzZtGDNmDEePHiUqKorhw4fj4OCQajmwnKZbt26ULFmSvn37EhERwcWLF9mxY4d+FQhvb2+2bdvGyZMnOXXqFEOGDMnUx306PD09CQsL48aNG9y7dy/NuA0bNsTHx4ehQ4fi6uqKv7+/kdzb25vQ0FAOHz7MmTNnGDNmDFeupO3T3aJFC3bv3s2uXbs4e/YsU6dONXrAcCriRNd+Xfnkk08IDg7mwoUL/PHHHyxfvpxVq1YB6mz06tWrOXXqFBcvXmTNmjUULFiQChUqZLo+JOaRBrBEIpGkwcGDV9i3L5YiRewYNqxO+gleMN5/vw5FitgRGhpLRETufKwlyR0cHR2pX78+S5cupWPHjjRt2pQZM2bQtWtX/cdU5li0aBFlypShU6dOvPXWW3Tv3p2SJUtSuHDhXNd3y5YtuLm50bt3b/3SazrDe+bMmZQqVYrXX3+d7t27U79+fZo0aZLpcqZPn054eDg1atSgZcuW6cbv27cv9+7do3fv3tja2hrJPvroI2rXrk337t0JCAigWLFidOnSJc38+vfvT8+ePXn//fd57bXXKF68OO3btzeKM2jUIMaMGcNXX31F48aN6dq1K9u2baNcuXKA+hHfypUree2112jevDk7duxg9erVuLu7Z7I2JJYQmfEHygoRERFKlSpVcrWM/IjLeOMv0HUuEJGRkXKplJcc2cZ5y5tvbiA0NJbx4xvy8ceN86zcvGznWbMO8dlnR2jZ0oMNG9K+ub+sPHjwIMc3PsgIcXFxODpadzOSO3fuUK1aNb755hu9+4Ek59C18e+XfzcKr+P+8j1Q5xVp9dfIyMjjfn5+9dPLQ84ASyQSiQVOnLhJaGgsTk4Fee+92tZWJ9d4773aODkVZN++WE6cuGltdSS5TFhYGNu3bycmJoajR4/yzjvvUKJECfz8/KytmkSSZ0gDWCKRSCyg2/I4MLAGxYvn7utha1K8eGECA2sAKecseXlJTExk5syZNG/enLfeeovChQuzdetWq89ESyR5iVwFQiKRSMxw7tw9Nm/+Bzs7G4YOfXlnf3UMHVqbZctOsHnzP5w7dw9v7+LWVkmSS/j5+cnZXkm+R84ASyQSiRkWLDiOokCPHlVwc0u9q9XLhpubEz17VkVR1HOXSCSSlxlpAEskEokJ1649Yu3aaISAESPqWVudPGPEiLoIAWvXRnPt2iNrqyORSCS5hjSAJRKJxIQlS06QkJBMx47eVKrkbG118oyKFZ0JCKhIQkIyS5bIHSYlEsnLizSAJRKJxICHDxNYseJPAEaNSnclnZeOkSPVGe9Vq/7i0aPUW+lKJBLJy4A0gCUSicSAtWujefgwgcaNy1C3ruXtTl9W6tZ1pWFDN+7ff8qPP562tjoSiUSSK0gDWCKRSDSSkxWWLTsJwJAhtaysjfXQrXm8bNlJkpNzd7MkiUQisQbSAJZIJBKNvXsvcu7cPdzdi/D6697WVsdqdOzoTdmyTpw9+y979160tjqSbBIQEMD48eNzPF8XFxc2bdqkP/77779p164dbm5u1KpVy2ycl5Hg4GA8PDyynU9utZPEPHId4BeU/g37W1sFieSlQ/fh16BBNSlQIP/ODxQoYMO779Zi6tQDLF16grZtvaytksQCt2/fZtasWezevZsbN25QrFgxqlatyqhRo2jdujUAq1atokCBnL/dR0dHU7x4ynrRM2fOxMHBgcOHD+Pg4GA2zovIpUuXqF27Nnv37qVOndTbF3fp0oW2bdumm08JhxJpyk3bqVatWgwaNIgRI0ZkXmlJukgD+AXly25fWlsFieSl4syZu/z22yXs7QvQr5+vtdWxOv36+TJ79mF+/fUSp0/foUqVtG/eEuswYMAAHj9+zPz586lQoQK3bt3iwIED3L17Vx/H2Tl3VjJxdTX2kb9w4QIdOnTA09PTYpyXEXt7e+zt7dON5+nimaY8t9pJYp78O8UhkUgkBuh8f3v2rIKz88u77XFGcXYuTM+eVYCUupE8X9y/f5+IiAgmT55My5Yt8fDwoG7duowYMYKuXbvq45m+Wr958yZvvfUWZcqUoWbNmgQFBdG0aVNmzZqlj+Pi4sKKFSsIDAzE3d2dOnXq8NNPPxmVb+je4OLiwp9//smcOXNwcXHR52XqAnHt2jUGDx6Mt7c3ZcuWpUWLFuzfvx9QDeg+ffpQpUoV3N3dadWqFbt27TIqs1atWsydO5fRo0fj6emJr68v8+fPN4rz4MEDxowZQ9WqVXFzc6NRo0asX79eLz98+DAdO3akbNmy+Pr6MmbMGB48eJClNoDULhCzZs2iadOmrFmzhlq1auHj48OwYcNISEhg+fLlVK9eHW9vbz755BOSk5P16QzbKSAggNjYWKZMmYKLiwsuLi76c3vvvfeoXLkybm5u1KlTh8WLF2dZ9/yMNIAlEkm+5969J/z4YzQAgwfn34/fTNHVxY8/nubevSdW1kZiiqOjI05OTuzcuZMnTzLePu+//z6xsbFs3LiRoKAgQkJCuHz5cqp4c+bM4bXXXiMsLIwuXbowYsQIYmNjzeYZHR1NpUqVGDZsGNHR0QwfPjxVnLi4OAICArh06RI//PAD4eHhjBs3zkju7+/P+vXrCQsLIyAggP79+/P3338b5bN48WKqVatGaGgoo0aNYurUqRw5cgQARVHo0aMHBw4cYMGCBURERDBjxgzs7OwAOHXqFN26daN9+/aEhYWxcuVKoqKictzNIDY2lh07drBmzRqWLl3K5s2b6dOnD5GRkaxbt46vvvqKb775hq1bt5pNv2rVKsqUKcO4ceOIjo4mOlodn2bOnMmpU6dYs2YNhw8fZsGCBbi5ueWo7vkF6QIhkUjyPcHB0cTHP6NlSw/5qt+AKlVK0LKlB/v2xfLjj6cZMqS2tVWSGFCgQAEWLlzIBx98wMqVK6lRowaNGjWic+fO1K9vfg3rs2fP8uuvv7Jr1y4aNGgAwMKFC6ldO3Xb9uzZkx49egAwceJEli5dSkREhNkPvlxdXbG1tcXR0dGi28PPP//MzZs32bVrFyVKqP2sfPnyenn16tWpXr26/njMmDHs3LmTzZs3M3bsWH1469ateffddwEYPHgwy5YtIywsjIYNGxIaGsrRo0c5ePAgPj4+AHh5eenTLliwgDfeeMPIQP/8889p2bIlt27dolSpUmZ1zyxJSUksXLiQokWLUq5cOdq0acPBgwf566+/sLOzw8fHh0aNGhEeHk6nTp1SpXd2dsbW1hYnJyej+oyNjaVGjRrUq6eu123obiLJHHIGWCKR5GsURWHlSnXji0GDalpZm+ePgQNrAPD991EoilwS7XmjU6dOnDp1iuDgYPz9/Tly5Ajt2rVj3rx5ZuOfPXsWGxsbo4+53N3dKV26dKq4vr4pvvAFChSgRIkS3Lp1K8u6RkVFUa1aNb3xa0pcXBxTpkyhcePGlC9fHg8PD06cOJFqdtpQL4DSpUtz+/ZtfRmlS5fWG7+mnDhxgpCQEDw8PPS/Dh06AKoLRk7h7u5O0aJF9cevvPIK3t7e+plogFKlSmW6PgcOHMimTZt49dVXmTRpEgcOHMgxnfMbcgZYIpHkaw4cuMLZs//i5ubI//1f+fQT5DM6dChP6dKO/P33v0REXKVp07LWVkliQuHChWndujWtW7dm/PjxjBw5ktmzZzN8+HAjgwvI1EOM6coRQggjn9XMkl7ZkydPZu/evUybNg1vb2/s7e0ZOnQoiYmJGdYrvTKSk5Pp168fQ4cOTSXLSVcCczoWLFgwVVhSUlKm8m3bti0nT55kz549hIWF0atXLzp16sSiRYuyrXN+Q84Av6C4jHcx+kkkkqzx/fdRAPTt65uvlz6zRMGCtvTtq864ffddlJW1kWQEHx8fnj17ZtYvuHLlyiQnJ3PixAl92JUrV7h+/Xqu61WzZk1OnTrFnTt3zMoPHTpEz5496dSpE76+vpQpU4aYmJhMl3H9+nXOnDljVl6rVi1Onz5NhQoVUv0yspJDdnj09BG/X/5d/0sPOzs7sw8cJUqUoGfPnixatIivvvqKtWvX8vTp09xQ+aVGjvYSiSTfcutWPFu3nsPGRsilz9Kgf39fbGwEW7b8w61b8dZWR6Jx9+5dOnfuzE8//cRff/3FxYsX2bhxIwsWLKBFixZGr+B1VKpUiTZt2jBmzBiOHj1KVFQUw4cPx8HBASFErurbrVs3SpYsSd++fYmIiODixYvs2LFDvwqEt7c327Zt4+TJk5w6dYohQ4Zk6uM+gJYtW1KvXj0GDBjA3r17uXjxIr/99hvbtm0DYNSoUURGRvLhhx/yxx9/cP78eXbt2sXo0aPTzfvcuXNERUUZ/TKrX2bw8PAgIiKCq1ev6h8a/vvf/7Jt2zbOnTvHmTNn2Lp1K15eXhQqVCjX9HhZkS4QEokk3xIcfIrExGTaty+Pu3sRa6vz3OLuXoR27bzYufMCwcGnGDXK/AdWkrzF0dGR+vXrs3TpUs6fP09CQgJubm507drV6KMxUxYtWsQHH3xAp06dKFmyJBMmTODixYsULpy7y/85OjqyZcsWJk2aRO/evUlMTKRixYrMnDkTUFc4GDlyJK+//jrFihXjvffey/TMpo2NDSEhIUyePJn33nuPR48eUa5cOT766CNA9R/eunUrM2fOpGPHjiQnJ1OuXDlef/31dPMePHhwqrDw8PBM6ZcZJkyYwIcffki9evV4+vQpd+/epVChQsyYMYNLly5RqFAh6tevT3BwcK7p8DIjcvujhoiICKVKlSq5WkZ+xNTt4e5n6qLnkZGR1K1b1xoqSfII2cY5Q3KyQv36K4mJecDatQG0a/d8+f8+b+28e3cMPXtuply5ohw/PgAbm9ydLcxrHjx4YHbGNLeJi4vD0dExz8s15M6dO1SrVo1vvvnG7IoEkuyha2NTt4c67ql3lZNkjLT6a2Rk5HE/P790n9LlDLBEIsmXhIZeIibmAe7uRfDzK2dtdZ572rTxxNOzKBcvPuDXXy/i7+9lbZUkWSQsLIxHjx5RrVo1bt26xcyZMylRogR+fn7WVk0iyTOkASyRSPIlK1aoS58NGFAdW1v5OUR62NraMGCAL9OnR7By5Z/5xgDO6kfGtcrW4rdRv5mVtf6qNSevmN9dT/c2LzdJTExk5syZXLx4EXt7e+rVq8fWrVutPhMtkeQl0gCWSCT5jmvXHrFjx3kKFLChT59q1lbnheGtt6rx3/8eYteuGG7ejOeVVxysrZIkC/j5+cnZXkm+R057SCSSfMePP54mKUmhfXt1jVtJxnB1daRdOy+ePUvmxx9PW1sdiUQiyTLSAJZIJPkKRVEIDj4FIGd/s4CuzoKC/pI7w71gXLp0CRcXF37/Pe01aAMCAhg/fnweaZU1tm/fTv369SlVqhTDhg2ztjoWCQgI4D//+Y+11cgzXoRrR4d0gZBIJPmKw4ev8c8/93B1dZAfv2WBtm29eJz8eN0AACAASURBVOUVB/7++1+OHbtOgwY5t3vW80hu+OQa+gbn5SoQZcuWJTo6Wr8VcXh4OJ06deLs2bMWtyfOKLq8dNjb21O+fHmGDx9Or169spW3OUaNGkXfvn159913cXJyyvH8c4pVq1aRkJCQZhxFUfjhhx9YvXo1p0+rb1YqVKjAm2++SWBgoFVWJ0mP4OBgPvroI2JjY43CV61alWoXvOcVOQMskUjyFbrZ3549q8qd37JAwYK29OypLm25evUpK2sjyQy2tra4urrmqoFy8OBBoqOjOXDgAN27d2fYsGEcPXo0R8u4f/8+d+7coU2bNpQpUybLBqIlwzQ9gzUzODs7p2ugv/fee3z88ce0bduWDRs2EB4ezsSJE9m/fz9bt27NctmmW0iDuhV0ZrdfzgzOzs4UKfJirKkuR3+JRJJviItLZOPGswC89VZVK2vz4qJzg9iw4W/i4lLfZCV5w+7du/H09OTZs2eAulOZi4sLY8aM0ceZMWMGXbp0AYxdIC5duqSfsa1UqRIuLi5GrgTJyclMnz6dihUrUrlyZSZNmmR2W15TSpUqhaurK15eXowaNQpnZ2eiolK20H769CkTJkzAx8cHNzc32rZty6FDh/Ty8PBwXFxc2LdvH/7+/pQtW5Y2bdpw8uRJvbx8eXXN7s6dO+Pi4qLfjOLw4cN07NiRsmXL4uvry5gxY3jw4IE+74CAAMaMGcOkSZOoVKkSHTp0AMDFxYVvv/2W/v374+7uzvTp0wHVmPf398fNzQ0fHx8mTpxoZBwHBAQwduzYNOspPReIDRs2EBISwtKlSxk3bhz16tXD09OTdu3aERISot+gIzk5mTlz5lC9enVKly5Ns2bN2L59uz4fXduuW7eOzp07U6ZMGVasWEFwcDAeHh7s3r2bpk2b4urqqt8mOigoiMaNG+Pm5kaDBg34+uuvjXR/8OABY8aMoWrVqri5udGoUSPWr19PeHg4w4cPJy4uDhcXF1xcXJg1a5b+fA1dIO7du8f7779P+fLlKVOmDF26dCE6Olov1+m3b98+mjZtiru7O506deLixYsW6yynkAawRCLJN2zZ8g+PHiXSoEFpKlfO2vJWEqhc2YWGDd149CiRzZv/sbY6+ZYmTZrw5MkTvU9veHg4JUqUMNqdLDw8nObNm6dKW7ZsWVauXAmkzNr+73//08tDQkKwtbVl586dzJ49myVLlrBhw4YM65acnMy2bdu4f/8+deqkbPgwZcoU/XbNoaGhVKtWje7du3P9+nWj9NOnT2fKlCn89ttvODs7M3jwYBRFoWHDhhw8eBCAlStXEh0dTcOGDTl16hTdunWjffv2hIWFsXLlSqKiohgxYoRRviEhISiKwrZt2/j666/14Z999hn+/v6Eh4czaNAgrl69So8ePahZsyahoaHMnz+f9evXM23atFT5Zaeefv75ZypWrEhAQIBZebFixQBYsmQJCxYsYOrUqYSHh/P666/Tv39/o4cLgGnTpvH2228TERHBa6+9BsCTJ0+YO3cu8+bNIyIiAg8PD1auXMmMGTOYMGEChw4dYvr06cyfP5/ly5cDqltGjx49OHDgAAsWLCAiIoIZM2ZgZ2dHw4YN+e9//4uDgwPR0dFER0czfPhws/oPGzaM48ePExQUxJ49e7C3t6d79+48fvxYH+fp06d8+eWXLFiwgF27dnH//n0+/PDDDNdhVpEGsEQiyTcEBamv7N96S378ll10s8CrV/9lZU3yL05OTtSsWVNv8B44cIBBgwYRGxvL9evXiY+P5/fff6dZs2ap0tra2uLs7AykzNoauhLoZjwrVqxIly5daN68Ofv27UtXp9q1a+Ph4UHp0qUJDAxk6tSpegM4Li6O77//nilTptCuXTt8fHyYN28epUqV0hteOiZMmMCrr75K5cqVGTduHGfPnuXq1avY2dlRqlQpQH3d7urqip2dHQsWLOCNN95g+PDheHt7U79+fT7//HO2bNnCrVu39Pl6enoyY8YMKleujI+Pjz68S5cu9O/fHy8vL8qVK8d3332Hq6src+fOxcfHh//7v/9j8uTJfPvtt8THx2e7nnScO3eOSpUqpRtv0aJFDB8+nG7dulGxYkUmTpxIkyZNWLhwoVG8d999l86dO1OuXDnKli0LQFJSErNnz6Zx48ZUrFiRIkWKMHfuXKZOnaqP2759e0aNGqVvh9DQUI4ePcrKlSvx9/fHy8uLtm3b0rFjR+zs7ChatChCCFxdXXF1dTXr5nHu3Dl27NjBF198QdOmTalWrRpLlizh4cOHhISE6OM9e/aMzz77jHr16uHr68vw4cMJDw/P0BuH7PBieCpLUlGrbC1rqyCRvFDExNznwIEr2NsXoEuX9G84krR5441KTJwYRkTEVf75518qVnS2tkr5kubNmxMeHs7o0aM5cOAAQ4YMYf/+/Rw4cAAXFxcKFixIvXr1Mp2vr6+v0bGbmxu3b99ON92mTZsoXrw4T58+JTIyko8//hgnJycGDBhATEwMiYmJNGrUSB/f1taWBg0a6F/LmyvfzU390PL27dt6o+7/2bvv+Kqr+4/jr3Pvzd4JI4QZNlgFURBxgODAiYCo4B5V2mq11tHaOmqrtjhqHfXX2ronqFjEhQMVJwqoyB5hh5AA2evm3vP743KvRNYFk3zvvXk/H4/7uPf7zc297yQ393zu+Z7vOT/2zTffUFBQwGuvvRbaF5ylpKCgIFQ0Dxiw+7Zz4MCBjbaXL1/O4MGDcbl+6CccOnQo9fX1FBQUhPLtz+8pOW7XebPDmUmlvLycwsLCRr+3YJ5333230b6de9uDPB4PBx98cGi7pKSEjRs3ct1113H99deH9jc0NITyLFy4kNzc3EYfEvbX8uXLcblcDB48OLQvPT2d/v37N/p7JyQkNPoQkJubi9frpaysLPQhrTmoAI5Se1phSER27/nnA+POzjijJ+npCQ6niX5pafGMGdOL559fzPPPL+bWW3ftZZTmd9RRR/Hf//6XpUuXUllZycCBAznqqKOYM2cObdq0YfDgwcTFxe334/74e4wxYfXIde3aNTSjRL9+/Zg3bx733HMPF110Uai4Msbs1/MH77+35/f7/VxwwQX84he/2OVrwQIa2OOMG8nJjYtTa+0ec+68f39+T33a71pM9uzZk+XLl+/2/nt73j3t+/HPAYEC0+12h7aD+e677z6GDBmy2+dqiikO9/YYO+f+8UmZ4fy9m4KGQIhIzPP5/LzwQnD4g05+ayrnnx8YBhFYWKR5GyvZvSOPPJK6ujoeeughhg4ditvtDvUK72n8b1CweGvOWQHcbndovGd+fj7x8fGNTnrz+Xx89dVXP6mnEQI9u0uXLqV79+67XJKSkvb78fr06cNXX33VqAj74osviI+Pp1u3bj8p687Gjx/PqlWreP3113f79bKyMtLT0+nQoUOj31swz4H83tq1a0eHDh0oKCjY7e8L4JBDDmHz5s279MwHxcfH7/N106dPH/x+f6NZQMrLy1m8ePFP/ns3BRXAIhLzPv54Axs3VtK1azpHHdXJ6Tgx44gjOtCtWzqFhVXMmbPB6TitUnAc8NSpU0PF7uDBg9m4cSPz5s3b7fjfoM6dO2OMYdasWZSUlFBZWfmT8xQXF1NUVMT69et57bXXeOmll0KzLaSkpHDJJZdwxx138O6777Js2TJ++9vfUlxczGWXXfaTnveaa65h/vz5XHfddXz33XesXr2ad955h9/85jcH9HiXXnopmzdv5vrrr2fZsmXMmjWLO+64g8svv3y3vawHauzYsYwbN44rr7ySe+65h/nz57N+/Xree+89zjnnHN544w0ArrrqKh5++GFeeeUVVq5cyV133cXnn39+wIuA3HTTTTz00EP885//ZMWKFSxevJgXX3yRv//97wAMHz6cww47jIsuuoj333+ftWvXMnv27FCezp07U1tby+zZs9m6dWujcdFBPXr04JRTTuE3v/kNn3/+OYsXL2by5MmkpaVx1llnHeBvrOloCISIxLwXXggMf5g4sR8u174Pv0p4jDGcfXZfpkyZy9SpSxkxoovTkVqlY445hvnz54cK4MTERA4//HAWLFiw1/G/eXl5/O53v+POO+/kmmuu4dxzz+WRRx75SVmGDRsGBA5rd+zYkYsvvpibbrop9PXbb78dCBR0ZWVlHHzwwUybNo3c3Nyf9LwHHXQQM2fO5M477+S0007D7/fTtWvX0DRi+ysvL4+pU6dy6623Mnz4cDIyMhg/fjy33HLLT8r5Y8YYHnvsMZ566imeffZZHnzwQYwx5OfnM378+NDsEFdeeSWVlZXcdtttFBcX07NnT5566qlGY3v3x4UXXkhKSgoPPfQQf/7zn0lMTKRv375cfvnlALhcLqZNm8att97K5MmTqayspGvXrqG/5RFHHMEll1zCz3/+c7Zt28aNN97I7373u12e5+GHH+bmm29m0qRJ1NXVccQRRzBt2rQD6pVvaqa5l7L8/PPPbd++fZv1OeQH8+fPZ9CgQU7HkGakv/H+qaiop2/f/1BT08A331xMly6Rt6rS7kTL33n16lIOP/xpUlLiWLr0clJS9n+8aSQoLy93ZMWtllwJTpyhv3HT29v/6/z58+eNGjXq8H09hoZAiEhMmzlzFTU1DRx5ZF7UFL/RpHv3TAYPzqWqyssbb6xyOo6ISFg0BCJKHfeP4xpta1YIkd2bOnUpABMm6EhUczn33H589dVmXnppKWefrd+zyO4sK2p8QtnuZoWQlqMCOEp9u/FbpyOIRLzCwko+/ng98fEuzjyzp9NxYtaZZ/bi97//iI8+Wk9hYSUdOuw6Kb5Ia1ft3fVEMXGOhkCISMx6+eXlWAsnnphPZmai03FiVlZWIieemI/fb3n55fDmNBURcZIKYBGJWdOmBYY/nHOODss3t+DvODjkREQkkqkAFpGYtHhxCd9/X0JmZgLHH9/V6Tgx74QTupGVlciiRSV8/32x03EOSHPPiiQiP521tkn+V1UAi0hMmjo1cMLJmWf2IiFBpzs0t/h4N2PH9gICK8NFm+TkZMrKylQEi0Qway1lZWVNMq2cWgURiTl+vw0Nf9CsBC3nnHP68vjjC3n55WXcdttReDzR08fi8XhITU2lvLwcY1pusZStW7eSk5PTYs8nLS/4N563cl6j/T3SeziUKHpZa0lNTcXj+enlqwpgEYk5n3yygcLCKrp0SeeIIzo4HafVOPzwXHr0yGTVqlI++mg9o0ZF19ATj8dDRkZGiz7nypUryc/Pb9HnlJYV/Bvf+OaNjfZfPuJyhxIJaAiEiMSg4PCHCRP6tGhvXmtnjGHChMDcpi+/vGwf9xYRcY4KYBGJKdXVXmbMWAFo+IMTxo8PFMBvvBFYgU9EJBKpABaRmPL22wVUVnoZNKg9vXplOR2n1enRI5NDD21HZaWXd94pcDqOiMhuqQAWkZiik9+cF+wFfvVVLYohIpFJBbCIxIxt22p4//11uN0mNCWXtLyxY3thDMyaVUBZWZ3TcUREdqECWERixuuvr6Khwc/w4Z1p2zbZ6TitVocOqRx9dCfq6/3MnLnK6TgiIrtQASwiMWP69MAh97FjezucRMaPD/wNNBuEiEQizQMcpbZN2eZ0BJGIUlRUxSefbCQuzsVpp2mCeaedcUZPbrjhQ+bM2UBRURXt2//0lZtEopna7ciiHmARiQkzZqzE77ccf3xXMjISnI7T6mVmJnL88V3x+y2vvbbC6TgiIo2oABaRmBCccUDDHyLHuHGB2SBeeUWzQYhIZFEBLCJRb8OGCr78spCkJA+jR2tZ2UgxenQ+KSlxfP31ZtasKXM6johIiApgEYl6wUPsJ5zQjdTUeIfTSFBKShwnn9wd0JzAIhJZVACLSNQLFsDjxmn4Q6Q56yzNBiEikUezQESpa1++ttH2A2c94FASEWcVFJQyf34RqalxnHBCN6fjyI+MGNGFrKxEli7dxuLFJfTv38bpSCKOULsdWVQAR6mn5z7daFv/SNJaTZ8e6P09+eTuJCXpLS3SxMe7GTOmJ08++T0vv7yMW29VASytk9rtyBLWEAhjzBnGGLUsIhJxggWwhj9ErrPO+mE2CGutw2lERMIfA/xnoNAY87Ax5ojmDCQiEq5ly7axaFEJGRkJHHdcF6fjyB4MHZpHXl4q69dXMHfuZqfjiIiEVwBbawcAxwM1wCvGmGXGmD8aY7o1YzYRkb0KLn182mk9iI93O5xG9sTlMpx5Zi/gh7+ZiIiTwp4Fwlr7rbX2BqAz8CtgArDKGPOxMeY8Y4xmlBCRFmOtDQ1/GDu2l8NpZF+Cf6MZM1bi8/kdTiMird1+Fa3GmB7ArcCjQOKO248BVwEvN3k6EZE9WLSohBUrtpOTk8ixx3Z2Oo7sw6BB7enaNZ3Nm6v44otNTscRkVYu3JPgfmWM+QL4EmgPXGCt7WOtvdNa+wwwCjixGXOKiDQSXFjhjDN64fHoAFSkM2bnYRArHE4jIq1duK3GycB9QJ619pfW2i92/qK1thoY19ThRER2Z+fhD+PGafhDtNh5GERDg4ZBiIhzwi2AP7TWTrPW1u+80xhzXfC2tXZWkyYTEdmD+fOLWLu2nA4dUhg6NM/pOBKmgw9uS48emZSU1PDJJxucjiMirVi4BfCte9j/x6YKIiISrp2HP7jdGv4QLTQMQkQixV5bDmPMSGPMSMBjjDkuuL3jcjlQ0TIxRUQC/H7La69p+EO0Gjs2sGDJzJmr8Hp9DqcRkdZqX6u7/XfHdQLw+E77LbAZuLo5QomI7MmXX26isLCKzp3TOPzwXKfjyH7q3z+Hvn2zWbp0Gx9+uJ4TTujmdCQRaYX22gNsrc231uYDzwVv77h0t9YOs9bOaKGcIiIAO8392xtjjMNp5EAEe4GDPfkiIi0t3JXgLmzuICIi+9LQ4Od//9PiF9Eu+Ld7441V1NU1OJxGRFqjPQ6BMMYssdb223F7PYFhD7uw1nZppmyyFx/8+gOnI4i0uE8/3UhxcQ09emRyyCFtnY4jB6hnzywOPrgNCxeW8MEH6zj55O5ORxJpdmq3I8vexgD/fKfb5zd3ENk/AzsNdDqCSIsLzv6g4Q/Rb+zY3ixcWML06StUAEuroHY7suyxALbWfrLT7Y9aJo6IyO7V1/t4/fWVgIY/xIIzz+zFHXd8xttvr6ampoGkpH2dky0i0nTCXQr5OmPMwB23hxpj1hljVhtjjmzeeCIiAR99tJ7S0jr69cuhX78cp+PIT9StWwaDBrWnstLLu++ucTqOiLQy4c4g/xugYMftu4H7gTuBB5ojlIjIj/0w/EG9v7Hih0UxljucRERam3AL4AxrbZkxJg0YADxkrf0v0Kf5oomIBNTWNvDGG6uAH6bQkugXLIBnzVpDZWW9w2lEpDUJd9DVemPMMOAg4GNrrc8Ykw5oGR+HPPnFk422Lx56sSM5RFrCe++tpbLSy4ABbenRI9PpONJEOnVKY8iQDsydW8g776xh/Hh9uJHYpXY7soRbAN8AvAzUA+N37DsNmNscoWTfrnv1ukbb+keSWBY8RK7e39gzdmwv5s4t5LXXlqsAlpimdjuyhLsQxpvW2jxrbTdr7bwdu6cBZzRfNBERqKry8s47gVMQNP439pxxRk+MCfTyl5fXOR1HRFqJcMcAY4zJMMYMMcaMNMaMBI7ZcRERaTZvv11AdXUDgwfn0rlzutNxpIl16JDKsGEdqavz8dZbq52OIyKtRLjToF0MbAJeB/670+U/zZZMRAR47bXA8Idx43R4PFYFe/anT1/hcBIRaS3C7QG+EzjLWtveWpu/00XL94hIsykvr+Pdd9dgTOBQucSm00/victlmD17HaWltU7HEZFWINwC2APMas4gIiI/9uabq6mv93PUUR3p0CHV6TjSTNq2TeaYYzrh9fqZOXOV03FEpBUItwD+G/BHY0zYY4ZFRH6qHxa/0PCHWKdhECLSkvZnJbg/AhU7lkEOXZoxm4i0Ytu21fDhh+txuw2nn97D6TjSzE4/vScej4uPP15PSUm103FEJMaFOw/w+c2aQkTkR15/fRUNDX5GjuxCmzbJTseRZpaVlciIEZ157721zJy5iosvPtjpSCISw8KdB/ijPV2aO6CItE6vvRY4FK7hD61H8G8dHPoiItJcwp0GLcEYc6cxZrUxpmzHvhONMVc1bzwRaY2KiqqYM2cDcXEuTjtNwx9ai1NP7U58vItPP93I5s1VTscRkRgW7hjgvwM/A84D7I59i4BfNEcoEWndZsxYid9vGTWqKxkZCU7HkRaSnp7ACSd0w9rAa0BEpLmEWwCPBSZZaz8H/ADW2o1Ax+YKJiKtV3AmAC1+0fqceaaGQYhI8wu3AK7nRyfMGWPaAlubPJGItGobNlTwxRebSEryMHp0vtNxpIWNHp1PcrKHuXMLWb++3Ok4IhKjwp0FYhrwlDHmNwDGmA7AA8CLzRVM9u7+cfc7HUGkWfzvf4He3xNO6EZqarzDaaSlpaTEcdJJ+UyfvoLp01fw618f5nQkkSahdjuyhFsA3wxMARYCycAK4DHgT82US/bh4qEXOx1BpFlo+IOMG9d7RwG8XAWwxAy125El3GnQ6q2111prU4H2QJq19jfW2vrmjScirUlBQSnz5xeRmhrHCSd0czqOOGTUqK6kpcXz7bfFrFy53ek4IhKDwp0Grb8x5kpjzO+BcUC/5o0lIq1RsPf3lFO6k5QU7gEqiTWJiZ7Q9HdaGllEmsNeC2AT8DiBoQ83A2cAfwC+M8Y8YYwxLZBRRFqJ4Jn/Gv4gY8f2AmD6dM0GISJNb189wFcAI4Ch1tqu1tojrbVdgCOBY4ArmzmfiLQSS5ZsZfHirWRmJjBiRBen44jDhg/vTHZ2IkuXbmPxYk04JCJNa18F8AXAr621X+28c8f2tTu+Lg74ZsM3jS4i0S54qPv003sSH+92OI04LS7Ozemn9wTUCyyxQe12ZNnXILv+wEd7+NpHwDNNG0fCNfLBkY22t03Z5lASkZ/OWhsqcoKHvkXGjevFU099z6uvLufmm4eiUXcSzdRuR5Z99QC7rbUVu/vCjv3hLqQhIrJHCxcWs2pVKW3bJnH00Z2cjiMRYtiwjrRvn0xBQRnffLPF6TgiEkP21QMcZ4w5DtjTx26dpi0iP9mrrwaGP4wZ0wuPR5+rJcDtdjFmTC/+/e9vefXV5Rx6aHunI4lIjNhXAbsFeHwfXxcROWDWWs3+IHs0blxv/v3vb5k+fQV/+tPRuFwaBiEiP91eC2BrbbcWyiEirdRXX21mw4YK8vJSGTKkg9NxJMIMHpxLp05pbNhQwdy5hQwdmud0JBGJATrWKCKOCvb+jh3bS717sgtjjOYEFpEmpwJYRBzj8/n53/8C4381/EH2JPjaeO21FTQ0+B1OIyKxQAWwiDjms882UlRUTX5+BgMHtnM6jkSoQw5pS48emRQX1/DppxudjiMiMUAFsIg4Jjj7w7hxvTXHq+xRYBhEoBc4OGRGROSnCKsANsbcb4wZ2NxhRKT18Hp9zJixEtDiF7Jv48YFXiOvv76S+nqfw2lEJNqF2wMcB7xjjPneGHOTMUYz1YvIT/Lhh+vZvr2WPn2y6dcvx+k4EuH69s2hf/8cSkvr+PDDdU7HEZEoF1YBbK29GsgDfgcMBJYYY94zxlxojEltzoAiEptee03DH2T/BIdBTJ++wuEkIhLtwh4DbK31WWtnWmsnAkOBtsCTwGZjzH+MMR2bKaOIxJja2gZmzlwFaPiDhC/4WnnjjVXU1DQ4nEZEolnYBbAxJt0Yc5kxZjbwMfAlcAzQD6gE3mqeiCISa95/fy0VFfUMGNCWnj2znI4jUaJ790wOPbQdlZVe3ntvjdNxRCSK7WspZACMMS8DJxEofP8PeM1aW7fT168DypoloezWhUMudDqCyAF75ZXg4hea+1f2z9ixvVmwYAuvvrqc00/v6XQckbCp3Y4sYRXAwBfAVdbazbv7orXWb4xp33SxZF8eOOsBpyOIHJDy8jrefns1oMUvZP+deWYvbr31E2bNWkNFRT1pafFORxIJi9rtyBLuEIhjdlf8GmNeDd621lY3WSoRiVlvvrma2lofw4bl0alTmtNxJMp06pTG0KF51NQ08NZbq52OIyJRKtwC+Lg97B/RRDlEpJWYNm0ZAGed1cfhJBKtzjorcOQg+FoSEdlfex0CYYy5Y8fN+J1uB3UH1jZLKhGJSUVFVXz00Xri4lyMGaPZH+TAjBnTi9/97mM+/HAdxcXVtG2b7HQkEYky++oB7rzj4trpdmegE7AemNCs6UQkpkyfvgK/33L88V3Jykp0Oo5EqZycJEaN6oLPZzUnsIgckL32AFtrLwEwxnxmrX2sZSKJSKx6+WUNf5CmMWFCX955Zw1Tpy7liisGOB1HRKLMHgtgY0w3a+2aHZvvG2O67+5+1lqdheCA7BuzG21vm7LNoSQi4Vm1qpT584tITY1j9Ojdvp2IhG306HxSU+OYP7+IVatK6dEj0+lIInuldjuy7G0IxMKdbq8EVuy43vmiY08iEpZp05YCcPrpPUlKCncGRpHdS06O47TTegA/HFkQEQnXHgtga23aTrdd1lr3juudL+6WiSki0cxaq+EP0uQmTOgLBD5cWWsdTiMi0STspZBFRA7UggVbWL26jPbtkzn22E5Ox5EYceyxnWjfPpnVq8uYP7/I6TgiEkX2NgZ4DrDPj9TW2mObNJGIxJzgfK1jx/bG7dbnbmkabreLceN68+ij3zBt2jIOOyzX6UgiEiX2NhDvPy2WQkRiVkODn+nTlwMwYYKGP0jTmjChL48++g3Tpy/nL385Bo9HH7BEZN/2WABba59qySAiEps+/ng9W7ZU06NHJgMHtnM6jsSYAQPa0qtXFitWbOfDD9dx/PHdnI4kIlFgjx+V/hSdFgAAIABJREFUjTEX7HT70j1dWiamiESrnU9+M8Y4nEZijTEmdGRBSyOLSLj2NgRiIvDMjtsX7OE+Fni8SROJSMyorvYyc+YqQMMfpPmcdVYf7rrrC958czVVVV5SUuKcjiQiEW5vQyBO2en2cS0TR0RiydtvF1BZ6WXQoPZ0766FCqR5dOuWweDBuXz11Wbeemu1ptoTkX0K+2wBY0ymMeY8Y8wNO67VmonIXr30UmDxi7PP7utwEol1wdfY1KlLHU4iItEgrALYGDMSWAP8GhgMXA2sMcaMar5oIhLNioqq+OCDtcTFuRg/vrfTcSTGnXlmLzweF7Nnr6O4uNrpOCIS4cLtAX4YuMJae4S19mxr7VDg58AjzRdNRKLZtGnL8PksJ57YjZycJKfjSIzLyUli1Kgu+HyWV15Z7nQcEYlw4RbAecArP9o3HdCs4yKyC2stL764BIBzzunncBppLYKvteBrT0RkT/Y2C8TOngZ+BTy4075f7NgvDhjQcYDTEUT2aOHCYhYv3kp2diInntjN6TjSSowenU9mZgLffVfMokUlHHRQG6cjiYSo3Y4s4S6F7AJ+YYy5EdgIdATaA180e0LZrdnXzHY6gsgevfBC4ESk8eN7Ex/vdjiNtBaJiR7GjevN448v5IUXlvCXvxzjdCSRELXbkWV/lkJ+rDmDiEhs8Hp9vPJKYEGCc8/V8AdpWRMn9uPxxxcybdpSbrttGHFx+gAmIrvSUsgi0qTef38tJSU19OmTraWPpcUNGtSe3r2zWL58Ox98sI6TTsp3OpKIRKD9mQe4vTHmdGPMJVoKWUT25IUXAicgnXtuXy19LC3OGMPEiYEjD88/v9jhNCISqcKdB/hMYBVwB/AvAvMA/4s9L5EsIq3Q9u21vPNOAS6XYcIELX4hzpgwoS8ul+HttwvYtq3G6TgiEoHC7QH+C3CJtfZQoGrH9RXAvGZLJiJR59VXl1Nf72f48M7k5aU6HUdaqby8VEaM6IzX6+fVVzUnsIjsKtxp0LpYa6f9aN9TwGbg+qaNJOE47h/HNdrW2aUSCYLDH4KHoEWcMnFifz74YB0vvLCEyy/X9FPiPLXbkSXcAniLMaa9tbaIwBLIRwIlgE6vdci3G791OoJII8uXb2P+/CJSU+M45ZTuTseRVu6UU7qTnh7PggVbWLJkK/365TgdSVo5tduRJdwhEI8BR++4/XdgNvAt8M/mCCUi0SfY+ztmTC+Sk+McTiOtXVKSh7FjewNaGU5EdhVWAWyt/Zu19pUdt58GegOHWWtvac5wIhIdvF5fqAA+//yDHE4jEhAcijN16lIaGvwOpxGRSLI/06C5jTFHGWMmAJ0BnVkgIgC8++4atmypplevLIYMyXU6jggAgwfn0qNHJkVF1cyevc7pOCISQcKdBu0QYAUwDbhhx/UKY4zOLBARnnlmERDo/dXcvxIpjDFMmhToBQ6+RkVEIPwe4MeBR4CO1tohQEfg4R37RaQV27SpknffXYvH4+LcczX3r0SWiRP743YH5gTesqXa6TgiEiHCLYB7Aw9Yay3Ajut/AL2aK5iIRIcXX1yC328ZPTqftm2TnY4j0khubgonnZRPQ4NfJ8OJSEi4BfCbwBk/2nc68EbTxhGRaOL3W557LrDc7AUX6OQ3iUzB1+YzzyxiRz+OiLRye5wH2BjzDBB8p3ADLxpj5gHrCZwEdxjwv2ZPKCIR67PPNlJQUEZeXiojR3ZxOo7Ibo0a1ZUOHVJYtaqUTz/dyNFHd3I6kog4bG89wCuBVTsu3wN3Ae8Ai3dc3wXorAKRVix4YtGkSf1xu8OeVEakRXk8LiZN6g/oZDgRCdhjD7C19k8tGUREoktpaS2vv74SgPPO09LHEtkuuOAg7r//K2bMWMlf/zqcrKxEpyOJiIP2Zx7g44wxjxtj3tlxPbI5g4lIZHv55WXU1voYPrwzXbtmOB1HZK+6dElnxIgu1NX5mDp1qdNxRMRh4c4DfDnwErAZeBUoBJ43xvy8GbOJSISy1vLss4GT37Tym0SL4MlwTz+tk+FEWrs9DoH4kRuBE6y13wZ3GGNeAl4BHmuOYCISuebPL+K774rJykrk1FO7Ox1HJCynnNKdNm2SWLJkK/PmFXH44Vq1UKS1CrcAziFw8tvOlgHZTRtHwrVtyjanI0gr9vjjCwGYNKkfiYnhvo2IOCs+3s255/bj4Yfn8/TT36sAlhaldjuyhDsG+BPgfmNMMoAxJgW4B/isuYKJSGQqLa1l+vTlAFx88cEOpxHZP+efH5gNYvr0FZSX1zmcRkScEm4BPBk4GCgzxhQBpcAA4MrmCiYikemFF5ZQW+tjxIjO9OiR6XQckf3Su3c2Rx3Vkaoqr06GE2nF9lkAG2MMkAQcD+QTWAEu31o73Fq7qZnziUgEsdby5JPfA3DJJer9leh06aWB1+5//7tQJ8OJtFL7LIBt4N1hIeC31m6w1s611m5o/mgiEmnmzNnAihXb6dAhhZNP1slvEp1OPbUH7dsns2zZNj79dKPTcUTEAeEOgVgA9G7OICIS+YInv11wwUF4PFr5TaJTfLybCy/8GQD//e93DqcRESeEe/r2h8DbxpgngfVA6JiRtfbxpo8l+3Lty9c22n7grAccSiKtxebNVbz55mrcbhMqHkSi1UUX/Yz77/+KN95YTWFhJR06pDodSWKc2u3IEm4BfBRQAAz/0X4LqAB2wNNzn260rX8kaW7PPruIhgY/p53Wg7w8FQsS3fLyUjnllO68/voqnn56ETfddITTkSTGqd2OLGEdw7TWHreHi5ZDFmkFfD4/Tz2lk98ktlx22SEAPP3093i9PofTiEhL2msBbIxJNsbcZYyZYYy53RiT0FLBRCRyvP12ARs3VpKfn8Hw4Z2djiPSJI45phO9emVRWBgY3iMirce+eoAfJjDt2VLgLODeZk8kIhHnX//6BoDLLz8El8s4nEakaRhjQr3AwRM8RaR12FcBfDJworX2xh23T2v+SCISSRYtKuGTTzaSmhrHeef1dzqOSJM699y+pKTEMWfOBpYu3ep0HBFpIfsqgFOstYUA1tr1QEbzRxKRSPJ//xfo/Z00qT/p6RoFJbElPT2BCRP6APDYY5oSTaS12FcB7DHGHGeMGWmMGfnj7R37RCRGlZRU8/LLyzAGfv7zAU7HEWkWV1wReG2/+OIStm2rcTiNiLSEfU2DtoXG05xt/dG2BbQclEiMevrpRdTV+TjxxG706JHpdByRZtG3bw4jR3bhgw/W8dRT3/Ob3wx2OpKINLO99gBba7tZa/P3clHxKxKjvF5faJWsK68c6HAakeb1i18cCsB//vMd9fWaEk0k1mktUxHZrRkzVlFYWEWfPtmMGKGpzyS2jRzZhT59siksrGLGjJVOxxGRZqYCWER2Kzj12ZVXDsAYTX0msc0Yw+TJgSMdjz66AGutw4lEpDmpABaRXXz99Wa+/nozGRkJTJjQ1+k4Ii3i7LP7kp2dyIIFW/jyy0Kn44hIM1IBLCK7eOiheQBcfPHPSEmJcziNSMtISvKElvr+5z8XOJxGRJqTCmARaWTVqlJmzlxFfLxLJ79Jq3PZZYcQF+fizTdXs3ZtmdNxRKSZ7GsaNIlQH/z6A6cjSIx65JH5WBs4HJybm+J0HJEWlZubwrhxvXnppaU8+ug3/PWvw52OJDFC7XZkUQEcpQZ2Us+cNL0tW6p54YUlAFx11SCH04g446qrBvHSS0t55plF3HDDEHJykpyOJDFA7XZk0RAIEQl57LFvqavzcfLJ+fTune10HBFHHHRQG044oRs1NQ38+9/fOh1HRJqBCmARAaCysj608MXVVx/mcBoRZ/3mN4cDgQ+FFRX1DqcRkaamAlhEAHj22cWUltYxZEgHhg7NczqOiKOGDs3jiCM6UFpax9NPf+90HBFpYiqARQSv1xea9unqqzX2VwTg2msDvcD//OcC6uoaHE4jIk1JJ8FFqSe/eLLR9sVDL3Ykh8SGqVOXsWFDBb16ZXHyyd2djiMSEU44oRv9+uWwZMlWpk5dxgUXHOR0JIliarcjiwrgKHXdq9c12tY/khyohgY/f//7VwBcd93huFxa9lgEwOUyXHPNYUyePIuHHprHpEn9cLt14FQOjNrtyKL/ZJFWbvr05axeXUZ+fgbjx/dxOo5IRBk3rjdduqSzcmUpr7++yuk4ItJEVACLtGI+n5977w30/l577eF4PHpLENmZx+MKjYu/5565+P3W4UQi0hTU2om0YjNmrGTFiu107pzGOef0dTqOSEQ6//z+5OWlsmTJVl5/faXTcUSkCagAFmml/H7Lfff90PsbH+92OJFIZEpI8HDddYEZIaZMUS+wSCxQASzSSr355moWL95Khw4pTJrUz+k4IhHtvPN+6AWeMUO9wCLRTgWwSCvk91vuuedLAK655nASEjQhjMjeNO4F/lK9wCJRTgWwSCv0v/+tYOHCEjp0SNHcpiJhOu+8/nTsmMrSpdvUCywS5VQAi7QyDQ1+7r77CwBuuGEISUnq/RUJR6AXeDCgXmCRaKcCWKSVeemlpaxcWUq3bumcd15/p+OIRJVJk/qFeoFfeWWZ03FE5ACpABZpRerqGvjb3wJjf3//+6HExWnmB5H9kZDg4aabjgDgrru+oL7e53AiETkQKoBFWpGnnlrEhg0V9OuXw7hxvZ2OIxKVzj23H336ZLN2bTlPPvm903FE5ACoABZpJaqqvNx331wA/vCHobjd+vcXORAej4tbbjkSgHvvnUtFRb3DiURkf6kFFGklHn10AcXFNQwa1J6TT+7udByRqHbyyd0ZMqQDJSU1PPLIfKfjiMh+0unfUer+cfc7HUGiyObNVfzjH/MAuP32ozDGOJxIJLoZY7j99qM45ZSXeeSRBVx66SG0a5fsdCyJYGq3I4sK4Ch18dCLnY4gUeTuu7+gqsrLKad05+ijOzkdRyQmDB2ax+jR+bz9dgH33juXKVNGOB1JIpja7ciiIRAiMW7RohKee24xHo+L228/yuk4IjHllluG4XIZnnhiIUuXbnU6joiESQWwSAyz1nLLLXPw+y2XXnowPXtmOR1JJKb065fDRRf9DJ/P8oc/zMFaLY4hEg1UAIvEsPfeW8uHH64nIyOBG28c4nQckZj0+98fQUZGArNnr2PWrDVOxxGRMKgAFolRXq+PW26ZA8BvfzuY7OwkhxOJxKY2bZJDi2P84Q8fa3EMkSigk+Ci1Dcbvmm0PbDTQIeSSKT617++Zfny7eTnZ/Dznx/idByRmHbZZQfzxBMLWbFiO//617dcffUgpyNJhFG7HVlUAEepkQ+ObLS9bco2h5JIJNq0qZIpUwJLHv/tb8NJSNC/ukhziotzc+edx3D22TO49965nHNOX02LJo2o3Y4sGgIhEoNuuWUOlZVeTjutB8cf383pOCKtwvHHd+PEE7tRUVHPbbd94nQcEdkLFcAiMeajj9YzffoKkpI83HnnMU7HEWlV7rrrWBIS3Lz00lI+/ni903FEZA9UAIvEkPp6Hzfe+CEQOPGtc+d0ZwOJtDLdu2dy/fWBGVeuv342tbUNDicSkd1RASwSQx5+eD4rVmynR49MfvWrQ52OI9IqXX31IHr3zmLlylIeeOBrp+OIyG6oABaJEcuXbwud+HbvvcfpxDcRh8THu7n//sAJTw888DUrVmx3OJGI/JgKYJEY4Pdbfv3r96mv93P++f0ZPryz05FEWrVhwzoyaVJ/6uv9XH/9bK0QJxJhVACLxID//vc75s4tJDc3hT//WSe+iUSCP/3pKHJyEpkzZwNPPLHQ6TgishMVwCJRbv36cu644zMA7rlnBBkZCQ4nEhGAnJwk7rnnOABuu+1T1qwpcziRiASpABaJYn6/5Zpr3qeqysuYMT059dQeTkcSkZ2ceWYvzjyzF1VVXq6++j38fg2FEIkEKoBFothjj33Lhx+uJzs7kb/9bYTTcURkN+65ZwRt2ybx6acb+c9/vnM6joigAlgkai1ZspU//elTAB54YJSWXRWJUDk5Sdx3X2BWiD/96VNWrSp1OJGIqAAWiUL19T4mT36H2lof553Xn9NO09AHkUh22mk9mDChDzU1DVxxxdvU1/ucjiTSqqkAFolCd9/9BQsXltCtWzp33XWs03FEJAxTpoygc+c0FizYwl/+8rnTcURaNc2UH6UuHHKh0xHEIfPmlfLgg9/jchkeffQk0tLinY4kImHIyEjgscdGc+qpL/Pww/M59thOZGc7nUpaitrtyKICOEo9cNYDTkcQB2zaVMnddy/DWrjhhsEccUQHpyOJyH4YMqQDN988lD//+XN+9at3eeSRg52OJC1E7XZk0RAIkSjh9fq4/PK3KS31Mnx4Z264YYjTkUTkAFxzzeEMH96Z4uIa7r57GQ0NfqcjibQ6KoBFosSdd37BF19sIicnnn//+yTcbv37ikSjwPClE2nbNokFC8o0HljEAWpBRaLAm2+u4sEH5+F2G/74xz60baspz0SiWW5uCo8/fjIuFzz44Dxee22F05FEWhUVwCIRbvHiEiZPngXALbcM4+CDMxxOJCJN4aijOjF5cj4AV1/9HosXb3U4kUjroQJYJIKVlFQzadJMKiu9jB/fm6uvHuR0JBFpQmPH5jFhQh+qqrxceOFMSktrnY4k0ipoFogolX1j47lztk3Z5lASaS719T4uuuhN1q0rZ9Cg9jz44PEYY5yOJSJNyBjD3/8+kqVLt7JwYQkXXvgG06aNISFBzXOsUbsdWdQDLBKBrLX89rez+fzzTXTokMIzz5xKUpIaRJFYlJwcx3PPnU5ubgqffLKRa6/9AGut07FEYpoKYJEIdPfdX/Lcc4tJTHTzzDOn0aFDqtORRKQZdeqUxgsvnE5KShwvvbSUv/3tS6cjicQ0FcAiEebxx7/j3nvn4nIZHn/8ZAYNau90JBFpAQMGtOM//xmNy2WYMmUuzz67yOlIIjFLBbBIBJkxYyU33PAhAPfffxyjR3d3NpCItKiTTsrnb38bDsC1136g6dFEmokKYJEI8d57a7jiirexFm6+eSgXXvgzpyOJiAMuu+wQbrrpCPx+yxVXvMO7765xOpJIzFEBLBIBZs9exwUXvEF9vZ8rrhjAb3872OlIIuKgG28cwq9+dSgNDX4uuugNPvlkg9ORRGKKCmARh3388XrOO+916up8XHrpwdx997Ga7kyklTPGcMcdR3PxxT+jttbHuefO4OOP1zsdSyRmqAAWcdCHH65j4sTXqa31ceGFBzFlyggVvyICBIrge+89jokT+1Fd3cC5587QcAiRJqICWMQh//vfCs45ZwY1NQ2cf35/7r9/JC6Xil8R+YHLZXjooeO55JKDqa31cf75M5k5c5XTsUSingpgEQc8+eRCLr30LbxeP5MnD+SBB0ap+BWR3XK5DPfeO4Jf/GIgXq+fSy55U1OkifxEKoBFWpC1lilTvuS662ZjLfzxj0dy553HqPgVkb0yxvCXvxzD9dcPweez/PrX73PnnZ9rxTiRA6QCWKSF1NQ0cMUV7/DXv36JMXDffcdx3XWDNeZXRMJijOHmm4dy//0jcbsN9933FZMnz6KursHpaCJRRwWwSAvYvLmK009/hVdeWU5qahzPPXcal1xysNOxRCQKXXzxz3j++cCyydOmLWPMmOkUFlY6HUskqnicDiAHZkDHAU5HkDB9+WUhl176JoWFVXTunMYLL5xO//5tnI4lIlHshBO68cYb45k48XXmzi1k5MgXeeKJUxg6NM/paLIHarcjiwrgKDX7mtlOR5B98PstDz88nz//+TN8PssRR3Tg6adPpW3bZKejiUgMOOSQdsyePZHLLnuLTz/dyBlnvModdxzNlVcO0NCqCKR2O7JoCIRIM9i2rYbzznud22//FJ/PctVVg5gxY5yKXxFpUu3aJTN9+lh++cvAqnE33/wxEye+zpYt1U5HE4loKoBFmtibb65i2LDneOedNWRmJvD886dxxx1HExfndjqaiMQgj8fFX/5yDE88cTIZGQnMmrWGo49+jnfeKXA6mkjEUgEs0kRKS2v5xS9mcf75b7BlSzXDhuXx0UcTGT26u9PRRKQVGDOmF598Moljj+1ESUkNEye+zlVXvce2bTVORxOJOCqARX4iay3Tpi3lyCOf5aWXlpKU5OGuu45lxozxdO6c7nQ8EWlFOnZM49VXx3LHHUcTH+/i+ecXM3Tos0ydulRzBovsRAWwyE+wePFWTj/9Fa68chZFRdUMGdKBjz6ayOTJA7W4hYg4wuUyXHXVIObMOY+jjupISUkNkyfPYty411iyZKvT8UQigmaBiFLH/eO4Rts6u7Rlbd5cxX33zeXJJ7/H57Pk5CRy++1HM3FiPxW+IhIRevXKYsaMcTz//BJuuWUOH320nmOOeZ4LLjiI3/3uCNq3T3E6YquidjuyqACOUt9u/NbpCK1SWVkdDz44j3/96xuqqxtwuQyXXXYwf/jDkWRmJjodT0SkEWMM553Xn5NO6saUKXN54omFPPXU97zyyjJ++ctDmTx5oN67Woja7ciiIRAiYSgpqeauu75g4MAn+fvfv6a6uoFTT+3OnDmTuOee49SAiEhEa9MmmSlTRvDpp+cxenQ+lZVepkyZy4ABT3LnnZ/rRDlpddQDLLIXa9aU8cgj83nuucXU1voAOProjtxyyzAGD+7gcDoRkf3Tu3c2zz9/Op9/vpEpU+by0Ufrue++r/i///uG887rz+WXH0LPnllOxxRpdiqARX6kocHPrFkFPPXUIt57bw3BE6dPOqkb11xzuJYaFZGod+SRHZk+fSxz5xZy771zee+9tfz739/y739/y6hRXfn5zw9h1KiuuN06UCyxSQWwyA7Ll2/j5ZeX89xziygsrAIgPt7FuHG9ueqqw+jfP8fhhCIiTWvIkA5MnTqG778v5rHHvmPatKW8//5a3n9/LXl5qZx9dh/OOacfffpkOx1VpEmpAJZWbe3aMqZPX8Grry7n++9LQvt79szkwgt/xsSJ/cjJSXIwoYhI8/vZz9ryj3+M4rbbhvHss4t54omFrF1bzgMPzOOBB+Zx6KHtOPPM3px6ane6d890Oq7IT6YCWFqVhgY/X31VyLvvrmHWrDUsXvzDnJjp6fGcdlpPzjmnD0cf3QljNJ2ZiLQu2dlJ/PrXh3H11YP48stCXnxxCdOnr2DBgi0sWLCF2277hL59szn11B6cdFI+Awe2w+PRMAmJPiqAJab5fH4WLSrh88838dlnG/n44w2UldWFvp6SEsfo0fmMG9ebkSO7kJCgfwkREWMMQ4fmMXRoHnffPZxZswp4443VzJpVwNKl21i6dBv33fcVaWnxHH10R449tjPHHNOZvn2zNRe6RAW19hIzrLVs3lzFd98V8+23W/j66818+WUhFRX1je7Xs2cmJ5zQjRNO6MaRR+ap6BUR2YukJA9jxvRizJhe1Nf7+Oyzjbz55mo++GAtq1eX8dZbBbz1VgEAGRkJDBrUnsMOy+Www9pz2GHtadMm2eGfQGRXavklKm3fXsuKFdtZuTJw+f77Er77rpgtW6p3uW/XrukceWSgJ+OYYzqRn6/xayIiByI+3s2IEV0YMaILABs2VPDxx+v5+OP1zJmzgcLCKmbPXsfs2etC39OhQwp9++bQr98Pl169skhLi3fqxxBRASyRqaKino0bK9iwoYKNGyt3XCooKChj5crtlJTsftL29PR4Bgxox8EHt+XQQ9sxdGgeHTumtXB6EZHWoVOnNCZN6s+kSf0B2Lixgnnzipg3bzPz5m3mm2+2UFhYFSqMd9amTRJdu6aTn59J167pdOuWQZcuabRvn0JubgppafE6F0OajQpg2S/bt29n9erVlJaWkpmZSffu3cnK2vOk6X6/parKS0VFfaPLtm01lJT8cNm69YfroqIqysvr9/iYEBi726NHJj17ZtGzZyb9+uUwcGA7unRJ3+Mb5vbt25kzZw5ffvkl1dXVdOnShfbt21NQUMDq1atJSEhgwIABnHjiieTn54e+55tvvmHp0qUYY8jLyyMtLY2ysjKKioooKSmhurqa7OxsBg0axMCBA8nKygr9njZs2EBlZSUpKSl07tyZ7Oxstm3bFtpvrcUYQ2pqKsYYCgsLKS4uxhhD165dSU9PZ9WqVaxbtw6v10tiYiI1NTUsWbKEYcOGhXICFBQU8Nlnn1FQUEBpaSk1NTUkJCTQuXNnevToQXl5OUuXLmXdunVs374day0JCQlYa4mLi6Njx460a9eOlStXUlBQQHl5OT6fD7fbTWZmJklJSRQXF1NcXIzX6yUpKYmcnBzi4uIoLi6mtLSU2tpaANxuNx6PB7fbjbUWn89HQ0MDPp8vlDchIQGPx4PP58Pv9+P3+/H5fNjgxMs/Eh8fT3JyMtZa6uvrqa2tDd03OTmZjIwMKisrqaio2Otrp6W43W6ARj/zT308ay1+v79JHi8hIQGAurq6fdwzPB6Ph4aGhv36nri4OFJTU2loaKC2thafz4fL5Qq9ZiDwd3e73TQ0NOD1ejHGkJiYGNqfkJCAz+ejpqYGay3JyckkJCTQ0NCAtTZ0X7/fH3rc6urqUNbg93q9XlwuFwkJCWRnZ5OVlUVubi6dO3fGWktZWRl1dXUkJiaGfnfx8fG4XC42btzIli1b8Hg8tG/fHq/Xy6ZNm/B6vSQnBw7/V1dXExcXR25uLr169SInJ4eqqiq+//574uPjiYuLIykpidraWmpra8nMzCQzM5Pq6mrKysqoqqqisrISgK5du3LUUUfRoUMH/H4/mZmZZGdns2bNmkbvVQCbNm3CWkvfvn1D70872/k93eVyUVFRwcqVK9mwYQPGGLKysmjXrh0pKSkApKam0qlTpz2+93fsmEbHjmmccUbPHb9fP+vWlbNkyVaWLNnGkiVbWbwBQNMEAAATOUlEQVS4hIKCstD7/7x5Rbt9fSQne8jNDRTD7dun0LZtMhkZCWRlJZKVFbjOzEzccZ1AWlo8CQnuAy6a97d9k+hm9tTYNJXPP//c9u3bt1mf48d8Pj8NDf7QAgY/XNtG27vb9+Pfx56+d9/buz5euN+7r2zWwqC/Dtrpjoavb/oaa2HRokUcdNBB+P2WhgZ/6Nrn23Xb5/vhuqHhh2u/f+dtP7W1PurqGigtrWTNmo2AG7/fRU2Nl+pqLykpGfj9LurqGqipaaCiop7KSu+O6/pdfq5wJCa66dgxjU6d0sjLS6VTpzQ6dkyla9cMevbMJC8vdb/e5LZv387MmTOZO3cuWVlZ+Hw+Fi9eTHl5OSkpKbRt2xZrLUlJSbRr144LL7yQzMxMPvzwQ9avX09GRgZVVVUsXryYNm3aEBcXR2FhIVu3bqVTp064XC6ysrLo06cPAwcOZOXKlfj9fgoKCnC5XPj9ftq2bcvKlSvp2bMnxcXF1NbWsnbtWrp27UpVVRVr1qwhISGBxMREPB4PxcXFVFdX4/P5Qt/r9/tp06YNRxxxBHV1dZx99tnk5+dTUFDA1KlTSUhIYNWqVSxZsgS3202vXr0oKioiLi4Ot9tNeXk5mzZtwhhDVVUVDQ0NeDwe8vLysNaydu1asrOzKS8vp6qqirq6ulBhuXMhF3wtBgtXESfsXHS7XC6MMaHXaVxcHC6Xi4aGBowxoSK4vr4eY8w+P0ykpaWRkZGBtZaUlBSSkpJwuVyUlZWRkJAQumzYsAGfz0dGRgYNDQ0UFRVhjKFt27ahwhUgIyMDn89HcnIySUlJpKen4/F42Lp1K+np6dTW1lJWVkZcXBzdunWjpKSE0tJSOnToQFVVFSUlJbhcLjp16oTX6yUuLo6hQ4dy3HHHUVtby+eff05SUhK5ublUVVWxYMECAA499FBSUlIoKyujc+fOjBgxIlTUbd++na+//pqUlBS8Xi9ffPEFK1euxO12U1VVhdfrxev1kpWVRVxcHH369CExMZH8/HxcLheHH374AReIfr+lsLCStWvLWbOmjDVrAtcbN1aweXMVRUVVVFfv3wcqALfbkJwcR2pqHCkpcSQnB66Dt5OTPcTHu0lIcBMXF7x24ffXU1S0ieTkRFJS4rG2AZ+vnj59epKenobbbXC7XbhcZsftH2+7QvtdLtdO9zEsWrSYgw8+iEP/eihgYUfT9d3N32EMobYseHvXfbu/T2C78e3d3Scc+9OehnvXuDiXIwupzJ8/f96oUaMO39f9YrIHeNq0Zfzyl+86HaOZXdNo6/BpT++0Na9lo1C+16+mpsaRlhZPamo8aWmBS1ZWIm3aJJGTkxS6bts2eJ1MdnZikx76Wr16NatWrSI7O5vU1FQ2btyIx+OhsrISt9tNRkYGXq8Xay1er5fPPvuMvn37hnoCEhMTKSwsJDs7m6KiQG+FMSbU6OXk5FBfX09paSmfffYZ+fn5rF27luTkZBITE6mtrWXZsmW0adOG5cuX07FjR8rKykINX7C3xVpLTU0N3bp1Y8WKFZSVldGrVy/WrFlDZmZg7PL27dupq6sjKysr9FyfffZZqOe5uLiY7OzApPUFBQV06dKFlStXYq0N9TR5vd5Qj2tKSgp1dXVs3bqVhISEUO+w2+0mPj6eqqoqjDGhHjSPx4PH46G+vl7Frzhq5x5nv98f6nWHwIc0ay0ulyt0lCH4Px7ctzf19fVUV1eHXuuBoibwP1FTU0NaWhpFRUWhoyFer5e6urrQh8OKiorQc/n9furq6sjMzAz1uAff34JHSyorK0OFenFxMQ0NDSQmJobeb4JHa2pra0lISMDr9bJlyxY2bdoU+l3U1taSlJREYWFh6PErKipo06YNxhhKS0tZvXo1hx12GBB4X0xJSSE5OZnFixdTV1eHMYbt27fTvn17iouLQ73CGRkZ1NXVhY5ide3atdFj7S+Xy4R6i4cN67jL1621VFTUs3lzVagg3rq1lu3baykrq2P79trQpbQ0sF1ZWU99vT90pLFpbGqix/kauLbRnkNefKKJHjsyPfroCZxzTj+nY+xRTBbALpchPj7wqWPnT0S7297TvnC+d9ftvd9vXzma4jnr6+tISEjAmMAnT4/HFbp2ucLdduHxBD7FejyGxEQPiYkeNm1aR0ZGCvHxrp0uhrq6agYPPpTERDcJCZ5QkZueHk9KSlxELKVZWlpKeXk57du3B3447NvQ0BBqRD0eDzU1NRhj2Lx5M7m5ubscwkxNTWXNmjUAJCUlkZKSQnV1NfHx8VRWVuL1etm2bRv9+/cPNRoQaLy2bt1K9+7dWb58Od27d6e6uprk5GSqqqqorq4mMTEx1FsFgQY4eMi1oqKCnJyc0OH/6upqunbtysaNGwHYvHkzHTt2ZOPGjdTU1ISed/PmzaECPNjj5Xa7qa4OnCwYPORcV1cX+r7y8nI8Hk+o2K2ursbtduPz+UKFcPBwvEgk2fk1ufOHs+BrP7gvnNdu4EhY4P8x2OMafO0HC+ngsI3gEI36+vrQ9waL5uBzeb3eRoVu8D0oOTk5dH+3243L5aK6unrH+6+H2tpa3G43cXFxoe9LTEzE7/dTVVUVGvIT7O2GwHtV8H0k+L+ekJBAWVkZpaWloZ+xtLQ09GE5+P5ljKG2tjbUux7czszMpKqqKvQ4SUlJbNu27QD+SuExxpCenkB6egK9e4e/Cp3X66O6uoGqKm/oUl0duK6srKempoH6ej/19b7Qpa7Ox4oVq3G7E2ho8OP1WrxeP16vn+rqOnJy2uw4Who4Mur320bbPp/dsc8f2r/zdl1dHXFx8Vjb+Ejuj48QBz607XqfnY8sN962e32ccO3Pe/n+PG6kT4cXkwXw2Wf35eyzW3bYRaSYP38+gwYN2vcdD8C8eYnU1dWFCkIIvLkmJLThsMO6NctzNpXMzEzS09NDRWxwDF+wNxMIDQew1pKbm0tmZiZxcXGhBic5OTk0nhcCb9DV1dUkJSVRX1+Px+MJjfEL9hAFv7euro6cnBzKy8vJyckJ/R6DRXBycjI+n4+4uLhQ5vj4eBISEqitrQ09VnB/cnIy5eXl5ObmApCbm0t5eXno8GpwLG5aWhq1tbUkJiaGeoBra2tDP5fb7cbv94e+LzjOONh7FvydwA9jUF2uxh8uRSLFzq9Jl8sVKkCNMaFtv98f1ms3WIAGi0+32x36IBgsRoOFqM/nw+PxEB8fHype4+PjQx8ag+PsGxoaiIuLw+PxhN6D/H5/6P7BIj05OTk0Zj4xMRH4oUc72BPscrlISUkhLS0t9DjB94/gWPngbWBHERYXOpIEgffFmpoakpP/v737j427ruM4/nzR3lZor+0qXTe2yVbiFgkQTXSBWyLoYKIMUBPjTwLBYEwwAcMPcSQa/AX+iBIC0RhIAFkkBlRwiZExYYkaCDIZBgHdFIXBZO21W9duLZO3f3y/d9xmy9rR63e97+uRXLb7fr933/f6Wb/3us99Pp87jra2NgqFQnXcdOV3v9ITXRkKMjo6SrFYZN++fQc919GiUGiio6OJjo65U3rck08WJnh9m3vEvdwV9XxdtiOXfdeczRq9vb3V3sqIYGRkhOHhYXp7e7Mu7bB6e3s56aSTKJfLDA0N0drayoEDB2hra6OlpaU6Vq/yQlUqlejt7aWzs7M6oay9vZ1yuUxPT0913O/u3btpbm5maGiIOXPm0NnZSalUYnh4mK6uLkZGRhgcHGRkZIQVK1YwMDDA8uXLqxfWPXv20NLSwuLFi6sv1Mceeyz9/f10dnbS3d1Nf38/S5cuZXBwkHK5zLx586pDFUqlEgClUomBgQHmzp1Ld3c35XKZ3bt3s2zZMvr6+uju7mbhwoXV4Q+VQNDU1MToaPKpwcKFC6tDKyqT08bGxmhtba2+4FZCcO1zmGWl8uYM+L//i8kYSFUD5jHHHFMNrrVv5CZSeaNZKBSqb0bnzJnD2NhY9U1vT09P9c1koVCgpaWlet5isVg9PyQ9sHv37q1OdisWi9X5CJVrUeUTnu7ubgqFAvv376enp4e2tjZGR0erYbQSpOfPn88JJ5xAV1cXzc3N1Umy7e3t1cBcCau1E7sqaq/pixYtqobrynCqSs95sVisTjisTPqdLdf+yZrNr292ZBpyElye1fud5myeJdtIq0CUSiWvAuFVIN4SrwLhVSBg+leBmM3q9frmHuCZNdlJcA7ADca/aI3PbZwPbud8cDs3PrfxzJpsAPbnl2ZmZmaWKw7AZmZmZpYrDsBmZmZmlisOwGZmZmaWKw7AZmZmZpYrDsBmZmZmlisOwGZmZmaWKw7AZmZmZpYrDsBmZmZmlisOwGZmZmaWK3X/KuRNmzbtAv5V15NYVblcPr6rq6sv6zqsftzG+eB2zge3c+NzG8+4E1evXt19uIPqHoBtZkn6U0Qc9juwbfZyG+eD2zkf3M6Nz218dPIQCDMzMzPLFQdgMzMzM8sVB+DG85OsC7C6cxvng9s5H9zOjc9tfBTyGGAzMzMzyxX3AJuZmZlZrjgANyhJV0sKScdnXYtNP0nfk/ScpKcl/VJSZ9Y12fSQdK6k5yVtk3Rd1vXY9JO0RNIjkp6V9IykK7KuyepHUpOkP0vakHUt9gYH4AYkaQlwDvDvrGuxutkInBIRpwF/A76ScT02DSQ1AbcBHwJOBj4l6eRsq7I6OABcFRHvBE4HLnc7N7QrgGezLsIO5gDcmH4IXAt4gHeDioiHIuJAevcxYHGW9di0WQlsi4h/RMQYcC9wYcY12TSLiFciYkv69yGScLQo26qsHiQtBs4Dbs+6FjuYA3CDkXQBsCMitmZdi82YS4HfZF2ETYtFwIs191/CwaihSVoKvBt4PNtKrE5uJumQej3rQuxgzVkXYFMn6WFgwTi7rgfWAWtmtiKrhzdr54h4ID3mepKPU9fPZG1WNxpnmz/JaVCS2oD7gSsjYk/W9dj0krQWeDUinpR0Vtb12MEcgGehiDh7vO2STgWWAVslQfKx+BZJKyNi5wyWaNNgonaukHQxsBZYHV7PsFG8BCypub8YeDmjWqyOJBVIwu/6iPhF1vVYXawCLpD0YaAFaJd0T0R8NuO6DK8D3NAkvQC8JyL6sq7Fppekc4EfAGdGxK6s67HpIamZZFLjamAH8ATw6Yh4JtPCbFop6aG4CyhHxJVZ12P1l/YAXx0Ra7OuxRIeA2w2O90KFIGNkp6S9OOsC7K3Lp3Y+EXgtyQTo37u8NuQVgEXAR9If3+fSnsJzWyGuAfYzMzMzHLFPcBmZmZmlisOwGZmZmaWKw7AZmZmZpYrDsBmZmZmlisOwGZmZmaWKw7AZmZmZpYrDsBmZmZmlisOwGZmDUzS3KxrMDM72jgAm5lNgqSzJMUEtzvf5HELJN0l6WVJo5JekbRB0vyaY1ol3SRpe3rMTkn3S+qpOWalpIcl7ZU0LGmTpJWHnOtOSS9JOkPSHyXtA75bs/8ySVsl7ZfUJ+kOSV3T+5MyMzv6NWddgJnZLLEFOOOQbe8Hvk3ytcUT+SlwInAN8CLQA6wGjgOQNAfYCLwLuBF4DOgAPgjMA/4j6TRgM/BX4BIggOuAzZJOj4itNefrAO4Fvg+sA/al57kJuAq4Ja1lEfBN4BRJpYj475R+GmZms5i/CtnM7AhIegfwOPA74OMxwcVU0l5gXUTcMsH+S4E7gAsj4sEJjrkPOBtYGhGD6bZ24AXg0Yj4WLrtTuBi4CMR8UDN45cC24EbIuLrNdtXAb8HPhoRv5rsv93MbLbzEAgzsymSNA/YAGwDLoqIkNQkqbnmpvTwJ4BrJF0h6dSa7RVrgJ0Thd/U+4ANlfALEBF7gAeBMw859kBaW61zSK7362trJAnwe9LnNzPLDQdgM7MpSIPjfUALcEFE7Et3bQdeq7ldnG7/BElQvRZ4Gtgh6auSKtfftwE7DnPaLuCVcbbvJBkmUevVcYYzVMYbbzukxteA9rQGM7Pc8BhgM7OpuQ14L7AqInbWbD8fqF1x4Z8AEfEqcDlwuaQVJMH4BmAX8COgDzjlMOcsAwvG2b4g3VdrvKEY/emfa4CBN9lvZpYLDsBmZpMk6UvA50h6fv9Su+/Q++OJiOeBdZK+wBuh9yHgk5LOj4hfT/DQzcB5kooRMZTWUiQJ3Y9OovSNwOvA2yNi4ySONzNraA7AZmaTIKlEsrLC3UBZ0uk1u3dFxPZxHtMBPAysB54jGXJwIcmwhYfSw+4BLgN+JulGknG5RZJVIG6OiOeAbwBrgU2SvkPSy/tlkpUkqpPaJhIR29PH3Zr2Qm8G9gNLSMYH3x4Rj0zhx2FmNqs5AJuZTc5yknkTl6S3WneNsw2SkLmFJOCeSNIL+zzwmcoqDRHxmqQ1wNeAz6d/9gN/IB3eEBFPSzoL+FZ6LpEsl3bmIUugTSgi1kl6lnQ4BkmIfhHYBPx9Ms9hZtYovAyamZmZmeWKV4EwMzMzs1xxADYzMzOzXHEANjMzM7NccQA2MzMzs1xxADYzMzOzXHEANjMzM7NccQA2MzMzs1xxADYzMzOzXHEANjMzM7Nc+R+h5VtaxgO0zwAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.plot(x_random_nonsig, y_constant_nonsig, 'ko', label = 'Non Significant Values', alpha = 0.2)\n", "plt.plot(x_random_sig, y_constant_sig, 'go', label = 'Significant Values')\n", "plt.plot(x, y, '-', color = 'navy'); plt.xlabel('z-score', size = 16); \n", "plt.ylabel('Probability Density'); plt.yticks([])\n", "plt.vlines([stats.norm.ppf(1 - (0.025/n_tests)), stats.norm.ppf(0.025/n_tests)],linestyles='--',\n", " color = 'darkgreen', linewidth=4,\n", " ymin = 0, ymax = y.max(), label='Significance Limits \\nwith Bonferroni Correction')\n", "plt.legend(prop={'size': 14});\n", "plt.title('Random Observations with Bonferroni Correction', size = 22);\n", "\n" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.6.4" } }, "nbformat": 4, "nbformat_minor": 2 }