{ "cells": [ { "cell_type": "code", "execution_count": 65, "id": "01455f76-175f-40bd-b060-db11958d75f5", "metadata": {}, "outputs": [], "source": [ "import pandas as pd\n", "import numpy as np\n", "import matplotlib.pyplot as plt\n", "from scipy import stats" ] }, { "cell_type": "code", "execution_count": 179, "id": "c60e9af9-428a-46d0-96c1-b8cb28c5621e", "metadata": {}, "outputs": [], "source": [ "class HashRate:\n", " def __init__(self):\n", " self.hr = pd.read_csv(\"../data/hash-rate.csv\")\n", " self.di = pd.read_csv(\"../data/difficulty.csv\")\n", " self.data = self.hr.merge(self.di, left_on=\"Timestamp\", right_on=\"Timestamp\")\n", " self.data.columns = [\"date\", \"hr\", \"nd\"]\n", " self.regression()\n", " def regression(self):\n", " res = stats.linregress(self.data.hr, self.data.nd)\n", " self.slope = res.slope\n", " self.intercept = res.intercept\n", " def graph(self):\n", " plt.plot(self.data.hr, self.data.nd, 'o', color = \"teal\", label = \"original data\")\n", " plt.plot(self.data.hr, self.intercept + self.slope*self.data.hr, \"r\", label = \"fitted line\")\n", " plt.title(\"Hash Rate vs. Network Difficulty\", fontsize = \"x-large\")\n", " plt.xlabel(\"Hash Rate\")\n", " plt.ylabel(\"Network Difficulty\")\n", " plt.legend()\n", " plt.show()" ] }, { "cell_type": "code", "execution_count": 180, "id": "4d061903-ab21-46ed-98c2-32641047ea0f", "metadata": {}, "outputs": [], "source": [ "hr = HashRate()" ] }, { "cell_type": "code", "execution_count": 181, "id": "6a9d05bf-2cd5-422b-9fbc-beaaee058bd5", "metadata": { "tags": [] }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAERCAYAAAB2CKBkAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8/fFQqAAAACXBIWXMAAAsTAAALEwEAmpwYAABE9UlEQVR4nO2deXiU1dXAfychgbAYWawbAloRkR3ivqHYVvZKtcWCAi4srljrikJwqbgWrKCmlYpKXYoUQaD9KmpVtCqhQMSlUDYRNxYjS5CQnO+PeWeYmczyzmRmMpM5v+eZh7n3vfe+Z94M98y959xzRFUxDMMwspecuhbAMAzDqFtMERiGYWQ5pggMwzCyHFMEhmEYWY4pAsMwjCzHFIFhGEaWk5GKQERmisg3IvKRi7ZnichyEdkvIhf61bd16leIyGoRGZtcqQ3DMNKTjFQEwNPA+S7bbgJGAn8Jqv8SOFVVuwMnA7eKyBEJks+oBSIyUkQ6O+8fitK2WESahqifJyITROQWt33ikLNYRO4KJ2c02cOMGWqc3iJSIiI3iMjtTt1FIlIoIk+JSDe/f7u6GT8e2Yz6S4O6FiAeVPUtEWnnXyciPwamA4cAe4ArVfVTVd3gXK8OGmOfX7EhmasU6yvDReQroLWI5ADXAkfiUein4vnurnTaXi4iRcBYVd3t1K1V1XsdZRDc/1hgmIj8Hbge2A68oapLnbZTVPVmEXkQ+DdwDPCxqi4MIad6v4sicghwm3c8oJOIXAocDjwJvAb8DLgU2Aq0Ag4DbgWWANOccU7G8+PkD3rgxOd8VX1VRK4QkU7O9UogHzjB79+dInIc0A5YAzQHljn3auqMf6wj2+nAIFW9RUQeVNWb3P1pjPpGfZr8SoBrVbUX8FtgRrQOInKUiKwCPgfuV9UtSZbRcM9zqjoV2Aw0xjPxbwHOBFbgmfiaOW2fwTPxtvPrf4yIPOC8D+6/FpgNnARUAF/gmexR1WpgmzNJfgKsdsYoDCPnY8DVzvveQeOtVtVngH/hURAvAzfhUQinqOrvOaAQVqjqK87nGkOgEgjGW78CeB94z/n3S6f+ZFV9yBlPAXFeOJ9xrSPbUmCdiAzFo/CMLKVeKAJnmX8a8FcRWYHn19fh0fqp6ueq2hXPL8QRInJoUgU14uVHQEs8k1oucDCwF+jkXK8EqvGb7IB1qnqz07dXUP/PgHHAB3gm3gJglV/f2Xh+WLyIZ5Leg+cXdBcROSVItgqgDM8v+38FjfeFiIwDPgT6A1OBfqpaBvxbRG5wxt/qyA+wD5gCFIuI/4p9kNO+pap+HOV5vS8iN4rIQGA5MAz4aVCbg0TkbOA5YDwwL8qYRj1GMjXWkLMcf1VVO4vIQcBnqhp28heRp532c8JcnwksCnfdMOojzmqgoarOqmtZjLqjXqwIVPV7YL2IXAQgHrpF6iMirUWkwHnfHDgDzy9Fw8gaVPUFUwJGRioCEXkez75oBxHZLCKX41n+Xi4iK/Hs6w522p4oIpuBi4AnRcS759sRzxJ6JZ4l/UPOkt0wDCOryNitIcMwDCMxZOSKwDAMw0gcGXeOoFWrVtquXbu6FsMwDCOjKC0t3aqqh4S6lnGKoF27dixbtqyuxTAMw8goRGRjuGu2NWQYhpHlmCIwDMPIckwRGIZhZDkZZyMIRWVlJZs3b2bv3r11LUpW06hRI1q3bk1eXl5di2IYRgzUC0WwefNmmjVrRrt27RCR6B2MhKOqbNu2jc2bN3P00UfXtTiGYcRAvdga2rt3Ly1btjQlUIeICC1btrRVmWFkIPVCEQCmBNIA+xsYRmZSbxSBYRhGvaWiAiZNgs8/T8rwpghSTL9+/fjuu+8itpk4cSKvvfZaXOO/+eabDBgwIGq73r17Rz2YN3XqVPbs2ROXHIZhJIiZM6FxY7jrLvjHP5Jyi3phLI6V2WVlTFiyhE3l5bQpLOTePn0Y1qVLUu+pqqgqixYtitr2rrvuSqosbpk6dSrDhw+ncePGdS2KYWQf330HzZsfKA8fDldckZRbZd2KYHZZGaMXLGBjeTkKbCwvZ/SCBcwuq10E6kceeYTOnTvTuXNnpk6dCsCGDRvo0KEDl156KZ07d+bzzz+nXbt2bN26FYC7776bDh06cMYZZ3DxxRfz0EOefOIjR45kzhxPfpx27doxadIkevbsSZcuXfj0008B+OCDDzj11FPp0aMHp512Gp99FjmVQkVFBUOHDqVjx45ccMEFVFRU+K6NGzeOoqIiOnXqxKRJkwB49NFH2bJlC+eccw7nnHNO2HaGYSSBKVMClcD//gfPPpu022XdimDCkiXsqawMqNtTWcmEJUviXhWUlpby5z//mffffx9V5eSTT+bss8+mefPmrFmzhlmzZnHKKYEZDj/88ENefvllVq5cSWVlJT179qRXr14hx2/VqhXLly9nxowZPPTQQ/zpT3/i+OOP5+2336ZBgwa89tpr3H777bz88sthZXz88cdp3Lgxn3zyCatWraJnz56+a/feey8tWrSgqqqKPn36sGrVKq677joeeeQR3njjDVq1ahW2XdeuXeN6ZoZhhGDLFjjyyAPlm26CBx4I3z5BZN2KYFN5eUz1bnjnnXe44IILaNKkCU2bNmXIkCG8/fbbALRt27aGEgBYunQpgwcPplGjRjRr1oyBAweGHX/IkCEA9OrViw0bNgBQXl7ORRddROfOnbnhhhtYvXp12P4Ab731FsOHDwega9euARP4Sy+9RM+ePenRowerV6/m449Dp8R1284wjDi4/vpAJfDVVylRApCFiqBNYWFM9bWlSZMmtR6jYcOGAOTm5rJ//34A7rzzTs455xw++ugjFixYELf//vr163nooYdYsmQJq1aton///iHHctvOMIwYWbMGRODRRz3lRx4BVTj00JSJkHWK4N4+fWgcFAKhcV4e9/bpE/eYZ555JvPmzWPPnj3s3r2bv/3tb5x55pkR+5x++um+CXzXrl28+uqrMd2zvLycI51fD08//XTU9meddRZ/+ctfAPjoo49YtWoVAN9//z1NmjShsLCQr7/+msWLF/v6NGvWjJ07d0ZtZxhGHKjCL38Jxx13oK68HG64IeWiZJ2NwGsHSKTXUM+ePRk5ciQnnXQSAFdccQU9evTwbeOE4sQTT2TQoEF07dqVQw89lC5dulAYw6rk5ptvZsSIEdxzzz30798/avtx48YxatQoOnbsSMeOHX32iG7dutGjRw+OP/54jjrqKE4//XRfn9GjR3P++edzxBFH8MYbb4RtZxhGjCxbBieeeKD87LMer6A6IuNyFhcVFWmw//snn3xCx44d60ii+Nm1axdNmzZlz549nHXWWZSUlAQYcTORTP1bGEZKqK6GM86A997zlA89FDZuBGf7N5mISKmqFoW6lnUrgnRi9OjRfPzxx+zdu5cRI0ZkvBIwDCMCr70GP/nJgfLixXD++XUnjx+mCOoQ7569YRj1mH374NhjD4SH6NEDPvwQcnPrVi4/ss5YbBiGkTJefNGz7eNVAu+9B8uXp5USAFsRGIZhJJ5du6Cw0GMTABg4EF55xeMmmoaYIjCMekwi42rFO1Ys/cK1rYv4YHEzfTpcc82B8scfQ5o7UCRNEYjIUcAzwKGAAiWqOi2oTW/gFWC9UzVXVdMj4pphZDjeuFrekCreuFpAzJNovGPF0i9c26WbNjFr5cqEfI6ksnUrHHLIgfLYsfD443UnTwwk00awH7hRVU8ATgGuFpETQrR7W1W7O6+MVQKPPvooHTt2ZNiwYcyfP58pU6YAMG/evIBQDE8//TRbtmyJaewNGzbQuXPniPXLli3juuuuq8UnMOobkeJqpWqsWPqFa1tSWpqwz5E0Jk4MVAKbNmWMEoAkrghU9UvgS+f9ThH5BDgSqJcBambMmMFrr71G69atARg0aBDgUQQDBgzghBM8OvDpp5+mc+fOHHHEEQm9f1FREUVFIV2EjSwlkXG14h0rln7h2laFOetUm/hgCWPjRmjX7kB58mSPUsgwUuI1JCLtgB7A+yEunyoiK0VksYh0CtN/tIgsE5Fl3377bTJFjYuxY8eybt06+vbty+9//3uefvpprrnmGt59913mz5/PTTfdRPfu3bn//vtZtmwZw4YNo3v37lRUVFBaWsrZZ59Nr169+NnPfsaXX34JeCKaduvWjW7dujF9+vSoMvgnpCkuLuayyy6jd+/eHHPMMTzqjWECPPfcc5x00kl0796dMWPGUFVVlZyHYtQ5iYyrFe9YsfQL1zY3jIE1WfHBXHPFFYFKYOvWjFQCkAJjsYg0BV4Gxqvq90GXlwNtVXWXiPQD5gHtg8dQ1RKgBDwniyPecPx4WLGi1nIH0L07ODkGQvHEE0/w97//3Rey2Rv757TTTmPQoEEMGDCACy+8EIDFixfz0EMPUVRURGVlJddeey2vvPIKhxxyCC+++CITJkxg5syZjBo1iscee4yzzjqLm266KWaRP/30U9544w127txJhw4dGDduHGvXruXFF19k6dKl5OXlcdVVVzF79mwuvfTSOB6Kke7c26dPwJ47xB9XK96xYukXru2Ibt0CbAS1+RwJ4aOPwN828fjjHntABpNURSAieXiUwGxVnRt83V8xqOoiEZkhIq1UdWsy5UoXPvvsMz766CN+4pw2rKqq4vDDD+e7777ju+++46yzzgLgkksuiTnIW//+/WnYsCENGzbkRz/6EV9//TVLliyhtLSUE50YJxUVFfzoRz9K7Icy0oZExtWKd6xY+kVqe3qbNnXvNaQK/frB3//uKTdsCNu2QQIiDNc1yfQaEuAp4BNVfSRMm8OAr1VVReQkPFtV22p14wi/3NMNVaVTp06854074hAtp7EbGvrFLvGGr1ZVRowYwX333Vfr8Y3MYFiXLgmbMOMdK5Z+4dom8nPExbvvgn+gxTlz4Be/qDt5EkwybQSnA5cA54rICufVT0TGioh3HXUh8JGIrAQeBYZqpkXBi4J/KOfgcocOHfj22299iqCyspLVq1dz8MEHc/DBB/POO+8AMHv27ITI0qdPH+bMmcM333wDwPbt29m4cWNCxjaMeklVFXTtekAJ/PjHnpAR9UgJQHK9ht4BIh6jU9XHgMeSJUM6MHToUK688koeffRR5syZw8iRIxk7diwFBQW89957zJkzh+uuu47y8nL279/P+PHj6dSpE3/+85+57LLLEBF++tOfJkSWE044gXvuuYef/vSnVFdXk5eXx/Tp02nbtm1CxjeMesWrr3pOBHt5/XVw8nfXNywMtZFQ7G9hZDx798IRR8COHZ7y2Wd7lEBOZodmszDUhlEPyYSwC+kgYygZlm7aRElpKVWq5IowulcvZvTvD7NmwciRBzovX+6JFlrPMUVgGBlIIsNHJIt0kDGUDCP+9reAQ2pVqsx+5x1mOOdwALj4YsiiMPGZvdbxI9O2uOoj9jdIHYkMH5Es0kHGUDIEn1T+7dKllDshYQBPMvksUgJQT1YEjRo1Ytu2bbRs2RJJ0zCv9R1VZdu2bTRq1KiuRckKEhk+Ilm4kTGZW0dXLVzIxgjP49CdO/nq4Yd95YdOPZXfvvuuR6apU9Nqyy3ZW2z1QhG0bt2azZs3k47hJ7KJRo0a+WItGcmlTWFhyEmuzsMu+BFNxmRuHV21cCGPBzmV+PPQP/7BjX7ndw678Ua2HnQQh6fBdlYwqdhiqxeKIC8vj6OPPrquxTCMlJHI8BHJIpqMkbaOajvBlZSWhqz/8bZtrP3DH3zlm37yEx5yzgiM69UrqTLFSypkqheKwDCyjUSGj0gW0WRM5vZWqIilf5kzh4s/+shXbnHbbexo2DDAayhn8uSkyRQvqdgGNEVgGBlKnYddcEE4GWeXlZEjEnLCrs32lncv3Z8eW7awvKTEVx7585/z5tln84cQijMdt9xSIVO98RoyDCMz8O55h1ICtdne8o7rnTSlupp/zZzpUwJbCwpoNGECs7p39+2zzy4rCxjj3j59aJyXF1AnQL/2NYIip4xQMiV6G9AUgWEYKSXUnjd48g6UDBwY9yrHf9ze69dTfdddnLVpEwD9f/1rDrnlFn7wm1BDubIO69KFEd26BcTGUWDWypU1lEaqGNalCyUDB9K2sBAB2hYW1uo5hcK2hgzDSCnh9rarVWs1uW0qL6dBVRWf/eEPHONE8F1x6KEUjRlDdZjwEKFkWbRmDcFrlbo2GCd7G9AUgWEYKSVZe95j1q/n8VmzfOXTLruM99q0oa0zrtt7ZsIZjURjisAwjJSScNfX3buheXMed8Zb2L49A379axAJ2N9/YtmygF/6+bm57Nq3j5zJk2lRUADA9oqKpBix0x1TBIZhuGZ2WRnXL17MtooKAFoWFPDLTp1YtGaNazdWf7fSjeXl5IoE7NeH6xvydO0778BVV/na3D19OpP8DpYq8KflyxGRGts9VdXVvs/h/RdCu56m2xmNRGOKwDAMV8wuK2PUvHlUVlf76rZVVASc4A0+9RouNIJ3snd7Yjb4dO3OL79kWNeuBxpccQX88Y88NXVqjQnfX15/Qk34/uSKUK2almc0Eo0pAsMwXDFhyZKwk6o//r/uI0304U7MDp87lwlLlgScQPbf35/45ptMfvPNA3LNnMm9o0YBid3Hr1aletKkhI2XzpgiMAzDFbFMspvKy6NO9JECwm0sL+eyV15BVX3Kp3V5OZ///ve+NpPPPpvic86BTZvYsXAhM/r3D2uIjof6bBMIxs4RGIbhCq9B1Q1tCgsjKo6N5eWR89gC+6qqfErg8QULApRAq5tu8igBhyed2EKhDl/FS10eIks1tiIwjDQgHTJ5JQqvYTXar37Fc2o30k59x2++4eMZM3zlq/r14/GTTqrRrtrZ7w+ObxTOA8gNJaWlnN6mTcb+HWLBFIFh1DHpkMnLDf6eNcG0dVYAwUos2E00GMVjlK0xWauy4C9/YcCaNQBU5uRw8K23sic/P6qc/sbo4GdLuPuFoEo1Lf8OycC2hgyjjkmHTF5uyA2T9CknRL13hbOnsjJsP/AokOqgSfnkzz9HJ0/2KYFfXngh+RMnRlQCTcJsB4UKzzDrggsYV1QUUS4veyorGfG3v9VZeIlUYSsCw6hjMuUka7hf0dWqvi2gUEbeKlVygFD+Rv3at2fRmjVsLC8np7qa0iefpPvXXwOwobCQ4669lsoGkaepHODJgQPDXg8VnmFYly6eZPWEXjX4kw0rA1MEhlHHpFPo40i2irYuPXL2VVXVqAvndFpSWsroXr3Y/OKLzH/2WV/9+SNHsqxjRyojbEd5GVNUVKsJ2ts3OKm9P96VgX/72th10s0mZFtDhlHHpCLMsBv8wzgr1AjVnEiPHC+5lZUU/+pXPiXwzlFHkTNxIq8fcwy/7NTJ1f0SERl0WJcujO7VK2Ib78pgdllZ1GcVidr0TRamCAyjjklFmGE3RLNVhJKzZQwupcHc+vbb/HDPPfxo924Aiq68kjMvvxzNyaGyuppFa9b47heJRNlTFjk2CTf3qo1dJx1tQrY1ZBhpQDpkG3NjqwiWM9T+en5uboCNACAvJwcRYV9VFa127+bbBx/0XXuhUycuvvBCCDLebiwv5/rFi9leUUHT/Hx27dsXs+yx4HaMSO3cjJGONqGkKQIROQp4BjgUj5dYiapOC2ojwDSgH7AHGKmqy5Mlk2GEYnZZGWMWLGC3M5nliDDGyWEb6zih9n1D1cMBX3f/yJeJ3i+OtBftveYN/BbOoVKBdlOnhpQrVF7ifu3b89Lq1QGB6ab17cvSTZs45+abuejjj339T77iCj5o3Tqs/N4xIikBSIw9xe2pZO+94rXrpJNNyEsyVwT7gRtVdbmINANKReSfqvqxX5u+QHvndTLwuPOvYaSE2WVlXDp3boAxs1rVF0jNrTIIdxZg6aZNzFq5MqA+2KvG3z8/kWcIIp1PgEAf/2h+9aHkClYyzw4ZUmNcgIr9+/n0gw+YccUVvrpPWrXihGuu8ZXznO2geEiUPSVUeOxI94o3lHbCw3AnANE4T93FfCORV4DHVPWffnVPAm+q6vNO+TOgt6p+GW6coqIiXeYX7dAwakO7qVPD/grMFWH/xIm1Gsft4aVg2hYWsmH8+Jj7uZEpUqIWt3KF2hJqnJdHQYMGNQ6eLX/iCXp89ZWvfNw117CmVauANuOKigKimMYiTzJXUF731kgrqkzxGhKRUlUtCnUtJTYCEWkH9ADeD7p0JPC5X3mzUxegCERkNDAaoE2bNkmT08g+Iu3LxjKBhxsn3vAGydzzrs3Y3r7hDJ7+dT23bKHUSRwPsOC44xj061+HHHfRmjW0LCiIeHo5mEQoy2BisdXUxq6TDjYhf5KuCESkKfAyMF5Vv49nDFUtAUrAsyJIoHhGlhNpX9jNydNo48S7IkjmnnekPW63ckVTJrvvuYfG+/f7yofdeCNfN2sWtv3G8nLywuQVDkVdb6XUN5LqPioieXiUwGxVnRuiyRfAUX7l1k6dYaSEe/v0CfufIJpfefA4oc4CjO7Vq0Z9fm5uxEkvkXve4c4n3NunT9Ton6HYWF5OqwceCBuJ9MJNm9DiYp8SePSkk5Di4ohKADwKM5SNQEK8ryv32vpM1BWBiFwLPKeqO2IZ2PEIegr4RFUfCdNsPnCNiLyAx0hcHsk+YBiJxjuZ1NZrKJT3jHff9/Q2berMa6igQQPfdo3Xe8c79tJNm2rk8XVDqO0bqa6m+q67Auqa3XYbuxo2jDpefm5uyNPIEBiZVDmgyEwJJJaoxmIRuQcYCiwHZgL/UBcWZhE5A3gbKOPACfPbgTYAqvqEoyweA87H4z46SlUjWozMWGwY0QlnzA3+JT27rCxiaAU3XLpiBbPmzfOVf/uTn/Dw6ae77p+Xk8NBDRu6tg8E2wbSLVxDuhLJWOzKa8iZsH8KjAKKgJeAp1T1f4kU1A2mCAwjOpE8hoINrDmTJ8e8KgDI37+fH+65J6Au78472Z+bG/NYLQsKqNi/P6LrphcBXwpJtwrPiKwIXNkInBXAV85rP9AcmCMiDyRMSsMwEkYsHkPxGKa94SG8/HrIEKS4OC4lAJ5tMbfhK/zlTcdwDZlIVEUgIteLSCnwALAU6KKq44BewC+SLJ9hZDWzy8poN3UqOZMn027qVNeBycJN7qHqYwkmd9DevWhxMff5TbQyaRLPd+3qqn842hQWMqxLFzaMH+87mLatoqKGQTvYkB5O4W0sL4/peWU7blYELYAhqvozVf2rqlYCqGo1MCCp0hlGFlObKJWhJnch/ARZECXmP3jyBpdPmeIr97n0UqS4uEaMoFjJy8nxTe7+nxkOpLOE0N5CkVYz6RDVM1Nwc47gGFXd6F8hIs+q6iWq+kmS5DKMrCfStke0/W9/LyZvonivHSBSqIlQHP7992x55IDj3/aCAlrdcgstCgrI37cvrMePWw5q2DBA3mBZlPCHx6KFhXD7vLIdN4qgk39BRHLxbAsZhpFEansy2Ht6NZTh2H8fPZIS+Mczz/DTdet85W5jx1LeoQPVzqTsH7guXrb7eQvF+pmDFV4sfY0DhN0aEpHbRGQn0FVEvndeO4FvgFdSJqFhZCmx7PNHItLkGu7a8d9+ixYX+5TAe61bI8XFrDrssBphqTeMHx/X4TQv/p8nns/slSFc3oK6jOqZKYRVBKp6n6o2Ax5U1YOcVzNVbamqt6VQRsPIShKVuSzS5Brq2rqpU/lk+nRfue348ZzmFzk0VJ9w92hZUEB+BE+i4M9Tm8+cLpneMpFIK4KeItIT+Kv3vf8rhTIaRlaSqMxl0UJNeK+d5oSHOPq77wB4rksXpLiYTQcfXKOf23tM69uXmYMHB7iCRjL+1uYzp0umt0wkko3g4QjXFDg3wbIYhhGEmyiVbhLitCgoIEfEl+DF30uoIDeX3RMmBA66bRvDW7Tg3YULebK0lGrn4Gm4LaBIITb8r4eSvd3UqSFDbTw7ZEjMk3i6RfXMFFKWjyBR2Mlio74SPHFD9PhD4U7Wntq6Na+vXx/xxHBeTg6DP/mEvz7/vK/u/t69af3ooz5Fctkrr9TwCsrLyeHPP/95whPnBGMnhBNLrfIRiMiloepV9ZnaCmYYhofgSTE4a9klc+eydNMmZvTvH6AwRMT3a93LnspKlqxfH/F+uVVV7CsuDqgrmDCBvXl5tHW8icLFIKqsrk6IS2YoV1F/9lRWMuJvfwNqn63NiIwb99ET/d43AvrgCUBnisAwEkS0SVGBJ5yVsH/qy3hW9OM++IAZixb5ymP79+fJEw/8N/eeM4gUiC5VyeKrVBOWutMIT1RFoKrX+pdF5GDghWQJZBjZiJtJUYGS0tK4I4UW7NvHnt/9LqAud+JEqoNyI+SKRA3+lspk8XYoLPnEk5hmN3B0ogUxjGzG7cQarxL43WuvBSiBwUOHIsXFNZSA23v0a98+Ljn8iSXGkR0KSy5ubAQLOHA6PQc4AU8YasMwEkS0UAnx0mLPHrY9EBgkWCZNqnV8oEVr1tSqPxzY6nGTD8EOhSUXNzaCh/ze7wc2qurmJMljGFlJsPtlk/x8n6tnvPxlzhwu/ugjX/nUyy/n30cdRaPcXPbWMj5QbUJK+OP93NG8h+xQWHJxowg2AV+q6l4AESkQkXaquiGpkhlGAkh29io3PvzR3D/DtfO/liPieluo7Y4dbJg2zVf+X/PmHHv99QCMKypi0Zo1tZ7Ic2u5ovAnWAkmM3WnERo3qSqXAaep6j6nnA8sVdUTI3ZMEnaOwHBLsrNXhRt/RLduAZ494e4bi3xus4h9UFLCiVu2+MrHX301nx1ySMD4idp+UidLmJEZ1DZDWQOvEgBw3ucnSjjDSBbJzl4VbvyS0lJX941Fvmh75N2//BItLvYpgcXHHosUFwcoAe/4ifg1Hy7Am5GZuNka+lZEBqnqfAARGQxsTa5YhlF74g1L7HZbJ9w44bZwgttHyq41u6ws4J6hjMneHAPfTZlC4d69vvrDb7yRr5o1C/fxapWoHmzPvj7iRhGMBWaLyGN4vnufAyFPGxtGujC7rCwgGYs/kX5dB2/X+CdxCVYG4fzgc8Ps5wffN5If/WWvvML1ixf79sn7tW9PQYMGPrlyRDhn7Vpee/ZZX5/pJ57INf37h/1s8dDWufeiNWuSZmcx6h43B8r+B5wiIk2d8q6kS2UYtWTCkiUhlYBAxF+zsWQFC/UrPZKNIPi+/dq35/Ew9q59VVW+MBMby8sD2kl1NVV33RXQ/qBbb2Vno0YBdfm5uagqldXVYT+vd4snlEIKlxXMqH+EVQQiMlxVnxOR3wTVA6Cqj4TsaBhpQLhtFyVyqIJYMmRFirh5eps2UbeX4vHFH7ZyJc858XcAbu3Th/vPPDNk22b5+QExi0Kxsbyc54YMCanQbPsne4i0Imjs/Bt+s9Ew0pRw2y7RjJzh+oXaTopkS/APh+xtd8ncuQHtYjktm7d/P/vuuSegLv+OO6iMkHQ+mhLwp2TgwKS62RrpTVj3URG5X1VvEZGLVPWvKZYrLOY+arghXtdRt/1q0w4Ia78Ixc3vvMP9r73mK19ywQU8162by97RaVlQwNabb07YeEZ6Em8Y6n4icitwG5A2isAw3BAtUUpt+7m1JYSLKupGCTT94Qd23ndfQF3OxIloiPhAtSGWlYNRP4m0IngQuBJoCuzxvwSoqh6UfPFqYisCI1EEb+2E8o4BQraJdDK3bWGhr328J3j/sHAh13z4oa/8k0su4bUf/ziusdxgh8PqP5FWBJEUQUNV/UFEXlHVwXHcdCYwAPhGVTuHuN4beAXwZtCYq6p3BbcLxhSBkQiiZccCTyYuEamRoSsSwVs+sWwBARy2cydfPnwgS+z3+fkU3n57DCPEjm0NZQfxbg29B/QEvo/zvk8DjxE5gc3bqjogzvENI26iJYIBIrpdhiLUpB+LElj43HP0W7vWV+4xZgwrDj88oE3jvDwKGjRI2HZOfm4u0/r2TchYRuYSSRHki8ivgdNEZEjwRVWdG2lgVX1LRNrVUj7DSAqJjm/fsqAg4uQcaWVw3NatfPbYY77yssMP58QxY3zlHCcdpfdw10urV8clox0OM8IRSRGMBYYBBwMDg64pEFERuORUEVkJbAF+q6ohv+EiMhoYDdCmTZsE3NbIdmqzfx+Kiv37wyqDts6E67U1NM7LY7ezGvns0Uc5bvt2X9ujr7+eDc2bB/SvViUvJ4d+7dvXOKgWTMuCArofdliNxPXecwE26RuhcBN99HJVfSquwT0rglfD2AgOAqpVdZeI9AOmqWrUtEdmIzASgRsbQazkiJArErClFNZl9b334LTTfMUXOnXi4osuijh+uNAVbtvZSeHsJi4bgYicq6qvAzvi2RqKhqp+7/d+kYjMEJFWqmoB7YykEyoG/g9VVb5kMNG2ekJRrUq1qm8bKMfJ/Tt87lyuX7yYaX37MqxzZwhy/2x1001sa9Ik6vhug8W5DXpnGF4iOSSf7fw7MMSr1gZeETlMnHgVInKSI8u22o5rGG4Z1qULG8aP59khQ6jYvz8gI5h3qycevNNwtd+EvK2igpenTAlQAlNOPx0pLnalBKD2yWAs3aMRjrArAlWd5Pw7Kp6BReR5oDfQSkQ2A5OAPGfMJ4ALgXEish+oAIZqtH0qw0gC4Q6HFTRokJBELjkhgsQ1vv12KvLdp/XIy8nhip49o9oIwmGxg4xIRIw+KiId8Bhpj3eqPgFKVPW/0QZW1YujXH8Mj3upYdQp4bZMtldU8OyQIQyfG/8u6Ohly3jy1Vd95av79WPGSSfFNEaOCH/++c8DgtnFauhOVFY2o34SyUZwKh7PoBLnJUAP4E0RGaKq/06NiIaRXMJ5ELUoKIiYzaxlQQEV+/eH/IXeqLKSinvvDahrcOedVOXmxiSbAGN69aqRB9nrCupmhdC2sNCUgBGRSCuCicDFqvqmX908EXkdzzaPnUIxMoJoGcf6tW/PE8uW1fDz31ZREdVgPKJbN15avTqg3V2vv86db73lKw/55S/52wknxCW7ArNWrgTnX/+EObNWrmREt24RQ15Ey79gGBA5xMR/VfW4MNc+U9UOSZUsDOY+asRCtCihtXUj9Y4F8MD8+ay8446A6zJpEiQgR3A0l9BIye0tjpAB8YeY2Bnh2u7aiWQYqSFalFA3oSYi4R1rw/LlDPNLG3nGqFEsbds27nGDieYSGm/+BcOAyIrgKBF5NES9AEcmSR7DSCjRMo7V1re+zXffsaG42FfeffjhNPULD5EoouVBDpc207aFDDdEOkdwE1Aa4rUMsFCFRkYQznfeW18b3/p3nnqKjVOnHqhYvZomW7ZEPH+Ql5MT8/mExnl5jO7Vi8Z5eTXqvRP9sC5dKBk4kLaFhQielYB5ChluiRpiIt0wG4ERDX/jsH9cHy/ek78tCwrYu39/jevR6PrVV6x84glf+V/HHsvmuXMZ1qULVy1cGNLwDNAkL49GDRqwvaKCNoWFVFZVsWXXrhrtcoDmBQW+dl7jdjSjt2FEIq58BOmKKQIjEsmIIeTPNw88wCF7DuRpOvI3v2HLQZ4cTZEijDbJy6Oyutp1bgOLC2QkmniNxYaRcdTW+BuMd3I/Z906Xn/mQGqNkp49GTNoUEDbSD+pYl11WFwgI5VEVQQi0kJVtwfVHa2q68P1MYy6ItETqKqikycH1BXeeivfN2qU0PsEY3GBjFTiJgv2AidkNAAicgKwIHkiGUb8JHICHVpWFqAEJpx7LlJcnHQlYIfAjFTjZmvod3iUQX+gA57Uk8OSKpVhxEkoN8pYaVBVReXddwfUNbzjDvY1SM1O6tiiIjMCGykl6jdbVReKSB7wf0Az4AI3QecMoy4IlWcAPAHkciIkdhHg3KOPpvfLL3PH/Pm++pGDBzOrR4+kyOr1IvKGp2hZUODJWWBKwEgxkYLO/YFA+1ch8D/gGhFBVa9LtnCGEQ/DunQJOZmG8yhqWVDA9LPP5lennBLYfsUKnl+wAEJ4+rQsKOCXnTr58v+2COGKGsmLqHFeHk+an7+RJkRaEQT7aJYmUxAje0mVf3zwaqGNE8Gz6IEH+NUtt/jajRg9miaDBrHojTfYV1XlSx4P7n61h1I4eTk5HNSwYY2zAYaRDkQ8RyAiucAzqpo2NgE7R1C/iBYULpm8/Oab/OKcc3zligYNaBwUNC4emdpNnRo27o+dDTDqikjnCCJ6DalqFdBWRNynUjKMGIgUFC6pDB4coASKrrwyrBKIVaZo8Y0MI91w4waxDlgqIvPxizqqqo8kTSoja0j5pLlmDRx3ILr6ikMPpce4ca66upUpXCRQOxtgpCtuzhH8D3jVadvM72UYtSZaULiE0qlTgBI4a8IE10ogFpnu7dMnYoA4w0g33LiPTgYQkaZOuWaULMOIk5SET37/ffD3CBoyBF5+mTFlZZS6PHMQi0yhjNJmHDbSGTchJjoDzwItnPJW4FJVXZ1k2YwsIKmTpirk5QW6f37zDRxySI17bywv98X89+YD9rqGxiNTOBdWw0hHokYfFZF3gQmq+oZT7g38TlVPS7p0ITCvIcMVCxfCgAEHyjfeCA89VHfyGEYdU9voo028SgBAVd8UkSYJk84wEsjslSsZ1r17YOWuXdAk/FfW7TmGUO3AtoCMzMeV15CI3IlnewhgOB5PIsNIK/5dXMwwvyBx159/Pn8680xK1q0LOzkHn2PYWF7O6AWemIr+fUK1GzVvHiLiyzEQrq9hpDtuvIYuAw4B5jqvVk6dYaQHe/eCCKf4KYEGd97Jo6ecEtX/3+05hlDtQiWaSckZCMNIMG5WBC0srpCRtkyeDH7J43954YX8tXPngCaR/P/dnmOI5VyDHRwzMg03imCmiLQGPgTeBt5S1bLkimVkG7PLyrh+8WL3kTh37IAWLQKq2j3yCBu//75GUwWuWriQGf3717jm9vBXuHahsINjRqYRdWtIVc8GOgJ/AA4GForI9oidABGZKSLfiMhHYa6LiDwqImtFZJWI9IxRdiODmF1WRrupU8mZPJl2U6cyu6ws4Nplr7ziUwIA2yoqGDVvXkA7H5ddFqgE3nwTVLn3vPNqHOTy8viyZVy1cGGNereHv0K1y8vJIT83N2pfw0h33JwjOAM403kdjOeU8dsuxn4aeAxPIptQ9AXaO6+Tgcedf416RjiD7NJNm1i0Zk3YX9qV1dW+/fYJS5ZQvWkTm37/+wMNDj8ctmzxePNMncqm8vKwigCgpLS0xqrA7TmGcO3c9DWMdMfNOYL9eEJQ3wcsUtV9rgcXaQe8qqqdQ1x7EnhTVZ93yp8BvVX1y0hj2jmCzCNcNM5I8fr9aZyXx8KSEnpv3Oire/XllxkwZEjYHAPh0EmTXEptGPWL2p4jaAWcDpwFXCci1cB7qnpnLeU6Evjcr7zZqauhCERkNDAaoE2bNrW8rZFqwhlP3SiBrt98w8oZM3zlN9u25ZxRo2i7aRMDCO3NE45cEVftDCPbcBNr6DsRWQccBbQGTgPCr7+TgKqWACXgWRGk8t5G7YnF0OrPFw8/zBE7d/rKrW+4gS8cQ6xXucTioTO6V6+YZTCMbCCqsdhRAg8DzfHs43dwDMi15Qs8ysVLa6fOqGeEMrRG+m1+1oYNaHGxTwnM7N4dKS72KQE44Jnj1kNnXFFRSK8hwzDcbQ2NVNW3/CtE5HRVXVrLe8/Hk//4BTxG4vJo9gEjMwmXInLWypWBUUcbNGB3UHKYl955h2v/9S8IE500VPRSf1KV7cwwMhk3imAqEOza+YcQdQGIyPNAb6CViGwGJuFsKanqE8AioB+wFtgDjIpBbqMOiBaTJ9L1YGWwaM0aRnTr5ovwOXb9embMmnXgZsXFzB4ypIYNIFeEEd26hR23RUEBgOUGNowYCKsIRORUPPaAQ0TkN36XDgJyQ/c6gKpeHOW6Ale7lNOoY6LF5Inn+qyVKxnVuTOPDR4ceLOKCmavWRPyl36VKrNWruR0x2nA/xBajgjbKipoW1gYNcG8YRgHiGQjyAea4lEW/pnJvgcuTL5oRjoRLSZPPNevePvtACVw+aBBNLn3XmavWRPRG2hPZSXXL17MqHnzAg6hVTuu0F4lFPIwmmEYNQi7IlDVfwH/EpGnVXWjiDRW1T0plM1II6LF5InleuN9+9j9u98FtMuZOBHNyQFHeUTzBvJXAKHwKiFbFRhGdNxEHz1CRD4GPgUQkW4iMiNKH6OeES23sNvrD/7jHwFKoP+vf40UF3uUgIPXxlBbLPibYbjDjSKYCvwM2AagqivxHC4zsohoMXmiXX+oRw+0uJjfvvceAPtFyJk0iUV+yeS9eI284cJFNM7Lo6VjFI6EBX8zDHe4UQSo6udBVVUhGxr1lmFdulAycCBtCwsRoG1hYYBbZsTrv/gFF/bu7Rvr5Cuu4NhHHmHsiSeGVR7+48GBU8Hecaf17UteTvivrwV/Mwz3uIk1NAd4BE8AuZOB64EiVR2afPFqYrGGMoj//Q+OPfZA+YQTYPXqgCZu00SGIjh0dY4I1U7yeXMbNYxAahtraCwwDU8coC+A/8PcPg2HcBP5jg4daP7f//razV+wgEEDBvjabywvJ1eEKmfifnbIkJgn7mFduthkbxgJIOqKIN2wFUH6ECry5+lff807jz/uK8/r0IELLr6Yxnl5jOjWrcZpYi92AtgwkktcKwIRmRhhTFXVu2stmZHRBPv67737bhr65fA99Le/5ZumTQGPO2dJaSlVYX54mLunYdQdkbaGdoeoawJcDrQETBFkOV73zPPXrGHx7Nm++qknn8wNffvWaB9OCQSPZxhGaol0oOxh73sRaYbHSDwKeAFPNFKjnuG/3+8fs6dFQQG79+1jb1Wgs5hUV6N33RVQ1/S229jdsGHI8b02gXCYu6dh1A0R3UdFpIWI3AOswqM0eqrqLar6TUqkM1KGd79/Y3k5iufk7raKCt/7YCUw4j//odpPCfzmpz9FiovDKoHGeXmM7tWrRo5fL/m5uebuaRh1RCQbwYPAEDwJYbqo6q6USWWkhKsWLuTJ0lJfjB435O/fzw/33BNQl3fnnewPM8HDgYihM/r356XVq0OGh2iWn2/2AcOoIyLZCG4EfgDuACbIgTR/gsdYfFCSZTOSyFULF/J4jN5Xt7/1Fve+/rqvPPQXv+BFF5O3f8TQ7WFiBIWrNwwj+USyEbg6dWxkJiWlpa7bHrR3L+VTpgTUyaRJEEMOYK9XULi0lWYfMIy6w82BMqMeEHwK1y0l8+dz5fLlvvK5l17KG8ccE5cMm8rLeXbIkBpnDywchGHULaYIMpBwp3NDhVWIVwEc8f33fPHII77yt40b86Obb66V3G0KC0OmrbRwEIZRt5giyDCCT/NWBSVjgQPpG0Od/HXDP2fN4rz1633lrmPHUnbYYbWS2/9Xv4WGMIz0wkJMZBjtpk4NucfujzdEc6yrgI7ffMPHMw6kmlh61FGccfnlsQvph4D96jeMNKC2QeeMNMLN6dtYFQDAht//nrZ+Y7cZP57PDz445nGCiSeYnGEYqcU8gzKMRHvXnLFxI1pc7FMCz3TtihQXJ0QJAL6cxYZhpC+2Isgw7u3TJ659/xqoopMnB1S1uPlmdjRuHPNQbcO4hILFDzKMTMAUQYbh3WYZPndu3GNc8PHHzH3pJV/57rPOYuK558Y8jgDVkyYB4W0Xdj7AMNIfUwQZSrQAbiH7VFWx/+7AoLEFEyawN0xu4Gj4T/KhVip2PsAwMgNTBBmG1yU0ViVw9fvv89jixb7y6AED+GNRSAcCVwRP8nY+wDAyF1MEGYB/eOicGFcCBfv2sed3vwuoy504keoIid+jEe7wmp0PMIzMxBRBmhPuAJkbpvzzn9yydKmvPGjoUBYcf3yt5GlbWMiG8eNrNYZhGOmFKYI05/rFi2P2EGq5ezdbH3wwoC7WIHHghJn1K9uev2HUT5J6jkBEzheRz0RkrYjcGuL6SBH5VkRWOK8rkilPpjG7rCzmw2Ev/PWvAUrg1MsvR4qLY1YCjfPyGFtURNvCQgTPSsCSyxtG/SRpKwIRyQWmAz8BNgMfish8Vf04qOmLqnpNsuTIZGI5jNVuxw7WT5vmK/+3RQs6XHddzPe0kBCGkX0kc2voJGCtqq4DEJEXgMFAsCIwwhDpMFaOiC+z2IdPPknRl1/6rnW45hr+26pVXPf0ngswDCN7SObW0JHA537lzU5dML8QkVUiMkdEjgo1kIiMFpFlIrLs22+/TYasaUm4w1gtCwqoVqXHli1ocbFPCSxs3x4pLo5bCXiD1RmGkV3UtbF4AfC8qv4gImOAWUCNI66qWoIndzJFRUWZFS61FoQ7pDWtb18Gn3giTX/4wVd/2I038nWzZnHfq0FODtP69q2VvIZhZCbJVARfAP6/8Fs7dT5UdZtf8U/AA0mUJ+MIdUjrqYIC+nTt6mvz2Ikncm3//q7Gy8vJIT83l91BXkgtCwqY1rev2QQMI0tJpiL4EGgvIkfjUQBDgV/7NxCRw1XVu7k9CPgkifJkLLv27YPqajbccENA/UG33srORo3C9ssBmhcUsL2iwgzAhmGEJWmKQFX3i8g1wD+AXGCmqq4WkbuAZao6H7hORAYB+4HtwMhkyZPO+J8cbuHs02+vqKBxXh67Kyu5ZMUKnpk3z9f+5vPO48Ezzog6bjXQND+frbVMMWkYRv0mqTYCVV0ELAqqm+j3/jbgtmTKkI4ET/w79+1jX1UVEJhUZl9FBXrPPQF98++4g8oG7v9sFgbaMIxo1LWxOOsIDhkR7sDYLW+/zRS/cwTDhgzhL362AbdYGGjDMKJhiiDFTFiyJGLIiGZ79/L9lCkBdfGEhwDP4TALCWEYRjRMEaSQ2WVlERPPT3/1Va5atsxXPu+SS1jy4x/HdS8BxhYVmXHYMIyomCJIEd4toVActnMnXz78sK/8XcOGNL8tdtNJrnPa2DyEDMOIBVMEKSLcltCi556j79q1vnKPMWNYcfjhMY/fOC/PgsIZhhEXpghSRLD3Todvv+XT6dN95fePPJJTrrwypjFtBWAYRiIwRZAkZpeVcf3ixT6vIH9T75pp0zh2xw5fud3117OxefOYxrcVgGEYicIUQRKYXVbGZa+84jsbAJ4EL6du2sS7M2f66v7SuTPDLrzQ1Zh9jj6atdu3Wz5gwzASjimCBOB/QKxNYSG7/A6IAaCKTp4c0KflzTezvXHjqGOHyw9sGIaRKEwR1JLgA2LB7qGDP/mEeS++6Cv/7owzmHDeeRHHzBVh1gUX2ORvGEZKMEVQS8J5A+VUV1N1110BdY1vv52K/PyI4wmYEjAMI6WYIqgloWL5jPnwQ55YuNBXvqpfPx4/6SRX4ymYEjAMI6WYIqglbQoLfdtBjSorqbj33oDruRMnUp3jPhFcW4sNZBhGijFFEAdXLVzIE8uW4Z8q7e4lS7jj7bd95Qt+9SvmdewYdoz8nBwa5ObWyD5msYEMw0g1pghi5KqFC3ncLx5Qiz172PZAYGK1aEHicoCZP/85EJh9zLyDDMOoC0wRxMiTpaW+98++/DLDy8p85TNGjWJp27YR+we7g9rEbxhGXWOKIATB5wKObdGC19ev920Ftd2xgw3Tpvnarz/4YI4ZPz7imPm5ucwcPNgmfsMw0g5TBA7eyX9jeTkCvkl/Y3l5wNmA9/74R0754gtfuePVV/PpIYdEHNuSwxuGkc6YIqDmoTAN0abrV1+x8oknfOX/O+YYfnbppRHHHVdUxIz+/RMpqmEYRsIxRUD0rGFb77+fln4pJY/4zW/48qCDIo5pSsAwjEwh6xXBVQsXhs0adu66dSx55hlf+YlevRg3cGDUMU0JGIaRSWS1Igh2BfURIkjcQbfeys5GjSKOlyvC6F69TAkYhpFRZLUiKPFzBfVy8apV/GXuXF/59nPP5b6zzoo4jq0ADMPIZLJWEVy1cCFVesAs3KCqisq77w5ok3/HHVQ2iPyI+hx9tCkBwzAymqxUBOc98wxL1q/3lW9cupSH/vlPX3nEz3/OM927RxzDXEINw6gvZJUi6DR9Oh9v3eorN/3hB3bed19Am5yJE9EIQeKa5ufzxIABpgAMw6g3ZI0ikCDj77RFi7jugw985Z8NH87/HXts2P5N8vJ40nIEG4ZRD0mqIhCR84FpQC7wJ1WdEnS9IfAM0AvYBvxKVTckXA4/JXDozp189fDDvvLuvDyaTpgQsb8Zgw3DqM8kTRGISC4wHfgJsBn4UETmq+rHfs0uB3ao6rEiMhS4H/hVQuXwUwJvP/UUZ3z+ua/ca/Rolh9xRNi+tgowDCMbSOaK4CRgraquAxCRF4DBgL8iGAwUO+/nAI+JiKhqqCgPtWL4ypU+JfCfww6j59ixYduaHcAwjGwimYrgSOBzv/Jm4ORwbVR1v4iUAy2Brf6NRGQ0MBqgTZs2cQmz+Nhj2dGoEeddemnYVUCeCPsmToxrfMMwjEwlI4zFqloClAAUFRXFtVrY1qQJLW69Nfw9Jk2KTzjDMIwMx30y3dj5AjjKr9zaqQvZRkQaAIV4jMYpxZSAYRjZTDIVwYdAexE5WkTygaHA/KA284ERzvsLgdcTbR+INMnrpEmmBAzDyHqStjXk7PlfA/wDj/voTFVdLSJ3ActUdT7wFPCsiKwFtuNRFomXxSZ7wzCMsCTVRqCqi4BFQXUT/d7vBS5KpgyGYRhGZJK5NWQYhmFkAKYIDMMwshxTBIZhGFmOKQLDMIwsR5IQzSGpiMi3wMY4u7ci6NRyGpGusplcsZOusplcsZOussUjV1tVPSTUhYxTBLVBRJapalFdyxGKdJXN5IqddJXN5IqddJUt0XLZ1pBhGEaWY4rAMAwjy8k2RVBS1wJEIF1lM7liJ11lM7liJ11lS6hcWWUjMAzDMGqSbSsCwzAMIwhTBIZhGFlOvVEEInK+iHwmImtFpEYGGhFpKCIvOtffF5F2ftduc+o/E5GfpViu34jIxyKySkSWiEhbv2tVIrLCeQWH8E6FbCNF5Fs/Ga7wuzZCRNY4rxHBfZMs1+/9ZPqviHzndy1pz0xEZorINyLyUZjrIiKPOnKvEpGefteS+byiyTXMkadMRN4VkW5+1zY49StEZFmK5eotIuV+f6+JftcifgdSINtNfnJ95HyvWjjXkvnMjhKRN5w5YbWIXB+iTeK/Z6qa8S88Ya7/BxwD5AMrgROC2lwFPOG8Hwq86Lw/wWnfEDjaGSc3hXKdAzR23o/zyuWUd9XxMxsJPBaibwtgnfNvc+d981TJFdT+WjwhzlPxzM4CegIfhbneD1gMCHAK8H6yn5dLuU7z3g/o65XLKW8AWtXR8+oNvFrb70AyZAtqOxBPrpRUPLPDgZ7O+2bAf0P8v0z496y+rAhOAtaq6jpV3Qe8AAwOajMYmOW8nwP0ERFx6l9Q1R9UdT2w1hkvJXKp6huquscp/htPJrdU4OaZheNnwD9Vdbuq7gD+CZxfR3JdDDyfoHtHRFXfwpM3IxyDgWfUw7+Bg0XkcJL7vKLKparvOveFFH7HXDyvcNTmu5kM2VL5HftSVZc773cCn+DJ7e5Pwr9n9UURHAl87lfeTM2H52ujqvuBcqCly77JlMufy/Foei+NRGSZiPxbRH6eIJlile0XzvJzjoh4U4+mxTNzttGOBl73q07mM4tGONmT+bxiJfg7psD/iUipiIyuA3lOFZGVIrJYRDo5dWnzvESkMZ7J9GW/6pQ8M/FsX/cA3g+6lPDvWUYkr88GRGQ4UASc7VfdVlW/EJFjgNdFpExV/5dCsRYAz6vqDyIyBs+K6twU3j8aQ4E5qlrlV1fXzyxtEZFz8CiCM/yqz3Ce14+Af4rIp86v5VSwHM/fa5eI9APmAe1TdG+3DASWqqr/6iHpz0xEmuJRPuNV9ftEjh2K+rIi+AI4yq/c2qkL2UZEGgCFwDaXfZMpFyJyHjABGKSqP3jrVfUL5991wJt4fh0kiqiyqeo2P3n+BPRy2zeZcvkxlKAle5KfWTTCyZ7M5+UKEemK5284WFW3eev9ntc3wN9I3LZoVFT1e1Xd5bxfBOSJSCvS4Hn5Eek7lpRnJiJ5eJTAbFWdG6JJ4r9nyTB4pPqFZ2WzDs82gde41CmozdUEGotfct53ItBYvI7EGYvdyNUDj2GsfVB9c6Ch874VsIYEGsxcyna43/sLgH/rAaPUekfG5s77FqmSy2l3PB6jnaTqmTnjtiO88bM/gUa8D5L9vFzK1QaP7eu0oPomQDO/9+8C56dQrsO8fz88k+km59m5+g4kUzbneiEeO0KTVD0z5/M/A0yN0Cbh37OEPti6fOGxpP8Xz6Q6wam7C8+vbIBGwF+d/xAfAMf49Z3g9PsM6JtiuV4DvgZWOK/5Tv1pQJnzn6AMuLwOntl9wGpHhjeA4/36XuY8y7XAqFTK5ZSLgSlB/ZL6zPD8MvwSqMSz/3o5MBYY61wXYLojdxlQlKLnFU2uPwE7/L5jy5z6Y5xntdL5O09IsVzX+H2//o2fogr1HUilbE6bkXgcSfz7JfuZnYHHBrHK7+/VL9nfMwsxYRiGkeXUFxuBYRiGESemCAzDMLIcUwSGYRhZjikCwzCMLMcUgWEYRhoTLUBeUNs2TtC6/zgRAfq5uYcpAiNrEJFdQeWRIvJYHOP0FpFXXbR7WkTWO1EqV4pIHxd9bo9VHqPe8zTuY1PdgeeMVA8856VmuOlkisAwkstNqtodGA884aK9KQIjAA0RIE9Efiwif3fiHb0tIsd7mwMHOe8LgS1u7mGKwDAAERkonjwV/xGR10TkUKf+bL+49P8RkWZOl6ZOIL5PRWS2E8k2Eu/hFwBMROY5/4lXewOXicgUoMC512ynbriIfODUPSkiuYn/9EYGUgJcq6q9gN9y4Jd/MTBcRDYDi/CEaY+KBZ0zsokCEVnhV24BeJPXvAOcoqoqngQ8NwM34vlPdrWqLnUCge112vfAE55kC7AUON0ZIxzn4wmq5uUyVd0uIgXAhyLysqreKiLXOCsIRKQj8CvgdFWtFJEZwDA8IQiMLMX5Hp4G/NXv90dD59+LgadV9WERORV4VkQ6q2p1pDFNERjZRIV3kgWPjQBPxFfwBOh60Ynrno8nTgt4JvlHnF/oc1V1s/Of7wNV3eyMswJP3JpQiuBBEfmdM/6pfvXXicgFzvuj8ETd3BbUtw+eQH8fOvcsAL6J6RMb9ZEc4Dv/77Ifl+PYE1T1PRFphCfuVsTvjW0NGYaHP+DJxtYFGIMnNhWqOgW4As8kvNRvL/YHv75VhP9RdZOqHgfcAswEj7EZOA84VVW7Af/x3i8IAWapanfn1UFVi+P+hEa9QD1hqdeLyEXgS13Zzbm8Cc8PCO+KshHwbbQxTREYhodCDoTsHeGtFJEfq2qZqt4PfIgn6mk8PAbkiCcndiGwQ1X3OIrlFL92lU4YYoAlwIVO3HtEpIX45bQ2sgMReR6PjamDiGwWkcvxbBFeLiLe4HfeDG43Alc69c8DI9VFQDnbGjIMD8V49lx34Ml4drRTP95J6FKN5z/cYgK3eFzh2B7uwWN76AeMFZFP8ES8/bdf0xJglYgsV9VhInIHnmxYOXgiZV4NbIznAxqZiapeHOZSDZdSVf0Yj70qJiz6qGEYRpZjW0OGYRhZjikCwzCMLMcUgWEYRpZjisAwDCPLMUVgGIaR5ZgiMAzDyHJMERiGYWQ5/w9tkYxsmRS4OgAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "hr.graph()" ] }, { "cell_type": "code", "execution_count": 79, "id": "63cade85-0222-43f9-adc5-6455cc7e98d3", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "LinregressResult(slope=137257.31246468204, intercept=26494280186.8291, rvalue=0.9904640902650609, pvalue=0.0, stderr=482.14977280148685, intercept_stderr=26760703665.486454)" ] }, "execution_count": 79, "metadata": {}, "output_type": "execute_result" } ], "source": [ "stats.linregress(np.array(hr.data.hr), np.array(hr.data.nd))" ] }, { "cell_type": "code", "execution_count": null, "id": "d6d2e165-eeb9-4539-a7ee-f40794bff54c", "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "id": "e6ac5bad-58ec-4670-bada-a10a9b501d79", "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "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.9.9" } }, "nbformat": 4, "nbformat_minor": 5 }