From 031e33f066f4977e7d85f7dc41fc703018f4bf90 Mon Sep 17 00:00:00 2001 From: 6a057d4c0b5da227e313d46089968294 <6a057d4c0b5da227e313d46089968294@app-learninglab.inria.fr> Date: Sat, 3 Sep 2022 17:15:22 +0000 Subject: [PATCH] fitting seasonal fluctuations --- .../exo3/carbon_dioxyde_concentration.ipynb | 186 +++++++++++++++++- 1 file changed, 183 insertions(+), 3 deletions(-) diff --git a/module3/exo3/carbon_dioxyde_concentration.ipynb b/module3/exo3/carbon_dioxyde_concentration.ipynb index a6ab66d..a93b4ba 100644 --- a/module3/exo3/carbon_dioxyde_concentration.ipynb +++ b/module3/exo3/carbon_dioxyde_concentration.ipynb @@ -2216,6 +2216,13 @@ "data" ] }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "Nous indexons les données par le mois durant lequel elles ont été relevées. Nous remarquons qu'elles sont déjà ordonnées par ordre chronologique." + ] + }, { "cell_type": "code", "execution_count": 146, @@ -3200,7 +3207,7 @@ "cell_type": "markdown", "metadata": {}, "source": [ - "## Quelques graphiques\n", + "## Premières observations\n", "\n", "Nous utilisons les données complétées pour la suite de cette analyse (colonnes `Filled`).\n", "\n", @@ -3243,7 +3250,7 @@ "cell_type": "markdown", "metadata": {}, "source": [ - "Nous observons une augmentation globale de la concentration de C02, ainsi qu'une oscillation périodique à plus petite échelle. Zoomons sur une partie du graphique pour observer ces oscillations plus précisément." + "Nous observons une augmentation systématique de la concentration de C02 couplée à une oscillation périodique à plus petite échelle. Zoomons sur une partie du graphique pour observer ces oscillations plus précisément." ] }, { @@ -3321,7 +3328,180 @@ "cell_type": "markdown", "metadata": {}, "source": [ - "Nous ponvons également isoler les fluctuations saisonales et les représenter." + "Superposons enfin ces deux types de données sur un même graphique, globalement puis en se restreignant à une période de quelques années." + ] + }, + { + "cell_type": "code", + "execution_count": 178, + "metadata": { + "scrolled": true + }, + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 178, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXoAAAD8CAYAAAB5Pm/hAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzs3Xd81eXZ+PHPdU422QsCARL2JmAAFVBcgOKeaFuttVU77XpqbZ9atfVX29r1VPtYn1q1FRVUHFVRAUFlCIY9Q0IIJGRvErLP/fvj/h7OQUECZud6v1553d/NdWK8cuf+3kOMMSillOq9XF0dgFJKqY6liV4ppXo5TfRKKdXLaaJXSqleThO9Ukr1cprolVKql9NEr5RSvZwmeqWU6uU00SulVC8X0NUBAMTHx5uUlJSuDkMppXqUTZs2lRljEk51XbdI9CkpKWRkZHR1GEop1aOIyMG2XKdNN0op1ctpoldKqV5OE71SSvVy3aKN/kSam5vJz8+noaGhq0PplUJCQkhOTiYwMLCrQ1FKdbBum+jz8/OJiIggJSUFEenqcHoVYwzl5eXk5+eTmpra1eEopTpYt226aWhoIC4uTpN8BxAR4uLi9K8lpfqINid6EXGLyBYRedPZ/72I7BWR7SLyqohE+117n4hki0imiMw70+A0yXcc/d4q1XecTo3+HmCP3/5yYIIxZhKwD7gPQETGAQuB8cB84G8i4m6fcJVSqpco2QPbX4JOWM61TYleRJKBBcA/vMeMMe8ZY1qc3Y+BZGf7KuBFY0yjMeYAkA1Mb7+QO9fDDz/M+PHjmTRpEmlpaWzYsAGA0tJSAgMD+fvf/37c9SkpKUycOJG0tDTS0tJYt24dBQUFXH/99QBs3bqVt99++9j1zzzzDAkJCceuv/XWWwG4//77WbFiBQBz5sw5NqAsJSWFsrKyNsf/zDPP8J3vfOfMvwFKqY7xytdh6dfh0PoO/6fa+jL2z8BPgIiTnP8asNjZHoRN/F75zrEeZ/369bz55pts3ryZ4OBgysrKaGpqAuCll17i7LPP5oUXXuCuu+467r5Vq1YRHx9/3LGXX34ZsIk+IyODyy677Ni5m266iccee+y46x966KGO+EhKqe6gpRGKd9rtihwYem6H/nOnrNGLyOVAiTFm00nO/xxoARZ5D53gss/8bSIid4pIhohklJaWnkbInaewsJD4+HiCg4MBiI+PZ+DAgQC88MIL/OEPfyA/P5/Dhw9/7nNyc3OZMGECTU1N3H///SxevJi0tDQWL1580nu++tWvHvvlcDLPPfcc06dPJy0tjbvuuovW1lYAnn76aUaNGsX555/P2rVrT+cjK6U6SmOtb/vwZt92TUGH/9NtqdHPBK4UkcuAECBSRJ4zxnxZRG4DLgcuMuZYQ1M+MNjv/mTgM5/EGPMk8CRAenr65zZSPfifXewuqGlDqG03bmAkv7xi/OdeM3fuXB566CFGjRrFxRdfzE033cT5559PXl4eRUVFTJ8+nRtvvJHFixfzwx/+8Nh9F1xwAW63m+Dg4GNNPQBBQUE89NBDZGRkHKvBP/PMMyxevJg1a9YAcM8993D77befMv49e/awePFi1q5dS2BgIN/61rdYtGgRl1xyCb/85S/ZtGkTUVFRXHDBBUyZMuVMvkVKqfay6zV46Ta45SUYNReenu871x0SvTHmPnwvWucAP3aS/HzgXuB8Y8xRv1veAJ4XkT8CA4GRwMb2DrwzhIeHs2nTJj766CNWrVrFTTfdxCOPPEJpaSk33ngjAAsXLuSOO+44LtGfqOnm85yo6eZUVq5cyaZNm5g2bRoA9fX1JCYmsmHDBubMmUNCQsKxZ+/bt++0nq2UamfbXrTlvmWYwdOONXuUho8h3AQQ2sH//BcZMPUYEAwsd7rqfWyMudsYs0tElgC7sU063zbGtH6RIE9V8+5IbrebOXPmMGfOHCZOnMizzz7L4cOHKS4uZtEi21pVUFBAVlYWI0eO7LS4jDHcdttt/OY3vznu+GuvvaZdJ5Xqbgq22LKmgK0Za5gC/Lz5aywqu5jJIVG83sH//GkNmDLGrDbGXO5sjzDGDDbGpDlfd/td97AxZrgxZrQxZll7B91ZMjMzycrKOra/detWWlpaqKur4/Dhw+Tm5pKbm8t9993Hiy++2KZnRkREcOTIkS8c20UXXcTLL79MSUkJABUVFRw8eJAZM2awevVqysvLaW5u5qWXXvrC/5ZS6jQ0HYVnLoc8pyGjpgBqiwCoKz3I8+98CMBHnokA/PGmtA4PqduOjO0Oamtrue222xg3bhyTJk1i9+7dDB8+nGuuuea466677jpeeOGFNj3zggsuYPfu3ad8GXsq48aN49e//jVz585l0qRJXHLJJRQWFpKUlMQDDzzAOeecw8UXX8zUqVPP+N9QSp2B/E8g9yN45Q67n/WeLQfPgJoChkoxrbg5bOJ587uzGJ4Q3uEhiemEzvqnkp6ebj698MiePXsYO3ZsF0XUN+j3WKkOsO4xeO/nEJEEP9pLy9K7IXsFz5v53Fr/b3aEn8vEwEJav7sFt+uLNbOKyCZjTPqprtMavVJKtafybFs2HaWusYWCrK1sOJLI1up+AIxv2ALxI79wkj8dmuiVUuqLyPkACrb69quc1f0aq5n3m9eJPXqALDOIc6dOBsDVUg9xIzo1xG47TbFSSvUI/7rSlg9Ug6cVSjPBFQCeFiY3bSU8qIEFc+eSMH4GOINhOzvRa41eKaXOVJPfEKLmeshaDjWHeTfKzm11Z+IuABJGTIPIgb5rNdErpVQPUbHft12dT876pdSaEH5eeD4AE4+ssbX7xLEQ6DcsKmF0p4apTTdKKXWmynzjbKqLc6nI2UopKRCegHFH4mqsgf4TICD4+PsiBnRqmFqj/xxut/vY9MFpaWnk5uaSkZHB9773PeD4KYAfeOABHn300dN6fnh4x/efVUp1oHJfjX7p6g0MklL6Dx3Da9+eiUQPtSdih/mu/+5m+P6OTg5Sa/SfKzQ0lK1btx53LCUlhfT0U3ZbVUr1RsbYF67uALu953Wa+yURWFdIbWEWSQEVMGI8xIRB9GAo3gHRQ3z3xw3vkrC1Rn+aVq9ezeWXX/651+zfv5/58+dz1llnMXv2bPbu3QvAgQMHOOecc5g2bRq/+MUvOiNcpVR7+vBReHgANNXRXLAdinbwm7oFlJsIzglwJg+MSbFlP2diw6jBJ3xUZ+oZNfplP4Widv5zZ8BEuPSRz72kvr6etDQ7D0Vqaiqvvvpqmx5955138sQTTzBy5Eg2bNjAt771Ld5//33uuecevvnNb3Lrrbfy+OOPf+GPoJTqZKt+DYDJfIdXPtjMQuDthslcG7SCdNltr/E22cz6AYREw/iruyZWPz0j0XeREzXdnEptbS3r1q3jhhtuOHassbERgLVr1/LKK68A8JWvfIV77723/YJVSrWvpqOw6RmY9nUICILW5mOn9mfvpq5oP0fdwaSkjmBs+BjIyrUnvTX62GEw91edHfUJ9YxEf4qad3fi8XiIjo4+6S8InUJYqR5i/eO2Bh8YAulfg8Ltx07lHcgiWUqRmKG8eNe5sMyZaDggFMITuyjgk9M2+nYWGRlJamrqsemBjTFs27YNgJkzZx6bztg7l71SqpuqzLVlsR30RJ5dCtsTGktrVT5TwisITXRernp71gSGQjeszGmi7wCLFi3iqaeeYvLkyYwfP57XX7e/7f/yl7/w+OOPM23aNKqrq7s4SqXU5/Im+qo8AEzuR3iihvDR0aGkUEhCw0EYMMFeE+8sOuRttulmdJriPky/x0o5WppgzxswbI6vt8yfJ9kJyvpPIOvifzJy0XSeaLmCCI7ypYCV9pobnrUvW42BjU/CoLMgufO6X+s0xUop1VbbnrcLhaz4pd33eI4t2u2pPsxfF78BwPutaQTF+fWLj3a6TorAjLs6NcmfDk30SilVU2jL/attWbQdPM3UhCbjaqhkYIMdAfunb13P1XNm+O6LSOrcOM9Qt0703aFZqbfS761Sfo6W27KuhPrGFoqW/4UGCeEPNRcBsDD+AASFMyh5KIHxfjNP9ut+PWxOpNsm+pCQEMrLyzUhdQBjDOXl5YSEhHR1KEp1D0fLbNnaxN1PvU9AzgreaplGZbjtVZNSs8lOXyBy/DQG7p7RQ73bRpmcnEx+fj6lpaVdHUqvFBISQnJycleHoVS3YOrK8XaKrMrbQ3xwDZmeZG655Fx4E/A0++aQD42BiTfAuK4f8dpW3TbRBwYGkpqa2tVhKKV6O2OoLdqHx4QRJUdJd2UC8LMvL4ARE2yiB1+iF4Hr/tE1sZ6hbtt0o5RSHWbpXbDMmYKkbB8RDUUsbZ0NwEVuZ1R74rjjFwuJ7ZqZJ9uDJnqlVN9SWwLbX4QNTwDQWGBHvq7y2AkMz3Xtsi9Z/eeRB4jv3OX/2pMmeqVU31KefWzT1Ffx6ooPAUgaNwsjTkr0vngFmPJliBoCSVM6O9J2o4leKdW7tbbAkSLffvXhY5ufbFxLaFUmRSaGh2+ehfRLsCfC+/uuv/IxuGcbuHpuumxz5CLiFpEtIvKmsx8rIstFJMspY/yuvU9EskUkU0TmdUTgSinVJu/8FP4wGuqcLpQ1+cdOZe3YyGzXDsLHXEiA2wWRA+0J/4FQIj06ycPp1ejvAfb47f8UWGmMGQmsdPYRkXHAQmA8MB/4m4i42ydcpZQ6TZ/8ny0z37ZlxQEIjcEEhBJetplYOUJ4ijN1QXCELaMGdX6cHahNiV5EkoEFgH+foquAZ53tZ4Gr/Y6/aIxpNMYcALKB6e0TrlJKnQZPq2+78iAATQU72FA/kPymcMZ6nMW9vXPWjLrUluOu6sQgO15ba/R/Bn4CePyO9TfGFAI4pXcs8CAgz++6fOeYUkp1rKY6ePOHUOdMaeDfNn+kkMqSw1C0g20tKZQTwSiX017vXdd1xt1wb+7xC3r3AqdM9CJyOVBijNnUxmeeaNb9z8xjICJ3ikiGiGTo6FelVLvY9RpkPAWr/5/dr/G9eG2uOszfn/wrQTTzXuAFDBrkt2i3N9G7XHbkay/Tlhr9TOBKEckFXgQuFJHngGIRSQJwyhLn+nzAf9nzZKDg0w81xjxpjEk3xqQnJCR8gY+glFKO5qO2dKYYpiIHgPLQFArzchjauI+jrnAevvMmEhKdhobAMAiL7YJgO88pE70x5j5jTLIxJgX7kvV9Y8yXgTeA25zLbgOcRRN5A1goIsEikgqMBDa2e+RKKfVp3lWh6qsA8JTspQU37x4ZRnRLGRfGlhI2OI3RSZEQ6fSsiRjQLZf/a09fZK6bR4AlInIHcAi4AcAYs0tElgC7gRbg28aY1pM/Riml2knJbltW5JBdUEbMmqcoMsnkmwQi5SiRR/fD8OvsNd6Rr32gU+BpdQ41xqw2xlzubJcbYy4yxox0ygq/6x42xgw3xow2xixr76CVUgqwS/h5eVqhwJmnpraI3z2zhDip4V+tc7nn2vPs8eY6iHEmS/T2mR/Yc0e8tlXPHgWglOq7mhvgd6mw9i92f+9bUF9BRpidnGx0ne0/cslF8wiO9etFE+sk+pTz4KL74bLfdWbUXUITvVKqZ8p6D+orYfn9ABzZv54mAvhjlU30X0rMBeDiWedClN/aC94avTsAZv+oV/ay+bRuOx+9Ukp9rqpDvm2Ph0N7NiGeQeQb24tvQGUGRCZDUD+I9BvKE9v31rnQGr1Sqmc6Unhs8/WPPiG6LptMk8zt887xXeOdWjgg2L58HTARQqI6OdCupzV6pVTP0FANLY0Q7gzC90v07737JlcFlbPgogsJmjMWNsRCfQXEj/Ldf/daCOib6yRrjV4p1TP881J4dCR47EwspiKH5mi76tP5rm0ABA0YZ691B9oybqTv/qCwHj8L5Znqm59aKdW9eVqhfL/fvgdK7EpQFGyGpjpM4U6eLhtDqxEujTpgz8Wk2LK1yZaJYzot5O5ME71SqvtZ+RD8darvhWvhVt+5gi1sWPJ7XKaZla1TORIYT0Sdc513MrL5v4X0r8HQWZ0bdzelbfRKqe5n+xJb5mfY5L33LRAXGA/vvb+CxroqDkkCA9MuJrpmGRwuhbB4CA63902+yX4pQGv0SqnuyDtrSvFOW+59C4bO5EDAcNx1xaSFlhIzdCK/v36SbyoDb7ON+gxN9Eqp7qWxFmqL7XZVHkcOZ0LpHrb2m8Whxn4kSiXJnsNEJI+zy/95e9b001lwT0YTvVKqeynNPLZZX5rLzx7/FwD3bo6ijCgmunKRlgaId3rUDL/AlkPO7uxIewxN9EqprtXaDO8/DDVOv/i8DQAcTZ5FXckBhrsK8OAiadh40sf79aLx1uST0+G+fDj3u50ceM+hiV4p1bX2vgUf/g7e/AEArQVbqXTH8Y/ceGJay7h2UA2umCE8843zGDp8rO8+/8FQwRHg6v3TDZ8pTfRKqa514ANbOouGFB7YzZ7m/oQlpOIWw+DqT3xJvf8E331hcZ0bZw+miV4p1bl2vQZ5n/j2i52BUNX5VNU1El57kKaIFL6+wM5CKQ3VvkSfNAkGTIKJN/b6VaHak/ajV0p1rpecFUh/aZf7o2SPLZuOcMOvn2F5cA0Jwyb6Bj+B78VrYCjc/VHnxdpLaI1eKdV5mup827XFUJ0HjTU0pV4IwAL3xwCMT59z/Bzy/u3x6rRpoldKdZ6qPN92ZS4epwnnkRw76Onr/bPtucRxtvbuFT+6syLslTTRK6U6zuZ/266TXmX7jm3+76srWbXkr+SbeF5unA5AePl2u+JTaLS9aPqdcN5/QT998fpFaBu9UqpjtDTCG9+x2xOvt1MGO10oAY6W5DA24CAbPGP56P5r4bH74Gi5b6k/gMt+38lB906a6JVSHaOmwLedtwFPQy2uo2Xs6TeDmNp9zAg7zMDmCi678AKCwwLtXDVHy/vkUn8dTZtulFLtI3sF/CoRqvPtvrcEPGXZvP76YgC+XP5Vil39menZAkBw0nh7kXdd1xhN9O1NE71Sqn289wtobYSs9+x+tX3xasTNli0ZNBftpcREU04UkUkjkdYGe513cZAEp4xM6uTAez9tulFKtQ9jbFnkTC18eBPN7jD2yjBMXSnTIoIgdAShpW6iBo6AAkDcEJ1irz//JxA50Lbnq3alNXqlVPuoOWxLZ1WoI5kfsK5pBLmN4cRSwxBzmMTUiex+aB6xyU6/eHegbx1XdyCk327nrVHtShO9UurMNNT4avH1VdBYA0BzeS63/OVtImqy2OAZQ7/YJIa5inA3VEL8SEQEBs+w9w09t4uC71s00SulTt+hDfD74fCGMzXw/vcBaE6YgKfyEOHFGwE475KruHDqON993hGuccPhzg/g6v/tzKj7rFMmehEJEZGNIrJNRHaJyIPO8TQR+VhEtopIhohM97vnPhHJFpFMEZnXkR9AKdUFst6F1ibY8TIYQ8XHz1FkYvlNwRSCaeTHQ7Mx7mDOnnUJRAzw3Rc3wrc9MO34c6rDtOVlbCNwoTGmVkQCgTUisgx4CHjQGLNMRC4DfgfMEZFxwEJgPDAQWCEio4zxLgKplOpxWptt6Q60pbePfEs9BYcPEZW3hjdbz6PYbRP3yJoNSNwICAiGgVN9z/GfqEx1mlPW6I1V6+wGOl/G+Yp0jkdh36EDXAW8aIxpNMYcALKB6SileqajFfDIEHj9275jfn3kX3h5Mf2kkZQxU/j9N64AQOpKIM5ZtDtxLKTMhrQv6+IgXaRNbfQi4haRrUAJsNwYswH4PvB7EckDHgXucy4fBPjNXES+c0wp1RMV7YDmo7DdDniitRlKM2mIsRONxZXZicnOO3sGYQkpvvtih9vS5YavvglXP96JQSt/bUr0xphWY0wakAxMF5EJwDeBHxhjBgM/AJ5yLj/RagDm0wdE5E6nbT+jtLT0zKJXSnW8I0W+7doSOLQe6kr4Zcl5AJwf7ExUFjsMQqJ81yZN6sQg1ec5rV43xpgqYDUwH7gNWOqceglf80w+MNjvtmR8zTr+z3rSGJNujElPSEg4zbCVUp3mSKFvu+oQ5Tl26oKVrVNpkUBSW3PBFQBRn2p/T57WeTGqz9WWXjcJIhLtbIcCFwN7scn7fOeyC4EsZ/sNYKGIBItIKjAS2NjegSulOsi2xfDEbKhxErx3IBRwpPgAy1e9T7mJIG3sKNze6Qqih4Lb6dtx47/h2n/oi9dupC29bpKAZ0XEjf3FsMQY86aIVAF/EZEAoAG4E8AYs0tElgC7gRbg29rjRqkeZNWv7ejW3a/D2XfTmLeFXDOY0ZLH22s2cpYri9q4SfzjtnR4KgmqD9lmG69xV3Zd7OqETpnojTHbgSknOL4GOOsk9zwMPHyic0qpbqy1xdd1svIANB4hoHgbH7ZeTJK7DFfpXkYGHIa02+014f1tGTe8a+JVbaIjY5VSPjX54GkB4EhBJk8/83fcnibeaZ1GmTuR+eFOC22cs1i3M+3BcQOhVLejiV4p5VNqe9C0BEVSeCib1oLtNBPIvXfcwrDhY4hocNrtY1JsOfVWW6ufeEPXxKvaRBO9Un1ZVR48vQBy19j9nS/jCYrkjYY0+kslX0qtI7D/aKaPGABRyb77vKtATbgOfrzPt8ar6pY00SvVl+1YAgfXwIoHaWxuoW7vSt6on0COZwBR1BJaut2ObAVfoncHH99fXnV7muiV6ssOrgfA1JXy5Ovv06+pjE88Yzh78gR7vq7Et/JTlDM8xqXrFfU0muiV6is8rfDrAbD6Ed+xskwATOVBDm95F4CvLryZWel+E5ElOtMMe3vWDJvT8bGqdqWJXqm+ouIAtNTD6t/Y/ZoCqDpEU7+BuPAw07UTjzuYkePTj+8X7226GTQV7lgBNzzT6aGrL0YTvVJ9ReFW37bHg8n4J0Zc/M/RiwGYHZaHK3qIXdovwm+B7uihvu3B0yAgqJMCVu1FG9uU6q28y/yJM89g1vJjp2rL8ylZ+wrFLWNY35IKwRDdeBiSL7QXuFzwrY/t4t0urQ/2dPpfUKneqKUJ/jkP/nY2NB6x+/uWQVg8AH979X0GteSxzQwjbcxI333+89MkjoWEUZ0cuOoImuiV6g3W/RU++qOvFl+4DfI2QOleKNhq9xuqqZn0VQDMwfUESzN3XTOXX9x8oe85UYM/+2zV42miV6o3eO+/YeWDULzL7vu1xx8tyeGN5bbZ5qYP4gC4NsqOgJW4ERDUz9cvXmec7JU00SvV09VX+bazbBdJCrdCcBRGXLz14XqO5GRwxISy1wymKTCSkUc32+u8c9Sc/S1b+r94Vb2GJnqlehpjIHuFb8Hu8v2+cwVbwBhM/iY8g86iKiABd00eV4TvpXXITF66eyZBcU4yD4qA8ES7ff69cOdqGDKjMz+J6iSa6JXqafIz4Lnr4K0f2f2cVbYcfDameDcvvfwCUrqHBzIHk9kYy6zwAiIbDhM9ehbpKbG+laDihvt65IjAwM/MRq56CU30SvU0hzfZcp9tpmnYv4Zc91CeOhBDbXkBB7baxP9q62zyPAkk1ufY6+OdHjTeOWt0auE+QxO9Uj1NgV2z1TRU89SHWVTnbmFT0xBKTRQRUs9ZrixM9FBW/OwK5s6c7rvPm+hjnKab+JGovkEHTCnV0xTYF6nSUs9/lr3NHcFVDB47ncoKN5TBnOBMJPE8+keGQJJfrd07h3zalyAoHCbd2Pmxqy6hNXqlujNj4M0f2vVbARpqoCyL8vh0AC4P2QbA9HPm8PV59kWqu+Wor1kmxq8XjTvQlqHRcNZtEBjaKR9BdT1N9Ep1Z/kZkPEUvPw1u1+4FTD8q2oyAHck7rXH+0/w9aABX6JPHGenGZ6rSzj3ZZrolepO6qug+rBvv3inLZ11XBty1uJBeLXxLACkeBdEJkNYrG+hbvC1v4dGw7c3wLnf6YzoVTeliV6p7mTR9fCncdB01O5X5x07daSikD3r3ma3Zyj/c+cCCLOjXBkw0Zb94n3P0R41yo8meqW6i+Z6yP/Ebh/40JZVh46d3rDuQwa05BOTmkba4OhjE5QxwFkNyuX2Pcu/dq/6PE30SnUVjwdKM337Oat926V7obUZk7Oa+gHTANi4fhVJUsHAYePtNUfLbTlgku++656Cr73nGwilFJroleo67/4MHp9uZ5YEOLQeXIEQGgulmVRmrUfqSvnxoXNoMS6mubwTkTlL+s35KQw5B0bN8z1z4vU6jYH6DO1Hr1RX2fumLXcuhaTJULgd038cRbWG/Zu3sTGjie8HCGs8EyglmnMD9oLBt8zf9G/YL6VOQWv0SnWVhhpblu4FY6jN3cTivBi2VIXQX6o4172Lon6jWPTd+cQnDaGfqbPX+6/nqlQbnDLRi0iIiGwUkW0isktEHvQ7910RyXSO/87v+H0iku2cm3fiJyvVxxz40L5wBWiohsZqu12ayb/eXUe4p4adJpX64ASGh9QwIySfgePPY8KgKAKjBtprw+J9c8cr1UZtabppBC40xtSKSCCwRkSWAaHAVcAkY0yjiCQCiMg4YCEwHhgIrBCRUcaY1o75CEr1APtXwb+vhpn3wCUPQXk2AC1RQ/FU5vHRhyu4NQhuv+5KBtdsxrX6LXuft/YeMcCW3vZ5pU7DKWv0xqp1dgOdLwN8E3jEGNPoXFfiXHMV8KIxptEYcwDIBqajVF+210ncB9cdV/6zfAJBtHC2aw9GXAyfMIOgmGTffd5EH+4k+ii/c0q1UZva6EXELSJbgRJguTFmAzAKmC0iG0TkAxGZ5lw+CMjzuz3fOaZU3/HRH2DbYt++d4Srs0jIoc3vsd+TxG6PnYvmq4lZSNxIu6yft/YOvhGu3j+I43TGSXX62pTojTGtxpg0IBmYLiITsM0+McDZwH8BS0REgBN14DWfPiAid4pIhohklJaWnvEHUKrL5X0C6x/37Tc3wMqH4NU77aRkxvjWcq2vYNGKT4gu/YQNnjHce90sANwV2ZDk9IePSPI9KybVllNvgylf1qkM1Bk5rV43xpgqYDUwH1tTX+o07WwEPEC8c9x/KflkoOAEz3rSGJNujElPSEg4w/CV6mQ5H8ADUVCW7Tv29KW2T3yWLBj9AAAgAElEQVS+syCId2EQgNpiO7q1sQZGXQpA5qpFREo911xzE0nJKb5rE8fZMirZvnAdewW4nP9FowbBVY9DcETHfTbVa7Wl102CiEQ726HAxcBe4DXgQuf4KCAIKAPeABaKSLCIpAIjgY0dE75SnWzd/9jSu3zf0QrwOGu3FjkDn4q2+66vOUzT/o8AeCDL1s6vCbMLh4QOO8dXYweIdpb4C+oH/7UfbnquQz6C6nva0usmCXhWRNzYXwxLjDFvikgQ8E8R2Qk0AbcZYwywS0SWALuBFuDb2uNG9RoNTpfIsixbOouAAHCkyJZFO3zHagrY9f4i4jwJLG2awgMhkNayHdzBdu1Wl19dy5vowTd3vFLt4JSJ3hizHfjMqsHGmCbgyye552FAJ8BWvY/zMvXYrJKHNwNiF/GocVooC7fjGTAZV9E2cvbvI7L2AHtkGBt/dR08+iOkscau9uRN8uOvgV2v6kAo1WF0ZKxSJ9PcYNvj1zrNNUcroL7CbnuT+uHNeOJGUBw8lKbKwzQeKcNTspenC1NoNAGs/DiDoa4SZk6fQUig29c90j+pX/80fH/n8dMMK9WONNErdTKFW225/Be2dAY5NYUmYI4UgjG05m/iP2UD2F4dxv6cbB565GFcppmXm86myMQyOzibAFoJHzjGPiPS6WnsP/BJBKL9+y8o1b400SvlVXUIqvN9+7lrbOmyLZymzM4e+XbtSKgt4d/vfIj7aAnbPcMpNtH0lwrGykGqTRi//dbNDEkdwZhWe8+xxB4SaUttplGdSGevVArsEn5/ngihMfCTA7aWnb3CnvO0YBqqyd2dwUATSIZnNFe717F/zcsQCJdfejlTmrfCqpXcPLgCV+BkJg2O8dXewbfi07RvQGCYnU5YqU6iNXqlAEp227K+0i4GUl8JeRsgYSwAL733AWWZa9lpUrhl3kwAvpG4B+MKYMq02cdGs7oLtyCJ9h4inYnIgiN9y/4NPQeuekwnJlOdShO96ptK9sKaP9lRqwBVfrN2FO/EHPoYjIftSdcCkLnhHaZKFs2DZzJulG1vH1T5CRI7HAJDINJvNGuCN9E7NfrQaF3xSXUpbbpRfdPSb9iBTQMmwYiLoOrgsVOe4j1sK6hhkhHu2jiA9SFwfUgG7lbD2ZffcXyTjHcumoiBvmPexbq967bq+q2qi2mNXvVNFTm2PPSxLfe/T2PMaCoDE3ntgw2UZWWw3wwkZkAKre5QxrY6a7vGjbDt+OIsxH0s0ftNRDZoqi1Tz4NJN8ENz3T4x1Hq82iiV31PXRk0OTNvl2eRveFtOLSe35Wkc6AxkkQqmNGviNTx03n7++fhjnb6vkcMhKAw2wzjHewdP8qWoTHQfwLM/jEEBNtjYbFw7ZM6tbDqcproVe9XXwWvftM3RYEzZbARN8UHdrH0jVcBWNx6AUExycyMriSyoYDApAn2eu9LVf++797pCuJH21IEvrkWLvpFR38apU6bJnrV++18BbY9DysesPvbX8IER7Im9AKC6gpIDyumLiSJxd+by4QxY5AjzqjX/uNteWw0q98EZNf8HS57FAZ+ZnYQpbodfRmrej/vbJLOYKjqnE/YdHQYGz3xzA6s5YKIPCR2AuMHRkGOX+8Z77TBoTG29F/0Y+i59kupHkBr9Kr3WfkruxC3l3cisqqDeFpaCKnOYZ9JJm2CbZqRihxIcJpgIv16z3ibZyYvhLNuh/TbOyF4pdqfJnrVu5Tvh48ehWev8B1zeth4qvKZ+4unCJZmJk+ZzvyZ03zXeAc5+a/u5O37PmAiXPFnXfRD9Via6FXPVpFz/Pzv3sW3ATyt0FADNQXU9RuCCw/nu2wzzpSp0yHKbyIx7yCngWkwegHc+nonBK9U59BEr3q2/5kCT8wCj8fuV+b6zlXnw+EMwPCP6nQAvp1sa/chSeOO7/ue4HSTDI6Am5+HYXM6OnKlOo0metVzVR3ybRfv+Myxv7+2gj3vPUWdCebD4DkAxBavg/ABdq4ZlxumfR0u+Lk2y6heTRO96jkOfAj/uhqajtr9g+t957wjXYt30Rpre8cczN5FXNEa3vVM47s3zIOAEHuNt/YOsOAPcP5POiF4pbqOJnrVcyy6wS7Kve8du++/CHfVIagphJJd/Kk0nUYTSLork0Sp4qIL5zFnzACIHmqvjR/12Wcr1Ytpolc9hzvIlvkZAJjyLDwJ4zAh0ezevZMnn/4/ADYEnEVr9FCuCbUrREWlpNn7AkNtqYle9TGa6FX31NIIy34Klc6sko1HoLHGbjsrPR3J282yonB2HY2m6NA+BpatpcjEcPW8SwjrPwLxzmfjnabAu9iHd8SrUn2EjoxV3VPOB7Dhf+HQerjrA1/CByjPJmNPDun1eezynMuAqCAmevKJpQYZfSlfOjsFKp3pCkKiIDzRbp/7XRg5zzfjpFJ9hNboVfd0yHnRWld63P4noTNpqMjn8X8/D0DaORdx1qTJJDTk4m6owJU0yV7vXZM1OPL4RT8SRukiIKrP0USvuoeP/giZy3z73oVAaouhtZn6vSs45EngnZoUQqSZK9zraXGHcvHcK31TFYCvtj5khi2bj3ZO/Ep1Y5roVderyoOVD8ILC48/BuBp4Vt/fYnm/R+wxjORS2bYF6tXhu0iYFAarqDQTyV650XrgEkw9Ta4/p+d9CGU6r400avOV5YF+9/37ees9m031dmyOo/KUNsdcmzZu0RKPbPm3cDZk+yMkgGNlXa1Jzg+0Uc6UwqLwJX/oyNclUITvepsxsBj6fDva6DR6RWTv9F3vjTT9oc/UsiLRyYDcGOwXe5vSNpFx0865m2m8ZaRg8ClP9JKfZr+X6E6V81h37Z3ArL8DIhJsdulmTTut1MMrwueiQkKp39rkR3VGp54/ELb3vnhA4Lhx1lwx/KOj1+pHuiUiV5EQkRko4hsE5FdIvLgp87/WESMiMT7HbtPRLJFJFNE5nVE4KqHKNkD/7rq2KIfFGz1nas6CA01mJI9bI++GIBtu3bw/jtLqTGhjJh0LuLtGhmVbJtjgsN99/t3kwxPhKhBHfxhlOqZ2lKjbwQuNMZMBtKA+SJyNoCIDAYuAY7NJCUi44CFwHhgPvA3EXG3d+Cqh/jkH7YNfs2f7H7hVsDp3lidT1X2BgTDo5lxlJsIdu3Zw6j6beSETuKnC8ZDvwR7rX87fHCULWP8lvZTSp3UKRO9sZzGVAKdL+Ps/wn4id8+wFXAi8aYRmPMASAbmN5+IasewxjIXmG3q22TjSnYSmPsaDzRKVQX57Jo6VIAZp8/j4DoZC5PLGO4q5C02QsIDnD7lvHzLusHcPeH8F/7wa3j/ZRqizb9n+LUyDcBI4DHjTEbRORK4LAxZpscPwBlEPCx336+c0z1BfWVtsbtctmmGe/88NV5YAx1uRksa5jIYFcproo9THBFUB85jG/MOwvKhkDWu/b6JPsiliCnqcZ/EW5ve75Sqk3a9DLWGNNqjEkDkoHpIjIJ+Dlw/wkuP9GwQ/OZi0TuFJEMEckoLS09nZhVd1VTCH8YA698ze47Sb41fgxN5Qf58VPLCG+ppDpmPK7owQwPrOScfocJHWoXBSHSr0dNwhhbnv8TuOxRGH9N530OpXqZ0+p1Y4ypAlZjm2dSgW0ikov9BbBZRAZga/B+a7SRDBSc4FlPGmPSjTHpCQkJZxa96l52vwYtDbDrVQCay3MBeK5oMEEtRzA5HwAw9+L5TJ88ibjWEoLqCn0Lc0c4C3MHhvl61ySMhunfsIuEKKXOSFt63SSISLSzHQpcDGwxxiQaY1KMMSnY5D7VGFMEvAEsFJFgEUkFRgIbT/J41ZNlLYd97/n28z/xbTceISfjPRpNAJs8NpH/bvxBjLgZMna67UXj5U303hp9YJjOR6NUO2pLjT4JWCUi24FPgOXGmDdPdrExZhewBNgNvAN82xjT2h7Bqm6k4gAsuh6evwFqnD/YCrYcO/3Rho0MKVrOa62zWDjvPADcWcuQhDEQFHb8wtzeaQsGnWXLCdd2xidQqs845ctYY8x2YMoprkn51P7DwMNfKDLVvbQ02het3gW1/WvvJXvsoh4VOeQPnEdywbusfncpswObiBp5DueeNRVWOdcOdBYB8e/z7u0m2X88/Gifrt+qVDvTkbGqbd76EfxhNBytsPtFO3znyvbR7Ixm/VXuWAAWhO4EYP6c86BfvO/aJG+id5puAkIgIMh3PqK/rfErpdqNdkRWbbPl37bc+yZMvRWKd2IGTKSpJJvXV6whpLGcc1xRrPCcRUtwNFMbnRGw8Z+a/91bow+OgDs/OH4glFKqQ2iNXp1aQ7VvuzQTPB5aCrbzZnEch1piiWgsZowcYn/QGHb9agEB0U77e2gMhMXZ7YsfhClf8bXDg036YbGd9zmU6qM00avPKsuCJbfBkSK7X7rv2KnKgzv4cNnzBNSXsbxxPCUSx7ykekYGFHP2jJmEBLpt8wscX5uf9X246jHtJqlUF9BEr2yN3ePx7b9zn+0Tv/H/7H7pXgAqw0dQnJ/D1vXLaTEuYtOvZ9L4CbhKdiKeFt+i294Xtro2q1Ldgib6vq6pDn6bCv/5ru9YaaYtvfPEl2XS4gpiWdUQ+ksl50SUQtxwHrh2KhGJKb77vIneOxFZnCZ6pboDTfR9Te6a41d02v8+mFbY8pzdry2Famcy0pI9AFTuW09mSxKh8UOIkVqmBR0ioL/tXXNcN0nvik8zvgkX3Q9pt3TsZ1FKtYkm+t6svur4JfuaG+CZBXZ+eC//+eHrK+HgWgDMmCugrpSNG9cSU5bBW60zuHymM5yiOs83m6T/CFd3oC0j+sPsH9k54pVSXU4TfW+24pd2yb5DzmSi/n3fa0ts6dTaAag8SEP2BzRICN/fbnvOrHjddqucfv4CAhNH+65NdCYdG3IOjJoPlzzUUZ9CKfUFaaLvzbxJfN87tiza7juXtxFamuDgWloSbO3894vfI2/r+2xsGUmesbXxa6Nse/2cmbN9bfDgq9EHBMMti2HmPR36UZRSZ04TfW/l8RxL9J4iO0qVQx9DcCTGFciR/evhcAY0VPHn2ksAkLJ9DDOHGD7lfBb9+HoAxhzdbPvC94uDkEjf82OHderHUUqdOU30vVXZPmisAeBw1lYeeGMXJvcjygddQIEk8sGGT3j8uRcBeKFyDEdd4fxw6H7ceBg0ZjqhMQPB5Qyc9s4ND3DPdvjSy772eKVUt6eJvrco2AKPDIX8TXbfmXRsXfAskijnjXXbkSOFPJHZjwNN0QyUcsZ7siiQRJb86ErCEofhKths7x0wyQ5sinCmDfbOLgkQMxRGXtKJH0wp9UVpou8t1v8NGqpgzxt2P38jjYGRvF47lgDxcFGQbb4ZPCqN2IHDmBJVx5x+hxg4bhbDE8J9c84ER0L0ULvtdiYb86/RK6V6HE30PdXzC2HNn3379ZW2LNkNQEXmOtY3phKaaNvSfzrA1vQvm3cp48aMQ44UQk0+JDvL+HnXYe0/wa73CpDo9JWPH9GRn0Qp1cE00fcEO5cet6gH9VWwb5ntPtnSZI85i3+05mVwyf97g+ja/VTFTOYH19tmlriS9ZA8jfgBQ47v+z7ISfTe5hn/hbeveQIu/T2kzumYz6WU6hQ6TXF31NJouy2CXXD75dvt9k8O2NkeC7f5ri3PgvABULILIy7cDRWkNX+EK9BwydzL6Tcg1Xft1FttGek3mjVpki0n3wz1Fccvwh0cATPubP/Pp5TqVFqj7262vgC/ToSqPLvvfUEKcHCdLXe85DtWsofmjf8E4M9NNkl/f4AdGNVv2HRw+/0uHzrTlt6ukTEpdmUosIt/zPrB8TV6pVSvoIm+u1n7F1tmL7elf5NNxX470+TOV2gaczUAj7y4nIz3l7LLM5RXPHZt1kHl6yFqiJ0PHuCrb8O0r/sSfGwqfPUt201SKdXraaLvbhqqbJm1wpYFW/AkjqclJJYVa9bzpz8/As1H+Wnh+dSYMFKliHR3Fms8E5gwZizG2/d9wETfM1NmwoI/HL/SU8osnUZYqT5CE31XaqqDl26H7JV2v7YUjhQC0Jy7jt2Hq2k6lMGSggS2HY2jX91BhtVtpdDEsrQ4kYbQ/twQsZNAmrnlxlv465enIf2cicQGTOiiD6WU6m400XemD34P2/3a13e9CruWwuKv2P0i+5I1NzKdwMZKvvPYEoKaqthhhpFrBnBWeAVXJhQRNOQsYvsFE5E4FFd9GQARQyYT6HZBnTNZmX+NXinVp2mi70itLb7tI8Ww6tew9Ou+Y7lrnesaoaUJT4FN9IvKbVfHC122ff6mq67k3GnTCTpahFTkEDfqXDb/4hJC45xuksFRvi6TF/zMts+P0NGrSilLE31Heffn8NepvoFMzpQEgG2yAV83SU8LDUWZ5K/8O1s9w6gMs10ifzY8F+MOYtLUc0lKHee7f9BUW3q7ScaP9LW/z/4RfH87BIZ00AdTSvU0mug7QksTrH8Mqg5CpjNFcJlvgW0qD0J1PpTsonTgBQA8+L9PM0SKWRV0Ib+9fT4ArkNrkf7jbZ/6GL/+8AOdBUC8a7P6r/IEx790VUr1eZroO0LFft+2N8H79YcvPbSXd156EoBfl18IwDVhdq74H9x6A27/KQe8ST1uuC2jh0JIlN0ecTFMvgXmP9L+n0Ep1Wtooj+VT56CTc/69lub4e3/sgt3eB1cD8vvh8Yjdt+ZbwaAsn1kb1wGe/7D/kF2Cb9//mcVCYfeZrdnKMuODMMjbqZ7nEVBEkZDUJjvfm+iD42G722FO97znYseAtf8L0QObMcPrJTqbU6Z6EUkREQ2isg2EdklIg86x38vIntFZLuIvCoi0X733Cci2SKSKSLzOvIDdChj4K0fwn++Z7cBdr8OG5+Ef1/ru27FL+1Apw9+a/cPfIQJCicnaga5+/ey4o1FeIxw9f4rqDFhTHbt5yxXFiGTrmbjf8/FFZFkX8iGRPsGOQ1wpibwJnqwA528zTVKKdVGbanRNwIXGmMmA2nAfBE5G1gOTDDGTAL2AfcBiMg4YCEwHpgP/E1E3B0RfLurOgTbl/j2K3I+u+0snk1LPXha7S+A0r0AmB2v0NTioSFnLe/Xj2RdeTiRTUWMkHzyA1O47txxtEQNYV6AbcYZNmUO0WFBvjb2WL92+JtfhIvuh0S/5fuUUuoMnDLRG6vW2Q10vowx5j1jjLf/4MeAd0rEq4AXjTGNxpgDQDYwvZ3j/uIaauA/34e6ct+xF26Gpd+wTTEAZVm+c1UHbXnYWdjD02JnjCzcBg3VNEUORY4UcO5/v0hgRRY7TCrJKaOIlVrODs1jyOg0HrhyPLHJo5GWevuMBGcaYG/vmf5+g5yiBtkeNC5tXVNKfTFtyiIi4haRrUAJsNwYs+FTl3wNWOZsDwLy/M7lO8e61upHIHeNb3/HEtj0NHz4O7vf0gjFztqqWe/a0v+latUhO5Nk4TY8Q2cD8OeXl1O36UU87mD+3mD7rZ/n2oZbDAvmXcqcabbZJbypDOKcF6zeScOCI33NMN429kFntecnVkopoI2J3hjTaoxJw9bap4vIsaqniPwcaAEWeQ+d6BGfPiAid4pIhohklJaWnn7kp6OuHFb/Bp5Z4DtWnW/LylxbFu04dqo1L4OWVg8U7aQlKBKPuGmtPMS+d/6GB+GO/bMAKDiwm92frGRL8xA+qY0D4L+GHwJg5ORZx8/7/ulEHxLt6wZ5wc/hpucg7Zb2/NRKKQWcZq8bY0wVsBrb9o6I3AZcDnzJGO/bSvKBwX63JQMFJ3jWk8aYdGNMekJCwhmEfhq8zS3g6xmT5wxg8nZ/zM8AoCQ2nfwDmTz0agbN+5bzXsM4CjwxvPHBxxTveJ+dnhR2BE6mVdxckVTDJNcBNntGUmrsu+ikvLdtF8iIpBMneu9L1vFX+c4FhcHYK3xz0CulVDtqS6+bBG+PGhEJBS4G9orIfOBe4EpjzFG/W94AFopIsIikAiOBjZ9+boc6vAl2v+HbL93j2y7Pts00zkhVU3GAJ1fu4Ej2Wird8bxekkiiVFG3+RUC60tZ1HIBzeHJTI44QlrQYUZMnEHGLy/FHT2E2Q2rCaaZb3zlNpbce4Pv3xh3pa2t+3d79E4RnJwOP9gNl/yqA78BSinl05YafRKwSkS2A59g2+jfBB4DIoDlIrJVRJ4AMMbsApYAu4F3gG8bY1o7JHr7D/qmFPDu/+tqWPIVKHLa3MuzfecrcuwvgtZGtkZegGB4fcWHuLLeZVnjREqIIVSauCFmH0dMKJdftZDUEWMZ1rCHiNZKwpIn2+dEJUNtMYgLhp5DRHQCx1qt4pzpf92BEGabdAiL9cUQNUhHryqlOs0plxI0xmwHppzg+ElXjDbGPAw8/MVCa6MlX4F978Ed79o+5xU50Fhjz5XutdP1Fm6jtN9IEuqyWLt5O+lJewhEeKxsKv8IWsVXIrfQr7GRqMlX8t2hIfD2Is5u3gip6dw8IwVWDQFPs32md/pfb209abJvpGrcCLu0X7Rfy9V3MuBoRad8K5RS6kR6ft+9nA/tYKMdzmpJBz48dur1D9bzyKK3oXAbT9VMp9aEsHdfJns+WspWz3DKI8YAsDDkYwAWLLiayEH2GE21NokDRPklbm8XyH7OewXv8nwAtyyGcVdDsl9v0rBYiD/p70SllOpwPTvRN9VBYzUAu3ZsZsvBCsyaP2IiB3HUHUFd0X6adr8NwDut0wiOHcTcuFImu3LoN34+r957HbgCoDrPJvCwWEgc63u+98Vp9BDfMW8TzNgrYfhFMN1v8ey44XDjsxAc3pGfWimlTkvPTvTeLpJARE0W333idaTqEL8sv5h9zYkMllJuid5FSUgq37v+EgKjkxlcnYFgGH32AnC5ISzePiDFdpk8tlg2+BbvGJhmy3F+PWWGzICvLIWYoR34AZVS6os7ZRt99ybsizmfXWUtXONey6Uu27lno2cs01yZXBa+D3d9DZzzHa6dmgwH/XrBeJtgaotsmXqe79z4a+3KT941VUOi4L7DtvavlFI9TI+u0ddGDmNu4V0833IRAN8LX0lLUCS/uvMGLp19Du76CjtVwci59gbvC9SowRASabfn/Mw20fivyHTN3+HH2bbXjFdwuC7moZTqkXp0ot9TaHvXjBpn29IjGosJGDiZaanxBCT7dRTyNsF4E71/zXzOvXD3RxAQ5DsWEAThHTyISymlOkmPTvTTUmLZ+eA8Hrz5AnA7iTphtC2HnOO70Ft7H3O5c+MdnRekUkp1sR7f6Bwe7P0IzgAkb6+Z8EQYOe/4RTwiBsDPi3WqAaVUn9LjE/0x0UPsYCX/mvyXlnz2Om1nV0r1Mb0n0S98Hna85JvjXSmlFNCbEn3CKLjw510dhVJKdTs9+mWsUkqpU9NEr5RSvZwmeqWU6uU00SulVC+niV4ppXo5TfRKKdXLaaJXSqleThO9Ukr1cmKM6eoYEJFS4GA7PW4IcKidnnUiUUB1Bz6/I+PvybFDx8bfk2MH/bn5PL3552aoMeaUU+12i0TfnkSktC0f/As8/0ljzJ2nvvKMn99h8ffk2J3nd1j8PTl25/n6c3Py5/f5n5ve2HRT1cHP/08HP78j4+/JsUPHxt+TYwf9ufk8ff7npjcm+o78ExNjTEf/0HdY/D05dujw+Hty7KA/NyelPze9M9E/2dUBfEE9OX6Nvev05Pg19g7W69rolVJKHa831uiVUkr56RGJXkT+KSIlIrLT79hkEVkvIjtE5D8iEukcTxGRehHZ6nw94XfPTSKyXUR2icjvulvszrlJzrldzvmQnhC7iHzJ73u+VUQ8IpLWVbGfQfyBIvKsc3yPiNznd093/94HicjTzvFtIjKni2MfLCKrnO/jLhG5xzkeKyLLRSTLKWP87rlPRLJFJFNE5nVV/Kcbu4jEOdfXishjn3pWl/zcn5Axptt/AecBU4Gdfsc+Ac53tr8G/MrZTvG/zu/6OGx/1wRn/1ngom4WewCwHZjsF7O7J8T+qfsmAjld+X0/g+/9LcCLznYYkOv8LHX77z3wbeBpZzsR2IStxHVV7EnAVGc7AtgHjAN+B/zUOf5T4LfO9jhgGxAMpAL7u+rn/gxi7wfMAu4GHvN7Tpf93J/oq0fU6I0xHwIVnzo8GvjQ2V4OXHeKxwwD9hljSp39FW245ws7zdjnAtuNMduce8uNMa30jNj93Qy84Gx3Sexw2vEboJ+IBAChQBNQQ8/43o8DVjr3lWC7/KXTdbEXGmM2O9tHgD3AIOAqbMLDKa92tq/C/pJtNMYcALKB6V0R/+nGboypM8asARo+9agu+7k/kR6R6E9iJ3Cls30DMNjvXKqIbBGRD0RktnMsGxjjNO0EYP9D+d/TmU4W+yjAiMi7IrJZRH7iHO8Jsfu7CV+i706xw8njfxmoAwqxNbFHjTEVdK/4Txb7NuAqEQkQkVTgLOdcl8cuIinAFGAD0N8YUwg2oWL/+gCbSPP8bst3jnVp/G2M/WS6/Hvvrycn+q8B3xaRTdg/sZqc44XAEGPMFOCHwPMiEmmMqQS+CSwGPsL+ad7S6VFbJ4s9APtn4Jec8hoRuaiHxA6AiMwAjhpjdgJ0s9jh5PFPB1qBgdjmgx+JyLBuFv/JYv8nNjlmAH8G1gEtXR27iIQDrwDfN8bUfN6lJzhmujL+04j9hLr6e/9pPXZxcGPMXmxTByIyCljgHG8EGp3tTSKyH1tTzjB28MF/nHvuxP6P3elOFjv2f9YPjDFlzrm3se20K3tA7F4L8dXmvfd0i9idWE4W/y3AO8aYZqBERNZimz9yukv8n/Mz3wL8wHudiKwDsv5/O/fvSlEYx3H8/SUpLBSlTMpiMBsUKQNlkjIog81iMvEHyGA2MIpSik2xyyCJyzVIEjHbDI/h+cp1csUdnB99XnW65z7nPqdP5z4998qnlXIAAAFrSURBVJzzPOf6tlSym1kDsaPcDCHsevGzmXWGEJ7MrBN48fIHvp7tdgGPaeX/Y/aqstJuIMdn9GbW4a91wBKw5u/bzaze17uBHuA2UacVmAPW/z959ezAAdBnZk1+uTcIlBJ1spr9o2wS2K5SJ9XsiSzJ/PfAsEXNQD9wnaiTyWPv7aXZ10eIZ/OptRszM2ADuAohrFZs2gdmfH0G2KsonzKzRr/11AOcpJG/huw/7SsT7QbIzaybLeItmTfir/8sME8cEb8Blvl8+GsCuCTetzwFxhP7KfkylbXs/vlpz38BrOQs+xBwXGU//5q9hnbTAuz4sS8BC3k59sTZQWXiwOEh8R8N08w+QBzcPgfOfBkjzkQ5Il5tHAFtFXUWibNtysBoWvlrzH5HHDh/9e+qN812/92iJ2NFRAout7duRETkd9TRi4gUnDp6EZGCU0cvIlJw6uhFRApOHb2ISMGpoxcRKTh19CIiBfcOJ8BLQ+bJ2FEAAAAASUVORK5CYII=\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "indexed_data.plot(y = [\"SAFitFilled\",\"Filled\"])" + ] + }, + { + "cell_type": "code", + "execution_count": 179, + "metadata": { + "scrolled": true + }, + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 179, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXoAAAEHCAYAAACgHI2PAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzt3Xd4VMX6wPHvpCcQQktoISSEXkNHmjSxISII2LH3a70Wrj+vetVrV6zXgoIiIgoiCBYUQZqUhE4ChBJIqIFACOllfn9MQk3IJtnds7t5P8+zzya7Z+e8m5O8ezJn5h2ltUYIIYTn8rI6ACGEEI4liV4IITycJHohhPBwkuiFEMLDSaIXQggPJ4leCCE8nCR6IYTwcJLohRDCw0miF0IID+djdQAA9evX15GRkVaHIYQQbiUuLu6I1jq0vO1cItFHRkYSGxtrdRhCCOFWlFJ7bNlOum6EEMLDSaIXQggPJ4leCCE8nEv00ZcmPz+flJQUcnJyrA7FIwUEBBAeHo6vr6/VoQghHMxlE31KSgrBwcFERkailLI6HI+itebo0aOkpKQQFRVldThCCAdz2a6bnJwc6tWrJ0neAZRS1KtXT/5bEqKacNlED0iSdyD52QpRQfk55uaGXLbrRgghXEZeJnxyMRzbDaFtoFHn07cGHcC/ptURXpBLn9G7gpdffpn27dvTqVMnYmJiWLVqFQCpqan4+vryySefnLV9ZGQkHTt2JCYmhpiYGFasWMH+/fu59tprAVi/fj0///zzqe2nTJlCaGjoqe1vueUWAP7973/zxx9/ADBw4MBTE8oiIyM5cuSIzfFPmTKFBx98sPI/ACEE/PECHE2EbrdCcEPY/hv88iR8cSm8Eg4f9IBZd8KK92H3EshJtzris8gZ/QX8/fffzJs3j7Vr1+Lv78+RI0fIy8sD4Pvvv6d3795Mnz6de+6556zXLVq0iPr165/12MyZMwGT6GNjY7niiitOPTdu3Dg++OCDs7b/z3/+44i3JISoqN1LYPUn0PMeuOJ185jWkHEADmw4fduzAjZ9b573CYAHY6F2U+viPoNbJPoXftpC/P4Tdm2zXeNaPHdV+wtuc+DAAerXr4+/vz/AWcl7+vTpvPXWW9xwww3s27ePJk2alNlOUlISw4cPZ+3atfz73/8mOzubZcuWMWHChDJfc+uttzJ8+PBT/wmU5uuvv+a9994jLy+PXr168dFHH+Ht7c3kyZN55ZVXaNSoEa1atToVvxCignIzYM4DULc5DH3u9ONKQa3G5tb68tOPn0yFXYvhhzth55/QbbzTQy6NdN1cwLBhw0hOTqZVq1bcf//9/PXXXwAkJydz8OBBevbsydixY5kxY8ZZrxs0aBAxMTH06tXrrMf9/Pz4z3/+w7hx41i/fj3jxo0DYMaMGae6biZPnmxTbAkJCcyYMYPly5ezfv16vL29mTZtGgcOHOC5555j+fLl/P7778THx9vhJyFENbXgWTieDCP/B341yt++Zih0vBZqNoCkpY6Pz0ZucUZf3pm3o9SsWZO4uDiWLl3KokWLGDduHK+++iqpqamMHTsWgOuuu4477riDxx577NTrSuu6uZDSum7Ks3DhQuLi4ujRowcA2dnZhIWFsWrVKgYOHEhoaOiptrdv316htoUQwI6FEDcZLnoQInrb/jqlILIf7F5qunhcYISbWyR6K3l7ezNw4EAGDhxIx44d+fLLL9m3bx+HDh1i2rRpAOzfv5/ExERatmzptLi01owfP55XXnnlrMd//PFHGTopRFVlH4e5/4D6rWDw/1X89ZH9YfMsOLoD6jsvL5RFum4uYNu2bSQmJp76fv369RQUFJCZmcm+fftISkoiKSmJCRMm8O2339rUZnBwMBkZGVWObciQIcycOZPDhw8DkJaWxp49e+jVqxeLFy/m6NGj5Ofn8/3331d5X0JUO789Yy62jvwYfAMr/vqoAeZ+9xL7xlVJkugv4OTJk4wfP5527drRqVMn4uPjiY6O5pprrjlru9GjRzN9+nSb2hw0aBDx8fHExMSc17dfEe3ateOll15i2LBhdOrUiUsuuYQDBw7QqFEjnn/+eS666CKGDh1K165dK70PIaqlbb/C+q+h36MQ3q1ybdRtDsGNXaafXmmtrY6B7t2763MXHklISKBt27YWRVQ9yM9YiHNkpcFHvSGoPty9CHyqMGLth7vNyJt/Jjqsn14pFae17l7ednJGL4QQJX55ErKOwjX/q1qSB9NPn5kKqVvtE1sVSKIXQgiA+LlmwtOAJ0xpg6qK6m/ud1vffSOJXgghMo/AvEdNgu//uH3arBMJIRGQZP0FWRleKYSo3rSG+Y9B7gkY+RN423Exnqj+sO1nKCoCL+vOq23es1LKWym1Tik1r/j7N5RSW5VSG5VSs5VStc/YdoJSaodSaptS6lJHBC6EEHax5QeInwMDJ0CDdvZtO7I/ZB+Dw1vs224FVeQj5mEg4Yzvfwc6aK07AduBCQBKqXbAdUB74DLgI6WUt33CFUIIOzp5GOY/Dk26Q5+H7N++i/TT25TolVLhwJXApJLHtNYLtNYFxd+uBMKLv74a+FZrnau13g3sAHraL2Tn8fb2PlWDJiYmhqSkJGJjY3noIfMLcWYJ4Oeff54333yzQu3XrOnaNayF8Hjrppoz7qs/BG8H9GSHhEOdKMvH09v6ziYCTwLBZTx/O1Ay+6cJJvGXSCl+7CxKqbuBuwEiIiJsDMO5AgMDWb9+/VmPRUZG0r17ucNWhRDuIH6uOZsPa+O4fUT1hy1zoKgQvKzp3Cg30SulhgOHtdZxSqmBpTz/DFAATCt5qJRmzpuVpbX+FPgUzISpCwbxy9NwcFN5oVZMw45w+asVftnixYt58803mTdvXpnb7Ny5kwceeIDU1FSCgoL47LPPaNOmDbt37+aGG26goKCAyy67rCrRCyGq6lgSHFgPlzh47YfIAbD2Kzi4ERp3cey+ymBL101fYIRSKgn4FhislPoaQCk1HhgO3KhPT7FNAc6sth8O7LdbxE6UnZ19qtvm3LIHF3L33Xfz/vvvExcXx5tvvsn9998PwMMPP8x9993HmjVraNiwoaPCFkLYIn6uuW87wrH7cYF++nLP6LXWEzh9oXUg8E+t9U1KqcuAp4CLtdZZZ7xkLvCNUuptoDHQElhdpSgrceZtD6V13ZTn5MmTrFixgjFjxpx6LDc3F4Dly5cza9YsAG6++Waeeuop+wUrhKiYhLlm3HzdKMfuJ7gh1Gtp+un7OuCCrw2qcvXhA8Af+L24LO5KrfW9WustSqnvgHhMl84DWuvCqofqHoqKiqhdu3aZHxBSQlgIF5CeAilrYPCzztlfVH/Y+D0UFjjmom85KjSCX2u9WGs9vPjrFlrrplrrmOLbvWds97LWOlpr3Vpr/Yu9g3ZltWrVIioq6lR5YK01GzZsAKBv376nyhmX1LIXQlgg4Sdz326kc/YX2R/yMsw1AQtICQQHmDZtGp9//jmdO3emffv2zJkzB4B3332XDz/8kB49epCe7lqrxAtRrcTPhbD2UL+Fc/YXWdJPb005BClTXI3Jz1hUSxmH4K3WMPBpc3OWD3tDrUZw82y7NSllioUQojRbfwI0tLvaufuN6g97V0JBnnP3iyR6IUR1Ez/HrAUb6sBJUqWJ7A/5WbB/rXP3i4snelfoVvJU8rMV1VLmEUhaZsbOO3sEXGQ/QFkynt5lE31AQABHjx6VhOQAWmuOHj1KQECA1aEI4Vxb54Mucn63DUBQXWjQwZL69C5bjz48PJyUlBRSU1OtDsUjBQQEEB4eXv6GQniS+DlmQZCGHa3Zf1R/iP0CCnKrvlRhBbhsovf19SUqysEz1oQQ1Uf2Mdj9F1z0gPO7bUpE9oeVH5nJWpH9nLZbl+26EUIIu9r2CxQVWNNtU6JZH1BeTu+nl0QvhKge4udASFNo3NW6GAJrQ8NOTq9PL4leCOH5ck7Azj+tGW1zrqj+pusmP9tpu5REL4TwfNt/g8I8a7ttSkQOMLEkr3LaLiXRCyE8X/yPENwIwntYHQk0uwiUt1P76SXRCyHcQ14WbPyu4iUEck/Cjj+g7VXg5QIpzz/YrDTlxH56F3jXQriZtV/B0rch+7jVkVQvf70KP9xlboUFtr9ux+9QkOMa3TYlovrDvjjzIeQEkuiFqIhdi2HuQ7DwBZjYEf58CbLSrI7K8504AKs+hbrRphtm7j+gqMi218bPgRqhEHGRY2OsiMj+Zqhn8kqn7E4SvRC2yjwCP9wD9VvC7Qug+UBY8ga80wEWPAsnD1sdoeda+iYU5cNNs2Dgv2DDN/DLk1BeiZT8bNi+ANoMBy9v58Rqi4je4OXrtH56l50ZK4RL0RrmPAjZaXDj99CoE0T0gsMJsPQt+PsDWP0ZdLvVrAtaq7HVEXuOtN0QNwW6jjfru178pFmtacX74F8Thj5f9mt3LIT8TNfqtgHwqwFNujmtn17O6IWwxZpJsP0XGPqCSfIlwtrC6EnwwBroMApWfwrvdoZ5j8LxvdbF60n+eg28fGDAE+Z7peCSF6H77bDsHVjyZtmvjZ8DgXWcWm7AZlH9Yf96M8bfwSTRC1GeQ1vgt2egxSXQ+77St6nfAkZ+BA+thZgbYe1UeK+L6c8vzHduvJ7kcAJs+BZ63m1WZyqhFFzxFnQaB3++CCs/Pv+1Bbmw/VdocyV4+zovZltF9gddCHv/dviuJNELcSH52TDzDggIgZH/K39WZZ1IuGoiPLyhOOF/aZKNqJw/XwK/mtDv0fOf8/KCqz8y/e+/PmU+XM+0azHknnDeAuAV1bQnePs5ZR1ZSfRCXMhvz0BqAlzzMdQMtf11IU3gijfAJxCSljsuPk+2Lw62zoM+/zC13Evj7QPXfgEthpqROJtnnX4ufi74h0DUxc6Jt6J8A81/G75BDt+VXIwVoiwJ8yD2c7joQWgxpOKv9/E3Z21Jy+wfW3Ww8EUIqgcX3X/h7Xz8YexUmHYt/HC3SZwthpoPidaXg4+fc+KtjDFTnLIbOaMXojTp+2Dug9CoMwx5rvLtRPaHQ5tlrH1F7V4CuxZBv8fMTNLy+AXB9d+aypDfjYdF/4Wc49BuhONjdQOS6IU4V1EhzL7HTLUf/UXVzggj+wEa9qywW3geT2tzNh/cGHrcYfvrAmqZcfb1omHZ26ZvP3qw4+J0I5LohTjXsnfM+OYrXjejaaqiSdfifnrpvrHZ9t8gZbUZL+8bWLHXBtWFm3+EsPbQ+fqKv95DSR+9EGdKiTX/9rcfZUbNVJX001dMUZEZLlknCrrcVLk2ghvAfcvLnzVbjcgZvRAlck7AzNuhVhMY/o79FqiQfnrbbfnB/KwGPVO1se9KuUalShchPwkhSsx/HNJTzEzXwNr2a1f66W1TmA+LXjbdLh1GWx2NR5FELwSYwlebvoOBT5saNvYk/fS2Wf8NpO2Cwf9X6tl4Tn4hG1OkNHRlSB+9EFrDktehdkTpMzCrysfffHhIoi9bfo6padOkuxn7foac/EKmr97LR4t3kpqRyzd39aJPdH2LAnVPNp/RK6W8lVLrlFLzir8fo5TaopQqUkp1P2M7X6XUl0qpTUqpBKXUBEcELoTdJC0zizX3echxNVEi+0k//YXEfgEn9sGQf5+6NpJbUMjUv5MY+MZiXvgpnqj6NQgJ9OWbVa5fLK6oyLUuBFek6+ZhIOGM7zcDo4BzCzWMAfy11h2BbsA9SqnIKsQohGMtextqhFV+lIctIvsj/fRlyM0wpZ6jLobmF5NXUMQ3q/Yy6I3FPDtnC+F1Avnmzl7MuLs313RpwoIth0jLrOBygk70w9oU2jz7K499t57thzKsDgewMdErpcKBK4FJJY9prRO01ttK2VwDNZRSPkAgkAc4vg6nEJWxfx3s/NNMs3fkmOvG0k9fppUfQ9YRCgb9H9+tSWbwW4v51+xNNAgJYOodPfn+3ovo06I+Simu7xlBXmERP6xNsTrqUu09msWzP26mUe0Aftl0kGHvLOHOL9cQt8fa/+Rs7aOfCDwJ2DAXmZnA1cABIAh4VGst/68K17T0bVP4qnsFZmBWho9fcT+98xaEdgtFheg1n3EwbADjvs1kb9pGOoeH8OLIDgxsFYo6Z4hr64bBdImozfTVe7mjX9R5z1upoLCIR79bj5eX4pu7ehPk682Xfyfx5YokRv/vb3pG1uXegc0Z1DrM6XGXe0avlBoOHNZax9nYZk+gEGgMRAGPK6Wal9Lu3UqpWKVUbGpqakViFsI+UrdDwk/Q8y4zfd7RpJ/+PNm7VqBOHuK/KZ0IDvDh8/Hd+fGBvhdMhtf3iGBnaiaxe445OdoL+9/incTtOcZLIzvQpHYgdWr48cjQVix/ejDPXdWOlGNZ3D4llsvfXcqP6/ZRUGjjmrd2YEvXTV9ghFIqCfgWGKyU+voC298A/Kq1ztdaHwaWA93P3Uhr/anWurvWuntoaAXKvwphL8sngk9A2YuJ2Ftkf3O/R8oWgzkDXjrnc3K1L5eMvIV5/+jHkLYNyj3bHd65ETX9fZi+2nUuym5IPs7EhYlc1bkxV8c0Oeu5ID8fbusbxV9PDuKtMZ0pLNI8MmM9A99czFd/J5GdV+jw+MpN9FrrCVrrcK11JHAd8KfW+kJXrfZiPgyUUqoG0BvYapdohbCX43th4wzoNh5qOGmonvTTn6K15pkfNtH+xBJSG/RlRM/WNndnBPn5MCKmMT9vOkB6tvWrd2XlFfDojPWEBfvz0tUdytzO19uL0d3C+e2RAUy6pTthwf78e84W7p9ma2dJ5VV6HL1S6hrgfSAUmK+UWq+1vhT4EJiMGZWjgMla6432CFYIu1nxvrnv8w/n7fNUP70k+vcW7mDr2sU08T8Kfa6r8Ouv7xHBN6v2Mmf9Pm65KNL+AVbAy/MT2H00k2l39iIkqPzhuV5eiqHtGjCkbRhrko7h5YTu+goleq31YmBx8dezgdmlbHMSM8RSCNd0MhXWfgWdr4OQcOfuO7KfWR4vK63sVZM83Hexybzzx3a+aJKAPuaDan1ZhdvoGB5C+8a1mL46mZt7N7PsouzChENMW7WXuwc0r/AkLqUUPaOc8zsgJRBE9bPyI7NwdN9HnL/vat5Pv3jbYSb8sIn+LeoxqGglKupiCKxTqbau6xlBwoETbExJt3OUtknNyOXJmRtp0zCYx4e1siQGW0miF9VLTjqsmQTtrob6LZ2//2rcT795Xzr3T1tL6wbBfDIsAHVsd5VWgLo6pjGBvt58u8b5F2W11jw9ayMZuQW8e10X/H28nR5DRUiiF9XLmkmQewL6P2bN/qtpP31yWha3Tl5DnSA/ptzWg6Ad80F5QZvhlW6zVoAvV3ZqxNz1+8nMLbBjtOX7ZvVeFm49zNOXtaF1Q1umF1lLEr2oPvKy4O+PzMLRjTpbF0c1G09/LDOP8ZNXk19YxJe39yCsVgDEz4Fmfas84un6nk3JzCvkpw377RRt+XamnuTFefH0a1GfW/tEOm2/VSGJXlQf676GrCNmwWkrVaN++pz8Qu78KpaUY9lMGt+dFmHBkLoNjmyDtlVfuLtrRB1ahtVk+ppkO0RbvvzCIh6dsR5/H2/eHNMZL2cMmbEDSfSieijMhxXvQdPe0KyPtbE07gq+QR7ffVNYpHnk2/Ws3XuMieNi6BFZPMIkfq65b3tVlfehlOK6nhFsSD5OwgHHl9R6f2EiG1PSeWVURxqGBDh8f/YiiV5UD5u+h/Rk0zdvdX0UHz9o6tn99FprXpwXz69bDvJ/V7bjio6NTj+ZMMe8/1qNym6gAkZ1aYKftxffOnimbNyeND5YtIPRXcPPfj9uQBYeEZ6vqAiWvQMNOkDLYVZHY0T2M4tge+B4+g3Jx3njt20s23GEO/tFcUe/qNNPpu2Gg5tg2Mt221+dGn5c1qEhs9ftY8IVbQnwtc8ImJO5BaxJSmPlzqP8vesom/el07h2IM+PaGeX9p1JEr3wfFvnwZHtcO0X1p/Nlzizn94OXRiuIPFQBm8u2MZvWw5Rt4Yfzw5vx23nXqxMsF+3zZmu69mUuRv28/OmA4zqWrlJcFl5BcQmHePvXUdZuesoG1PSKSzS+HorYprW5sFBLRjTvSnBAQ5anMaBJNELz6a1WdSibnNoN9LqaE5r3OV0P72bJ/rktCwm/pHI7HUpBPn58OjQVtzRP4qa/qWkl/i50CgG6jSzawwXNa9HZL0gvl2dXKFEn5qRy9SVe1ix4wgbUo6TX6jx8VJ0Cg/h3oubc1Hz+nRrVodAP9ceJ18eSfTCs+38Ew6sh6veAy8X+mN1oX56rTXxB07QsFYAdWv42VxOIDUjlw8X7WDaqj0opbijXxT3DWxB3Rp+pb8gfR/sizXLBdqZUopxPSJ47det7Dh8khZhNS+4vdaaWWv38eK8eDJy8unYJITb+0VxUfN69IisS43SPqTcmGe9G2f7+0OzDNrAp62ORJSmMN/UlQluZOrauBoX6af/fNluXppvVgmtHeRLdGhNmtevQXRYTaJDaxIdWoOIukH4eJuxG+nZ+Xy6ZCdfLEsir7CIsd3DeWhISxqFlLNCV8JP5r7t1Q55H9d2C+etBduYsWYvz1xZdj96cloW/5q9iaWJR+jerA6vju5U7geDu5NEX1nH98Lvz0FRvpn4EdXf6ojEuf58EfavhWsng4+/1dGczwX66XcczuD137YxoFUoA1uFsjP1JDtTT7J4eyrfx51ers/XW9GsXg0i69VgTVIa6dn5XNW5MY9d0oqo+jVs21nCXAhrB/VbOOS9hAb7M7RtA2at3cc/L219XlmCwiLN5OW7eWvBdrwUvHh1e27s1cxtxsJXhST6ylo20dzXagK/PAn3LAFv97tI47ESf4fl70L326HDKKujKZ3F/fQFhUU8/t0Gavh589aYzoQGn/1hmJ6dz67Uk+xMzTQfAIfNh0CPyLo8eklL2jcOsX1nGYfMwugXP2Xnd3G263o25dctB/k9/hDDOzU+9XjCgRM8PWsjG1LSGdwmjJdGdqBxbQeuEexiJNFXRvo+WDcVut4M0UNgxo2w+jOzwLSw3on9MPseM5zy0v9aHU3ZLO6n/2TJLjakpPP+9V3OS/IAIYG+dImoQ5eIylWXPMvWeYA2xeQcqH/LUJrUDuTb1ckM79SYnPxCPly0g/8t3klIoC/vXd+Fqzo1cqm1Zp1BJkxVxvKJoIug36PQ5kpTO2XxK+asRVirsABm3Qn5OTBmCvi6+FmbRXVvth48wcQ/tnNlx0Zc1blx+S+oqoS5UK8FhLV16G68vRRjuzdl2Y4j/LhuH1e+t5T3/9zBiM6N+eOxixnRuXG1S/Igib7iThyAuC8h5gaoHWHGZV/2GuRnwx/PWR2d+Os10+c9/B1ryhBXlAV1b/KLu2xCAn15cWTZS9/ZTVYa7F5qats4IcmO7RGOl4JHZqwnJ7+IKbf14O1xMdQpazRQNSBdNxW14j0oKji7MFb9FmZJumVvQ7dbIaK3ZeFVa7sWw5I3IOYm6DzO6mhsY0E//Qd/7mDL/hN8fFO3sodC2tO2n0EXVqn2fEU0CgnkoSEtyc4v5KHBLT1uqGRlyBl9RWQcgtgvzFC9ulFnPzfgn+bC7Px/mu4D4VwZh2DWXVC/FVzxutXR2M7J/fSb96Xz4aIdjIxpzGUdGjpln8TPNf/9Nopxzv6AR4a2YsLlbSXJF5NEXxEr3oPCPOj/+PnP+dWAS1+GQ5sgbrLzY6vOigrhh7vMnIYxU8yxcCdO6qfPLSjk8e82ULeGHy+McEKXDUDOCdi1yGndNqJ0kuhtdTLVnM13HAv1okvfpt1IiLrYjN/OPOLc+KqzZW/D7r/MmXwD9ys45ax++nf/SGTboQxeG92JkCAnDQXe/ps5ObJD7XlReZLobfX3B+aC64B/lr2NUnDFG5CXKRdmnWXPClj0X+g4BrrcbHU0lVPST797qcN2sW7vMT7+aydju4czqE2Yw/ZznoQ5ZmZyeA/n7VOcRxK9LTKPmnHyHUaXP5IjtDX0vs+sZpQS65z4qqvMozDzDqgTaUbZuGvXgI8fNO0Je/92SPM5+YX88/sNNKwVwP8Nd+J/PHmZkPiHWRfWS1KNleSnb4uVH0J+Fgx4wrbtL37KnMXMf9z0Hwv7Kyoyk6Kyjph+eX/XX6D5gsJ7wqEtJjna2VsLtrEzNZPXru1ELWeW2E38HQqyHT5JSpRPEn15stJg1afQfiSEtbHtNf7BMOwlUzVx7ZeOja+6+vsD2PG7mflq5ULf9tK0pxmCuG+tXZtdk5TGpGW7ubFXBP1bhtq17XIlzIWg+tYv3Sgk0Zdr1ceQl2H72XyJDqOhWT9Y+B+nz3r0eIcTYOEL5gJfjzutjsY+mnQz9ymr7dZkVl4B//x+A01qB/KvKxw7I/U8+TnmQmybK12rPHQ1JYn+QrKPw8qPzUSWBu0r9tqSC7M5J0yyF/az5nNQ3jB8ovv2y58rqC7UawnJa+zSnFmzNYE9R7N449rOzh9Pvmsx5J102iQpcWGS6C9k1SeQmw4Dnqzc6xu0g173QNwUu/9LXm3lZ8Om70wCqVHP6mjsq2lPSFljVsWqAq01r/26jemr93LPxc25KNqCn1PCT+AfApEDnL9vcR5J9GXJOWEuwra+Ahp1qnw7A5+GGqHw8xPmAqKomoSfICfdfYdSXkh4D3Nx+djuSjehteatBdv5+K+d3Ngrgqcvs/G6kj0VFpiyB60vMyOKhOUk0Zdl9acmoVxcybP5EgEhcMkLZgm1nQvtE1t1tvYrM5wy0gMXemna09xXoftm4h+JfLBoB9f3bMqLV3ewplLj3hWQnWaGVQqXIIm+NLkZZlRHy0vNZJaqaj8KfAJhxx9Vb6s6S9sFSUuhy02eOS47tA34BVf6gux7CxN5d2EiY7uH8/LIjtatnJTwk/l9bzHEmv2L83jgX4sdrJkE2ceqfjZfwjcAIvuahapF5a37GpQXxNxodSSO4eUNTbpCcsUT/YeLdvD279sZ3TWcV0d1si7JFxVBwjyT5N2t5pAHsznRK6W8lVLrlFLzir8fo5TaopQqUkp1P2fbTkqpv4uf36SUCrB34A7ggC7OAAAgAElEQVSTlwkr3jcrR4V3L397W0UPgSPb4Xiy/dqsTgoLYN00aHEJ1HLCQhl2kHgog/unxbFo62G0rRdYm1Z84tTHf+3kjd+2cU2XJrx+rYVJHmD/OsjYb9kauKJ0FTmjfxhIOOP7zcAoYMmZGymlfICvgXu11u2BgUB+1cJ0om2/QNZR6PeIfduNHmzu5ay+cnb8AScPQtdbrI7EJgfSs7nli9X8vOkgt01Zw/WfrWRjyvHyXxhesYlTny3Zxau/bGVE58a8OaYz3lYvdL31J/DygVaXWhuHOItNiV4pFQ5cCUwqeUxrnaC13lbK5sOAjVrrDcXbHdVau08dgK3zi2fz9bVvu6GtIbixJPrKWvsV1AhziwSSnpXP+C9Wk5FTwJwH+vLCiPYkHjrJiA+W84/p69h7NKvsF5f8F5lS/gXZz5ft5uWfE7iyUyPeHusCSV5rU3s+sj8E2mGdWWE3tp7RTwSeBGwZH9gK0Eqp35RSa5VSpXZ0K6XuVkrFKqViU1NTbQzDwQryzJlj68vtP5tPKXNWv2ux1L+pqIxDsP1XiLkevJ1Yq6UScvILufOrNSQdyeLTm7vRuWltxveJZPETA/nH4Bb8Hn+QIW8v5oWftpCWmXd+AyUTp8pJ9F+uSOLFefFc3qEhE8fF4OPtApfbUrdC2k5oK6NtXE25vx1KqeHAYa11nI1t+gD9gBuL769RSp13+V1r/anWurvWuntoqJNrcJQlaSnknjDTth0hehDkHIf96x3Tvqfa8I3pznDxsfOFRZqHpq8jds8x3h7XmT4t6p96LjjAl8eHteavJwYxums4X65I4uLXF/Hhoh1k553zwd+0p7kgW0q/vtaaqSv38NzcLQxr14D3ru+CryskeTAXYVEyrNIF2TIvui8wQil1BRAA1FJKfa21vqmM7VOAv7TWRwCUUj8DXQHXH0S+7WdTF7z5QMe033wQoMx4+vBujtmHp9Ea1k6FiD4uvdi31ppn52xmQfwhnruqHcM7lX7BuEGtAF4d3Yk7+kXx2q9beeO3bUz9ew+PDWvFkDZhpJ7Mxdu3DS2zpvH1r4vZkR9GakYuhzNyOJyRy+ETuWTnFzK0bRgf3NDVdZI8mCJm4T0g2ElLFAqblZvotdYTgAkASqmBwD8vkOQBfgOeVEoFAXnAxcA7VQ/VwbQ2F2KjB4NvoGP2UaMeNI4x/fT2Grrp6fasMN0BF1rwxQW8/+cOvlm1l3svjua2vlHlbt+yQTCTxvdg5a6jvPJzAk/O3HjquTbKl1/9IXbpAv7wHUhYsD+hwf50Cq9NWLA/zeoFMa5HU/x8XCjJH9sDBzfCJS9aHYkoRaUrHSmlrgHeB0KB+Uqp9VrrS7XWx5RSbwNrAA38rLWeb59wHejAejixDwY949j9RA+GZRNNiYWAWo7dlydYNxX8a7l0TfPpq/fy9u/bGdW1CU9d1rpCr+3dvB4/PtCX3+MPse94NmHBAYTV7EnR9Jd4o2sOviNc/+IzAFvnmXvpn3dJFUr0WuvFwOLir2cDs8vY7mvMEEv3sfVnMxmn1WWO3U/0YFj6lrke4KhrAZ4iJx22/Aidr3PZyTe/xx/imdmbuLhVKK+N7lSpkgNKKYa1P6e7I7wbXgdsvSzmAhLmQYMOULe51ZGIUji5dqkL2zofmvZ2fEXE8J7gVxN2LJREX55NM80KRV3Pvgi7fMcRnv5hIzX8fIioG0SzekFE1Kthvq4bRJM6gU7pu47bk8aD36ylY5MQPrrRzv3lTXvC0rfNxCkX/ZA75eRhswziwKetjkSUQRI9wLEkOLzFrArlaD5+ZpyxjKcv37qp5iyxcddTD+1MPcl9X8dRr6Y/4XUC2X0kk7+2p5JbcHrkr7eXonHtAJrVrUFEvSA6NglhSNswwoLtN0E78VAGt0+JpXHtQL64tYf9672XTJzavw4i+9m3bXvbOh/QMtrGhUmiB9NtA6YksTNED4btv5giXfKvbukObjJJ7rLXTi0uciwzjzumrMHX24uvbu9J07pBABQVaQ5n5LLnaCZ70rLYezSr+D6T+RsP8M2qvSgFXSPqMKxdA4a1b0hU/cqfJR9Iz2b8F6vx8zFx1Kvpb5e3fJaSiVPJq90g0c+DOlEVX5xHOI0kejDDKkPbQr1o5+zvVDmERZLoy7J2Knj7Q6exAOQVFHHftDj2H89h+t29TiV5AC8vRcOQABqGBNCr+dldb1prth7MYMGWQyyIP8grv2zllV+20jKsJsPaN2BYu4Z0bBJSan2YoiJN8rEsth7MYFvxbevBEyQdzSLQ15tv7+59Vhx2FVQX6rWwaYaspXLSYddf0Ptez1ntywNJos9KM0P47F3b5kLqRUNIhOm+6XGH8/brLvJzYOMMM4IjqK4Zo/7jZlbuSmPiuBi6Natrc1NKKdo2qkXbRrV4eGhLUo5l8Uf8IRbEH+Ljv3bx4aKdNKwVwNB2YfSJrs+B9By2HTzBtoMZbD90kuz8wuJ2IKJuEK0bBHNlx0Zc3rERbRs5eNRUeE9IXGCG/rpqEt2+AIryzfq9wmVJok9cYPpCWzvxwqhS0GIwbP4BCvNdflq/022dZ2YQFxcwm7R0NzNik/nH4BaM7NKkSk2H1wni1r5R3No3iuNZefy59TALthxiVtw+vl65F4B6Nfxo3TCY63o2pU3DYFo3rEWrBjUJ8nPyn0vTHmZW8LEkqFv+2HxLJMyFmg2hiR0rvQq7k0S/db75RbXHAiMVET24eC3ZOIjo7dx9u7q1X0LtZhA5gD/iD/HfXxK4omNDHh3ayq67qR3kx6iu4YzqGk5OfiEJB04QXieI0GAH9LlXRnjxilMpa1wz0ednm9pQna/3zIVgPEj1Pjr5OWaYY+vLnf+LGjXAjNuX0TdnS9sNu5dAl5tJOHSSh79dR8cmIbw1JsahddYDfL3pElHHdZI8QFhbMxS3EguROMXOPyE/SyZJuYHqneh3L4H8TGvGswfWgSbdJNGfq3gVqSMtR3Pnl7EEB/jy2S3dCfSzczVRd1Cy4lQllxZ0uIR5Zk1kT1y/18NU70S/bb45Y4oaYM3+o4eYrpvsY9bs39UUFsD6byhsPoQ7Zx8gLTOPSeO706CW+yxQZnfhPeHg5gqtOOUUhflmtFqry+Uakxuovom+qMgUMWsxBHws+nc9ejDoIjM8TZiqnhn7mZzdn/XJx3lnXGc6NAmxOiprNT1j4pQrSVpmLpjLkoFuofom+n1xcPKQtbP5mnQzBbvcrPvmRE4+D3yzlru+imX/8Wz7Nbz2K7J86/LarkievKw1l3VoZL+23VV4D3Pvav30W+eBT+DpOSHCpVXfRL9tPihvaHmJdTF4+5huo52LSl1kwhXtPpLJNR8u57fNB1mWeIRh7yzh29V7bV/8ugw6fR9F237hq+w+XNW1Gfdd7KTJa67u1MSpWKsjOa2oyPTPtxwKfg6aMCbsqvom+q0/Q2Rf69e2jB4M6Xvh6E5r47DB0sRUrv5gGWmZeXx9Zy9+e2QAHZuE8PQPm7jli9WkHLvAWqgXsGLHEaZ/8l+8dCFbGo3ilVEdK1UF0mOF9zAXZF3lZGBfnFmovY1027iL6pnoj+6EI9ucO0mqLC2KV1nc6boLcGmt+XzZbsZ/sZrGtQOZ+2A/ejevR0S9IKbd2YuXRnZg7Z5jXPrOEqat2mPz2f2G5OPcNGkVN09awdCsnzlQvw/v3HsN/j7VcITNhYT3gMxUM3HKFSTMBS8ft1ioXRjVM9FvLV4HpY2TiphdSJ1IU+/GRfvpcwsKeWLmRl6cF88l7Row674+59WZual3M359ZABdIurwzOzN3DhpFclpZZ/dJx7K4J6psVz94XLiD5zgfz1TCSONRkMecI1Frl1N0zMmTllNa9M/HzUAAmtbHY2wUfX8q9r2MzToCLUjrI7EiB4Mu5dCQZ7VkZzlcEYO13+6kplxKTw8pCX/u7FbmeV4m9YNYuodPXllVEc2pqRz6cQlTP07iaKi02f3yWlZPP7dBi6duITlO47y6NBWLHlyEMOyfobgRo5f9MVdhbVznYlTh+NN1VUZbeNWql8JhMwjkLwKBjxhdSSnRQ+GNZNMP6yLlKTdmHKcu7+KIz07n49u7MoVHcsfAaOU4vqeEQxoFcrTszby7JwtzN90gCcva8Ocdfv4ZvVelFLc0S+K+wa2oG4NPzMTdsdCuPgpc3FanM+VJk4l/AQo1+j2FDarfn9Z2381Y9edVXveFpH9TZ/njoUukejnrN/HkzM3Ur+mPzPvu4j2jSs2lr1J7UC+ur0n38em8OK8eEZ9tAJvL8XY7k15aEgLGoWcsfj62i9NkbfiAmaiDOE9Ydk7kJdl3UgXrWHLbNOVFNzAmhhEpVS/RL/1Z6gVDo06Wx3JaQG1zB/yzj9h6HNO3XVhkeZEdj5pWXkcy8xjQfwhPl2yi56Rdfnopq7Ur+SiGkopxvZoSv9W9Zm34QBD2zU4f7GPglxTd77V5RBStaqUHi+8xxkrTvW1JobkVZC6Fa5615r9i0qrXok+L8sk0643u1597+jBsOhl07VUo75dm45NSmPh1sMcy8wjLTOPY1kl9/kcz8qj6JxBMtf3jOCFEe3x86n6JZxGIYHcNaCMxVUSfoKsI9Dj9irvx+OVTJxKWW1doo+dDH7B0OFaa/YvKq16Jfpdi8xi067UbVMiejAsegl2LYaO9vtDituTxg2TVlFUpKlbw4+6NfyoE+RHm4a1qFPDl7pBftQ54/GGIQG0DKvpnHHssZNNOeLmMruyXDXqQd1oSLZo5E1Wmum26XIT+Ne0JgZRadUr0W/9GfxDXKIf/DyNY8zkrZ2L7Jbodx/J5M4vY2lSO5BZ9/UxFz9dReo22LMMhj4vtcxt1bSnqf9uxYpTG6ZDYS50v825+xV2UX3+wooKzYXYlpe4ZrU9L29oPtBMnLLDDMi0zDxum2xGaUy+tYdrJXkwZ/NevhBzk9WRuA+rJk5pbY5XeA9o2NG5+xZ2UX0SffJq0x/sCpOkyhI9GDIOmAteVZCTX2gKjqXnMGl8dyLPvQhqtbwss0ReuxFQM9TqaNyHVROnkpbB0UToLtdS3FX1SfSbZ5kzyBYWFjErT/NB5r4Ks2SLijSPf7eBuD3HKryQttNs+QFy0qG7LIxeIaFtwbeG8xN93GSzwEj7a5y7X2E3np/oC/Jg/uOw5jPoMMoMZXRVtZtC/VZmPH0lvfbrVuZvOsC/rmhj0yQnS8R+AfVbQ7M+VkfiXrx9zMQpZ86QPZkK8XPNurC+geVvL1ySZyf69H0w5Qoz67TPP+Dqj6yOqHzRg2HPCrOebQVNXbmHT5bs4ubezbirfxlDGq22f72pftj9dtcb4uoOmvaEQ5tN95czrJ8GRfnQTS7CujPPTfS7l8AnA+BwAoz5Eoa95B5T7JsPMkNAk1dV6GV/bj3Ec3M2M7hNGM9d1c51y/zGTTYLVnQeZ3Uk7iniIigqgD3LHb+voiKImwIRfSCsjeP3JxzG8xK91rBsInx1NQTVg7sWQfuRVkdlu8i+phzCrkU2v2RTSjoPfrOOdo1r8f71XVy3AmTOCdj4PXQYbf06AO4qaoDpL9/0veP3tXsxHNstF2E9gItmhErKOQEzboI/noO2I+CuhRDayuqoKsY/uLgcgm2Jft/xbG7/cg11gvz4YnyPMqtLuoSNMyA/U2bCVoWPP7QbaVZ4cvSC4bGTIbCuGR0l3JrnJPrDCfDZILPg97CXYcwUkzTdUfQgOLDBzEa8gPTsfG6bvJqc/EIm39aDsFoBTgqwEkrGYjfqDI27Wh2Ne+s01nxgbvvFcfvIOGjWbYi5wXy4CLdmc6JXSnkrpdYppeYVfz9GKbVFKVWklOpeyvYRSqmTSql/2jPgUm2aCZ8NNmf043+CPg+694W+5oMAbcohlCG3oJD7p8WxKzWTT27qRqsGLv6hlrwaDm+Ri7D2ENHHFObb+J3j9rFuqimiJhdhPUJFzugfBhLO+H4zMApYUsb27wAOPOUACvPh1wkw6w4zY++eJdYVfLKnxl1MqYYyxtOfyMnn1i/WsHzHUV4d3Yk+LexbBM0hYr+Qglj24uVlymTs+MMUwbO3okKI+8pcD6jfwv7tC6ezKdErpcKBK4FJJY9prRO01tvK2H4ksAvYYo8gy7R3Jaz8CHrdC+PnQS0XHTdeUd4+ENXfnNGfUw7h8Ikcxn2ykjVJabwzrjPXdgu3JsaKKCmI1XmcFMSyl05jzRn3ltn2b3vHQrNgvVyE9Ri2ntFPBJ4EisrbUClVA3gKeKGc7e5WSsUqpWJTU1NtDOMcUf3h3uVw+Wvg42K1XKoqehCkJ5uFzIvtTD3JNR+tYM/RTL64tQfXdHGDJA9mLHZhriQOe2rQHsLaO6b7Jm4y1AiVVaQ8SLmJXik1HDistY6zsc0XgHe01icvtJHW+lOtdXetdffQ0CrUO2nYofKvdWUl5RCKh1mu3XuMa/+3gtyCQr69uzcDWrlJjZiiInMRtmlvk5yE/XQaY+rTp+22X5vpKab4X5ebPe/kqRqz5Yy+LzBCKZUEfAsMVkp9fYHtewGvF2//CPAvpdSDVQ202qnb3NRq37mIhQmHuOGzldQK9GXWfX3oFF7b6uhst3MhpO2U8raOUHK9Y9NM+7W5dqrpLuw23n5tCsuVm+i11hO01uFa60jgOuBPrXWZtWW11v211pHF208E/qu1/sBeAVcbSkH0IPJ2/sV9U1fTqkEws+7rQ7N6LlaJ8kK0hkX/hZAIaD/K6mg8T+2m0KyfmZ9gh9LWFBbA2q9MGY46kVVvT7iMSs+uUUpdA7wPhALzlVLrtdaX2i0yN7Ah+ThrktII9PMm0NfcAs74uuTxAF9vavh7E+Rn+49ba80vmW24ouAkNzU9wuO3X+7ak6FKs/032L8WRrwv3QCO0mkM/PQwHFhvRmtVReJvkLEfrnjDPrEJl1GhzKG1XgwsLv56NnDBS/5a6+crGZfLmxmXwtOzNlJw7oKrFxBeJ5AOjUPo0KQW7ZuE0L5xLcKCz5/kVFik+feczcxbX5fLAhTPtD6It7slea3NGrh1Ik3lQ+EY7a6Gn58wpSWqmuhjJ0NwI2h1mX1iEy7DzbKH9bTWvP/nDt7+fTt9W9TjrTExKGUW+8jOLyQ7z9zn5BeSnVdkHssv5ER2PgkHTrBl/wl+3XLwVHthwf50aBJCh8Ym+bduEMx/f05gQfwh7hvYGbW3C167FwMTLHvPlbJ1PhzcCCP/55orenmKwDrQcphZb2HYi2alsso4tseMy7/4Sfco/icqRI5oBeQXFvF/szczIzaZUV2b8OqoTvj5VLyKREZOPgkHMti8L53N+9PZsu8Ei7cdpuSfA6XguavacVvfKFg4yBRpy0k3xazcQVERLH7FLGbdcazV0Xi+jmNg6zxTsTV6UOXaWPul+cXreot9YxMuQRK9jU7mFvDAtLX8tT2Vh4a05NGhLStdCjg4wJeeUXXpGXV69afsvEK2HjzB5v0naBFak4ui65knmg+CpW+Z5dzauMm45oS5pmb6qM/k7NAZWl0G/rXMmPrKJPrCfDPapuUwCHGTuRmiQjynqJkDHTqRw9iP/2bZjiO8Nrojj13Syu713gP9vOkSUYebezc7neTBLDThG2RzNUvLFRWas/n6rUw5YuF4vgGmWmvCT5CfXfHXr/oEMg/LhDYPJom+HNsPZTCqeDbq5+O7M65HhHMD8PGHZn0rVJ/eUltmm8XNBz5d+f5iUXGdxkJeRsUrWiYtg9//DW2GmzN64ZEk0V/A3zuPMvp/K8grLGLGPRcxsHWYNYFED4KjO+B4sjX7t1VRISx+1Sxi3U4WknaqyH5mxExFFiRJ3wff3wp1o8xFc6kq6rEk0Zdhzvp9jP9iNQ1rBTD7/j50aGLhhdBzyiG4rE0z4Whi8dm8/Go5lZe36SpL/L3cdQwAKMiF724xXT3XfQMBtRwfo7CM/DWeQ2vNR4t38PC36+kSUZuZ9/YhvE6QtUGFtYWaDV27n76wAP56FRp0MP3Fwvk6jTULedtS0fKXp2BfrDmTD23t+NiEpdw60ecXFvHuH4mczC2wS3t5BUU8NWsjr/+6jRGdG/PVHT0JCXKBMeDF5RDYtdgMXXRFG2dA2i4YOEHO5q3SsBPUb11+983ar0yFyn6PyjKB1YRb/0WuSUrj3YXbGfHBMrYePFGlto5l5nHz56v4LjaFhwa3YOK4GPx9XOhiYvNBkJ0GBzdYHcn5CvPhr9dMonGXIaCeSClTEmHv33B8b+nb7FsL8/8JzQfC4GedGZ2wkFsn+j7R9fnmrt5k5BQw8sPlfB9buYuVOw5nMPKj5axLPs6718Xw2LDWeHm52IWp5gPNvSt236z/Bo7vgUHPyAU9q3UcY+5LO6vPPAIzboaaDWD0FzIqqhpx60QP0Lt5PeY/1I8uTevwxMyNPDlzAzn5hTa/fsn2VK75aAWZuabO+9UxTRwYbRUENzALTbjaBdmCPFjyhlnwu1W1qmnnmupEmtr/G787u6JlYQHMvB2yjsC4qVCjXplNCM/j9okeICw4gK/v7MU/Brfgu9gURn64nF2pF1z3BIAvVyRx25Q1hNcJYs6DfekaUccJ0VZB9CCzfGJeltWRnLZuqlkJS87mXUenMWYuw8FNpx/78z+w+y8Y/g40jrEuNmEJj0j0AN5eiseHtWbKbT04dCKHER8sZ/7GA6Vum19YxLM/bua5uVsY1DqMmfdeRJPagU6OuBKaD4LCPNi7wupIjPwcU54hvCe0GGJ1NKJEu2vAywc2FS8zuGU2LH8XetwJMTdYG5uwhMck+hIDW4cx/6H+tGpQkwe+Wcvzc7eQV3B6pEp6dj63TV7D1JV7uGdAcz65uZv71Hlv1ge8/Vynn37tV3BiHwz6l5zNu5Ia9aDFUNg0Cw7Fw48PmA/jS1+xOjJhETfJcBXTuHYg3959Ea/9upXPl+1mXfJxPryhCwWFmtu/XENyWhavX9uJsd2bWh1qxfgFQdNeZpil1fKzzdl8RJ/TF4qF6+g4xqz9OuUK8KsBY7+SxV+qMY9M9AB+Pl48O7wdPSLr8sT3G7jyvWUoBQqYdmfvsypHupXoQbDwP3DyMNS0qCQDmEUqTh6E0ZPkbN4Vtb4C/GpCbgaM/wlqNbI6ImEhj+u6OddlHRoy76F+RNQNokFwAHMe6Oe+SR7OKIew2NIwiJsMERdBVH9r4xCl8wsySziO+9p0+YlqzWPP6M/UrF4N5j7YF61xvfHxFdWos1lVaOefZsq7FdJ2w5HtcNmr1uxf2KaDLMgujGqR6AGUUp7Rw+DlbfrEdy4y46SteFOJC8y9lLUVwi14fNeNR2o+yPSPp261Zv+JC8wygfWirdm/EKJCJNG7o5Ll4qwYZpmXCbuXyixYIdyIJHp3VDvCnFFbUQ5h91IozIWWlzh/30KISpFE766iB0HSclNrxpkSfwPfGmZ5QyGEW5BE766aD4L8TNiz3Hn71Bq2LzAfMj7+ztuvEKJKJNG7qxZDwD8ENkx33j4PJ8CJFOm2EcLNSKJ3V76B0PFaiJ8DOenO2Wfib+ZehlUK4VYk0buzLjdBQQ5snuWc/SX+Dg06Qq3GztmfEMIuJNG7s8ZdIKwdrPva8fvKPmZq4beSs3kh3I0kenemlDmr3xdn+s8daeci0IXQUsbPC+FuJNG7u07jzCITjj6rT1xgauyEd3fsfoQQdieJ3t3VqA+tL4cN30JhvmP2UVRk+udbDJUFpYVwQzYneqWUt1JqnVJqXvH3Y5RSW5RSRUqp7mdsd4lSKk4ptan4frAjAhdn6HKzWfR5+2+OaX//OtO+jLYRwi1V5Iz+YeDMjuDNwChgyTnbHQGu0lp3BMYDU6sUoShf9BCo2dBx3TeJv4HyMmf0Qgi3Y1OiV0qFA1cCk0oe01onaK23nbut1nqd1np/8bdbgACllEyjdCRvH4i53vSjZxy0f/uJCyC8BwS58YItQlRjtp7RTwSeBIrK2/Aco4F1Wuvcc59QSt2tlIpVSsWmpqZWsFlxnpibzKiYDd/at92MQ6brRmbDCuG2yk30SqnhwGGtdVxFGlZKtQdeA+4p7Xmt9ada6+5a6+6hoaEVaVqUpn4Ls7Tfuq9NTRp72fG7uZdhlUK4LVvO6PsCI5RSScC3wGCl1AU7g4u7emYDt2itd1Y5SmGbLjfB0URIXm2/NhMXQHAjaNjRfm0KIZyq3ESvtZ6gtQ7XWkcC1wF/aq1vKmt7pVRtYD4wQWvtxNKKgnYjTQnhdXa6/l2YbyZKtbzEmiULhRB2Uelx9Eqpa5RSKcBFwHylVMnYvgeBFsCzSqn1xbcwO8QqyuNfEzpcA1tmQ+7Jqre3dyXknpBuGyHcXIUSvdZ6sdZ6ePHXs4vP9P211g201pcWP/6S1rqG1jrmjNthRwQvStHlZsg7aapaVlXib+DlC80vrnpbQgjLyMxYT9O0F9RrYZ8x9dsXQGRf8A+ueltCCMtIovc0JYXO9q6AIzsq386xJDiyTWbDCuEBJNF7os7Xg/KG9dMq30aiDKsUwlNIovdEwQ3NSJkN06GwoHJtJC6Aus3N+HwhhFuTRO+putwEGQdg558Vf21eFuxeIt02QngISfSequWlEFS/cmPqk5aaJQol0QvhESTReyofP+h8HWz7BTKPVOy1iQvANwgi+zkmNiGEU0mi92QxN0JRPmz8zvbXaG2GVTYfCD5SdFQITyCJ3pM1aAdNulWs0FnqNkjfK902QngQSfSerstNcHiLKTVsi8TiShZSllgIj+FjdQDCwTqMhl8nwM9PmOQdEl58awq1moBvwNnbJ/4ODTqYbYQQHkESvacLCIF+j0LsZFj8yvnP1wg9nfxrhcPev6HPQ86PUwjhMJLoq4OBT5tbQS6c2A/pKWfcks196nbYsdD05be72uqIhRB2JIm+OvHxh7pR5kwLUaUAAAUTSURBVFYaraEwT0bbCOFh5GKsOE0pSfJCeCBJ9EII4eEk0QshhIeTRC+EEB5OEr0QQng4SfRCCOHhJNELIYSHk0QvhBAeTmlbqxo6MgilUoE95zwcAey1sYkQIN2O20mb9m2zIvu29bh72s+ourcpx71y2zXTWoeW+0qttUvegNQKbPupPbeTNu3bZgX3bdNx97SfkbQpx92ebZ57c+Wum+MV2PYnO28nbVq3b1uPu6f9jKp7m3Lc7dvmWVyi66Y0SqlYrXV3q+MQziXHvXqS4+5YrnxG/6nVAQhLyHGvnuS4O5DLntELIYSwD1c+oxdCCGEHkugtppQ6Wc7zi5VS0nfpYeS4V09WHXfLE315b1x4Jjnu1ZMcd2tYnugFKKUGKqXmnfH9B0qpWy0MSTiBHPfqyYrj7hKJXilVUym1UCm1Vim1SSl1dfHjkUqpBKXUZ0qpLUqpBUqpQKvjFfYhx716kuPufC6R6IEc4BqtdVdgEPCWUkoVP9cS+FBr3R4zqWK0RTEK+5PjXj3JcXcyV1kcXAH/VUoNAIqAJkCD4ud2a63XF38dB0Q6PzyHK+DsD90AqwJxMjnuctzluDvhuLvKGf2NQCjQTWsdAxzi9JvPPWO7Qlznw8me9gDtlFL+SqkQYIjVATmJHHc57nLcnXDcXeWHGAIc1lrnK6UGAc2sDsgZlFI+QK7WOlkp9R2wEUgE1lkbmdPIcZfjLsfdCcfd0kRf8saBacBPSqlYYD2w1cq4nKg9sBNAa/0k8OS5G2itBzo5JoeT4y7HHTnuTj3ulpZAUEp1Bj7TWve0LAiLKKXuBR4CHtFaL7A6HmeS4y7H3epYnM3q425Zorf6jQtryHGvnuS4W0uKmgkhhIdzlVE3QgghHMRpiV4p1VQptah45tsWpdTDxY/XVUr9rpRKLL6vU/z4JUqpuOKZc3FKqcFntNWt+PEdSqn3zphsIVyMnY/7y0qpZKmX4vrsddyVUkFKqflKqa3F7bxq5ftyW5Vdg7CiN6AR0LX462BgO9AOeB14uvjxp4HXir/uAjQu/roDsO+MtlYDF2EmXvwCXO6s9yE3S4977+L2Tlr9vuTmnOMOBAGDir/2A5bK33sljoeFvwhzgEuAbUCjM345tpWyrQKOAv7F22w947nrgU+s/kHKzbHH/ZzHJdG72c0ex734uXeBu6x+P+52s6SPXikVifkEXwU00FofACi+DyvlJaOBdVrrXMx06ZQznkspfky4uCoed+Gm7HXclVK1gauAhY6M1xM5fcKUUqomMAszzOpEed3rSqn2wGvAsJKHStlMhg65ODscd+GG7HXciydbTQfe01rvclC4HsupZ/RKKV/MQZ+mtf6h+OFDSqlGxc83Ag6fsX04MBu4RWu9s/jhFCD8jGbDgf2Ojl1Unp2Ou3Azdj7unwKJWuuJjo/c8zhz1I0CPgcStNZvn/HUXGB88dfjMX15Jf+mzQcmaK2Xl2xc/O9ehlKqd3Gbt5S8Rrgeex134V7sedyVUi9h6uM84ui4PZXTJkwppfphrphvwpQmBfgXpt/uOyAC2AuM0VqnKaX+D5iAKfpTYpjW+rAyaypOAQIxo27+oZ31RkSF2Pm4vw7cADTG/Bc3SWv9vFPeiKgQex13zEibZEw9nJI++w+01pMc/iY8iMyMFUIIDyczY4UQwsNJohdCCA8niV4IITycJHohhPBwkuiFEMLDSaIXQggPJ4leCCE83P8D19l9YQHSM0MAAAAASUVORK5CYII=\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "indexed_data[-30:].plot(y = [\"SAFitFilled\",\"Filled\"])" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Fluctuations saisonales\n", + "\n", + "Nous ponvons isoler les fluctuations saisonales afin de les étudier indépendemment. En voici une représentation (sur une échelle de quelques années par souci de lisibilité)." + ] + }, + { + "cell_type": "code", + "execution_count": 165, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 165, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXYAAAEHCAYAAAC6IG0BAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzt3Xd8nFeZL/DfmaIZSdNkaUbdktUsW7Ylx05x4pLYToGQZEPgLhA6Kexd2n5gdym7S7mwd4G7LASSgElY9i7ZzYWYkISEBFsuspPYjuS4SxoVS1adot5G0879Y2ZsxZaikeadt83z/Xz0wZbG7xzl8D5z3nOecx7GOQchhBD10EjdAEIIIcKiwE4IISpDgZ0QQlSGAjshhKgMBXZCCFEZCuyEEKIyFNgJIURlKLATQojKUGAnhBCV0Unxpjk5Oby0tFSKtyaEEMVqamrycs7ti71OksBeWlqKxsZGKd6aEEIUizHWHc/raCqGEEJUhgI7IYSoDAV2QghRGQrshBCiMhTYCSFEZSiwE0KIylBgJ4SQefgCIambsGwU2Jegb3QGVEowtQRCYTx1pBMjU36pm0JEdKxzCOu/9RoOOz1SN2VZKLDH6c2OIWz9/gE8ebhD6qYQEf3q6EV89+Vm/Oxgu9RNISIJhML4pxfOIRDieEKh/U6BPQ6+QAhff/4sOAf+/fUuzAaV+4hG4jcwNoOf1LdBq2F49sQljM0EpG4SEcF/vtkNp2sS2ypzcPziMM70jkrdpCWjwB6Hxw+246J3Co9uL4NnYhYvnR6QuklEBN/9YzNCYY7HP3IdpvwhPHM8rt3cRME8E7P4t31O7Kiy4/EHr4PJoMMvj1yUullLRoF9Ea2DE3jyUAfev7EQX31PNapyTXjqSCfNtavckTYPXj47gM/dVoG71uVha0UOfk1Pa6r3/Vdb4AuG8M171sJi1OPDNxTjlbMD6BudkbppS0KB/V2Ewxxf+/0ZmI06fOPuNWCM4aGtZWgZnMAbHUNSN48kyWwwhG++cB6l2Rl4eHsZAOCR7WVwT8zihVP9EreOJEtT9wiea+rFQ9vKUGY3AQA+ecsqAMC/H1XWqJ0C+7t45ng3Tl4axT++by2yTQYAwL11BcgxpeGpI50St44ky1NHLqLTO4Vv3VsDo14LANhWmYPqPDN+2dCJcJie1tQmFOb45ovnkGcx4nO3VVz+fqEtHXevz8ezb/Vg3KecNRYK7AsYHPPh+6+2YmtFDu7fWHj5+0a9Fh+7qRQHWz1od09I2EKSDH2jM/jpgTbcVZOHW1c7Ln+fMYZHtpehzT2p2BQ4srBn37qEc33j+Prda5BpeOdp5g9vK8PkbBDPnrgkUeuWjgL7Ar754jkEQmF87/51YIy942cfvWkl0nQa/Or1LmkaR5LmOy+dBwD84z1rr/nZPbUFyLca8YsGSnlVk5EpP374WituXLUC92zIv+bn64usuKlsBf799S4EQmEJWrh0CQd2xpiRMXaCMXaaMXaeMfZtIRompVfPDeK18y58aXcVSrIzr/l5tsmAB64rxN6mXgzTxhXVONjqxmvnXfj8zkoU2tKv+bleq8Gnb1mFY53KTIEj8/vXfa2Y8AXx7ftqrhnExTy8rQwDYz68clYZGXFCjNhnAezknNcCqANwF2PsJgGuK4lxXwDffPEc1uRb8NC2VQu+7tO3rMJsMIxnjlEKnBr4AiF868XzKLNn4uFtZQu+7kM3FMNs0GFPA62xqMG5vjE8c/wSPr6lBNV5lgVfd9tqB8rtmfilQjLiEg7sPGIy+ld99Ev+v/kCfvhqK9wTs/iX96+HXrvwf57KXDNuXW3Hf7zZTSlwKrCnoRPdQ9P49r01SNMt3O9mox4fuXElXjk7gJ7haRFbSIQWDnP80wvnkJ2Zhi/trnrX12o0DA9tK8O5vnG82Sn/jDhB5tgZY1rG2CkAbgD7OOfH53nNI4yxRsZYo8cjz8Wnpu5h/OZ4Nz55cylqi22Lvv6hrWXwTs7iRUqBU7Se4Wk8frAdd6/Px7bKResE41O3rIKGMTytsBQ48k7Pv92Hk5dG8fd3VcOarl/09fdvLER2ZhqeUsCGJUECO+c8xDmvA1AE4AbG2Lp5XrOHc76Zc77Zbl/85hGbPxjG135/FvkWI758x+q4/s0tFdmozjPj6aMXFfF4Rub37ZcuQKth+If3rYnr9XlWI+6tK8D/e6sHo9O0xqJE474A/vefWrBxpQ0PXFcU178x6rX42JYSHGhxyz4jTtCsGM75KIBDAO4S8rpi2NPQAadrEt+5bx1MV6U7LYQxhs9sXYWWwQm83i7/xzNyrfpmF/Y3u/DFXZXIt167YLqQR7aXYSYQwm9ojUWRfrK/DUNTs/j2vTXQaOZfMJ3Px24qgUGnkf3TmhBZMXbGmC3653QAuwG0JHpdMXV6JvHYgcij+O61uUv6t5ENSwY8dZQW05TGFwjhWy+dR4XDhE/dsvBC+Xyq8yzYUWXHr9/oVvS53anI6ZrAr9/owoeuX4kNRYtPuc6VbTLggU1F2HuyD56J2SS1MHFCjNjzARxkjJ0B8BYic+x/FOC6ouCc4+vPn4VBp8E37702d3kxBp0WH99SgkOtHrS55P14Rt7pyUMd6BmewXfue/cF04U8uj2yxvKHt/uS0DqSDJxzfOvF8zAZdPjbO+Obcr3aZ7augj8Yxn/K+GlNiKyYM5zzjZzzDZzzdZzz7wjRMLHsPdmHY53D+Pp718BhNi7rGg/euBIGnQa/el3ej2fkiu6hKTx5uAP31hbg5vKcZV1jS3k2agos2HOEjhlQilfODuKNjiF85c7VWJGZtqxrlNtN2L0mF785Jt+ntZTfefr/3rqE1blm/OXm4mVfI9tkwPuvizyeDU3K9/GMXPHj/W3Qaxi+cXd8C6bziR0z0OmZQn2LW8DWkWTZ09CBqlwTPnLDyoSu8/C2VRie8mPvyV6BWiaslA7sQ5OzaOoewV3r8pa0gDKfz2wthT8Yxm+OKec8iVQVCIWxv9mF967PR65leU9pMXevz0ehLR2/pA1Lsuca9+F07xjuqyuENsH7/YZVK7ChyIqnj1yU5dNaSgf2g60ehDmwe83SFkznU+Ew47bVdvznsS7ZPp6RiLe6hjHhC2KXAP2u02rw6a2rcKJrGG9fGhGgdSRZDkSfqnatcSzyysUxFtmw1OmV59NaSgf2/RdcyLUYsK5w4a3ES/HQtjJ4J/20YUnm6pvdSNNqsK1yeXPrV/vQ9cWwGOmYAbmrb3ah0JaO1blmQa733nV5kac1GR7hnbKB3RcIoaHNg91rchc8+Gepbi6nDUtyxznH/mYXtpRnX3M863JlGnT46E0lePX8ILqHpgS5JhGWLxDC0XYvdq9xCHa/67QafOqWUpy4OIzTPfI6FC5lA/uxziFM+0NLzlt/N7HHs1bXBI62ewW7LhFOh2cS3UPT2C3A4/hcn7y5FHqNRhHbzVPR6+1e+AJhQabf5vrL6yOHwsltw1LKBvb9zS5kpGmxpSxb0OveU5sPu9lAN7hM7W+OzbMKe4M7LEbctS4PL58doKc1Gdrf7IbJoMONZSsEva7ZqMed6/JwtN0rq35PycDOOcf+C25sq8y5XPpMKAadFu/fWIg3Ory0iCpD9c0urM23oGCe89YTdf2qFRie8qN3RFmFj9UuHOaob3Zhe1UODDph73cAqCu2ya7fUzKwn+8fx+C4T5BsmPlsXJmFQIijeWA8KdcnyzMy5UdT94jg0zAxG6Mngp6S2XxrqjvXPwb3xCx2VSfnfq+L9vtpGRVfScnAvu+CC4wBO6uTc4Nf7mi6wWXlYKsbYS78NEzM6jwz0nQa6neZ2d/shoYBtyXpfpdjv6dkYK9vcWHTyixkmwxJuX6e1YhciwGne8eScn2yPPXNbtjNBqwvtCbl+nqtBusKLLIauZHI9Nt1K7OWfYTAYvRaDWoKLDjdI5/7PeUC+8DYDM71jSdt1BZTW2ST1Sd4qvMHwzjs9GBXtSPhXcbvprbYhrN9YwgqpOix2g2MzeB8vzj3u5z6PeUCeywr4va1yXksi6kttqHTO4Wx6UBS34fE5/jFIUzOBpO2rhJTV2yDLxCG0zW5+ItJ0tWLdL/XFdswEwih3SOPfk+9wH7BhdLsDJTbTUl9n9g8+5k+GrXLQX2zGwadBrdUCLPbdCG1RfJbSEtl+5tdKBHhfq+V2bpaSgX2ydkg3uwYEnS36ULWF0XmceXS0aksttt0a0UO0tOET3ebqyQ7A9Z0PfW7DEz7g3ijYwi7qpN/v5dmZ8Bi1OGUTObZUyqwH23zwB8SfvfZfCxGPcrtmZT6JgNO1yR6R2ZE6XfGGGqLbdTvMnCkzQt/MJy09Na5Yv0ulw/0lArs+y64YU3XY3NplijvF7nBx2S1Iy0V7W92ARDmVL941BVZ4XRNYNofFOX9yPzqm10wG3W4fpWwu00XUltkQ6trQhYbE1MmsIfCHAdaXLhttR16rTi/9sZiG7yTs+gf84nyfmR++5tdWF9oTfjs9XjVFtsQ5sC5PtqgJpVwmONAiwc7qsS732uLbQiFOc73Sz8dkzKB/eSlEYxMBwQ99GsxcltQSUXeyVmc6hlNejbMXLECydTv0jndOwrv5CxuF/N+j66ryWGePeHAzhgrZowdZIw1M8bOM8a+KETDhLb/ggt6LcP2Krto71mdZ0GaVl470lLNgRY3OBdvGgYA7GYDCm3pOEWZMZLZ3+yCVsNwa5V4/e6wGJFvNcrifhfiQOoggC9zzk8yxswAmhhj+zjnFwS4tmD2Nbtw46psWIx60d4zTafB2gILLaRJqL7ZhXyrETUFwhRTiVedjBbSUlF9sxubS7JgzRDvfgci8+xnZPCBnvCInXM+wDk/Gf3zBIBmAIWJXldInZ5JdHqmRFkdv1pddCdiSIZ1EdXOFwjhSJsXO6uFK64Qr9piK3pHZuCl4uai6x2ZRsvghKjTbzG1xTZ0DU1jdNov+nvPJegcO2OsFMBGAMeFvG6i6pN0Bnc8aoutmPaH0O6Wx460VHK5mIoU/R6dZ5fD6C3VXLnfxR/I1RZH969IfE6UYIGdMWYCsBfAlzjn16QDMMYeYYw1MsYaPR6PUG8bl33NLlTnmVG8IkPU9wXm7ESkx3LR1Te7ka7XYku5sMVU4rG+yAoNk8dCWqrZ3+xCWU4mypK823Q+6wutYEz6+12QwM4Y0yMS1J/hnP9+vtdwzvdwzjdzzjfb7eItYI5M+dHYNSzq6vhcpdmZkR1pNHITFeeR4grJKKYSj4w0HapyzZLf4KlmwhfAsc4hUbPf5jIb9Si3myTvdyGyYhiApwE0c85/lHiThJXsM7gXo9HIa0daqrgwMI7+seQVU4lHXbENp3tHaYOaiI60eREIcexK0tnr8agtsuF0r7QbE4UYsd8C4GMAdjLGTkW/3ivAdQURO4N7Q5LO4I5HbZENLYPy2JGWKuqb3WBJLK4Qj9piG0anA7g0PC1ZG1LN/mYXrOl6bCoRZ3f5fOqKrZJvTBQiK+Yo55xxzjdwzuuiX68I0bhEzQZDOOz0YPea5J7BvRg57UhLFfXNLtQW2WA3J6eYSjxi6yuU7iqOUJjjUKsHt622QyfSbtP5yGGDmqp3nh7vHBblDO7FyGlHWipwj/twundMkvTWuapyTTDqNRTYRfL2pREMT/klm3aNqc43S74xUdWBfX+zC0Z98s/gXozDYkSB1Ug3uEgOtETS3aRaQIvRaTVYX2il9RWR7G92Q6dh2LFavOSM+Rh0WqyRuESiagM75xz7L7iwtcIuSVbE1WgBVTz7m10otKVjda5Z6qagtsiGc/3jCMikZJqa1Te7cGPZClF3ly+krsiKs73SbUxUbWBvHphA/5gv6SWx4lVbbMOl4WkMT0m7I03tfIEQjrZ7sXuN+LtN51NbbIM/GEbr4ITUTVG17qEptLknsata2qe0mA1FNkz5Q+iQqFSeagP7/mYXGAN2yqSjqWSaOF5v98IXEKeYSjxiJRJpGi65YrWMpV5Pi6mVuN9VHdjriqXNiphrfZE8dqSp3f5mNzLTtLixTJziCospykrHisw06vckq292odJhwsps8XeXz6csJxNmg06yIyVUGdi7vFM40zsmm09vADAZdKh0SL8jTc04jxRT2V5lh0En/boKEC2ZVmSlJ7UkGpsJ4MTFYdk8pQGRjYkbiq04LVEmnCoD+y8aOpCm0+CDm4ukbso7yGFHmpodbffCNS5ucYV41Bbb0OaexOQslcpLhmdPXEIwzHH3+nypm/IOG4psaB4Yl2RjouoC++CYD3ub+vA/NhfBYRanFFq86lbaMDzlR+/IjNRNUR3OOR6rb0OexYi7N8jrBq8ttoFz4KzEJ/6p0Yw/hD0NndheZcf6Iul2l8+ntsiGYJjjwoD4JRJVF9ifOtKJEOd4dHu51E25Bu1ETJ5jncN4q2sEn91RJptpmJg6WjhPmmeOd2Noyo8v7qqQuinXqJOwNKaqAvvIlB//deIS7tmQL8kRvYtZnWeGQUel8pLhsfo22M0GfOiGlVI35RpZmWkoyc6gfheYLxDCLxo6cXN5NjaVyGOxfK48qxG5FgPOSPCkpqrA/h9vdmHaH8Jf3Sq/T28A0Gs1WFdIC2lCe6trGG92DuHR7WWy2Iw2n9oi2qAmtGdPXIJnYhZf2FUpdVMWtEGifldNYJ+aDeLXb3Rh9xoHVudJv+NwIbVFkVJ5QdqJKJjH6tuQY0rDgzeWSN2UBdUW29A/5oN7XLoT/9TEFwjhycMduGHVCtxUJn4hlXjVFdvQ6Z3C2HRA1PdVTWD/7xOXMDodkO1oPaa22ApfIAyni0rlCeHkpREcafPi4W1lSE+T52gdiBzlCkhfMk0tftfUC9f4LL4o49E6MKdEYp+4o3ZVBPbZYAhPHbmIG1etkPQc5nhcXlCh6RhB/LS+DVkZenz0JvmO1gGgpsAKrYbRdIwA/MEwnjzYjk0lWbhZgrKHSxHL1BF7nl0Vgf0Pb/dhcNyHv75N3qN1AFi5IgO2DD3d4AI40zuKg60ePLStDJkGndTNeVdGvRbVeWb6QBfA3pO96B/z4Qu7KmVxHtC7sabrUZaTKXomnOIDeyjM8fPDnVhXaMG2SmmP541HZCeijVIeBfBYfTus6Xp8fIu8R+sxsRM+wxKd+KcGgVAYTxxqR22xDdsVcL8DkX4/1SNuiUTFB/ZXzw3ioncK//PWCtl/esfUFtvgdE1ginYiLtu5vjHsb3bhM1tXwSyDY1rjUVdkw7gviK6hKambolh/eLsPPcMz+OIuBd3vRVZ4JmYxKOLCuaIDO+ccjx9sR1lOJu6syZO6OXGrK7YizCPBiSzPzw60w2zU4RM3l0rdlLjV0vpKQoKhMB4/2I51hRbctloex3HH43K/i3hujKID+2GnBxcGxvHZHeXQSljTdKk20E7EhLQMjuPV84P41M2lsKYrY7QOABUOEzLStJIdDKV0L53pR9fQND6/U/5z63OtybdAr2Wi3u+CBHbG2K8YY27G2DkhrhevJw51IN9qxF9sLBTzbROWYzKgKCudbvBl+umBdmSmafHpraukbsqSaDUM6wuttL6yDKEwx08PtKM6z4zbZXSKYzwiC+cWURMmhBqx/xrAXQJdKy5N3cM4cXEYD20rQ5pOeQ8esQUVsjRtrgm8cnYAn7i5FLaMNKmbs2R1xTZc6B+HP0gb1Jbi5bMD6PRM4Qu7KqFR0NN5TG2xFWd6x0RbOBckInLOGwAMC3GteD1xsANZGXp8+IZiMd9WMHVFNvSNzsAzMSt1UxTlZwfbka7X4qFtZVI3ZVlqi23wh8JoGRT/xD+lCoc5flrfhkqHCXcpaC1trtoiGyZng+j0irMxUbShLmPsEcZYI2Os0ePxJHStlsFx1Le48albViEjTd75ywuJLahIVWFFiTo8k3jpdD8+dlMJVmQqb7QOUKm85Xj1/CDa3JP4vEJH68Dckx7FmX4VLbBzzvdwzjdzzjfb7faErvXkoQ5kpmnxiS2lwjROAusKLdBQqbwlefxgO9J0GsWO1gEg32qE3WygwB6ncDhyzn6ZPVN2hTSWosxuQmaaVrQFVMVNTl8amsZLp/vx4E0lsGYoJyPiahlpOlTlmnGKzg6JS5d3Ci+c6seDN5bIpo7tcsQ2qNEHenz2NbvQMjiBz++sUFTm29W0Gob1RVbR+l1xgf3nDR3QaTT4jMIyIuZTF92JSKXyFvfEoXZoNQyPblfuaD2mrtiKDs8Uxn3invinNLGqWCXZGbhnQ4HUzUlYbbENFwbGMRtMfqk8odId/xvAmwBWM8Z6GWOfEeK6V3OP+/BcYy8e2FSEXIu8yt4tR22xDWMzAXQPTUvdFFnrGZ7G70/24SM3rIRDJf0OUKm8xRxsdeN8/zj++rYK6LSKG4Neo67IhkCIo2VgIunvJVRWzIc55/mccz3nvIhz/rQQ173a00cvIhgO47M7lD9qA64c6Ukbld7dE4faoWEMj6qk39fmWwAArYPJv8GVinOOH+9vQ1FWOu5X2D6VhdSttGFXtTg7ZhWVUrK1Mgdmow4l2ZlSN0UQFQ4TtBoGp4tu8IVc9E7ht429+OiNK5FvTZe6OYJYkZkGW4Ye7R46k38hr5134UzvGH7wgQ3Qq2C0DgD51nQ8/cnrRXkvRQX2bZV2bKtMLKNGTtJ0GpRkZ6DDTYdCLeRH+5xI02rw1zvlfyRzvBhjqLCb0O6mwD6fUJjjR/taUWbPxPtVMloXmzo+ChWs3G6ikdsCzveP4aXT/fjULaVwmJU/tz5Xud2ETur3eb10uh9O1yT+ZneVKubWpUD/1SRW4TChyzuFANVAvca//tkJi1GHR7eXS90UwVU4TPBO+jE67Ze6KbISCIXxb/udWJNvUXTeutQosEuswm5CMMxxaZgyY+Zq7BrGgRY3PntruaL3KyykwmECENlNS654rqkX3UPT+PLtVYrdZSoHFNglVh69wWm+9QrOOX7waivsZgM+qaDz1pei3E79fjVfIITH6ttQV2zDrjXKOW9djiiwS6zcHsnwoZHbFYecHpzoGsbnd1Yo9iygxRRmpcOg06DDQwvnMf91/BIGxnz42ztXK+q8dTmiwC4xs1GPXIuBRm5R4TDH/3mtFcUr0vGh61dK3Zyk0WoYVuVkUr9HTfuDeOJQO7aUZeOWCmXUMpUzCuwyUOEw0cgt6pVzAzjfP46/2V2lyHP2lyLS7xTYAeDXb3TBO+nHV+5cLXVTVEHdd45ClNtN6HBPpvyZMcFQGD/6sxNVuSbcV6f+/OVyuwk9w9PwBZJ/doicjc0E8IvDndhZ7cCmkiypm6MKFNhloMJhwuRsEO4UL7rxXFMvOr1T+ModqxV9kl+8KhwmhHlkd20qe/pIJ8ZmAvjyHVVSN0U1KLDLAGVIRDIifhLNiLh9rbJqWi5XrN9TeTpmaHIWTx+9iLvX56OmwCp1c1SDArsMUE4z8Jtj3RgY8+HvUigjosyeCcZS+wP954c7MBMI4W9ur5S6KapCgV0GHGYDTAZdyt7gk7NBPHGoA1srcnBzCmVEGPVaFGWlp+zCuWvch//7Zjf+YmMhKhxmqZujKhTYZYAxhnJH6h4K9dSRTgxP+fG3KZgRkcqHgf30QBtCYY4v7aK5daFRYJeJCntqpr4NT/nx1JGLuLMm93IBilQSOwwsHE6tjKie4Wk8e6IHf3l9MVZmZ0jdHNWhwC4T5Y5MuMZnU65c2pOH2jHtD+Ird6TeaB2IrK/MBsPoG52Ruimi+vH+Nmg1DJ/fSXPryUCBXSYqohkSnSk03zowNoP/eLMb928sQmVuas6xVqTgWUHt7gk8/3YvPnZTCfKs6jqOWS4osMtEKh4G9lh9Gzjn+NLu1B21pWLK47/ta0O6Xou/ulV9xzHLhVDFrO9ijLUyxtoZY18V4pqpZuWKDOi1LGVu8L7RGfy2sRcP3liC4hWpO8ealZmG7My0lPlAb3dP4OWzA/j01lXINhmkbo5qJRzYGWNaAI8DeA+AtQA+zBhbm+h1U41eq0FJduocCrW3qRehMMdntq6SuimSK0+hhfPfNfZCp2H4hEqPY5YLIUbsNwBo55x3cs79AJ4FcJ8A1005qZIZEw5zPNfUi5vLs1N6tB6TKqmuwVAYz7/dh9uqHcih0XpSCRHYCwH0zPl7b/R778AYe4Qx1sgYa/R4PAK8rfqUOzLRPTQNf1DdZfLe6hrGpeFpfGBTkdRNkYVyeyZGpgMYnlJ3mbwj7V64J2bxwHXU78kmRGCfb//3NUm5nPM9nPPNnPPNdrtdgLdVnwqHCaEwx6VhdWfG/K6pFyaDDu9ZRzUtgdTJjHmuqRdZGXrsrKbqSMkmRGDvBVA85+9FAPoFuG7KqbBHUv7UfINPzQbxytkBvG9DPtLTtFI3RxZSITNmbDqAfRdcuK+uUPXn7MuBEP+F3wJQyRhbxRhLA/AhAC8KcN2UUxYtk6fmwP7K2QFM+0M0DTNHoS0dRr1G1f3+0pl++INh6neRJFxQknMeZIx9DsBrALQAfsU5P59wy1JQpkGHAqtR1YdC/a6pF6tyMqmgwhwaDUNZjroXUJ9r6kV1nhk1BRapm5ISBHkm4py/wjmv4pyXc86/J8Q1U5WaMyS6h6Zw4uIwPrCpKGWO5o2XmsvktbsncapnlPpdRDTZJTOxnGY1lsnb29QLDQPef536y94tVbndhL7RGcz41Vcmb+/JXmg1LCXKHcoFBXaZKXeYMO0PYWDMJ3VTBBUOc+w92YetlXbkW9Olbo7sVDhM4Bzo9Kpr1B4Kc/z+ZC9urbLDbqbcdbFQYJeZCpVmSLzRMYS+0RlaPFuAWlMej7Z74RqfpX4XGQV2mSl3qDMz5rmmHliMOtyRIvVMl6o0JwMaBtUtnO9t6oUtQ4+dayh3XUwU2GXGbjLAYtSpasQ+7gvgT+cGcW9dAYx6yl2fj0GnxcoVGehQ0Qf62EwAr50fxL21BTDoqN/FRIFdZtRYJu+PpwcwGwzjA5uKF39xClPbYWAvn4n1O03DiI0CuwxFDgNTzyP5c009qHSYUFtklbopslbhMKHTO4WQSsrkPdfUg6pcE9YXUr+LjQK7DFXz63GOAAAWDUlEQVQ4TPBMzGJsRvll8trdkzh5aRQf3Ew5zIspt5vgD4bRMzwtdVMS1umJ9PsD11G/S4ECuwzFzg5Rw3TMc02RHOa/2Eg5zIuJVdFSw3TM3pORPQv3U79LggK7DFWo5AYPhTmefzuSw+wwU23LxVSo5AM9krvehx1Vdjgs1O9SoMAuQ0VZ6UjTahSfIdHQ5qEc5iWwZuiRYzIo/gP9zY4hDIz58AD1u2QosMuQTqvBqpxMxd/gzzVGzt/etYZy1+NVbld+ecTYnoXd1O+SocAuU+UOZd/go9N+On97GSKHgU0p9qygcV8Ar56nPQtSoztOpirsJlwansZsUJmHQr14uh/+EOUwL1W53YSxmQC8k8osk/fKmQH4ArRnQWoU2GWq3GFCmANdXmWmvv2usRdr8i1YRznMS6L0M2P2nuxFuT2T9ixIjAK7TCm5XFrL4DjO9o3hgzRaXzIlZ0R1eafwVtcIPrCpmHLXJUaBXaaUnMv+XGMvdBqG++oKpG6K4uRbjchI0yqy3yl3XT4osMtUepoWhbZ0xY3cAqEw/nCqD7vWOJBtovO3l4oxpsgzY8LR3PWtlXbkWSl3XWoU2GWsQoGHgR1q9cA76ccHafFs2crtmYrbw3Csk87bl5OEAjtj7IOMsfOMsTBjbLNQjSIRsZFbWEGHQr1wqg/ZmWnYsdoudVMUq8JhQv+YD1OzQambEre9J/tgpvP2ZSPREfs5AO8H0CBAW8hVKhwm+AJh9I/NSN2UuITCHEfavLh1tQN6LT0MLldsfaVTISd8hsMch1rd2FntoNx1mUjo7uOcN3POW4VqDHmncruyqimd6R3F2EwA26typG6KoiktM+bCwDiGpvzYUUVPaXIh2rCKMfYIY6yRMdbo8XjEeltFu3KDK2PkdqTNC8aAbZV0gyeiJDsTWg1TzAf6YWfkfqZ+lw/dYi9gjO0HkDfPj77BOX8h3jfinO8BsAcANm/erJxJYwmtyEyDLUOvmBu8wenB+kIrVmSmSd0URUvTaVCyIkMx/X7Y6UFNgQV2M2VBycWigZ1zvluMhpBrMcai1ZTkf4OP+wJ4u2cUf7WjXOqmqEKZQvp9whfAye4RPLy9TOqmkDlohUvmKhwmRaS+vdHuRSjMsZ3mWQVR4TCha2gKwVBY6qa8qzc6hhAMc5pfl5lE0x3vZ4z1AtgC4GXG2GvCNIvElNtNGJryY2RK3odCHXZ6YTLosHGlTeqmqEK5PROBEMclmZfJa3B6kJmmxXUrs6RuCpkj0ayY5znnRZxzA+c8l3N+p1ANIxFKyJDgnKPB6cGW8mxKcxSIEg4D45zjsNODLeU5dDSzzFBvyJwSDgPr9E6hb3SGpmEEVK6AjKiL3in0jszQZjQZosAuc4VZ6TDoNLIeuTVE0912ULqbYCxGPRxmA/U7WRYK7DKn1bBomTz5jtwanB6UZmdgZXaG1E1RlUg1JfkG9sNOD1blZFK/yxAFdgWQ82Fgs8EQjnUO0zRMEpTbIxlRciyTd7nfK2mXsRxRYFeAcrsJPSPT8AXkVyavqWsEM4EQttPjuOAqHCZMzAbhnpiVuinXaIz1O32gyxIFdgWocJjAeWSxSm4Ot3mg1zJsKc+Wuimqc3nhXIZPa4edHqRpNbipjPpdjiiwK4Ccqyk1OL3YVJKFTMOim5jJEl1OeZThPHuD04PNpdTvckWBXQHK7JlgTH4pj+4JH5oHxulxPElyLQaYDDrZjdgHx3xoGZyg3aYyRoFdAYx6LYqz5Hco1BGnFwBofj1JImXyMmU3Ym9oi6Q50ge6fFFgV4hye6bsAntDmwfZmWlYm2+RuimqVeEwo3VQZv3u9MBhNqA6zyx1U8gCKLArRFWeGZ2eKfiD8jgUKhytlrStMgcaDZO6OapVU2CBd3IW7nGf1E0BcKVK1vYqOxijfpcrCuwKUVNghT8URpt7QuqmAIhUzRme8tPjeJKtLYg8DZ3vH5e4JRFXqmRRv8sZBXaFWCezG5yq5ojjSmAfk7glEYednkiVrAramCRnFNgVojQ7E5lpWpzvk8cN3uD0YG0+Vc1JNotRj5UrMnBhQB4f6A1ODzYU2ZBFVbJkjQK7Qmg0DGsLLLIYsU/OBtHUPUKP4yKpkUm/j00HcKpnFDvoGAHZo8CuIDUFVlwYGEcoLO3ZIW9Gq+Zsr6IbXAw1BRZ0D01j3BeQtB1H270Ic9AxvQpAgV1BagosmPaHJD9aoMHpQUaaFptLVkjajlRRU2AFADRLPGpvcHpgNupQW0RVsuSOAruCrCuM3OBSL6Q1tHmwpSybquaIRA6ZMbFqSdsqc6CjKlmyRz2kIBUOE9J0Gklv8O6hKXQPTWMbzbOKxmE2IMeUJukCapt7EoPjPtplrBCJFrP+IWOshTF2hjH2PGOMntGSSK/VoDrPLOmIPVY1hxZOxcMYw9oCq6Qf6NTvypLoiH0fgHWc8w0AnAC+lniTyLupKbDgXN+4ZMUXGtq8KMpKx6qcTEneP1XVFFjQ5prAbFCaM/kPOz2odJhQYEuX5P3J0iQU2Dnnf+acB6N/PQagKPEmkXdTU2DF2EwAfaMzor93IBTGmx1DtJ1cAjUFFgTDHG0u8c+NmfGHcPwiVclSEiHn2D8N4E8L/ZAx9ghjrJEx1ujxeAR829QSW0A91yf+Y/nJ7hFMzgZpnlUCscwYKabhjl0cgj8YpmN6FWTRwM4Y288YOzfP131zXvMNAEEAzyx0Hc75Hs75Zs75Zrud/g+yXNV5Zmg1TJIbvKHNA62G4eYKqpojtpIVGZGdxxLMszc4PTDoNLhhFaW3KsWi5U8457vf7eeMsU8AeB+AXVyOVXdVxqjXosJukugG9+K6lTZYjHrR3zvVaTQMa/ItuCBRYL+xLBtGvVb09ybLk2hWzF0A/h7AvZzzaWGaRBYTWUAVd8Q+NDmLc/1jNA0joZoCC5oHxhEWcedx78g0OjxTNA2jMInOsf8MgBnAPsbYKcbYzwVoE1lETaEV7olZuCfEO6P7aLsXnFO6m5RqCqyY8ofQNSTezuOGaJWsHXR8hKIkVImWc14hVENI/Grm7ER0rDaK8p6HnR7YMvSXF2+J+ObuQC2LFjhPtsNONwpt6ZcLqhNloJ2nChS7wcWab+U8UjVna0UOtFQtSTJVuWbotUy09ZVAKIw32oewvSqH0lsVhgK7AlmMepRmZ4g2z948MAHPxCxNw0gsTadBhcMs2tECp3pGMUHprYpEgV2hagqsOCdSymOsKj0toEmvpsCCC/1jouw8bnDG0ltpfl1pKLArVE2hBT3DMxibTv4Z3Q1OD6rzzMi1iDOfTxYWKW7th3tiNunv1eD0oK7YBms6pbcqDQV2hbq8E3EguaP2qdkg3uoaptG6TIi1A3V4yo8zfZTeqlQU2BWqRqQF1GOdQwiEOM2vy8SafDMA4HySj5S4kt5K0zBKRIFdoXJMBuRZjElfQG1wepCu12JzaVZS34fExxxdOE92ZkxDNL11A1VLUiQK7Aq2rjD5RY4POz3YUp4Ng462k8vF2gJLUjNjIumtHtxC6a2KRYFdwWoKrOjwTGLaH1z8xcvQPTSFrqFpbKdqSbJSU2DFpeHkFbdudU3ANT6LHTS/rlgU2BWspsCCMI/kmSdDrGrOjtWOpFyfLE+yN6jF+n0bza8rFgV2BYtt77+QpAyJw04vilekozQ7IynXJ8tTk+Ti1kfavKh0mJBvpWpJSkWBXcHyrUZkZeiTUnTDHwzjzQ4vtldStSS5cZiNsJsNSUl5pGpJ6kCBXcEYY1hXaE1KLntT9wim/CHKX5epyA5U4T/Qj0erJVFgVzYK7ApXU2BF6+AE/MGwoNdtaPNAp2HYUk7VkuRobb4F7e5JwYtbNzi9MOg0uJGqJSkaBXaFqymwIBDiaHMLu4B6uNWDTSVZMFO1JFmqKbAiGOZwDgpb3LqhzYMbVq2gakkKR4Fd4WILqELuRPRMzOLCwDg9jsvYlQVU4abh+kdn0O6epOk3FaDArnAlKzJgMugEPenxCJ3mKHsro/0uZGZMLM2RPtCVjwK7wmk0DGvzhd2BetjpQY4pDWvzLYJdkwjrSr8L94He0OZBnsWISgdVS1I6CuwqsDaaIRESoMhxOByplrSt0g4NbSeXtbUFFrQMTgjS78FQGEfbvNhWSdWS1CChwM4Y+1+MsTPRQtZ/ZowVCNUwEr91hVbMBEK46E28yPH5/nEMT/lpGkYB1hZYMC1QceszfWMY9wVpGkYlEh2x/5BzvoFzXgfgjwD+SYA2kSVaVyjcQtphpxsAsJXOh5E9IXegNjg9YAzYStWSVCGhwM45n/v/qEwAya/XRa5RbjchTacR6Ab3Yn2hFTkmgwAtI8lU6YgVt078A73B6cGGIhuyMtMEaBmRWsJz7Iyx7zHGegA8iHcZsTPGHmGMNTLGGj0eT6JvS+bQazVYk2dO+Gz2cV8ATZdGqLiCQqTpNKjKNSe8A3VsOoBTPaPYQU9pqrFoYGeM7WeMnZvn6z4A4Jx/g3NeDOAZAJ9b6Dqc8z2c882c8812O83jCW1tgRXn+hIrcvxG+xBCYU7l0BQkdrRAIv3+eocXYU5pjmqyaGDnnO/mnK+b5+uFq176XwAeSE4zyWLWFVow7guid2Rm2ddoaPPAZNDhuhKqlqQUNQVWDE354RpffnHrBqcHZqMOdcVULUktEs2KqZzz13sBtCTWHLJciRY55pzjcKsHN5dnQ6+lLFilWJvgDlTOORqcHtxSngMd9btqJNqT/xKdljkD4A4AXxSgTWQZqvPM0GrYshdQO71T6BudwY7V9DiuJGvyLWBs+ZkxHZ5J9I/5qKiGyugS+cecc5p6kQmjXotKh2nZC6iHW6PbyWl+XVFMBh1KszOXPWJvcHoBUL+rDT17qcjaAgvOLXPk1tDmQVlOJopXULUkpVlbsPwjJajf1YkCu4qsK7DCMzEL97hvSf/OFwjhWOcQZUUoVE2BBb0jMxibWVpxa+p39aLAriLL3Yn4VtcwfIEwHSOgULGF86Xmszd2jcAXCNO+BRWiwK4iy82QaHB6kKbT4MYyqpqjRLFTOJfc720epGk1uKmMqmSpDQV2FTEb9SjNzlhycevDTg9uKF2BjLSE1tKJROxmAxxmw5JH7A1ODzaXZlG/qxAFdpWpWWJx64GxGThdk/Q4rnA1S1xAdY370DI4QfPrKkWBXWXWFVjRMzyDsen4FtKORNPddlQ5ktkskmQ1BVa0eybhC8RX3PpytSRKc1QlCuwqE1tA/V1TT1xZEoedkao5VblUNUfJagosCIU5XjjVF1dwb2jzIsdkQHWeWYTWEbHR5JrK1BbbkJ2Zhu++3Ix/fqUZ64ts2FaRg1sqcnBdiQ0G3ZXq86Ewx9F2L+6syaWqOQq3qSQLFqMOf7/3LL754nncUp6DnWsc2FntQL41/R2vDYc5jrZ5cNtqB1XJUikK7CpjTdfjza/twuneURxp8+L1di+ePNyBnx1sR7peixtWrcC2ykign/aHMDYToHlWFXBYjDjxjd04fnEYB5pdONDqRn1LpGjKmnwLdlbbsbPagbriLJzvH8PINPW7mrFEjvtcrs2bN/PGxkbR3zdVjfsCON45jKNtHhxt96LDEymllqbVIBgO4+Q/3g5bBhVYUBPOOTo8k6hvduNAixuN3SMIhTmyMvTItRjRMjiBxn/YTQVVFIYx1sQ537zY62jEngIsRj1uX5uL29fmAgD6R2fwersXR9u9yLemU1BXIcYYKhxmVDjMeHRHOcamA2ho8+BgixsHW924uTybgrqK0YidkBQTDnMwBlpXUSAasRNC5kULpupH6Y6EEKIyFNgJIURlKLATQojKUGAnhBCVocBOCCEqQ4GdEEJUhgI7IYSojCQblBhjHgDdV317JYBLcfxzK4B4DxyP97V0TemuSf2emtekfl/e60o454sf8sM5l8UXAE+cr9uzhGvG9Vq6pqTXpH5PzWtSvwt4zau/5DQVMxrn615awjXjfS1dU7prUr+n5jWp34W95jtIMhUzH8ZYI4/jDASiLtTvqYn6PbnkNGLfI3UDiCSo31MT9XsSyWbETgghRBhyGrETQggRAAV2kTHGJhf5+SHGGM09qgz1e2qSqt9FDeyL/ZJEnajfUxP1u3RoxC4BxtitjLE/zvn7zxhjn5SwSUQE1O+pSYp+Fz2wM8ZMjLF6xthJxthZxth90e+XMsaaGWO/ZIydZ4z9mTGWLnb7SHJQv6cm6ndpSDFi9wG4n3N+HYDbAPwru1J8sRLA45zzGkQ2MDwgQftIclC/pybqdwlIUfOUAfhnxth2AGEAhQByoz+7yDk/Ff1zE4BS8ZsniiDe+aFqlKohIqJ+p36nfo9Ier9LMWJ/EIAdwCbOeR0AF678orNzXheCeottdwNYyxgzMMasAHZJ3SARUL9Tv1O/i9TvUvyHtAJwc84DjLHbAJRI0AZJMMZ0AGY55z2Msd8COAOgDcDb0rZMFNTv1O/U7yL1u2iBPfZLAngGwEuMsUYApwC0iNUGGagB0AEAnPO/A/B3V7+Ac36ryG1KKup3ANTv1O8i97toRwowxmoB/JJzfoMobygzjLHPAvgCgC9xzv8sdXvEQv1O/S51W6Qgdb+LEtil/iWJNKjfUxP1u/ToEDBCCFEZ2nlKCCEqk7TAzhgrZowdjO4uO88Y+2L0+ysYY/sYY23R/82Kfv92xlhTdHdaE2Ns55xrbYp+v50x9ticDQ5EZgTu9+8xxnrozBH5E6rfGWMZjLGXGWMt0ev8i5S/l2Itt6ZeHPX68gFcF/2zGYATwFoAPwDw1ej3vwrg+9E/bwRQEP3zOgB9c651AsAWRDY7/AnAe5LVbvqSVb/fFL3epNS/F32J0+8AMgDcFv1zGoAjdL8voz9E7PgXANwOoBVA/pz/M7TO81oGYAiAIfqaljk/+zCAX0j9H46+ktvvV32fArvCvoTo9+jPfgLgYal/H6V9iTLHzhgrReQT+jiAXM75AABE/9cxzz95AMDbnPNZRLYg9875WW/0e0TmEux3olBC9TtjzAbgHgD1yWyvGiV9gxJjzARgLyKpT+OLTY8zxmoAfB/AHbFvzfMySuWROQH6nSiQUP0e3eD03wAe45x3Jqm5qpXUETtjTI9IJz/DOf999Nsuxlh+9Of5ANxzXl8E4HkAH+ecd0S/3QugaM5liwD0J7PdJDEC9TtRGIH7fQ+ANs75j5PfcvVJZlYMA/A0gGbO+Y/m/OhFAJ+I/vkTiMzFxR67XgbwNc7567EXRx/fJhhjN0Wv+fHYvyHyI1S/E2URst8ZY99F5IyZLyW73WqVtA1KjLGtiKxon0XkuE4A+Doi826/BbASwCUAH+ScDzPG/gHA1xA5JCfmDs65m0VqAv4aQDoiWTGf58lqOEmIwP3+AwAfAVCAyFPaU5zzb4nyi5AlEarfEcmE6UHkTJnYnPvPOOdPJf2XUBHaeUoIISpDO08JIURlKLATQojKUGAnhBCVocBOCCEqQ4GdEEJUhgI7IYSoDAV2QghRmf8P7OUOscBxXskAAAAASUVORK5CYII=\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "seasonal_fluctuations = [data.at[index,'Filled'] - data.at[index,'SAFitFilled'] for index in data.index]\n", + "indexed_data['SeasonalOnly'] = seasonal_fluctuations\n", + "indexed_data['SeasonalOnly'][-30:].plot()" + ] + }, + { + "cell_type": "code", + "execution_count": 276, + "metadata": {}, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "[13.96250778 11.99932062 -3.48316557 21.26924621]\n" + ] + }, + { + "data": { + "text/plain": [ + "[]" + ] + }, + "execution_count": 276, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD8CAYAAACMwORRAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzt3Xl8VOd56PHfqxUBArGITSxiFSMbDLaMF4xtFgM2Hpm4t0ua5tPettdN7k2a3Fw7sZNbN22S2i1J2t40ba/TNOme1DcO1rCYxThesGMswAbMSOybxCIWgQAhtLz3j3dGCDFa58x5zznzfD8fPjaj0TnP0cM8Ouc973kfpbVGCCFEsGTYDkAIIYTzpLgLIUQASXEXQogAkuIuhBABJMVdCCECSIq7EEIEkBR3IYQIICnuQggRQFLchRAigLJs7HTkyJG6uLjYxq6FEMK3tm/fflZrXdib91op7sXFxVRWVtrYtRBC+JZS6mhv3yvDMkIIEUBS3IUQIoCkuAshRABJcRdCiACS4i6EEAFkZbaMEEJ42eqdNazaUE1tfSPjCvJ4ZlkJK+cW2Q6rT6S491IQki36TvKeflbvrOG5V3bT2NwKQE19I8+9shvAV7mXYZleiCe7pr4RzY1kr95ZYzs0kUKS9/S0akN1e2GPa2xuZdWGaksR9Y8U914ISrJF30je01NtfWOfXvcqKe69EJRki76RvKencQV5fXrdq6S490JQki36RvKenp5ZVkJeduZNr+VlZ/LMshJLEfWPFPdeCEqyRd9I3tPTyrlFvPDkLIoK8lBAUUEeLzw5y1c3U0Fmy0DTZWi6BEPGdfmWeFJl1kSANNZD81XJe7q5XAet12Fo9zlcObfI93lWWmvXd1pWVqY9syrkq5+DaAV8rhIGj7IdjXDLTz4Fx96Dz2+HvGG2oxFu+fHjcO4gfO4DyB1sO5o+U0pt11qX9ea96T0s09psCvu1i7D5T2xHI9zS1AD7N8HVc/DGn9mORrjl8hk48g401MLb37EdTcqld3E/8rYp7OPmwof/Cic8cjUhUmv/JmhtgqIy+OAf4NQe2xEJN1StBTRMuAfe+xtzBh9gSRd3pdQApdQ2pdRHSqmPlVL+OQXeWwHZg+BT/w8Gj4F1T0Nbm+2oRKpFIzCoEH7zpzBgKKz/ClgYnhQui1bA8Cnwa/8Mmbmw4au2I0opJ87cm4BFWus7gDnAcqXUvQ5sN7XaWs1v8umPwKCRsPQbULsTdv6L7chEKjVfg/0bYeYKk/dFfwRH34GPX7EdmUilxgtw+C0IhSF/DDz0Zdj3GuzbaDuylEm6uGvjcuyv2bE/3j8NOr4NrpyB0nLz91m/ChPvg9f/xPxDEMF06A24ftl8yAHu+h0YMxs2/hFcv2I1NJFC+zZAWwuEnjB/v+czMGI6vPYstDTZjS1FHBlzV0plKqU+BM4Am7TW7zux3ZSKRiAzB6YvNX9XCh79C1PY33jBbmwidaIRyB0KxQ+av2dkwmOr4FJNWtxkS1vRCAwpMvfXALJyYPmLcP4g/PJv7caWIo4Ud611q9Z6DjAemKeUur3ze5RSTymlKpVSlXV1dU7stv+0Nsmeughy82+8PnY2lP0ufPADuckWRK3NUL0OSh41H+64iffC7F+Hd78X+Jtsaen6FTiwGWY+DhkdSt70JVDyGLy5Ci6dtBdfijg6W0ZrXQ/8Alie4Gsvaa3LtNZlhYWFTu62705+CBeP3bg072jh12BAAaz/stxkC5oj75grs0R5f+RPzZVcwG+ypaX9m6DlWuK8L/szM1yz6Xn340oxJ2bLFCqlCmL/nwcsAaqS3W5KRSOgMs1v7c4GDofFfwRHt8Ken7kfm0idaASyB5orts7S5CZbWopGYOBImHT/rV8bPhnu/zzs/k849kv3Y0shJ87cxwJvKKV2AR9gxtzXOLDd1IlGoPgBVlc3Mv/FLUx+di3zX9xyY53uO38bxt5hbrI1Xe5+W8If2tqgag1MW8Lqjy8kzvs9nw38Tba009JkbqbOfIzVH51KnPcFXzLj8eueNrPoAsKJ2TK7tNZztdaztda3a63/1InAUuZMFZzdx0f5D3bdiCEjEx5dFXuS7du2IxZOOPEBXD7NBwMXdJ33NLjJlnYOvQnXG3g35/6u854zyEyFPrUbdvyT7Ygdk35PqEYjADxfXdx9I4aJ98Adn4R3g/8kW1qIVkBmDl/dM677vN90k63WQqDCUdFXIXcIX905vPu83/YkTHoAXv8GXD1vIVDnpWFxr4AJ97DrYuI1uW9qxLDkTyBrgDzB6Hdam7xPeZgDF1XCt9yU9wDfZEsrrS1QtQ5mLOPoxZaEb2nPu1Lw6J/DtfrArDeUXsX9whE4tQtC4d41YsgfDQ9/BQ5sMjfahD+d2gX1x3qf9+GTYf4fwu6X4ei7LgUpHHfsXWg8D6Hy3uV9zO1w9+9D5Q/NEI3PpVdxjw3JMPPx3jdiuOczMLLE3GRrvuZSoMJR0QioDChZ0fu8P/AlGDIe1n05UDfZ0ko0All5MG1x7/P+8HOxqdD+v1pPv+I+ZhYMn9z7biuZ2eZy7cIReO97NqIWyYpGYNJ8GDSi93nPGQjLvgmnd8P2H/W4i9U7axLPxBB2tLWZvE9bDDmDep/3gcNh8fO9ngrt5bynT7OOhlPwnRJY+L/hoWf6/v0//S3Yvxk+XwlDx3f5ttU7a6Rzj5fU7YPv321mP93zVN++V2v4pzCc3gOf32E++Ams3lnDc6/svumGXV52pi9bswXG8Q/gh0vgEy/BHb/et+9ta4UfLDRdm7pp6mEj79KsI5Gq2NT7RE+p9cbSb5k1wCu7PouLJzvhdCthR7TC/Df0eN+/N36TrfEC7PzXLt+2akN19zMxhPuiFZCRDTOW9f17MzJh2QtmKnQ3q4V6Pe/pU9z3VpgHVAr72dx42CRzaR8ft0/A68lOS9EIjL+7216p3Rp9m1k1spu83zTTphevixSLrx015SHIK+jfNibdDwWTfJ339CjuV8+bdUVCYXM21l+hcjhbDXWJi7XXk5126o+ZdYT6e7UWFyqHE9u6XFyqVzMxhHtO74ELh5PLu1Lm+w/9wnRrS8DreU+P4l69HnTrjbXb+yt+ad/Fb3OvJzvtRGNDcTP7MSTTUfzfTVXiVTV6PRNDuCMaARSUrEhuO6FyaL1uFh5LwOt5T4/iHo3A0Akwdk5y2xkyzlzid1HcvZ7stBONwOjbYcTU5LZTWAIjZ9wYv++k1zMxhDuiETOsMjjJ1WfH323ab+59NeGXvZ73LNsBpFxTAxzcAnf/XnJDMnGhsHly8cJRMw7fQTypMlvGAxpOw7H3zLxlJ4TC8M5fwZVzMGjELV9eObdI8uwFZw/Amb2w/M+T31ZGhrla//Df4fpVMz22Ey/nPfhn7vs3mlkuyY67xsW308Ul+sq5RWx9dhGHX1zB1mcXeTbxgVcd63TvZN51K+xb78z2RGokMzsqkVAYmq+aE0SfCX5xj0Zg0CiYcI8z2xs+BUbP6vYuuvCAaASGT4VRIWe2N3YODJ0oefe6aATG3dntsyh9Mmk+5A3zZd6DXdybr5nGCzNXmLmrTgmFzcL+Daed26ZwTsdO904MxcGN2RMHt5ihPuE99cehdodzV2tgnlAvecxMymi57tx2XRDs4n5wCzRfcTbZENue7nJoRlhW/ZpZ1THZ2VGdhcJm9sS+Dc5uVzijaq35b8jpvJdD00U48paz202xYBf3aAQGDIXiBc5ud1TIXPL78FItLbR3ur/T2e1OmGeG+CTv3hStgFGlMHKas9ud8jDkDPZd3oNb3OOd7md06nTvhPgl+pG3A7Owf2A0XYaDrzs7JBOXkWlu1O3fBM3yYJqnXD5jlmd2+iodIHsATF9qrgx8tEJocIv7kXfMwvtOX5rHlZabS3+5RPeWA5u77nTvhFDYDPUdfCM12xf9U70OR2dHdRYKw5U6XzXRTrq4K6UmKKXeUEpFlVIfK6W+4ERgSeuu070Txt1pLv27eLBFWBLvdD/xvtRsv3iBGeqTvHtLNALDis1Da6kwfSlk5vpqaMaJM/cW4H9prUPAvcD/UEqVOrDd/ot3up/+CGSn6NH/+NDMgdfNUICwr73TvcOzozpqnz2xzgz9Cfsa600j7FQMxcXlDjZrw0cjvmnikXRx11qf1FrviP1/AxAF7D65c2IbXD7t/F3zzkJh84DUgcRrTwiXHfoFXG9I3aV5XKjcLCZ15O3U7kf0zr4N0Nbszuf90gmo3Zna/TjE0TF3pVQxMBd438nt9lk0Apk55lIqlSbeZ4YAfHSpFmjRCsgdApMfSu1+pi6E7EGSd6+IVkD+WCjqVQ+L/puxHDKyfDMk51hxV0oNBn4GfFFrfSnB159SSlUqpSrr6uqc2u2t2jvdL4QBQ1K3HzCX/jNXmDMH6a9qV3un++XOz47qLDvPDPlF1/hq9kQgXb9ihkZnPm7WgkmlgcPNPZe9Fb4YmnHkp6GUysYU9n/TWidsXaK1fklrXaa1LissTHK1tu6c/Ki9070rQuVw/bIZEhD2HN1qOt2nanZUZ6EwXDkDx7e5sz+R2IHXoaXRxc97GM4fhLoqd/aXBCdmyyjgh0BUa/3d5ENKUjQCKtPc9HLD5AfNUIBcotsV73Q/dbE7+5u+1Az9Sd7tilZA3nCzBowbZj4OKF/k3Ykz9/nAp4FFSqkPY39cqqwJRCNQPD/hsqwpkZVjhgKq15qhAeG+9tlRSxIuy5oSA4aYabY+mj0ROO2zox6DTJdWL88fDRPvNUMzHufEbJl3tNZKaz1baz0n9medE8H1WV21aYOX6rvmnZWWm8Wqjr7j7n6FUbMdGk66n/dQGC7GWvkJ9x1+C5ou2cn76d1w/pC7++2jYD2hGr9Umplke62+mrrYDAn44FItkKKv9r/TfTJmPGqGACXvdkQrICffrP3ipnjbxqi3Fw4MWHGvSK7TfX/lDDRDAtE1ZohAuKe90/3D5slRNw0aYYYApbi7r7XFrPUyYxlk5bq772GTYOwdns97cIr7haNmpozbl2hxoXK4fApqKu3sP12d3gMXjrg3W6KzUDmc3WeGBIV7jr0HV8/ZzfuJbXDppJ3990Jwint8bXWn2mv11YxlZmigi2a6IkWiEVAZ7g/FxbVfonv/BlugRCOQNQCmLbGz//hJpId7OgSnuEcjpv3d8Cms3lnD/Be3MPnZtcx/cQurd9akfv8DhpqhAZk94a69FWYa3KCRdvY/ZCyMn+eL2ROB0dZmPmdTF5s1X2wonAEjSzz9Sz0Yxb3htFmKMxRm9c4anntlNzX1jWigpr6R517Z7U6BD4Wh/iic2p36fQk4ux/qou2X5lZ+qYPZ/6ldZnhIpF7tDmiobX9gzWrej2yFK+fc2V8fBaO4V60BNJSWs2pDNY3NNz8S3tjcyqoNLoyJljxmhgg8fqMlMDrMjrL+Sx08P3siMKIVZo2XGcvs5r20HHQr7Fuf+n31QzCKezQCI6ZB4Uxq6xN3yOnqdUcNLoSJ93v6Ui1QohEouguGjrf7S334ZBgzS36puyE+O2ryg5A3zG7ex8yGgomeHZLzf3G/et4svRpby3lcQeL127t63XGl5Wbdibp97uwvXbV3ujeX5lZ/qYOJ4/j70HDKnf2lqzN7zcNDsaslq3lXyuT90Btw7Za1Eq3zf3Hft8G0u4sl+5llJeRl39yoIS87k2eWlbgTT3zWRpWcxaVU++wok3frv9RDYUB7evZEIOytABSUmM+ZJ/Leeh32b3Rnf33g/+IerYAh49s73a+cW8QLT86iqCAPBRQV5PHCk7NYOdel/iFDx5uhAo9eqgVGNAKjboMRUwEP/FIvnGmGBmVoJrWiEdNHIX804IG8j58Hg0d7cijWpdV2UqTpslnys+x3b2qvtXJukXvFPJFQGDZ/3QwdFEywF0dQxTvdP/SV9pfi+V61oZra+kbGFeTxzLIS9/4dxC/Rt/61GSocONyd/aaTcwfhzMew7IX2l6znPSPDPOvw0U+guTF1bT37wd/F/cAm0+bO1lNqXQmVm+IejcB9/912NMHTRad7T/xSf+e7UL0e5n7KXhxBFb8q6vSgoifyXvlDOLjF3sN0Cfh7WKa90/29tiO52YipZshALtFTY28FDJ8Co2+zHcnNxs01Q4SS99SIVsDYOWaGipcUPwADCjyXd/8W9+Zrqe90n4xQ2Kx/cfmM7UiCpbEeDqe4031/KWXiOrgFmhpsRxMsF0+YpZ3d6rTVF5nZ5hmX6nXQ2mw7mnb+Le6HfmHa23kx2dBh9sRa25EES/vsKI/mvbTcDBXu32Q7kmCJf468nPdrF80a8x7h3+IejUDuUCh+0HYkiY2+zQwdePAuuq9FKyB/XPvsKM+ZcA8MKvTcJbrvRSNmRtLI6bYjSWzKQsge5Km8+7O4t7aYtnYlLnS676/4Jfrht0yXJpG8eKf7UDj1ne77KyPTDBXu32iGDkXyrpw1DdC9NnGio+wBMGOpucJoa+35/S7w6CekB0e3moLp5WSDuYRsazFDCSJ5Bza72+m+v0JhM2R46A3bkQRD1VrQbd4dkokLheHKGTi+zXYkgF+Le7QCsge61+m+v8bdaYYQPHSp5mvRCAwcYR5i8bLiB80S0JJ3Z0QjUDDJrN/jZdOXQmauZ4ZiHSnuSql/VEqdUUrtcWJ73WprM6vvTXOx031/ZWSY3+YHNpsHrkT/xTvdl7jY6b6/snJMf1WPzZ7wpWsXzeQJL86O6iw3H6Yu8kxPB6fO3H8MLHdoW92rqTTt7Lx+iRYXCkPLNVPgRf8detN0ui99wnYkvRMKm6HDI+/YjsTf9m2AtmZ/fd4vHoeTH9qOxJnirrV+CzjvxLZ6FK2Idbpf6srukjbxPjOUIJfoyYlWQO4Qs9SrH0xdZIYOJe/JiVbA4DGm8b0flDwKKtMTeXdtzF0p9ZRSqlIpVVlXV9e/jcTXcp660P1O9/2VmWWGEvZtMEMLou9aW8wQh41O9/2VM9AMHVatMUOJou+uX4X9m81yA16dHdXZwOHmidW9FdaHZlz7iWmtX9Jal2mtywoLC/u3kVO77Xa676/SJ+B6gxlaEH1nu9N9f5U+AZdPw4kPbEfiTwdf98fsqM5Ky+HcfqhzoWFIN3zy6zAm3um+5DHbkfTN5AfNkEL0VduR+FO0ArLy7HW676/pSyEzxzOzJ3wnGoG8YaYBup/MfBxQ1odmfFPcV++s4fn3WvmH5uXM/94u95rgOiEr1wwpVK0zQwwx1hr7+sjqHcc5s+1nbGi6jfnffd9fP6MBQ2DKw6a4d7hEl7z37NXth2nYFeHlhlnMX/W2v35G+WNgwrxbfqm7nXenpkL+B/AeUKKUOqGU+j0nthsXb4L7zw138c2W33K3Ca5TQmFoPA/H3gWw29jXJ1bvrOE/fr6aUZxjfes8f/6MQmGoPwandgGS995YvbOGNat/Sj5XWd/m47yf2gXnDwN28u7UbJlPaq3Haq2ztdbjtdY/dGK7cVab4Dpl2hIztBC7VAvEMaXYqg3VLNTv06wz2dI2F/Dhz6hkhRlKlLz3msn7L7msB7C17XbAhz+j+H2CWNtFG3n3xbCM9ebHTsgZBNMWmw95W1swjinFauuvsjxjG++23cYlBnV43Uc/o0EjzJhxrLhL3nt2qv4KSzO380bbHJq4sXaUr35Gw4phzGyrefdFcbfeBNcpoXJoOAk124NzTCm0YMgZijNOs75t3k2v++5nFCqHuiqo2yd574XlQw4zUl1ifWsA8n78fWg4ZSXvviju1pvgOmXGMvMAVrQiOMeUQl+dvJ82rdjUelf7a778GcXbwknee+VL46tp0tn8om1O+2u+/BnFh2aiESt590VxXzm3iBeenEVRQR4KKCrI44UnZ9ntm9gfeQUw5SGIRlg5Z1wwjimFZl54k/Mj72JAwRh//4yGjDNPWEYjwfm3nCpaM/XsG5wf+wDDCob5+2dUWAIjplvLu9IWnqIqKyvTlZWVru/XE7b/GCJfgM9shTG3247Gu84dhO/dCctfhHs/azua5G39a9j0PHxxt/d6gHpJzXb4wSJY+Xcw5zdtR5O81/8U3vkreOaAeXo1SUqp7Vrrst681xdn7oFSsgLzgIM82NKt+M/Hb08ndmVmfGhmjd04vC4agYwsmOHOOoQpFwqDboXq9a7vWoq72wYXwqT7rT+95nnRiFkPf+h425E4Y8RUGH27/FLvjtZmTZbiBY6c5XrC2DkwdIKVvEtxtyFUDmf2wtkDtiPxpnin+6CctceFyuHYL6HhtO1IvOlMFM4fDFbe4+02D26BpgZXdy3F3Yb47IkqOXtPKD504Zc1vHsrFAa06f8rbhWNAOrGEFZQhMqh9brpq+siKe42DB1vhhz2yiV6QtEIjCqFkdNsR+KsUSEYPlWG5LoSjcDEeyF/tO1InDVhHgwa5XrepbjbEgpD7Q6oP247Em+5XGfW3wnSpXlc/BL98FumS5O44fwhOL07mHnPyISZK2DfRmh27ylbKe62xIccquQS/SbV62Kd7gP4IQeT97YWqH7NdiTeEj+rDdqQTFwoDM1X4OAbru1SirstI6eZoQe5RL9ZNGLW5Rgd0GcAiu6EIUWS986iERh7BwybZDuS1CheYLrHuZh3Ke42hcJmCOJyP9sOBk1jfazTfbn3O933V/vsideh6bLtaLzhUq3pVhXUqzWArBzTZKh6HbQ2u7JLKe42hcJmCEJmTxj7N/qr031/hcLQcg0ObLIdiTfEhyZDT9iNI9VCYbhWD0fedmV3UtxtGn07DJssl+hx0QrIHwtFd/X8Xj+beB8MHCl5j4tWwMgSKJxhO5LUmroIHvkGFM50ZXdS3G2KX6IfetMMSaSzeKf7mT7qdN9f7bMnNkDzNdvR2HXlHBzZGuwhmbjsPJj/h2YhORdkubIX0bVQObz7f8yQxOxfsx2Na1bvrGHVhmpq6xsZV5DHd2cd4x4/drrvr1A57PgnOPymWQo6TXTO+/dm7uFO3Zo+eXdRwE+RfKDoLjMUsfdV25G4JlE/ydPvv0xTToH/Ot331+QHIXdIWq01kyjvl3b+nCsDi8xMGeEopxpkL1dKVSulDiilnnVim2kjI8OctRx4Ha5fsR2NKzr3k8ymhYfVdja23gWZaXIxmZVjVj6sWgetLbajcUXnvA/mKvexi4qmu4I7O8qipIu7UioT+D7wKFAKfFIpVZrsdtNKKAwtjabAp4HOfSPvz/iYIaqRnzfeaSkiS0rLofE8HN1qOxJXdM77wowPyVUt/OzqXEsRBZsTZ+7zgANa60Na6+vAT4CAz2ly2MT7IW942sye6Nw3cnnGNhp0Hofy77YUkSVTF0NWXvrmPXMbZ3QBp4bMthRRsDlR3IuAjguknIi9dhOl1FNKqUqlVGVdnTy0c5PMLJj5GOx7DVqu244m5Tr2k8ygjaWZlbyp7+SLywP6VGpXcgbC9CVQtQba2mxHk3Id857LdRZmfMTr+m6eXh6yHFkwOVHcEw2W3dK7T2v9kta6TGtdVlhY6MBuAyZUDk2XzOyJgOvYT3KeqmaEamDUPf/Ffz0ynRAqh4aTUBP8tpMd8/5Qxi4GqiYmzv/19My7C5wo7ieACR3+Ph6odWC76WXKw5CTnzazJ1bOLWLrs4v4yYN1kDWAeUvSZxroTWYsg4zstMv7S3efhAEFzF+y0nZIgeVEcf8AmK6UmqyUygF+A0iPf6lOyso1H/SqtdDW2vP7g6CtzYw3T10MuYNtR2PHgKHmF3s0YtrMpYPWZrPGSsljkJltO5rASrq4a61bgM8BG4Ao8J9a64+T3W5aCoXh6jk49p7tSNxRuxMu1ZhZI+ksFIYLR+D0HtuRuOPwW3DtouQ9xRyZ5661Xqe1nqG1nqq1/pYT20xL05ZA1oD06dAUrYh1uk+fJzQTmrkCVEbazJohGoHsQTBloe1IAk2eUPWS3MGmwEcjwZ89obUp7sULIG+Y7WjsGjTSPJmbDr/U21rN0OOMpZA9wHY0gSbF3WtCYWioNUMWQXZmr2mtJpfmRigMdVE4u992JKl1fBtcOSNrybhAirvXzFhmhiqCPnsiqJ3u+2vmCvPfoA/NRCOQmQvTl9qOJPCkuHtN3jCzqFS0ItizJ6IRs6754FG2I/GGoePNInJBLu5ax2ZHLYLcfNvRBJ4Udy8KlZshizN7bUeSGucOmpkhcml+s1A51O6A+uM9v9ePTn4IF49J3l0ixd2LZq4AVHDP4uLHFZIhmZvEi17VGrtxpEo0AioTSh61HUlakOLuRYNHmSGLIBf3sXOgYKLtSLxlxFQYdVuw8z55AQwcbjuStCDF3atCYTN0ce6g7UicdbHGrKMis2QSC4Xh6Ltw+YztSJx1pgrO7pMhGRdJcfeq+IcgaGdx7Z3upbgnVFoO6Bs/p6CQ2VGuk+LuVQUTYNzc4BX3aIXp/j5yuu1IvGlUKQyfEsy8T5gH+WNsR5I2pLh7WShshjAu1tiOxBlXzpquQ3Jp3jWlzM/n8JvQWG87GmdcOAKndsnVmsukuHtZ/MMQlEv06nWg2+RD3pNQObS1wL4NtiNxhsyOskKKu5eNnG6GMILytGo0AgWTYMws25F427g7IX9csPI+ZjYMK7YdSVqR4u51oXIzlHHlrO1IknPtIhz6hRlykE733cvIMD+nA5vh+hXb0SSn4RQcf1+u1iyQ4u51obAZyqheZzuS5OzbCK3X5UPeW6EwtFwzBd7P4g9kyX0W10lx97oxs8xQht9nT0QrYPAYGH+37Uj8YeJ9MHCE//O+twJGzoBRM21HknakuHtdfPbEoV+YoQ0/un7VnIGGHjdDDqJnmVmmDd2+DdDSZDua/rl6Ho68I2ftlsgnzQ9KnzBDGvs22o6kW6t31jD/xS1MfnYt81/cwuqdsSmcB7dA81X5kPdV6RPQdAkOvWk7km51mffq9aBbJe+WSHH3g6IyM6QRfdV2JF1avbOG517ZTU19IxqoqW/kuVd2mw96NGKWMp4033aY/jL5Qcgd4ulZMz3mfehEs46QcF1SxV0p9atKqY+VUm1KqTKnghKdZGSYIY0Dr5shDg9ataGaxubWm14p5sFxAAATCUlEQVRrbG7lL1/bY87gpNN932XlmuYtVWuhtcV2NAl1lffvv/ahuWKT2VHWJHvmvgd4EnjLgVhEd0JhM7RxcIvtSBKqrW9M+Hpxw3ZouiizZPorFIbG83DsXduRJNRV3ksa3oPWJhmSsSip4q61jmqtq50KRnRj0gNmaMOjl+jjCvISvv6JvJ2QMximPOxqPIExbQlk5Xl21kyXeR+wAwaNMuvJCCtcG3NXSj2llKpUSlXW1dW5tdvgyMyCkhVQ/Rq0XLcdzS2eWVZCXnbmTa8NylYsz9pu+mVKp/v+yRkE0xZDdA20tdmO5haJ8j40u5WH1A7TdCYjs4vvFKnWY3FXSm1WSu1J8OeJvuxIa/2S1rpMa11WWFjY/4jTWShshjiOeG8UbOXcIl54chZFBXkooKggj79/qJkBTefk0jxZoXJoqDUt+DwmUd7/730XyWptlDX7Lcvq6Q1a6yVuBCJ6YcrDZogjGjGX6x6zcm4RK+cW3Xhh/bOxTveP2AsqCGYsg4ws2PsqjPfevIVb8v7zz8KAoVC8wF5QQqZC+kr2gBuzJ9pae36/TfFO99MWS6f7ZOUVwOSHzM9Ta9vRdK+12SyVIbOjrEt2KuQnlFIngPuAtUqpgKxR6mGhMFypg2O/tB1J92p3wqUTMiTjlNJyuHAYTn9sO5LuHXkHrtVL3j0g2dkyP9daj9da52qtR2utlzkVmOjCtEfMUIdHZ0+0i3e6n7HcdiTBULICUP7Ie/YgmLrIdiRpT4Zl/CZ3cGz2hIcv0bU2Uzal071zBhfCpPu9Xdzb2swqkNMfgezEUySFe6S4+1EobIY8anfajiSxuio4d0AuzZ0WCsOZj+HcQduRJHZiG1w+LXn3CCnufjRjuZk94dEHmqTTfYrEi6aX856ZY55rENZJcfejgcPNNLO9Fd4cmolWwIR7pNO904aONy34vDg0Ex+Km7oIBgyxHY1Airt/hcJw/qAZAvGS84fh1G65NE+VUBhqtsPFE7YjudnJj6D+mOTdQ6S4+9XMx/Hk7AnpdJ9a8QXYomvsxtFZ++yoR21HImKkuPtV/miYeK8ZmvGSaATG3iGd7lNl5DQoDHnzl3rxfBg0wnYkIkaKu5+FwnB6N5w/ZDsS49JJM2NCLs0d0WWHo9JyswTwZY8swFdXDWerZVlnj5Hi7gNdfsjjs1G8cone3ulePuTJ6rbDUSgMus085u8F8dk7MjvKU6S4e1y3H/Jhk8wQiFcu0aOxTveFJbYj8b2uOhyt2lANo283w15emRIZjcD4eTBkrO1IRAdS3D2u2w85mLPkE9vgUq2F6Dq4cg6ObJWzdod01eGotr7RtK0LhU3j7MZ6lyPr5MJRM1NGhuI8R4q7x3X7IYcbxbRqrUsRdWGfdLp3UlcdjtpfD5VDWzPs3+hiVAm0D8VJ3r1GirvH9fghL5wBI0vsX6K3d7q/w24cAZGow1FedibPLIsNeRWVQf5Y+3nfWwFjZsHwyXbjELeQ4u5xPX7IwZw1HdlqhkZsuHZJOt07LFGHoxeenHWjKUZGhrmBuX8zXL9iJ8iGU3D8fRmK86geOzEJu+If5lUbqqmtb2RcQR7PLCu5ufNNaTm8/W0ze+LOT7sf5P6N0HpdLs0ddkuHo85CYfjgB3DgdTst7arWAlry7lFS3H2gxw/5mNlQMNEMjdgo7tEIDB5t1pMR7pk0H/KGmZ+/jeIejcCIaVA40/19ix7JsEwQKGUujQ+9YYZI3NTcCPs3xTrdyz8nV2VmmSYe+16Dluvu7vvqeTjytgzFeZh8GoMiFDZDI27Pnji4BZqvyKW5LaXl0HQJDr/p7n73vQZtLTLe7mFS3INi/DwzNOL2A03RCAwokE73tkx+CHLy3Z81E43AkPEwbq67+xW9lmyD7FVKqSql1C6l1M+VUgVOBSb6qH32xCYzVOIG6XRvX/YAmLHU3Nxsa+35/U5oumxu4sqQjKcle+a+Cbhdaz0b2Ac8l3xIot9CYTNEcnCLO/s78jZcuyhDMraFwnD1HBx7z5397d8IrU2Sd49LqrhrrTdqrVtif/0lMD75kES/FT9ghkjcGprZWxHrdL/Qnf2JxKY9AlkD3Fv+ORqBQYVmyWnhWU6Ouf8usL6rLyqlnlJKVSqlKuvqPLJUadBkZptZK9XrzJBJKrW1mqEA6XRvX+5gmLrYFN22ttTuq/maOXOfuQIyMnt+v7Cmx+KulNqslNqT4M8THd7zNaAF+LeutqO1fklrXaa1LissLHQmenGrUNgMlRx+K7X7Ob4NrpyxM79a3CoUhoZaqN2R2v0cegOuX5YhGR/o8SEmrfWS7r6ulPpt4HFgsdZe7NacZqYshJzB5ixu2uLU7Uc63XtLyXLIyIotv1uWuv1EI5A7FIofTN0+hCOSnS2zHPgKUK61vupMSCIp2QPMUEkqZ09obT7kUxdBbn5q9iH6Jm+YmY4arTD5SYX22VHLISsnNfsQjkl2zP1vgHxgk1LqQ6XU3zsQk0hWqNwMmRzflprtn/wQLkqne88pLTctF8/sTc32j26Fxgvy4JJPJLW2jNZ6mlOBCAdNfwQyc81Z3KT7kt7c6p01Ny1c9o8T1lOiMs38duEdJStgzZfMVdXo25LeXOe8/8uYnzIle6C5YhOeJ0+oBlFuvvkARiNJX6InavOXtW8NZ0bcDQOHOxOvcEb+aDM90YGpsJ3zXlt/hfzDr1Ez8gHIGZh8rCLlpLgHVSgMF4+bIZQkdG7zN02dYKqq5Z/rZycboUiFUBhO74FzB5PaTOe8z1UHKFT1/OBs8lcEwh1S3IOq5FFQmUmfxXVu87c84wMAXr4sHZc8KX4fxOm8Z37AdZ3JKw1S3P1CintQDRwOkxeYpxaTGJrp3OZveeYHVLbNIKtgXLIRilQomAhj5yRd3G/Ou+bRjG280zaL/IIRycUnXCPFPchCYTi3H+qq+72Jjm3+xqsz3J5xhNeZd3ObP+EtoTDUVMLFmn5vomPeb1NHmZBRxxZ1j+TdR6S4B9nMxwGV1Flcx16ej8aGZOYu/XT3naGEXfGpilVr+72JjnlfnrmNVjK4/zHJu58oGw+VlpWV6crKStf3m5Z+uNQsAfyZtx3a1lX4zDvJb0uk1vfvMYt7/c4ab21LJEUptV1r3atHkOXMPehC5XBqF1w4ktx2pNO9v4TC5qGjK2eT207dPqirkgfWfEiKe9CFHjf/TXbuc1XsrE2Kuz+EwqDbzHIByYh3eJr5ePIxCVdJcQ+6YcUwZnbyxT0agRHToVBuqPnCmNlm5owTeS8qg6Ey1u43UtzTQajcDKk0nOrf9186CYel072vKGXyfvANswR0f5zdbx6CkyEZX5Ling7ia65X9fOG2OY/No0Z7vy0czGJ1AuVQ1sz7NvYv+/f8DXTfPuOTzobl3CFFPd0UFgCI2f0rw3bsV/Crp/CfZ+D4VOcj02kzvi7YfDoG+PmfbFvA+zfAA9/xaxZI3xHinu6CIXhyDtw9Xzvv6etFdY9DUOK4MGnUxebSI2MDHMj9MBmuN6HdgstTfDas+Yey7w/SF18IqWkuKeLUBh0K1R32eb2Vtt/DKd2w9JvQM6glIUmUqi03DybcPD13n/Pe98368I/+ufSlMPHpLini7FzYGgfZk9cPQ9bvmG6+9z2ZGpjE6kzab7p0tTbvF+qhbe+bc74U9mmUaScFPd0oZQ5ez+4BZoaen7/lm/CtUvm7E1myPhXZrZpqlL9GrRc7/n9m56HthZY+s3UxyZSSop7OgmFobUJ9vcwe+LkLtj+I7j79x3p6CMsC4Wh6SIceav79x19F3a/DPO/AMMnuxObSBkp7ulkwjwYNKr7S3StYf2XzaX8wufci02kzpSFkDO4+7y3tcK6L8OQ8fDA/3QvNpEySRV3pdQ3lFK7Ys2xNyqlZJFvL8vINMsR7NsIzdcSv2f3y3DsPVj8x6bAC//LHgDTl5pVIttaE79n+4/g9G5Y9k1poxcQyZ65r9Jaz9ZazwHWAM87EJNIpVAYmq+YsffOmhpg4x/BuLkwVx5YCpRQGK7UmecWOrt63txjKV4ApSvdj02kRFLFXWt9qcNfBwHurx8s+qZ4AQwYmvgS/a1VcPkUPPZtM0daBMf0RyAzN3Het3wjdvP8L+TmeYAk/QlWSn1LKXUc+BRy5u597bMn1kFr843Xz+6H9/4W5nwKxvdquWjhJ7n5MHWRKe4dezic/AgqfwTz/huMLrUXn3Bcj8VdKbVZKbUnwZ8nALTWX9NaTwD+DfhcN9t5SilVqZSqrKurc+4IRN+FyuFavXliFWI3Ub8C2Xmw5Os2IxOpVFoOl05A7Q7z93jeB46Ah+XmedD0WNy11ku01rcn+PNqp7f+O/Ar3WznJa11mda6rLCwMNm4RTKmLoTsQTcu0avXmScYH34WBo+yG5tInRnLQWXeyHv85vmSP4a8AruxCcclO1tmeoe/lgNVyYUjXJGdZ8Zgq9bA9Svw2nNQOBPmPWU7MpFKA4fD5AVmAbmON8/n/JbtyEQKJDvm/mJsiGYXsBT4ggMxCTeEwnD5NLz8X6H+qHkSNTPbdlQi1UJhOH8QXvkDuXkecFnJfLPWusthGOFxM5ZBZo5Z1rX0CZjysO2IhBtmPg5rn4bqtXLzPODkV3a6ys2HqYshKw+Wfst2NMIt+WNg4r2QO0RungdcUmfuwuce/y5cPgMFE2xHIty08m/NmLvcPA80Ke7pbMg480ekF+molRZkWEYIIQJIirsQQgSQFHchhAggGXMPmNU7a1i1oZra+kbGFeTxzLISVs4tsh2WSDHJu+hMinuArN5Zw3Ov7Kax2azZXVPfyHOv7AaQD3qASd5FIjIsEyCrNlS3f8DjGptbWbWh2lJEwg2Sd5GInLkHSG19Y69fl8v44JC8i0TkzD1AxhXk9er1+GV8TX0jmhuX8at31rgQpXCa5F0kIsU9QJ5ZVkJeduZNr+VlZ/LMspKbXpPL+GCRvItEZFgmQOKX1z1ddvflMl54n+RdJCLFPWBWzi3qcQx1XEEeNQk+0F1d3gvvk7yLzmRYJg319jJeBIvkPb3ImXsa6u1lvAgWyXt6UbpjJ3SXlJWV6crKStf3K4QQfqaU2q617lWHFRmWEUKIAJLiLoQQASTFXQghAkiKuxBCBJAUdyGECCArs2WUUnXA0X5++0jgrIPheEHQjiloxwPBO6agHQ8E75gSHc8krXVhb77ZSnFPhlKqsrdTgfwiaMcUtOOB4B1T0I4HgndMyR6PDMsIIUQASXEXQogA8mNxf8l2ACkQtGMK2vFA8I4paMcDwTumpI7Hd2PuQggheubHM3chhBA98FRxV0qVKKU+7PDnklLqix2+/rRSSiulRsb+XqyUauzw/r+3F31iXR2TUurrSqmaDq8/1uF7nlNKHVBKVSulltmMv7O+Ho+fcxT72udjefhYKfUXHb7HdzmKfe2W4/FzjpRSP+3w2hGl1Icdvsd3OerqePqVI621J/8AmcApzLxOgAnABsz8+JGx14qBPbZj7c8xAV8Hnk7wnlLgIyAXmAwcBDJtx57E8fg5RwuBzUBu7GujfJ6jro7Htznq9Pp3gOf9nKNujqfPOfLUmXsni4GDWuv4w05/CXwZ8PNNgs7HlMgTwE+01k1a68PAAWCeK9H1XW+Ox286HtNngRe11k0AWuszsff4NUddHY/f3PLvTimlgF8D/iP2kl9zBCQ8nj7zcnH/DWIHppQqB2q01h8leN9kpdROpdSbSqkFrkbYd+3HFPM5pdQupdQ/KqWGxV4rAo53eM+J2Gte1JvjAf/maAawQCn1fiz2u2Ov+zVHXR0P+DdHcQuA01rr/bG/+zVHcZ2PB/qYI08Wd6VUDlAOvKyUGgh8DXg+wVtPAhO11nOBLwH/rpQa4l6kvdfxmGIv/R0wFZiDOY7vxN+a4Ns9d7XSh+Pxc46ygGHAvcAzwH/Gzqj8mqOujsfPOYr7JDcXSL/mKK7z8fQ5R54s7sCjwA6t9WlMwZgMfKSUOgKMB3YopcbELrnOAWitt2PG1WZYirknHY8JrfVprXWr1roN+AE3LhlPYO4vxI0Hal2NtHd6dTx+zhEmF69oYxvQhlnvw5c5oovj8XmOUEplAU8CP+3wPr/mKOHx9CdHXi3u7b+1tNa7tdajtNbFWutiTNLu1FqfUkoVKqUyAZRSU4DpwCFbQffgpt/ESqmxHb72CWBP7P8rgN9QSuUqpSZjjmmba1H2Xq+Ox885AlYDiwCUUjOAHMxCTr7MEV0cj89zBLAEqNJan+jwml9zBAmOp185sn2nOMGd44HAOWBoF18/wo3ZMr8CfIy5K74DCNuOv7fHBPwLsBvYhfmHOLbD176G+c1cDTxqO/5kjsfnOcoB/hXzi2oHsMjnOUp4PH7OUez1HwOfSfB+3+Woq+PpT47kCVUhhAggrw7LCCGESIIUdyGECCAp7kIIEUBS3IUQIoCkuAshRABJcRdCiACS4i6EEAEkxV0IIQLo/wNnpHkQXIIGTwAAAABJRU5ErkJggg==\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "import numpy as np\n", + "from scipy.optimize import curve_fit\n", + "import math\n", + "\n", + "def signal_triangle(x, a, b, c, h):\n", + " v = np.divide(x - h,b)\n", + " return a * np.abs(v - np.floor(v + 0.5)) + c\n", + "\n", + "x_data = np.array(data.index)\n", + "y_data = np.array(seasonal_fluctuations)\n", + "\n", + "params, _ = curve_fit(signal_triangle, x_data, y_data, p0 = [12, 12, -3, 21])\n", + "\n", + "fitA = params[0]\n", + "fitB = params[1]\n", + "fitC = params[2]\n", + "fitH = params[3]\n", + "\n", + "print(params)\n", + "\n", + "model_seasonal_fluctuations = signal_triangle(x_data, fitA, fitB, fitC, fitH)\n", + "\n", + "plt.plot(x_data[-30:], seasonal_fluctuations[-30:], 'o')\n", + "plt.plot(x_data[-30:], model_seasonal_fluctuations[-30:], '-')" ] } ], -- 2.18.1