{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Test of Resampling Methods" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Utility Functions" ] }, { "cell_type": "code", "execution_count": 38, "metadata": {}, "outputs": [], "source": [ "import pandas as pd\n", "import numpy as np\n", "import neurokit2 as nk\n", "import matplotlib.pyplot as plt\n", "\n", "%matplotlib inline" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Original Signal" ] }, { "cell_type": "code", "execution_count": 21, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX8AAAD8CAYAAACfF6SlAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJztvXuwdclVH/br87yP7zXvl2Y0EhpjD8RB+IssBXAoA0bYKQ2hMBHGMBDwVJkidsCEiJAQgVMVTB44lOXAWAQLgxFIAjMhUilCFnlgaaQR6IGQRjNIMDOa9+t73HvPY5/T+aN379279+rutfqce+/33bNX1dR8997TZ+3du3ut3/qttXorrTU66aSTTjrZLOkd9wV00kknnXRy9NIZ/0466aSTDZTO+HfSSSedbKB0xr+TTjrpZAOlM/6ddNJJJxsonfHvpJNOOtlA6Yx/J5100skGSmf8O+mkk042UDrj30knnXSygTI47gsIyfXXX6/vvPPO476MTjrppJOrSj7+8Y8/r7W+IfW5K9b433nnnXjooYeO+zI66aSTTq4qUUr9OedzHe3TSSeddLKB0hn/TjrppJMNlM74d9JJJ51soHTGv5NOOulkA6Uz/p100kknGyid8e+kk0462UDpjH8nnXTSyQbKiTf+n3riZXz4T18QjVksNd798SdwYX8uGvfZpy7i//n8c6Ixy6XGb//RE3j+8lQ07tFnL+ODn31GNEZrjQc++SSevjARjXvshX2879NPicYAwHs//RQef3FfNOapCwf4nU98SazrA3/yDL74/J5ozPOXp3j3x5+A9FWmH3r4WXz+mUuiMRcO5viNjz2G5VKm6w8efR6fefKCaMzetMCvf/QxzBdL0biPfvFFfOLxl0VjJvMF/vWDj2FaLETjPvH4y/jYn70oGjNfLPHrH30M+7NCNO6Pv3QB/+7R50VjFkuN3/zY47hwILMBjzxzCR96+FnRGK013vNxuQ1YVU688f87/+JBfOe/+AgKwUb4wJ88jR991yfxc7/3eZGuH3jHQ/ie//2josX54S+8gB/+jU/iv//dPxHp+s9//Y/w/e94CC/uzdhjPvXEBfyDX/8j/Nh7PiXS9WPv+ST+/q/9ociQ/+lzl/GDv/aH+MFf+0ORrv/udz6Df/jOT+CzT11kj3n6wgR/71cewnf/0oMiXT/zvs/hR9/1STz05y+xx1zYn+P7fvlj+LZ//u9Euv7X33sE/9V7Po3f/zzfMEzmC3zX2x/E3/r5/0+k6+3/7xfx47/1afwfn3ySPWax1PiOX/wwvvVtfyDS9c6PPob/+rc/jV/9yGOicd/6tj/A3/6FD4uc4e9+6kn8+G99Gv/8Q38q0vUdv/hh/J23P4hZwbcBv//ws/ix93wKP/O+z4l0ffcvfRTf98sfw6UJ32l89Isv4h+965P4b//NH4t0rSon2vhP5gtcnhpD/OTLfLT7+WcuAzAolCtaa3zpZfP5P3uebyQ/97RBkFLUao3jo89eZo95uNT18NN8wwoAH/mCQWiPPMtHu597ynz201+SoVaLmiTI+uHys0+8xH9eAPDgF18Q63r0OTPfdl1x5Y/LeXjkGf7zesxxtguBkbT3I1lTz1ys94cEvHy+XH8SYOAaxuf3+GjX7pUnBfuyWCyxPzNRyRMv8a/R7i+JDQCAp8t5lMy91SW1AavKiTb+z12qF9bjggdvH/jTF/kL0w0PHxNshKdLXZcExsRFMBJdT5V0z3zBNyQuJfL4i/yN8PRFGbVkZTzoA5AZ8mccGksS4SkoALL7co3kRQG625+b5ytZhy4994xgPl8oDarEILvPSzL3z5Z7RGJYn3LuSzL3VteXBNf3nEOliNZUqcu1ISlx14NoTZU6rJM6KjnRxt/dMBJ6xC7OFwWo5BnHUby0z9f1ZKWLP8Zd0C+J7sssyIsHczbPfXFSO6UXBLqedhDTZM5b1PPFskLUL1yWPy8AeJnJ0Wqtq40tm0NnTQmu0a6Pl/b4DsM1/pL1YY3ki4KclatLMvfPXpKv3ydfrtfGCwKe2+5nyf5yn9cLov1c6hKtefd5He59rUNOuPGvH8DLgom1RuFlwUZ1HY3kIdrFf+FgzuY/c3U9X27qYqmxx0QZzzq6JHP47CV37nnz6Ca8RM/rsvwaL00LHJRO6cWMtQHw575YLKt7y3bygmvMMVzPXcqb+9pw8feKOwdcZ52rq/G8JPv5knWguWuDr8s660uTQpykX0VOtPF3qZgXBQ/eos9LU/7DcHVJNt2lEllrzacR3CokiVFwuVbuNb7cmEO+rstOxMC9RtdJSDbdpYaujDkUPa96HNepXZ4WsIGW5Hm564F7X8ViWTl2yfO6lKELqNe9TFf9vCQVdRcOzLiX9mbsyNXVJXFqF8v7msyXOGACpaYu+RxKx60qJ9r4X56aiVRKtulcw8Ut9XITgJLN446TGBMrIoM8XfW+BBs8Y9xeOUYpmUFuOBrmuL1ZrUviaHLm47JzX7kOlGu49qaLSpfE2OU8r2KxxGS+RE+Z9cRNSueuqcvTOXrKRK7cHNnlTKd2aVJAKdk1NveybE1Jda1DTrjxNxvhptNbouqMS9MCt57dAsA3yNZw3XJ2q7Fpk9c4cXQJDfItZ7eqzc4dJ70vey+3nN2q5pM77rZz2wBqFJUSu6FvOSN/XlYX26lNal17Al3ufUkdqFiX87y4CPlSCXhuPrOFvdmCb5AnBa7bHWE06LGfl117N58x18itEro4mWM06OH6UyP2mreOxurizqOd+1vPbolq9vemBW4pdXHXonU0t57dEtFZlx1dlwS2Y1U52cZ/UmB31MeprQF7YU6LBWbFEjedlS1o+9BuPD2uUCV3XKWLucjswr9RaiQdXdxrdHVJDdeNZ8YA+FUMTV18R7Pn6DpgJpcvZ+q6tOp9zRbs3M7lSYFrdkcY9hX2mfdl7+UmoUG+PC1wamuA3VGffV+XZ/V9ubqT4yYFTo8H2B0P+Gu+vKYbhcb/0rTAeNDD2Z0R+74WS42D+UKsy66pm85uicplL0+KSpe0gW0VOdHGf88u6PGAjVotIrzptGxB700L7Iz6OL01ZC+WabHAbOGgGSG3eOPpsXiRSVFataBPj9kb1Y6r5pBrTDIdaEMX9znb+zpjdHE5ZIuQ+z0lBgY3CR3UpWmBU+MBdkZ8I+neF8B3UJcntS7u3F/27os9blrg9NYA28M+f2149yXZz6crp5Y3h9w1dWliHM257aHY0dS6OuO/Frk8K7A7HuDUuC/ePDcLkf/laqMKkFO1eWS6akczYC9M62hukqK0CrWORVHG5UkhXtD1pjNRhsQgX396hJ6SzaHVpbUsYji9NcTuqC8ABk00zp6P0nCdGg8yjKSMsrCO5tR4gH2xAxWi8UkNyrjJ1D1Pl2Q/nxoPsCOYQ18X16ldKp3aznjAngv73dJ9uQ452cbfQTPsTTCpjR0gQK3O5pEgoIYuwabblery7ovt1CYFtod9nNkyaIZjkK2juf7UGErJNipgkP9SA1NmO74xXEPsZjxn6abbK5/z7phPJdrCg5uEEV4TUMidGgC+IZ9Yw9UXU4JSR2NonyF2Rnxd/vOS6Dq1lQcAc5z1qbGQOvPXYUf7rEf2XDQjRU6nMxBGiWYkYWJDlyRstiG61KlVupjGblbfV7HULINs7//M9hA7gtB+b1o6mu0hAN4Gtzkaa7i4xs5e4w2nZdHJpfI57wiQv6UobhTqsmtqR7CmXOoMkFExxnDJ94rVxZ37i5N5mV/gRxl7vi4uRepSZ5n3JXLWW7J96d9Xh/zXJBYhm42ai5D5i3N3VCInKdKV8rPTms6aLzTrwCqr6/TWAFvDHn/zOGiGe412Dm24LXFquyWqBnjGxG4Wa7j4xm6OrWEPZ7aMo+GMs47GIn951CVHkrvjXL5anqwU75VMemR/tsDuqH9kUYZdv7m6RM9rNMDuWID8y+++/tQYPdVx/muTHCrGcr/X7crokUvlRj01GmBaLFlnzNhjD3bHA4wGPdHi3B33KyPJWTD2yN3tYd8YSQEy2R33sSPQZe9jZ9QXceOXpwucGvcrR8PZ4PZ6tktjwt90i9KI98vvSY+z/PT2sG+oGC5qnRWloymdGuMai8USs8XSAIqRoGDBpxKZ87E/qx2vlLKw+THuc57MF9ge9UVRxiXPIHP3pdFVghBh9CSlLSeFua+dkYmSJaDMRkId7bMmOZjVD2MyX7Jqni2tcXprgGFfsTfPwXzhGcn0uMnc6NoaGIPHXZwGOZnFAvCMZKVrKDOS+7MFdkbGgXJ12TncGvbLcJvpeGcFdkYO8meh8aYurlPzdXE2nTuHEkCxX+qSzOGkvK/xoGf4agF4GfV7OFtSZ5z50FpjMl9iPOiJ8gs1UBoBkBnk8cCuQ/4YALi21MV1hpP5AluDHnZHfcwWS5ZBPqhAWR87w75QlwEGAG8+9l1AIWAN1iFrMf5KqTcqpR5WSj2qlHpL5HPfrpTSSqnz69CbkmmxxNagj61hr/yZY5DNZ8bDnqjEblpuHoku+xmri2tMpsUS42EP40oXP8rYGvZEyH9a+PfF12Xui4/87X1tDfvl9wjuqzJcAl2DHrbKU0SnjGof+7y2hj1sC1BrvTbsfTF0Vc+rj+0Rn/M3hrW+L46u2aJ2aqYsmrvmF1AKOF1SZ9wD/Cblc94Z8kGZ/e6dUR+jQY/98phJsayAAcAzyK6T3xr2MeHqmi+r/QXwoi4XvGwN++wih3XIysZfKdUH8DYA3wLgbgDfqZS6m/jcaQD/AIDsjRsryLRY5BsTG9ozN/ikMGim3nRyXdyyN4swJMbEXs940Mf2qM8ubZzMF2ZhCoxJZfwH0s2zaDhrlpGsHKi5RrZRsPdV6eI8L2ejDnoshwHUBmg86FU/c8YYXcbxcu/LONB+BQwk92WcRg/TYsms6jJObdhX6CmermVJh0gR8tS/RoYuwHGG5V4RgZdynEyX2V+AiS5FugZ9tgNdh6wD+b8OwKNa6y9orWcA3gngHuJz/xjAzwLIO+hdKIulxnyhM4xJvcjG5UbgyLT0+pVBZkUZlvYxKJ6tyyL/gTyi2SrHcRf0rNzgY8Ec1kayNFxMXfV9STZqPYdbwx7LAFW6Bj2MK6MgcWq90qlx18aiqSvDgYqN3YB/X26UMR72oXUdDaR0bQ37UEqZ+RDsL9fxsp5zscCgpzDo99i67Hc3nbzsGsfDngy8OPuSsxb9+eCu33XIOoz/bQAed35+ovxdJUqp1wK4XWv9u2vQxxKXUhkLUatSwEiwyLTWJsoYOOhOsMG3MlDr2EH+HMMwKRxdQzlClhjkadE0JkeB/O01ipF/5kYdC5H/WGjspo4DHQ96mC2WrGMhjLHroddTGPV5xsRF/jWg4M2H/TyXsmiCEMGamte6xsMea23MF4ZSkuqy3z3qGycqjfC2BIDCfre1U1cb8lfE76pVqpTqAfg5AP8o+UVK3aeUekgp9dBzz8lehO6LjwgB/oPfGhg0w0X+84XGUlukm4FaLcIQIuS8KMMYLilClvHVXvR0FPclRE459+UarhzkP+r3oBRTl0tnCdbUtAQGAMrnLHOg0vmwn99i6po0gEGvumbOuFqX1NHIkf9oYBwod19qbeissUPvcdZ9BSgGNsq4upD/EwBud35+BQD3zdGnAXwlgN9XSv0ZgNcDeIBK+mqt79dan9dan7/hhhtWuiiXCx6LFvSyWihcJFnpEqLWSbHAsK/Q7ym2rmqRDaR8dRNhyJE/f0GvFGU0oic5P8tOOvr3JeKCzbjFUvNKektEKAEUdSLbvUae4arWLxuNNyMagG+4ajTOc4aTRkTDz4+5yJ/7nJsRjQz5V7rYNsC5L2F+QSlg2FclvXd1If+PAbhLKfUqpdQIwJsBPGD/qLW+oLW+Xmt9p9b6TgAfAfAmrfVDa9AdlIkXNru/i4+rEQbX6+cuaGvs7FgZ9eAuaF54qZQ1kvkImafLp0fykL9s0/UrwypJVg76PQx6iums3VxGv/G7qC7XmHC5cWu4MgoWxs6a4qJqAOLn3NgrXOTvOjVJFV7pQAH+XqkBjxz55+rKyTFapuGqq/bRWhcAfgjA+wF8FsBvaq0/o5T6aaXUm1b9/lxponE+jeA++DEzvMxF/rYyo9YloVSEyL80dgZ98nQViyWKpc6oivETo3loXFYxJUVc9dwbg8wz4kCJ/AWURcOYMHndiopZAflz59531gA/ByJH426ORvacRwPhfTWAAT8/luOsm9GTMMqwz4vpQNclg3V8idb6vQDe6/3uJwOf/fp16ExJo+JE8OAbIR8XOREoTRxeMqtiakpFVss99REhM6EHNBEhl87q9xSGfRN1zRcai6VGv0elh4yYpHntoEbMhKpbwupSOPZ6w/dWz/14wKvoqBqv3LkXIn9uVVezYioP+XOjLtqB8gy5bZLjril6r0hBWQ/PX5ag8TrvJweAzLXhAB5J0cd0vqzWEjeRvS45sR2+JEpjbvDxCshf/OCd/IIsmereF2/T+bpS9Ihr/Ac9U8vN1uUgJ/Nd8XubLZbQGg5ClnHjDcqCPfc18uciQvt5Uemri8aZz5lCyIeF/Bs5mlzOn1kV00j4CvMtbn5BkouTIv+Jh/xz8lxGP2/c2FkbVxXtc6UKyc9yKxhcNC6qOOkJk8tNzpS7WKwuUeNV0dS11OZdqKnrM7qktdwLx4HyjInraKxOcWUG05gslxqzhVc6KKERBjk0QiZfPZDnkeS6nMSoqNJq4VBnzJxVI5EtiFw96oxnxN3INQ/5S3M0bo6RS+PWeT8eKFuXnFjj7yJ/US13wyAzuWCC8+ca8jq53GclK13kL6vlXjQWmf1dVJdDc9hx3PJLH/mnxtUJTtmm81Ga/V1Ul8MFAwLU6lRmjAXGxEX+7AiPbIYSVvsIq9UaupjGteGsBUBpPJSVYDd5eHkJay7y55ZfZucyXOQ/NKBsvuiM/0qSj/zdUjke8ne9vq3l5vLV9UblbQQX+dtxOffF0uWgNKuLu3m2hMjfRbqVLibl1tKVGFc7azfC4ztQW7IJpO+rKBuNxg6vKy+/FCJ/caVKO6LhIuQqypAmshsGWV6IIU7CZnP+fcyKdIOda29sBZkU+Uvmfh1yYo1/g/OXcIteM4lNVkZ1OcjfGgYuWvDReNJwOcgfkNRyO0aBWY1AI38mFzz0oowU8vfQuCTcdjluIO3kXaMASHh4QlfivlyjYMdJe0C4yN8mzauoS5hfaKBxdpScmcgeyA8mzHHW9vOiPhVPF5A+7sLtATH/ZzpeD/lzr3EdcmKNv1uZIarldgwyN1Hscv7m//xN10bIKV1N5M+v5XbQDDNZSSF/KRWTi/zZqLVoVrcAecifS++1dKVyGYRR4B9n4CHChC6bNB8PZVFG8zwrGQ8/HsicdaPrVnRkyLIBQgpGg51bCVaDMnl+wb3usC5/X/JLt10gx9G1Ljmxxt/fdNxablP/Xdfdmu9iGpMMJNkyklw07t4Xs8lrS2iQfeQ/PkQ07id8uXw1icZzOP9MXUlgQCB//nEGHiJkRk/uHPLLgHtlUp+3Dhdl0nzL2SuTeTpn5UZdw77iH3cxbzZEAuky2ylBkeZU+7jXHdblReRcXQ3kzy99XYecXOPvbTqDuHicpJuEBRihvRPKmnF8vrpFjwgRhmhBi+/Lc6CCiiQxD+8gQvN/eSJbmstoHhR2OLkMCvlzKUEp8nc7Ws3/c3I0zCjDM3a2ASuVrHTpLMlxFy7y5+YKWvQeu0R0SazfVMFCc/1yj7tocv78RPE65OQa/3l9OidgFuc8ESb6nCkfITc33bCf1mW/V5rw9RHGsN9jHb/brMmWIX97jZw5tN+bW4Fjx7Hn0KUDmCG6j/xHEl3SKiZqDrm0j5egT96XD0LKtcFB4/b6THNeOlnpgxBr/FNr0Z5LZYUTCVWnczpr3vyeV6rcaLATIv/qvpKAwouSmU2KUwf51w60M/4ridstCpRGMvEAi6WG1mi0kQNy5D9i6ALMRvF15Ww6jq75CroqdMe8L5eH5xqFiYecRgOeU3ORP3+j1gnOSpcwaS7W5VzjlOVoFi1EmNLlgxAuGndzGYB5zjkgBEDSsc2KJYb9utObkxj16Ufu3NseEGsDOE2bbqc5UN9X0qktPMotA/lboMpZi+uQE2v83bAZ4BkTO+n2gduHwUH+/fJFE1YXZ5HNFstGZGKuW84tcozkrHB09YW6hMh/5mwe7oImkT/HgRa1U+M6mhohy6IntzGsWhvCiikbZaTQ+Gyhq/tRyqDxtANtgpAaIacNl9UF8J5zLvKfL5bVdQG8fAsFDIyudMQwdnSNGM+5fl4yx1vZjl4ddXFzOy2n1iH/1cT13oBF/vENZxf80DPIqY3g0hwAj0ZYlFGGryu9OJuhLJcecTe4HVssZXw1J3oCyg1ur485hz5fzUX+rjEZMh2NjeTc+WDfl+dA2fflOA3N6K6eFYtKB8BzhiHknzZcumGQOc/ZByFcJ+87mmFfMSKTJjAYlZEDZz8PXV2M50yBEPNdaV2DnkKvPL+K40Dt2wbbujrjv5LMHC4YYCL/8u8j72Fwxo095JRaZHYxjQaeMWGGzdUiYxtkXYXbXCPph7JceoSKMtLz4W9w/ot0fKfGQZ8NXQMeNz4vdHU/vZ7CoKfYz7lFgzHGDQc1PcKLXMs1JQQv88WyMqhcXXNvr0g4f9epjRhUjJ2rHF0NiokBlGYe8uc/L7lTq9dhR/usVWZ+eNnvYZY6XMwustbmYSCMftP4p8bYhduOMtava7E0jWqjvs8FJxZ0aUxcZD1jtJ67Brla0Mn5aOpiU0wE8p+nEKG9L8fJa41kM18LSTKMiR9NctHdfLFsI3+2QW46+ZRx9Q0XB1DU67cGIQDPSDbWb1+J59CuY841NiKaAcdZWxtQI3jzXWnb4e9L8Rwyiz7WJSfW+Lc2z4DjiXX1WUAQyrbC5vQiqx2Nh8aTPKY8RK82z0CG/OeLJfo9VR3FbOiRNI9pNri3edibrnYa7ER2qctea3KjBjYdhw8e+giZO/fC+WgZE0Yk5BtJSSTUBhTc9euhcUZEI80v+Lk4O5dpp6FbTk3urHkU04zYl2lw5UU0TIppXXKCjb9ubR6+ARIi5FbCrJ8dNqdQq49mOMnlmW9YK6PAiTJkdMByqVEstXyjFm0judRIdnH6RlISbvvUVDJiINB4MqIJFBFwqMRhy0jKdLETvq05lACKHF31mpLoEjuaov28DtOB+tRvWlc72gU62mdloQwXH2Goxv95oay7oPnhpW8kU2WAZkE798UJmwMIQxw2MwyrT2cN+j30FG8OlUIVZUgSbUMP3bErM+zc23CbUT3iOnlOoriV25FwyC3wwo1opM9ZtwAF14FKDVc74ctB42360X5XSpc8ovEMMpdi8qPCjCijq/Nfk+SgmVkAzcjL1/LLSlMJ35axE3GL0rDZN0D9Kn8QGwOgVWmVrJgq59DWZHOMidamWmLkbXBu1FWhOzYN1qb30sZk0dDFLoksdDZtKadi/Mo4Dm2pq8+6upLgJSNy9akzbvllCwD2JYnsEoQM+HulxfkzC0ykYHNdcmKNf5uKYSwyj8dk1/hmUEyhhG96cbZRGseQuLosGuckfH1d5hrC43yUZsel594z4nYjMHTJk5VNw1VvcDnlxkX+fqUVh2JqVPusktthGS43ShbQlkLw0qr2YQADH5RJ8iat/IK03FuQH2tFNLl0Vof8VxOK8+eGl36pJ4dWkVdm2BJAT5eYiukl0bhfwmrHsZDToElnAfFqBD+iAfghsE/TuddO6vKQE2AMBFdXFWUwQnutdWm4mnw123C1NriseoTTiJZNxbSSy3LaUtLk1U74pqqzmkZSUjElncNWlMykYrKoM7+artdx/muRnGSl3ZBSDq6F0jhUjGckBz1zwmFOTXbqGv3yteoahZwpJxLyQ/RKl7QEkKMr4Gg4iUCpLtuUJY66gkgyHWW06CzhmhpXa0OYN2E6UPtZQFIZ184j8aMnmaOh5jCNxn3Ov6SzkvfVbspL9Y74EU2vx+vkXpecWOPfCtFZiUA6RM8xXMVSR9/+4xtJ08KfPvelXSfNQMi5RnLRjmjcayd1EVEGK9FGoLSULt8A2XHS6qwhi2LKj57cclnLJYvLL1eIJln0iBghe6iVDZQohCzj/PnVPu3yS05BBdCms6RNeRzHSwIlxppal5xY49+u8+cvsnrz2I0qr/E1v5cZyXGfwUkGeHgOGs+pssjWJaYsNPwjMlK6/BDdXmO69DXjvoo28uck9lsJaQbFtFhqLHXbgeZTMenekVY9vDA/ln+8Q0Z+TJLLIEBIDI37eRPbyS0tjhABCqGdWpecXOPvG8m+4eB4aNyMU0plIWRObTCFxoeDHqtByS9tdL8vNMbXxadiZDw8ZSTNHKZ5XdIgM+ZQmkAMPa/YOFsGKu3wbdW1M6pH/HXI1TVfLNFrlMuWuli9I828CZvOGnhonNGLMMqswMk5LNAvIkidq5S7V6imPCCew6Oec4f81yAUDw8A88hhZqRBZraf+40rAO/Bt84DYTUaEbqkjoZTs08kHd3vI8dQRpJpTPzqlpSuUGWRvC+Dr6tluIRVIJz7mlbPS5hHCuRNYlSHLRbwDZC0GSr7eIdBmhv3czsiNE7ZgKjjbZbLutcYEx+UDVngJS+PtC45scaf4vwBbrJSaEwyUCsV8nETiHIqhlrQvHI+n3oAuAvaO1Ar03DJ6Sxeh6/4eRWELk5Fx9ruK91N7JfLjkuKSYw+WTkacy2DnpCHJ2if1LlK9bEr8nyLFLysEnWNM3WNhCW965ITa/xbDUqMxTnzSj0BbjmfFlMx/uFiVldOnT8Qr0agjAmvnK+pyy5ujjFpNA0N0hUMoY3KSS6LETKRCEzpyg3RqZJj99pjunxqKtXh63PcQ0ZymXpeEqrTlsvapHZsDpflEcZ07wiHipFWkAVyO6wouWmQOSePkpVx0jXFuK91yQk2/qENnub8c9C4n7Rxv4+SKbWgOQg5xC1KF7SgHt4KJ5QNGUlOclnqrP2ko9XLK81tz2Fsg9dUjHuNaafm8+m5yWV+RLMTJ8HIAAAgAElEQVQ6dTbsm3OVomjcy2XYa4yujSUd0aSvsZz7niwS8pu8RLqkvQhetU8uoOCs33XJiTT+JI/JQf4Br5/TRg4k0DiVrOQ6GoLHlCZ8uVUWUuqMzi8cTnLZr5O241hHGHuvEnSvPTQGaCeXcynBeETTzpvk9GX0Gb0jZBkwcz7c6wPS8xHi0+21x3S5L0oB0ucq2TflUWicE2VIz1UyjZ7NV2ECPKah4/zXKH4lAsAre7PVEgPvYfAQsmyR0VSMXBcnbCapGKau3Dp/6YIORTTiHE3GHPJQGo2Q5T0FHAfaTi4P++lObj964lSrUbXmnDJFv1wWSD9nGlzxjvGQ6rIVPVRxhDTqYifbB+0oObqmAtVqHe2zgtAoLX0ut48UgDRqtYeLiatishFyqKcg7NRC6I5z3o48v9BsNAIE0VNmcllembGeXokc3tl2cksbyliVKkV7/aaeM9lvwpyPka+L62gIZxijO/1yWTNOZZVRun8LjXPLZStdEXDlv4/b1RWzNxRINZVxV9F5/kqpNyqlHlZKPaqUegvx9x9RSv2JUupTSqkPKqVeuQ69IQmVAALpZI+/oFMPnjxcTGBMBj3Zg1+piiknyiC6FnMqpsScvyC57FdLcBqU5FQMhZDTndw+5287uXPyJu51hMa5z8uO4/DO0gjPp+mAtDMk+zKYtI+7vypdjCjDr84C0sURZJQRGeO/j9uM4Sfbr9o6f6VUH8DbAHwLgLsBfKdS6m7vY38E4LzW+i8DeDeAn11Vb0zoSe03/hYa1+IxEw8+FjbHjYnpqrTVEuYa44bLRBn+2T7pLuQQ15pb5x/lTInQPofz59VJE9ETkxt3HYZ1wNJIjWu4fGMyTsxHLg9PRa5cg5yTH6OiZI4DzUnCtu6LqUtcHEE6mn60V4I+XTZdZuuXywIGYF1NCd/XAXhUa/0FrfUMwDsB3ON+QGv9Ia31fvnjRwC8Yg16g0J6fcZZ2X6dtBmXgdIYCIOKMlLVIyTCYBwXQDUNsc5s90PZyijE8yaAvFpivZy/Trbwu2OUUoYekVIxTITcfs4JgxxI0HN0UcafVXFCUTEZunKiwrSuNuc/7MfpLBrwMGwASWfFKSYahPCQv1sua74jvS/XJesw/rcBeNz5+YnydyH5fgDvW4PeoFBonF/B0AxlU56Y2qjcGt9W2JwM0du9AbzN0zzC2I47jOQyVS2RU1nEOVcpxo2nGrZIJMngZ+XlqLTh4jxn8VlHJG2ZokfCFTjxXIGmo+TovBO6OPeVQTFVXcENSpAR/RPPK0kxVWujXfGXztHI7mudMljDdyjid+QuUkr9XQDnAfxHgb/fB+A+ALjjjjuyLyhkgNy/UTIl0EzKcFEIg5PEonjMlC4aEXIimoBRSPCYSx3QxdjgVLWP1rrhgBrjvPngVKqQpXIONTUOrO5wpUokogk05Vld4WtchUps89VphNwGFDldpkldRZN+tNcrPfqDVa0WctYsg1yXXw6ZyD87b0KsjWR+gbABVxPt8wSA252fXwHgSf9DSqlvBPATAN6ktZ5SX6S1vl9rfV5rff6GG27IvqDQCYyA3Egmw2YKYTBr1KUhep0wy4lo2puHRWcRRiG1oN0jjN3rTSXO15FAZDnDRfO1hfY7YsifbMpjc+Pt+8o56heQJytTRQQr1fn7uYxkErZ99AenWi2Y8GVRMUTTm3BfJp1aZA5TeT8pdbZOWYfx/xiAu5RSr1JKjQC8GcAD7geUUq8F8Iswhv/ZNeiMSqj5B2DwfUI0Hq3MyAjROcc00KWeckcTO1CLip5MpUqckwyF6O71+2Lr112UBpT8Z8Io2M/VupiJfR+1JvItdDcxn9dt6Eo851A3sfk+WZRhksucOZRFNKG9wkHI/vuCgVRJpJyKieUyxJx/Ti5jBV2pdw6sS1Y2/lrrAsAPAXg/gM8C+E2t9WeUUj+tlHpT+bH/EcApAO9SSn1CKfVA4OvWInQFDscTB7w+69wRaYhOoxlWck5awRBIZMcO1KLekmX1FYnohJpD+zdSFxFlAGY+ioQB8pvyuMifui9p7T239DU3CUtVj+QYE1aOhsybRHpHKIOciJ4ohMw5TpsCFMN+vDKOomN5zZd0jobF+RNlpWIQ0ovvr3XKOjh/aK3fC+C93u9+0vn3N65DD1eya5eLNvWQOlCLrF1moHESzfSZdICzyOzxttKIxuVaB31qTJs6s7pTRrJFqSQ2AvW8rO6ciMb+jRKqKc/qkjYN8apiqOccPwAtFmWk1i/Vp3JxIs1ZcdA4EdFkgBcuGj/lJXCSoIys+LO64pEQTdPJ5pCL/Kl1aM9VcqnTw5B10D5XnMRKylIofjRoUw+cJi/SKCRQEIVmYi38VNu/HZfKFVC67N8ooaInjq7QggbCfHXI0SSTlYSxS206qibb/pw6zgBo1mRzcztyKjFc1SVNVubQlrnVasO+EncTc6vV/OecqsKL2wBZ5JoCZXS3frpazT9dFqij36NI+p5I40++JIGF0toVDHwqph7HOd42hlpTCJmmEfh17a6u0EagNipgDXLKqbU3qnv91PUBofuKzeGCrJbg6ZI70HZNtnXysSiDrorhlQ9TiVEhbZlMwkYo0kRxBPWceQULspwVXS6b3l/2c76unFxGDJRRa4p7rhJViOF+52HKiTT+VPkahwsO8dWxFn4quWz1STl/e41hhEwbSQ490jbi8WRlrkGmqnZSm44yQHZcvASQrhByr799feFcRqrOP6QrNIdUU54dx6rAcY4wzkbIzCSstMmLqvPn9i9kVRZlFmLIK8HkaDwElFLFERQwsJFlqilyHXKijb/7EDnH2+YgZIqftbrlHb7xTUcll43uRAVOJAkbQndUuaz5OV3t06bOmJz/WhyoNSYJZ52Ry6CiIPO3ECJs0zd2HCcR2GuUy+blF4apJCxhJHMdDdsgE81rySomAhhEQRlBxdhqNWkSNpXDo3QBZcFC5LWxodJcAEeS9N0Y41+FYYnNE97gcr5a7GiYukY+r8ugpnIdTU45X5vmSEU0IUeTPpGyveHiKI1Kzpmf0ygt6ECFzpqDCENGIWTIqXdYmHGp5HI4Z5XqrpZGalGKKXf9BoxrCFBwmsOCdGwiIpfuFapiikODrUtOpPGnOj8BRiladIPz+T6AR4+EUGvoGmPll7m60rQP0Xi1dqeWyfkHeiXs30hdRFOeHRfTNSWjpziNUFMPOQn6gK7UHAoTvrPFgmjK4xlkms5i9I44a7HfU+gpORWTihhCNiD13uVQRBPTFa1Wi1GJkefc0T6ZEqRiWMYkxFdL0V0aBYXQnZSHTy0ysst0hQUtNlxsOkveVi9F46vkMvIdqNCpRYxd2tG0x6VLjulILaQrVC5r906RSoz2ZGsqZiRTTl5+1pG8T6WKaDL6dlpRPMPxrktOpvEPoiB5PTwb3QmTPf6xwkC9ecSOJgONZxvJJJIkDvxK6Yps1HgSlqDpEg6U6p6tdCXOl6HKDc11SPMmaW7c75VIP69AtJsRPaXyY6FyWY6R9F/HCKQdlP+WLEAAKIhzepJRBmE33O9s62ofW2F/XjeVuE452cZfiFpXSfhKK3Bi4WWqRr19Lk3iyNkMp0aVy9pxUiPJTsL615isvV9vRJM+hCuA0kJzGLgvHnW2niqmYT/+MnYK8KTyYyFKkAMofF0AEyFLI6EgoEg32FElrPY6YrrGreNJOJVxHe2zVqEacoD4plsuw52fAAdxyY6CpkPZw6Is6G5MIE1ZtKKTjFA2nYTNizLyIppQPkjel1Gjz5SjaSNCaV9G9cKZjHJZ9+/UOF+XueZw1BVbh+Yaw3NP6YpFeNTpsq6u2JpSCq0u2XQSNsb5Z+RbpGuqQ/6riTV2/vHBMQ7ZVg0Eq30SqFWe6adrl93vbF1jBHFJefj6vjLyC5kLOh1ltJ2GPOHLdTQyp7ZSErY19/1EJ3dbl0XjqfsKRniRcaRBjsxHsLQxEXVRSXOjK/ycY+sQiFMxQ8oGDBJJ2Kz1uwoo6zj/tQr1kgQgzrWGSgDZSVgiiSVOmCWQZCixlKYswtx4uoqpvcEPK6KRJpdjJZHBjZrprEkevpei6UI5mlRup503Md8T5pBDEU26u7q9NlK6qN4AoL6v2Hz4c2ivOdVpHqSYQqCMAAZ2XGguQlEGF5S1mAZGHqmjfdYscwJVA/Fwex4Mm9NIkkpixc5RD3P3KTQeqlGX1y6nq5jykpWzqKPJiDKkTi21URMJ31DTEOWs7aF6YYQcLjlOXeN4TWicU60WpGJS1FkrmRo/Tnu+WGJAgLIYvVdx98LEfjC/EAWAcV2xstLRoB1lpBvK6C5pIGwD1ikn0vhTpXJAnPMPJucYHHKYMxWimVxuPFK7rLUOlLAm6IBAwoxzjlAw4ZvixslSOWG5bAIRxg52AxBtGqKiyVhoX0cZskQxlVwG4k4+hZBjBs9fG1aXvOQ4jfylaDzcAJjeK2T0H6GYYiXHQCSxH4kyYtF/DJTFOoPXJSfS+OcsMurQKTvG/Ts1jjYKER4zgdJim65HJrHCumzNdZiflSexUi+zDlfgZEQZi3jTkG9YU4fqxegs83c+D2+/Jx1NSpFkWFeIrw5SZ6mSyFDCNyfKSOoKJXwjEXniecWoKXFEE4gyUi9jDzuasAO1+zLE+ce6q9clJ9f4Z4Z8OXw1zZnKkdMqCzr9ohSZU4vx1fKeghKlCcvyOE1DYWMiRcjxa6SqQMy4HIS8foMces71mpLOYcSphSI1xpoKRU/ZnH9kDqnnZSKaxFlMGTYgGP1nRk8d7ZMpQU+cEV6OKq+fsXlCm2BJc/f1Igshyfb59XZcrLzO/e5qTLJJJlAuG0Hjy6VGsWzz8LZpKOdoAiBeyx2aj9TJqDmJfRpQhHssYmf7APHTW0NJWDklmE4uy3UFck+MIgJ5fiFMP7p/p8aR98VxoGsDZSpyVEvc0XS0T6aEwssYhxwsAbRJrEgSNpxfSNABRA29ey3UNVKJ7FhPQb15QonsSPNPny6XDb3+0fLloTLFnGMQ3L+3x9FVMVHElaLchM85FtrHEvTm7+FxUg65ftGPkEqM0Fk5R4u7f/eliOXHgtcXN5IxFB9KLos5/4RTKwIOlJMPCpd7d7RPlsQQhpgeYZzlTYayDDQzaJWHZiaxGCH6IGCAYqg6xGOa7yWMf8DY2d+FeyXoKCNmuEIJMyBFj9DHLKdKIqPNUCnHm1H95M8FwKOYspLLmYnsEHUWchqzBV3tM4yujRSdJaNi4uWyqbwf/frH2L5MPi8h/bhOObHGnzRAEb66MgpEvb79TnJcEDnJa7JTKC3G+YeahoqArtTZ5qEoI3aNdsGSGzwRnYSiDPP3sKMZBoxkLJlqPhPikGNJ2FDviCyazM0jsTj/QPSUUxIpjZ44RQRkRBONngJzmAJl0QS9NMpIR8k+kANyOf+O9llJQqVyHG4xWCoXo0cCRjKNCGWNK0E6ILIRQsYudY3BctlIOV8I6Vr9cQdKl+WFdBWBjmyrK6cpL6QLMC/XkHLIsaY8IIKQo1FGXlVM7IUz4ZxVZrmskPbJQeMpR1MsM6izUJSRcKDhpjz580qBkHXKiTT+YcMVScJmh5dhzj/0pqEQ7ZMsU4wYBSCAxgOLzP4uxfm3dEU2eCjBaX8XDdFJ9BluGgqhNPu72HMmm/IifQ+hzk9zjTmJ0RSSbB+zbK8xlaCX5jKCEU1GlMGqipFGGSFdiU7uEBUzzsj7WaAQdmp0r0QsIo8dnGj+3iH/LAlSMdEHH0qY5XH+1UYgwre4kYyjVrL5xxpkCo0njGQsqRe9L9LR0Gf0APGmoZCzjm2EWEQTLR2MlFGaewg70BBfLS2XTRcRZOSsIqdYutfC16Wi3L373e71AQlHQ9B08Z6CvC7pWbFs5bns90gri2wnd+w5U7RPLCIPRf9KxbvG1ykn0/iHuMV+uEwxtMEH/R56Sr55YhFDUYXNMmMS4/zNtcQqcCiDnOD8o7qo+wo7tVTTUIgLDumKRTRxJEmj6hiNEDIKQJqHJ5vyIkZBa1MuK68sspErXakSbDbKqcAJUTGVrsPoKZBSpGEbkIrIpbmCEO2zSkQe6m1Zp5xY4x8ydubv8gcfTZhFUKtUV6psM2aQSR4+gAiBFF8d5oIBmq9O0T7x+u+wQY45UClfTb3bwP0eOpcR1xU8Zjn1vCLRE3WN4wx6JN1dnVExlZnwLRbtF6XY75GWAacoUuNAI3kkiraMRMmxMvEQ7ROLyGN7ZZB4Ccy65IQa/zAPb/4e2XTBXEG4WYPWFeGrA5w/kApLQ0ms8IK290WFwCl0F+opcL+X0iXlq4PRU0RXKDlnf5dT126vhbo+9zONa4zNYaApL+bU6rUhLYlMlMsSuszpsjnlw3S5LKdLmq6KMWsjFpGHqKnc/JiEirH6xcdixyLyCKBIHdS4LjmRxp86Kx/I3+Cp5FfoBFEgjiRDi0yKkGM16jHaJ1WmGEpiBXWloqcMjjupi+KQE44mVEYJ0JSFvW6a85eXy8acWjSiSZx7T5XLxqizxVJDBxLZMWMX4vw5ndyhqEtr+hiPUCLb/i6cl9BBzt9cS/g5hxL7saKPuC5hRN7vVevgMOVEGv8Y3wckEJe4yiLcJGP/HtIlTfjGapfd723oiizopFOL0VmROQyheHmpXDiiSdE+4uqsikMmdAWO47C6pA1UsecVjWhi6DN0wGDE2IUO/TPjjLGLonEPxac6uYtQFVMMvESMZAwhF8sEFbNGABg63DEekdM5Gvu7DvlnSmoj0HxfPmqNOZp49UiGMVkjXx1DyMEKnFL/NEIxyekR+vz6GD0SM5Kp7tQQqna/t3F9CaMgb8pL90rQEY28XDa3YspeN4XGTXVLu1zWXiNFTdlzn6LghRgXM5KpvIQ4P5ZJ7xlHk7l+qWvsxd+Wty45mcY/yfmHUVD4JMBcJBkLZWVlirEzX8z3xha0sEwxxZlGkstBDjl2smSshFVI+6SoMymdFaN94gZ5lecldWr0OrSd3LGO7JiTD80HNcaOi9OPMcpNTiXG8hI0uIokfCN5vzhFSp8jFAeAcV1XDe2jlHqjUuphpdSjSqm3EH8fK6V+o/z7g0qpO9ehl5LYmS88emRNCCOyoG1JpPzFFnGKKaeKKc6Nr1FXgmIKJb/d7yV1ZRynTSa/I01DUWAQ6x3JSjrSyVTA3NcycKheqFzWjovRWVLwEiqXNd9Fc+PxHpBYhzqdyI7psjYg9JKakK44FZNRRJBNMV0ltI9Sqg/gbQC+BcDdAL5TKXW397HvB/CS1vo1AH4OwD9ZVW9IYqVynHDbr8k24+JnxdDla4zQPhReCs9ej6Lx6HzEzjbPcKARYxLNZUSOGLB/9yXG+UcTvokKHLkDjSdhKUQYe/1jVXGSkWyndAFh8DKLRWoR8BICVzFdsR6Q1H1RiWw7LtSRndRFObVEwjd8RHuK9omtqfZ9Da4i2ud1AB7VWn9Baz0D8E4A93ifuQfAO8p/vxvANyjqaa5BojRHLLwszwKhF1m6ysIXDroLHYAmTvgy+GoKOcV7CuRdt7nJuWCUweJMQxRTmIePheix+wrpCjUNhXoK7Lh4qaeMsgjxzpUuoVOLofHQOUe1LiHHnYjIQ1FGaP1GoycGxUTtldRJsTSdFYloAuc+AWbdXy20z20AHnd+fqL8HfkZrXUB4AKA69aguyWxBT2uvL5skYU2auzMFw6H7FdLAOmzYmINSnJ6JFXFlMfPktz4IMyNFwldOZUZ0qa8fk8FO7ln1X1F5p44xiMUPZlx9KsLY88r2jQU6DcBwnmJVPOa+V56XExX1CBHHC9VthlKZNtx0ecVSEjb723pigJAWteyPLtnrZz/1UL7AKAspr+qOZ+BUuo+pdRDSqmHnnvuuayL0Rr48ptO47pT49bf4mfR01wwEE7CppJz7mdcKZbhaokszj+xUd3PtHRFq5gyOdOM5Nxay0pLXVSZYug4AyBGWYTn0PZYhOY+hMZHg360KS9mJENrMWwk4xRTTjlqCCiFnHz86I9wN3yIErTfJT1ddhQ5KTYFAEnnFElks0AZMe6Oa3dx67lt8jrWKYM1fMcTAG53fn4FgCcDn3lCKTUAcBbAi/4Xaa3vB3A/AJw/fz4r7rlmd4T3//BfI/9Wb572SxlSiyxakx3dqHTCLMTPhoxkrFSOk4QlE2aD+NvG5As6US0hjDJs05AUSbrz4dNJKSQprslO5ApOb9FbLMQh15y/NAlLH5oGhJ0aJ8oIVQnF94q88MD9jD8uSGcNetg/aO/lVA9ITFdobaQimnhVFz0f1LlPAPA/fNu/R17DumUdyP9jAO5SSr1KKTUC8GYAD3ifeQDAveW/vx3Av9UULDtkqflqChHSSNeOkyeIInx1YvOQqISBMKg3Dc0XS6jAIjN0QOjtRPGuxWgim+Tv4y38lBE3ZYq9QE9B2EjGTrKMotbA3Ofy1amSyDh1JuerpeAl9AIjIJ0Aj9FZ8XNzItFkYE0FE9mByDXX0cwWmsy1AOGqrtgLjFI5qxDTcFSysvaSw/8hAO8H8FkAv6m1/oxS6qeVUm8qP/ZLAK5TSj0K4EcAtMpBj0KixyDEvP6gH/f6ZGIpkjCLJOdCbxuLI4xw2GwplXAiuz2mKpWLdt3Gogx60wXf/RsxkuM+Xf1Ub7oMdBfY4CEjWUSMZCwBHnqhSExX3EiG11TofcZWV4hPD+mKzWHo0DQ7TlrhlnJqsb2y3og8AQDFuafI8yo02dh4lLIO2gda6/cCeK/3u590/j0B8LfXoWsVSVEWFDVixoUMchr5hx68OGGW6EMA5JsnaOyqunYajWcl9Zxr9A32PEBn2XG5CcSQE6UiEyCC7mIRTQpQRKMMqtkwFtGk6JEIQhbSPtHO4JizHvTwYgDwuPfgX5+5Hqp/gY5A7XdJj/6oAGAgLxGlBMmIPE37SMHmUcnxaj9iiVEW0drlHM4/mggMh7KhJGx0o0beNBQNmwd0mWLMqZnfh/nqYV+RUUaIr7YnS4oRF6djNHiNMsqCQyOIufEUD59BucUT2eGTJWPvUggCCmkFTuKlQvZ7SV25e0X6vFLJ5RgoiwCD8PM6lGp3tmyk8adD+9SCltZJhxfZLJHEih5hTIyLvWkoxQUDaJUpxt7+Za8xZBRCnGnImMROlrTXQG86RndqMLSXobvYyZKpYzyCFWShTtgi3AMSj/BiCDlOWcT46lDjYDhKTu0VeUNkzKmRp7BG7iudD4o50MgcCvMmoeOtj1I2yvjXBjnMjVMSDPk4aIZMLqcRhp8YjdEcdlyIHglypoHFGTs0zY4L8s5B6oHedLH3DQDh0sFYtUSK3pOXejIivBzUGokm4zQCfbzDOnn4VAWZdK/EaJ8YGo+W5iYo0mhEI0Tjo6BTy41o6J6do5TNMv6xw6oWyzCaCSRhY7RP7E1D8QoGOnmbomJifHXsHBagje4qXRF0F2pZj0VP7ndbic0hENngyzh9A7Q3+CLSkBPVFaUR4p2w+Y4mFmXQxjUWTcqLCPL46lBPTOotWSFdoXdzmO+ioyfWmUUZTo16GXtVdUbQPrHGwaKjfY5WopUZEcNljYKPxmMHtFl90vAytOnSPHyEiokYBXs9El2xEFh6X9UcSim3wBk9QH1fOXMYO4VV3DGaqMAJRaDuPTTGRF8IHkbIodc/cpyauHkto6EsVeqZGz1Jq5iiB9YF5j5G+1h9UjrrqGSjjH+qaSgWXlJlirGNCpgHT5XYFZHNEzb+YZRmfh8+TTFlkP1x6fsK64oll4F2C38sbLa65NETfV+x0zkBc7+hnoJQIjv0vGKny9proHosYkcTZCeXA06Nw/mHewrWmITNpJiCBQuR/ptURB57XkB7PmK5ODsudHx0Z/yPUGzTUDhZGfL69OJMIuQI4oqVGwJyNB6rHomVALrf7esKjgudEsnYPDl0VvDMl8gYUlekXNZcI13FlKJvgJijiVUxSU8rDfPVsea1IG0ZO88/QftIj0IpYhFN5MA6DnhpFSwwAEUIjYd7CujopDr6I7Cfg/sycgT3UclGGX8g3DQU5xZTXl/eVh/sJAwZyaThiuiSGslMiqmIbJ5czj+WyM6mziJ8da4un8KxRiE296HnZdEpNcZeDzUujvzD3HjoPQX2e31JvjuAOFcpSsVESpVnRU6El15ToQqyWAm2ucYMABg6jiPgMI5KNs74hzxxjIoJV6qkOP9Q2WasqYlOwqYqcIJJ2CLejQlkcuNiI5lILmdET+E8AY0kqzkUdviyIprQ8xLPYSQCDTyv2OmyQBiNp46qNrqoRrRIVVcKUBBznypVjuXUKF2xiAZIRORiijTB+Qe6kGMA8Khk84x/JORLGcnggxeGlya5LDXI4SMGgEhFB6cqxjeSiTr/8BHBkea1AJ0Vozns70PIKWQkw3SWpR5kJZGxyoxUIjCWrKTorGgPSAKEBKnERMKXLpfNq70PXWOM9jHjwt3mYafWJ3WlaZ/Amoq+E4HeK/MU7dNx/leORBGXlEaInMNix+UcwgVQidHUBg933Sbr/INoXHYMQqztP9lTEDsiOHMOQ5VFcTpL3gNC6cqtmCoSFUJAOB8Ui0CXur7/alzk/PrQu39NR3a8OAJoO400vRd4v0E0kU07qCTtE0LjkWNXQueC5UauMQB4VLJxxj9ch8xI9rQ2T3zTxd40lOb812NMojxmbnI5UIETO8gslBiNNeTYawhWMUW4e0pXiooJPa/U+fUA4awZ0ZO0LyPUTVw1yiXponYeKVZrTuUKks2GKSMZPDWTfr9B6DWpZkwgFxd5SxZg5jFY1SWOyFPRf1fqecVI0BMnTioE2kdBx0rK7O/DfLU0sZQoUwzpitI+CXpEmPDlUGchXeHEaJgzjVEqQDsJm57DECLMqGJaWqMQvkaqaSjW9l9x4977KDiJbPvd/rjYsfIkd4cAACAASURBVMIUoODkg0K6+oEXGAF5kWs2RRpIgMfmPkT9zlLOMBTRdLTP0UuOJw69i5NT559zXgmtK73pQiWRSR6zhSRLNC7kkE1lhpQzTdEBK1RMBXM0Qs4/dupoAhGGewrCayrW9k8VEaSpB1pXinemaLBUjibkDItIAxVAF2IkE9nBSIhDW4ZAmcyppRs9w+807oz/EQvliS2PKeX8kzxm4OCpGO2TayTD55yEaZ9sHjNw3k6U9smlswIGeRbjxkOGlRGphU45Dc17qIWfQ/uY+2gbkzgV0y4i4ESFAMGNR6In+33BNZ9BuaUcTcipyYsIwi8wCuky41aoYpLSsYm5PwrZOONPcf7phpwAX53iMSPhZZJb9DdPiseMnu2z3mRlaPOwaJ+Q4Qoa8tBGTdM+7RxNwiDbDb5sz30aIecZSSo6ielahYqhxsXOlKcdjUW6cm48NYd+5Bo75M5eH0BTMcNAItuO859X7DWp7jWE1q/4vKhIlHxUsnHGn6IRDo3HDPDVsWaoMMJI8ZjylvVwRMMplZMt6OzO5X7olNPwfQ3KYzxanH8iuRykLJZxKoZq5OHQWZSu2DHQdlzreUVeJQjEqMTw0cwhXRUaF0auKdqHOhOITzERjiZyX9Sb+WLvhwDqaDK4fqP5hY7zvyKEohHqM9SFSb0Uj0kgwhSPWTuaQFJPfGxv2HAFE6OJbuLwKaccbpwu9QzX7JvfF8S5SqHnVR3jISyXDSFJk1+IGZNe63kVmVSMaftP6Aqi8ZSu9txnRxmHQvvI6Uf3c+64WERDnYPFOTI9pGsQAYBULqN6TWpH+xytkMgpEcpGF1lqQYuTjrYTto0kUzym9MyiUE12atONA2jcHCssN6xAmPYJlogmuPEYX5164QyFrNPPOcRXCxFypDqr0kVEJvZvlIwCqDV2kq39vnDFlLzBTupoZikQkqmLdDSMrmD3mlxdob1sr92nsxaJ53VUsnHGn+LGYy+aML8Pd1bGDmei3jSUQmmhRGDsRexWF4UwopUqgcRoypgM+4FTTiNG0jYNSUP7GOKKzT2lK30yasTRxIwkQe+ligFCyXZeBU7IcK3XSFI5q+z8AqOnILhXEpEr1RAZXxtEcjlhA0J0bAoYUCAkVSV4VLJxxp+kA1KvLQyG6IzNE3rwGaFs1NFQCzqRTA2eS5OkfcI0mJSv5nLjlDOMI662rqrDN9E0RDnDGIdMGmQmFSM1khQ3znnzmv1uf1yU8yfoPXZ+rGgDLHmUnFdWOl+Ejxmx35dfMeXRe5EjIey4MDDojP+RyohY0LF3fpox4UUWO5kvVi2RPj1QRnMM++2moTo5Fy5TpN5vYHWFqyXaBpnDY1JRF+c8dHNN/CjDfl+of0Gc2M9Cd1wjmUElhgCF8P3Jsbr2kK5ZlR8Tdt0maR8qUkvkg0LOOkdXVU2XoGNb65BDZ4Wi/47zP1KJoU+p189JYqVOHAwlRpNGgaBwUgaofr+BTBfFw3N4TKqcL1kPH+L8E4hrTJT0cqqY3Gtyx0UTiOQGTzcAmu8mooyUMQk5tUhTHkDV+TMoJmFyOVgSyaJ9ZBFNiI5NJrIjkVqwcTDaJZ0XFXalnkcsFI+Z5p0DXn+xDKItAOSbhpJvkwrVLjMoJqC5ODncYoiaiuuiHE2urvCL2IEVk3oZB9YBdBI2FdqHe0Bk0SQnueyfS8Pn4fnNa/Ya2xRYQlcwj5RB+6TOfcosBqB0JY/gDrxzYB5pbDTXbqIMtzgiBXiOSjbO+FNNQ5yz8t3PWTGHi8UXGYBG0xAPjdNca2qjAk3ExQkvqcRo7CAz99rdcamSzZCu2NlDdgzQjjJiL2IH6CRsusM33Bmcuq9wElaIkDMMVzKRHUwuJ872IfJIqeO+Q8466yiJVBVTcF+momSi6CNxNHPonQMc+tEvjkgVAxyVbJ7x77fLFJMoLcrDM9C4s2BSDTn2GqWdnxS6S4WydlzOgna/3/23lEOOvYgdAFl+mXKg9m8to7CMn34ZK+eLd8LKKzpCTUNp2ifMjaecGsn5JxxvkCJNJM3b9F7cqVEnqnLe8ma+O6//hrQB4nxLgvYhIrwicRrwUcnGGf8RQeGkXjRhvb68zp/QlaB9zDg5FUMlEGsuWGgkmZy/u8FT1Fmli+CrY9c3Ju4rfw6XwfPr3Wt3x6VexG6vI9twCZPLsZyV+HWHHCMZKliQvoaUsVekkVroZexpUNZuHEz1gADhggXxvuxon+MRKgTme/02wojX+bffNMTRFeoPSJ3DAtCcf4z2GVPJysjx1ub7YlFGqhOWQE7RckNCV6IUFQjcV6qMklgbi8S5T3acjz4rZ5gwknQnbKqEtZ0nsH8LjQHkkWssWZmiSMly2QzAA6SiyRBFKo1c47SPHUcdW5Hi/IHmvkzRzEclG2f8V+Gr5YmlNl89Z9A+ocRoyojbz7ljgDTtQ1dmyDZP1SWdPAOHMnZ5FFOqp4A8poGD0hpOLX70R6Ur4NTC5z7lIeT4wYRCzj/xnCkqJnUMQuiUU1aX9EKTVEzq/CG6/ybt5Kl9mYoY6AZAGaBIFR4clWys8afokajhokoHU/XfRPUIp8EjxLXmGklpd2q6/tseQSGkfQhdReTsIVcXNYdSvjpZa07laLj5BdLRpJPmbvlw6twne42ht6FJD5GLvcDIjgvXqCf2ivQ4jsy9QjooRmWRrytVxQTYRDHl5DPzYx3tc7RCVQhw+er8JGxbl7wrMIfzj6M0Oy50lERIKs4/izqTGYUxgdLqHI00RI9HT1SZYv2yjpjhClVMyfJBtbMWVkyVP4fLZc3vqWMQUs8r1DgopYvStE94PsRFBIlENkVbpnpA7LjcXJw0cj0KWUm7UupapdQHlFKPlP+/hvjMVymlPqyU+oxS6lNKqf90FZ2rChVuc0qvwklYRklkId3goQocBpIUIgxKV7Ltv7ovgjrLOL8le/PEknMhR5Ogpdq6mMBAaoAIpFtXncmBQSyRTZ2rlHqBEUA7w9T59UBg7rlUYlYRQZv2SZUB+7o4eaRQwUIqTwA0j7uoy2WvbtrnLQA+qLW+C8AHy5992QfwPVrrrwDwRgD/VCl1bkW92RKlRxKbNYefBehQNhoqhhK+DF1kElaMZuK0D6kr0ZBjx1GGi8PDU/mF5PMiSwDj8w7QlGD6zKJ2wjeqi2ga4kZqFDeeMiR+voVTDEBFyVVyOVEVI01k53Lj1Jv5VqJIE+uDzu3E5oKIJk9Iqec9AN5R/vsdAL7V/4DW+vNa60fKfz8J4FkAN6yoN1tiSVh55UOq/juMnJJdgX5ZXqr+m0LjDISxrs2Tasix1yg9v55KwqYacuw46riAaDUSZewYBojMByUSgVTTUH3wnMzxFomIxn4nZYBYgKJojoslsu135q4pKec/GvSzq32mhJNPJXypBjvpfXF0HYWsqv0mrfVTAFD+/8bYh5VSrwMwAvCnK+rNFuqQq9pwxQxDIDEqTCCmTpYEAgg5FTYTTUMchEG28CfoLKpagkedEQnfRHKOrMzg5heoo34T/DHQpOm4OZpZ4TcN6Sg6rq5RauxII5l+K5TvDDm15lSimPMGKr/8kpvIbuliXCO9L3nrVwzKAglfDu3TpJjSR6EchQxSH1BK/R6Am4k//YREkVLqFgD/CsC9Wutl4DP3AbgPAO644w7J17NlTDwMbkWHNFmZy1cP+z1cnhQtXaxFlpGEJVErIznXNCZ5czhbaOyMZDw8jzpr13+nqBiqaYhHndVNQ24lVMpI+lQij/Zx+OpxfY2xI6fNNTb5am6ey35/NS6Re7LjqGNGWEeheHNvn0lsnPRk1FifSuoaqX2ZH9EcL+efNP5a628M/U0p9YxS6hat9VOlcX828LkzAP5PAP+N1vojEV33A7gfAM6fP69Dn1tF6KahzMQSm4dvojSAgSSlnD8RZcwTx+/acfnVEvLEqPwoCSI5x6B9Ql23acOlSEfDNZKNf0euz46T0gGh4y5SKNJHrTxg0K4SSq1DoH2iKi+RTdCxCZrOjGveV+oFRnaMe13uv1OgrF0xpaNjxpFc3NXO+T8A4N7y3/cC+B3/A0qpEYDfBvArWut3rahvZak9cV1fPV/EX5EItJOwrGoJgq+ukaTUIPPKL+kqphQ94nfCyjdPzSELu24TBpl650DBvK+ld6AWx3D5G5xL+wBNumieMEDmO1UWnQW0OWQOFUOVHPMQslSX72gYEU2gOCJlIP1ueG5XcFsXh/pdc/R/lRv/nwHwTUqpRwB8U/kzlFLnlVJvLz/zHQD+GoDvVUp9ovzvq1bUmy01wmgaheQia4WygmoJIrmcbFmnyi+FTUP2O+IHp9Ht8VHnRHH+mbRPypgopTDq98iegtQ5LO5nObqA9gZn0T5U2SaD9vGTsBwjSSNJJuefUZprvr/pNDgRzbxR2mjXfD84hipYKBIFFbUuam1w9oq/NsIvMDLX3wQv1QuMMos+rnjaJyZa6xcAfAPx+4cA/ED5718F8Kur6FmnUGhmWsQfINBu5OGF6O0HPyujjNTilFIxKzXJCI9csNc+I6KMlC77fgNbLcIxXCPPmExZDrRXfXZrWJ6xlGj7t9/pc9wpXfS7FGq9sXGkrn54XAhJpt4K5VNMVtc4dl8kbclzagcH85Yucfkwg6bz38zH0UUdhVLkRE9LDa15USEZ4SUKAg5bjlf7MUhoQcc2AdA2yLlh86yInywJtBHholxkrKqYDMMlfb2fReO04Uo7Q99Ixrh7M47WFTVcFOJKVBYB7TLFqeC+fAQaQ9VW14xEyHJjkqIQWk6N6ayNLnmU3IwyFmxdYkfjr0PBffloXBw9CZyav+ZT5bJHIRtn/KnNYw1yahzJp0fPzWnrYkUZ3iJjRRlE0xAnysgp9TTXoryeAj5ClnPIdBJWvsE5z1mR9IjccDGMiZcY5RlJqtKKEdH465cRZVANSibPlUbjTWedpkhDp9JK6VgeCGnnaHjFAHT0xCsfljmao5Djv4IjllD5Wsog+4nAOmMvbxpKRRk+t8ipOKGahlhRhpcYXS7T1RL2GrM3gsdXx7j7Spd0g1NJ2ETbvx0npUeoJGzBoSxCDpQTTbprkRXRtNcGkBllZPYviIEBM28yoyI1sQNN35d/iBwPGLTpWE657FHIxhn/2iA3H0YajWdw/oQB4kUZqrUJgPgis/qkuQzfGc4ZnZ+ULs6x2CTlxrxGMZ0Vqh5J6GptcFFE03TYnEoV/3mldFGlyvzIdXUqhhu5SqOMMJ3F6V9oFznEKcEA9TsURv/MI03s91tJ0apHJcd/BUcsoSoQqWHlVIFQTUNcY+e++L06TIuBkH0DlM5lNFEQtxKB4qtHg3SUYT8L8KolgHAFjrQc1RjJeBK2ZSRLw8LZ4A0jOZevqVzDNS3ShisYqWXkkZJratCMkllRRm5y2TueRJbLkFO/7imnLGdNnuGUjjKOQo7/Co5YqKYhziILojQpX80JZT0agbN57Dg/lOWgT1cHxyjYcS0DxJgLwLmvRdrYmXE04oq38LebhthRBtX9LTUmTEAxpXQxnJpvXFncuJeTADIiGobxb+uS5DKEyeVQpBbNZWRGoF4OhLNXqFfAzop0A+BRyMYZf9s0NPMMFzfha89v4aA0O87PFXCQrv0swNs8Zlw7Wcmpk27oYtRk22v0Q3TpfXEdaMvRLPhRhh23XPKiDJ9D5uUXiKahYokxYw7FFBNVPpzj1CR5EyF1RgGelK5x35bjyiIaP0rmAINQqTJ3/dZRsoSO9SPy+No4Ctk442/ONm+HwOk6/yYK4qLxNgribR5XFzvKIFAQx6kBdV6CHWW0KjoE+QWprr5qId10lOFFT8woI1Q9wntvg3xN5fRluJ8FmGtqoMjoiYXGpVFGKzHK6LoNOTXxc04DJapUmQMA/VNOORSTuUa6EOO45fiv4Bhk3G93IHK58bnUmBD0SHqRNV/8bg2fONwWOLVZpWvB0tVyoBzk5OmSRE/ZjsY6a/bz8stKdTLK8IFBlcvIoASBDDTOyGWM+n3POTGomBUS9PNFG41n3Zc0Si7Sjsboa1dacddUi44VVsZxkuZHIcd/BccgLYTM5GeB+oFP55khH2uRNRGXxEj6JzdyqCKAoGJyHA03yhDmMnKcWtXF6ecyMiqLUlFGm6YT6BJGGbUDFa4pvxOWcY1UqfK0WPAd75I/9xQVM+XkF/yInJtHIhLgKSomP3Lttemszvgfj+SEYS3DtULIx1/QQsNFIYxDui+q9j49xnNqAl1+Q85hoTR/o84WDGPXip4Ezjo3l1HIchm2k9vPWUl7Jbi5DMClSPlUTE6C3ujyn7OsqmtaLNj70n/OPHpPdl9HIcd/BccgVNmmNNkjMSZSeqRFWbARsmodMZC7eTicequuPXujytE4h1IhdXGc2poitRyjwI0yciqmtNPMx4kyyFJlCb2XwY3n5sfEc7+GKHmVarUO+R+TUDXPuQnfceLwLqppiJ3EykHIrVCWp8uiGC6a8TdqTpTB15WPCMVzmJM3yXRqVDMUuxigkBkgf/1yogyrz+pYLs0x5tn0XlbBQiLKaOWRmJVx69iXgijDp866ap9jEpqvZvJ9XmJU3FlZ8BpXACoJK21Q4lMxLaeWCu29d6fKEmbCiKaVhF2hOosx9616/UPKZfhNQ5w59N9vINEF1GtqXuhklGHH5SB499o4Z+UD7X05LRbp5rUV6Nj8IoLm3Kdq9v3TATjr9yjk+K/gGMRHktMF561LzQWdi1qz6vwFuvymIW5DmTwx2t48KfTZSsJKyi+FKK3dU8BDhHkO1HM01X0lnLXfNMRYG3WpssypjVoGOZ3LsN+bU+EGNNcvL8qoqRhxlOEnYTNAmbRarT7fKyNy7Uo9j0fclnCtNa9u3F/QoqoC4TlCoYQvi0YQGsnWfXEbytqdsFLkNJ3b6IkRNkt1DXzjz6fO3KYhTlleVamSwTsDzTwSxyiMXYMsRP6+QU6JS49wK9yo++JEGa4uNoKvDPKi0gWkowwKKHEpUmkuLsfRHIUc/xUcg7icP+e1b4CzoFt8H8Mge1GG2EiK0LiUQ25vVJ6uNl8tpbMqB5rq4hzIE2b++2e5dBZVpph6xlWlijDKqByvc42pOQSa3LjEqQFy4+8+Z0mCE3DQ+GKR7DQ345QcXHlv5uPmMlxHU0UZUhsgcFDSEtajkOO/gmMQ1xNLWrTN52uvP+ynX8jg8tXsKKNaZM0OX2lDmaRl3e8mZukShrL++w2k1T5umWLS0fTaxg7gUUyNa2Tys27TEP/cJ8/xCnS1e0BkOSsOJVjpkgKDVhKWF9G4EZ6NMthUolDXiHBqnF4JVxf7eJI+sVc643884lIWkkQg0Hzw7AUtjTK8N16JjGRmVUwOjSDu8M0Mm+0GL5b1Zk2N8Q/UkiYr3WtkPefBClRMwaeY7DjJcQauLjcSkhpkbnWLn2/hoOpKl/h55UY0KsuIt3T1OVGGakQZxZI3H4ctx38FxyAu8q8XmbSkTBA2i6slPDpgwYsy3M1jjxhgn7QpRK2jEhFaND7PcTS51SNcdOfQRWyazjtaQ/Kc/SSs/DhtHh3golb+4Xjt86KknL+kOsv9PDu/sEIuw35+zqjOsuPa0ZNMl4w6k635o5Djv4JjENIoZFTgcGp13aYhCYL3dfGNnRlTvVxamBiVXKPbNMTJZfhNQ4edrCRRq7hMMcNwMTd4NmXhUG789Zvn1Jr3JYsypFGyuy/rUmoenSUGZQN5lNECgMyKKTdvUtNZXZ3/sUiDW8w1QDlGwVYiHFZlRr/t1FK8rl+ZMWXnMmoKh5vLMNej2nMv5caz6JG8Dc55J4Idl+NA3c9nIUlhlGGdBjfKMJGr8HTZzL3SSPgy72vsrEOJrlWijGaCPn0uv7sOp0wHehRy/FdwDDIaEEmsVGlYVZlRIxNpco5bF1wboFoXBym4TUO5vDM3onENMjeXYcflnIoI5GxwlXWwm9ElN1y+U+OcRZ+tS1oJVv596upi5jJyzrIxOuqzfQ4vUsuMnhynxk+aE/kg5jpsUUyMazxsOf4rOAbJ4eCq0sGMCoacEzOBDJrDqUbgVjG1qBhuKOtU7kh4zJE3H6wow0HI3FyGvcbspF5GEjYb+a9QFSOtYnKdYbbhklZMHWrSvJ2z4u4VcR+NHxUymtDMNcrv6yjk+K/gGMQtvarROLOCwVlkKWQHmAdvm4ZyQ3RJKGs/zzV2Vp8cOdWOhrt5AMtXCytpHGfIzWXYcTkvxAFyKSbpsRVeRJOThBXSWVJdeRRTE7xwq32avRJMXVTe5LAqpnpNp2YcDSNKJkBIV+d/TOI2DVUJGMYLsIG85Bxgmoa4Z/T4r5oU6yqWDvXAo3Ckxs6dj9rY8XS5UQbv+mp0J0FOI/+++j1GxVStSxJlNJKwXIOcyflTuZ3DKiLIQa3+XpnM08clG12qnV/gzqHjDPkVU7KI3C8f5hwDXeuqc2McXUchx38FxyBu09CkPGJgK6NJRoTGiyUmTEfjv2pSQj2Ya3Tvi2kk3YhGpMt1ajK+WuzUFsvqvlh5CY9GkDyv+WJZrhGeAyWNSTLKaPaOTIoFtpjOelV6ZFIsWboalXGZpbmTYoFt6X0thFFG6TQmc37OSlr04V/jdL7EFiP6d6OurtTzmGXkJG4mpeFKPUT/TUNT9iKjdDGNiYOCJAbZOBq+Lt+YSBzNtFjioNS1PZJtOqlBnhVLTMqxHGPS4Ku5B5k1HI0Zy5vDNnXGOcjM6lqURwykQAiQWaboIeTJfMEyXDkJej+XMZkveQ50QOlKR8n9nqqom6ngvpZlqbKYtixtAN9Z106+Oh+pS/gej7gb/GDGM5L+S5+54aXLSU5mvCjDXmNO16K9NpHhcl7uLY1oTJRhdcn4arZTG7i6SuTP3OANlCaNaOY8YGDH1e9E4NMBgO+smdFTURtxpSQVU3XEywMhTV1AOhLyk+ZcR9OMniTPuV6/7Pty5sM+N+lznswXPGfdiCZ5c3gUspHG330YFklyF9nM2eBy1Crb4FIjOW6gVr6uRqJtzuRMKV1szt9s1IP5gh0tWF1cZ23HzRxd3DFAM6LhonGX495i3JdbPSKJ1FwjaQ2QJMqYLZZYanlUaJ18KupqceNzHsVENUNxyzbtvjxgI38HKM1s5Dpg6KrvS7p+Z8UJQv5KqWuVUh9QSj1S/v+ayGfPKKW+pJT6Z6voXIc0H4Zg0w16FSI5mC1ZD765wflo3OiqNzhnYTadmuS++tWiPJjz+VnAJJe5iBBoOpqDGRd91nM4Fd1XnYTlI10qokmPG7tObSafwwNhlJHrQBu5J2bU5Rq7Ub+HfiJpbvXNSjprtliyo11Df9Xrl7vHXAfF1QWU67fg58eGvi4B5z91dO0w7uuwZVX38xYAH9Ra3wXgg+XPIfnHAP7vFfWtRVxjYg0XZ7NuD/vVxplKjWSGLotyjTGRcv58bnx72KuujW9M6r4HEWXhbh520rFOIErua9TP36iN5zWSUUyT+ZK5Npx8kBAY1E5Npmu20BXgYRlWp1KFS99YfVI6y97HwXyByczQWdxcl63OmhRCZ+hEk/z5cKuYeHvZfv5gtmTrOmxZ1fjfA+Ad5b/fAeBbqQ8ppf4KgJsA/F8r6luLDJ1Drg7mC/R7KnkMAmAeWGWQBYklwPCfkg2+PRpUaFBSLQGUC1rAjW+P+rUuIT0yFxrk8aBXfX4idGpSOms86FVIS+rUpHTW2HOgEkcjpn28++I8YzdnJaGzGqiVuTbsuLnwvuzzmcwWZg4ZdBZQO15bncWK1Mr7msyF9J5bJcSksyzK3y/vC+Bd42HLqsb/Jq31UwBQ/v9G/wNKqR6A/xnAf7mirrWJe8jVhJkIBGojaZuNWOiT2OAcNLM97DWQv4RDnjvojrXphgPsz2rjzzXiVpdkg++M+jiYFUYXs1qikaAX6NpynDX3vtzO5dqB8tDdfgMYMHQ5TUMSOmtn2Md8oasiAs59AbaOXkY/+oaLi1ito5kIkqn2PsS6ygocCZ1l7+tgtsBkboocUj0gAKoXs1R0ljSiEYCXw5YkkayU+j0ANxN/+gmmjh8E8F6t9eMpL66Uug/AfQBwxx13ML9eLm4dsgTNmA1eiAxQQ1eZJOYssp3RAM9dmmK51JgWTB6Tqh7hhKWjfo1ahXz1bLGQoTvXSEp1FUv0e6UxYdzXjmeQuTkJwCJdgeEa9TEtDMc9nS9wbmeUHGMTo418EAcYOAaZ60CBukRUYoBq41+wq1sqXcKS48ogl/QI36n598WLrI2ugg0MAFv9tKicNYtWdRzNwWyBHqM66ygkafy11t8Y+ptS6hml1C1a66eUUrcAeJb42BsAfJ1S6gcBnAIwUkpd1lq38gNa6/sB3A8A58+f1/7f1yVVeDlflFww0/iPBrhwMHeqAxjo0/J9xYJdbgjUjkaS+LKb2d5XT9XOJyY7pS6tNRu1Vvc1X+JAaCQPhAjZ1aWUzHAdzBeGC2Y6GlOfb3I6IjTuGq75ArcIAIXhgiW6SsNVGpMdRjGA1WWQrtxITuYLHMyXrAgUqA2yvS8ON77Vcmq8vWKoRNl9NSKaGV/X1rCPS5NCSOE2kf/2kEdnHbas6n4eAHBv+e97AfyO/wGt9Xdpre/QWt8J4EcB/Apl+I9S7IPfKxcZhzMFysSow9txjMluuXn2pwt2dQtg0biTjOIs6HGpy6I05iKzBtmWAHIczc64nMNpIYoydkZ97JcGmcuZ2nvfmzm6GNe4PRpAa+M0THI5/Zx7PYWdYR97s4Vo7rddgyygLHbGfexPFyIn76LxAwF42Rn1xbzzjkPFcHM0VtfBzEHInPtyEqOTmWAOR/2KvgH4BRVA7Wi4yN/qklRn7QztvpStjcOWVY3/zwD4JqXUIwC+qfwZSqnzSqm3r3pxhyU7lUEuMJnxQ9md0QD7cxntUzuaQhSi+SDJFAAAENVJREFU+xQTz9FYg7xgo2qgzmVMZgIu2Ns8EjprsTT87GzBC+37PVXRRVJHY66xYFNMgDHk+0KEbOfDGiEukjRrauHQPnwkuT9bsDtaK12zQkRn7bTQON9I7s0W2XQWt+TY6BqUuuR01sFMBsp2R4MmCGHlnsz12KjwSkj2AgzaJyZa6xcAfAPx+4cA/ADx+38J4F+uonMdYtG4Rf5cT7w17ONgJkPju2PH0Yg2ar+ZIGIhQstjLmSJ7DKBeHEyr35OyaDfw2jQKw2QTBcAvLg/A8BPfBnUWpjD2Zh0Vq4x2R1byi3DSM5LQCFB4270JDFcc6mRNAZZQmf5fPX2Ob6RfObSJJ/zny9waswzT9XzyqgssmuD/bzGTRDCobPsvpzMatrnSpDjzzocg1jK4qBEQSKDLEz41kZBll/YKnsK9qZ8RzMa9DDoKUPFFLwKIfcaXyoNMqeuHTCRhkXVEmMHAC9etrpk9IjdqBw6y+q6cGCcGrelfmc0MNGT4DiO3CTsrhdlcCuLrC6J4dodD1Y2kmxHMx4YOkvAjdvPHAgjNfu8RM562KRIJbrc6ImzfhvPS0BnHbZspvF3OP8DAe2zPWyicc5DHA8MUt3PoGIAxyALeV0JnWWv6cU9qS676eT3VSN/ebgtdjTi++qbKpBiwToG2v3uS5MC8wWvDBgw89E0yHx65MAaLqYx2bacf+XUpEaSd0AbYIDB3qzAflnWy8pZOXmTKbMB0OqyDgPgoXE3opECwMl8Wd0XZw77PYXRoHfF0T4bafxN84ihYvZmBU5tMaslRn0sdY0kOQtaKVUZrr1pgdPMUDYbjY8H2JsWuDzl35fV9fK+uS9ZArHA3nRR0VvpMeZzL+3JjL81XPuzBXbHXGdtdL0o1LUz6mNvusD+lK+rfV/c52WomMvTBYZ9JaIR9qYyem93ZKKnvdJIcp6ZX6kiyi9MF9ibFqUuAUIuq59EUcZ8gculLg5dZKPkfaFBtpTxsxen5mfmujeVVoUo+j9s2Ujjbys69mdmcUoeIAC8cFnIV49NhYBE15ani2MUgLqaZm9WsDnT6r6kRnJsKIvL0zlOsY2k+dxzl8zm4Y6z9MilSVFtQKkuriHfHQ1wMFuIHOh2SxffGdq1wX1e9r5e2JPr2psZYGAcDT+/sDeVrSnj1IrK0XDmsSpVni2wPyvYEc3uqA+tgecvT9m6gLrKbV+wLy1l/Gz5nE8LANZB5dSuDLN7ZVzFMcjO2FQIXJoIjKTd4HaRcTdCWY1wWbLIPF3cRbY7HmB/WuDyRODURnZBTwBI7qtG/twx1rHUuoascQaNl9GT0CBbXZKNao2k1NFI59DqEoEQq+uizNjZMkW7NlhlwBUImUJryX0NsCwNMjeiUcpUde3NTHQidYbPXJTP/f6swKUMx2ufswQ47ltAwVzzhy2ba/xHfVyczDEtluIH/9xFa5B5D3G7rOi4LDFcw0xdw5pGkCJ/q+sMU5elR3Kc2rOXZA7U5jIuCzbq9tDXxbyvsqLj8oT/vPz74o8z9IjIALXuiw8MiqXGi/sz9pheT2Fr2Kt1sUFIfY3ctQGYeXxpb4bFUguc2qDSNejxIho7zq4pyfOyugD+3Nsu+kuTOVvXYcsGG/9BZeykfPXTFycYD3rsV7HtjgZVVyAfSZrPPXVBiMbHLo3ArwJxdUk2wsFcunmaqJU9rqSYJAjZ1yWJ1PZnJh+Ui8a5znp31MdsscSF/Tn7+gb9Hkb9XlZEAxgnz9Vlxg0yIhrLjU/Yax4wjle6Dq2jee7ilB3RAAYcvHB5Jopo7L08d3GKvsDRWBsg2SuHLRtr/HdHfWdB84zk2W2zoR9/aV/0AHfG/Yq+4fLO53aMride3hc5mp1RH5cmcxzM+UlYe19PvLwPQIbu9kqKibt57OeeFobou06Izp17ex9Wl4QumsyXuHjAN8ij0iA/deHA6M6gEiUI+dTWwAEG/EgNMJSFxPjvjvt4WgwMLBUzFe2Vc9sjfOll2RzWaFx2X6fGg+p5SaJ4q2t3xD+m4ezOEE9fnGCp+XN42LK5xn88cAwQ78FfYw3ySwfsxWJ1SY3CuUxdpxr3Jde1NeyxjrcGDJq5cJDnaB57UepoDPK/cDBnI8nxoI+dUR+Pl7rYTsNxUNw5VErh3M4Qj78ke86VrgsT9lwA5pnZ+5JEhYC5L4mjuWZn5NyX3NFIdJk5NPfFpR8toHpGGNE0nlfG2pDsy3Pb7vPqOP9jlWt3R1WjBvfBny2NpCRMBIBrd2pd3AVzbntU6TojMArX7Lq6+Au631PQmn99Vpd925iEsjizZY54kDiaa3fNfMyKpchIXrMzQrE0ZwRyjZDVNZnz80FW16LUxZ37a0pdB/MFuwzY6ipVsXVdu+Pcl3AOxfe1kzeHZ7eH0OV9ca/R6jqYL/Lvi3mN1+ya/SG9r3M7Q/HzOmzZaONv5brd9PG7QG2QAeDMtsD4O99/LVPX9qhf8YmSxXJdQ9eYNUYphXMlIs/XxbsvoDZ4NgrgSM7zAuqoZnfUFzsaALj2lFyXUnxneP2pzDncqefuDHMerztVrwfJHObpWu15AXzk795XzjoE+Pv52p285+Ue8d0Z/2OWHIM8GvSqw9OuP8UzrIC3ETKMiURXzn25um7I1HWdYJzdCPn3xR9nUeH1pwXPazfXSJrPXrszwoDtaBxdGXM4GvTYkWHTgcp1AU1nFRP3+0XrcFuu68zWoDrriTsGaDq1G05tscYM+r1qr0iAgQt0JOv+MGVjjX8uarUbIddI5myEGwSGKx8hy3XlortrM5ya+1mJA7XoTuRoTuUZSUsJ5Dq1HDR+w6kxO+noGjuJ4bJObWfUZzcbbo9MFz0gdWr1NV7DeCEOYPMtpePNcNYAcP1pwfq1ayozorlRsMcOUzbW+LtGTtJubUO2m87wkALQdBTcBQ0AN54Zl/8X6HLuS3KNt53bNrpOS+6r/uzNZ/njbj5rdeUZf8l93VJel8Swup8V3dcZc1+SsN5F7ZL7snPYE+xgNxq5OWMO7YvLuWK5+5vP8p/zK67Zqf7NOVOppUs0h/VnuS/EAer+kZsEa+N2574kDuowZWON/923ngEgD8Fec+Mp8/+bTrHH/MVbjK5Rv4e+YEG/8jqzYO66ka/rNc5nJacHWl1fduMue8wrrtmu/i0xrvZ+7rh2J/HJWtwNc7ujNyV/4abTAJpGJSXm1FDz71ddz5+Pv3iL0XW74L6UUhUif43gOf+lm42uL7uBPwao70eiy97X3beeFel67R3njK4bTrPH3HWTfG0AwF9+hbm2L8u4L0nuCQC+vJz71wjm3p1vLiV42NJ/61vfetzXQMr999//1vvuu+/Qvv/01hCjQQ/3vuFO3HEdf6H9hZtOY9jv4c2vux0DJuzaGvZxemuI7/qrd+DVggVz142noTXw3W94JTtZOej1cMPpMb7tq2+rDB9HvuyGUygWGve+4U52JKSUwu3X7OCbv+JmkWF45XW7mBZLfN/X3ikqA7zrxtP4mtdch6+6/Rr2mNuv3cZkvsD3f+2r2MlKAPjKW8/i37/9HP7qq69jj7nt3DYOZgv8va97VSOZyNH1mptO4eu//Eb2mJvPbmFaLPEDX/dqEVX3Fbeewa1nt/HGr7yZTRfdcGqM+cLokkRCX3nrWZzZHuI/ee1tImpKKeD7v/bVuE3g5P/SLWcwHvTwHf/B7WyAtTseYDzo43u/5k688jqBky+N/3e9/pVsQ25yM0P83de/UuywpfJTP/VTT731rW+9P/U5pfWhvSp3JTl//rx+6KGHjvsyOumkk06uKlFKfVxrfT71uSsj/uikk0466eRIpTP+nXTSSScbKJ3x76STTjrZQOmMfyeddNLJBkpn/DvppJNONlA6499JJ510soHSGf9OOumkkw2Uzvh30kknnWygXLFNXkqp5wD8+QpfcT2A59d0OVe7dHPRlG4+mtLNRy0nYS5eqbW+IfWhK9b4rypKqYc4XW6bIN1cNKWbj6Z081HLJs1FR/t00kknnWygdMa/k0466WQD5SQb/+Spdhsk3Vw0pZuPpnTzUcvGzMWJ5fw76aSTTjoJy0lG/p100kknnQTkxBl/pdQblVIPK6UeVUq95bivZ12ilLpdKfUhpdRnlVKfUUr9w/L31yqlPqCUeqT8/zXl75VS6ufLefiUUuqrne+6t/z8I0qpe53f/xWl1KfLMT+vuG/gOEZRSvWVUn+klPrd8udXKaUeLO/tN5RSo/L34/LnR8u/3+l8x4+Xv39YKfXNzu+vqrWklDqnlHq3Uupz5Tp5w6auD6XUD5f75I+VUr+ulNra5LVBitb6xPwHoA/gTwG8GsAIwCcB3H3c17Wme7sFwFeX/z4N4PMA7gbwswDeUv7+LQD+SfnvvwngfQAUgNcDeLD8/bUAvlD+/5ry39eUf/sogDeUY94H4FuO+74Z8/IjAP41gN8tf/5NAG8u//0LAP5++e8fBPAL5b/fDOA3yn/fXa6TMYBXleunfzWuJQDvAPAD5b9HAM5t4voAcBuALwLYdtbE927y2qD+O2nI/3UAHtVaf0FrPQPwTgD3HPM1rUW01k9prf+w/PclAJ+FWeT3wGx6lP//1vLf9wD4FW3kIwDOKaVuAfDNAD6gtX5Ra/0SgA8AeGP5tzNa6w9rs/J/xfmuK1KUUq8A8LcAvL38WQH46wDeXX7Enw87T+8G8A3l5+8B8E6t9VRr/UUAj8Kso6tqLSmlzgD4awB+CQC01jOt9cvY3PUxALCtlBoA2AHwFDZ0bYTkpBn/2wA87vz8RPm7EyVlWPpaAA8CuElr/RRgHAQA+yLY0FzEfv8E8fsrWf4pgB8DsCx/vg7Ay1rrovzZvYfqvsu/Xyg/L52nK1VeDeA5AL9c0mBvV0rtYgPXh9b6SwD+JwCPwRj9CwA+js1dG6ScNONPcZAnqpxJKXUKwHsA/Bda64uxjxK/0xm/vyJFKfUfA3hWa/1x99fER3XibydiPmCQ7lcD+N+01q8FsAdD84TkxM5Hmde4B4aquRXALoBvIT66KWuDlJNm/J8AcLvz8ysAPHlM17J2UUoNYQz/r2mtf6v89TNlSI7y/8+Wvw/NRez3ryB+f6XK1wB4k1Lqz2DC7r8OEwmcK0N9oHkP1X2Xfz8L4EXI5+lKlScAPKG1frD8+d0wzmAT18c3Avii1vo5rfUcwG8B+A+xuWuDlJNm/D8G4K4yqz+CSd48cMzXtBYpOchfAvBZrfX/4vzpAQC2IuNeAL/j/P57yqqO1wO4UIb97wfwN5RS15QI6W8AeH/5t0tKqdeXur7H+a4rTrTWP661foXW+k6Y5/xvtdbfBeBDAL69/Jg/H3aevr38vC5//+ay4uNVAO6CSWxeVWtJa/00gMeVUl9e/uobAPwJNnN9PAbg9UqpnfJa7Vxs5NoIynFnnNf9H0wVw+dhsvE/cdzXs8b7+lqY0PJTAD5R/vc3YbjJDwJ4pPz/teXnFYC3lfPwaQDnne/6z2CSV48C+D7n9+cB/HE55p+hbAK80v8D8PWoq31eDbNBHwXwLgDj8vdb5c+Pln9/tTP+J8p7fhhOBcvVtpYAfBWAh8o18m9gqnU2cn0A+CkAnyuv91/BVOxs7Nqg/us6fDvppJNONlBOGu3TSSeddNIJQzrj30knnXSygdIZ/0466aSTDZTO+HfSSSedbKB0xr+TTjrpZAOlM/6ddNJJJxsonfHvpJNOOtlA6Yx/J5100skGyv8PUO7vmWpEPAIAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "rsp = nk.rsp_simulate(duration=90, sampling_rate=1000, respiratory_rate=15, noise=0, method=\"sinusoidal\")\n", "nk.signal_plot(rsp)" ] }, { "cell_type": "code", "execution_count": 25, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 25, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX4AAAD8CAYAAABw1c+bAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAHERJREFUeJzt3X1wXXd95/H39z7oXj1cWbKlOH6I48QQ00wYlCDcUBgw4SlRWYJ3urPJtA0L3XWWrXeAZYaFbXeSFmYfGNKUh9aMY6eBlpq2CbRMcEtpQnF3lwbkxHUc7MSO4ySyHUuOn2TJerr3u3/cI/tKOldXvpJ15XM+r5k7Oud3fufc3z0+/uin3z0P5u6IiEh8JGrdABERmV8KfhGRmFHwi4jEjIJfRCRmFPwiIjGj4BcRiRkFv4hIzCj4RURiRsEvIhIzqVo3IExbW5uvXr261s0QEbli7Nq164S7t8+k7oIM/tWrV9Pd3V3rZoiIXDHM7OWZ1tVQj4hIzCj4RURiRsEvIhIzCn4RkZhR8IuIxIyCX0QkZhT8IiIxE6ng/+oTB/jJC321boaIyIIWqeDfsvMQOxX8IiLTilTwN2VSnBsaq3UzREQWtGgFfzbFuWEFv4jIdKIV/JkUZ4dGa90MEZEFLVLBn1OPX0SkokgFv8b4RUQqi17wq8cvIjKtaAV/Vj1+EZFKIhX8uUyKcyNjFApe66aIiCxYFYPfzB42s14z21tSdr+ZHTGz3cGrK2S9tSXLd5vZWTP71Fx/gFK5bBp3GBzNX863ERG5os2kx/8IcHtI+YPu3hG8dkxe6O7Pjy8H3goMAt+bVWsraMoWnySp4R4RkfIqBr+77wROzvJ93gu86O4zfiZkNZoyQfAP61x+EZFyZjPGv8nM9gRDQa0V6t4FbJ/Fe83IeI+/Xz1+EZGyqg3+zcAaoAM4BjxQrqKZ1QEfBv5qug2a2UYz6zaz7r6+6m60lrvQ41fwi4iUU1Xwu/txd8+7ewF4CFg3TfU7gKfd/XiFbW5x905372xvb6+mWRrjFxGZgaqC38yWlcxuAPaWqwvczTwM88DFMf5+9fhFRMpKVapgZtuB9UCbmfUA9wHrzawDcOAwcG9Qdzmw1d27gvkG4P3jyy+3XCYNaIxfRGQ6FYPf3e8OKd5Wpu5RoKtkfhBYUnXrLlFjJgloqEdEZDqRunI3lUxQn07qdE4RkWlEKvhBD2MREakkcsGfy6Q0xi8iMo3IBb96/CIi04te8OthLCIi04pm8KvHLyJSVvSCP6sxfhGR6UQu+HPq8YuITCtywT/+5a67nsIlIhImesGfSZMvOEOjhVo3RURkQYpc8Ocu3JNfV++KiISJbvBrnF9EJFTkgv/C4xd1Zo+ISKjoBr96/CIioaIX/HrurojItCoGf/Aw9V4z21tSdr+ZHTGz3cGrq8y6LWb2qJntN7N9Zvb2uWx8mPGHsajHLyISbiY9/keA20PKH3T3juC1o8y6XwH+zt3fBLwF2FddM2fu4nN3dVaPiEiYisHv7juBk5e6YTNrBt5F8LQudx9x99OX3MJLdOEpXOrxi4iEms0Y/yYz2xMMBbWGLL8e6AP+xMyeMbOtZtY4i/ebkUwqSV0qodM5RUTKqDb4NwNrgA7gGPBASJ0UcAuw2d1vBgaAz5XboJltNLNuM+vu6+ursllFOd2aWUSkrKqC392Pu3ve3QvAQ8C6kGo9QI+7PxXMP0rxF0G5bW5x905372xvb6+mWRfoYSwiIuVVFfxmtqxkdgOwd3Idd38NeNXM1gZF7wV+Uc37XapcNsXZ8/pyV0QkTKpSBTPbDqwH2sysB7gPWG9mHYADh4F7g7rLga3uPn56538Gvm1mdcAh4GNz/QHC5DJpnccvIlJGxeB397tDireVqXsU6CqZ3w10Vt26KjXXp3jpxMB8v62IyBUhclfuAjRn1eMXESknmsFfn9YYv4hIGZEM/lw2xcBInrG8HsYiIjJZJIO/OVu8X4+Ge0REpopm8Ncr+EVEyolm8Ac3ajurG7WJiEwRyeDPBUM9+oJXRGSqSAZ/c/14j19DPSIik0Uz+Md7/BrqERGZIprBX6+hHhGRciIZ/OMPXNdQj4jIVJEM/mTCyGVS9GuoR0RkikgGP4zftkE9fhGRySIb/LlsSl/uioiEiGzwN2d1ozYRkTDRDf76lG7ZICISomLwm9nDZtZrZntLyu43syNmtjt4dZVZ97CZPRvU6Z7LhlfSnE1rqEdEJETFJ3ABjwBfB741qfxBd//yDNZ/j7ufuNSGzZbuyS8iEq5ij9/ddwIn56EtcyqXTdE/PEah4LVuiojIgjKbMf5NZrYnGApqLVPHgb83s11mtnEW73XJmrNp3GFgROP8IiKlqg3+zcAaoAM4BjxQpt473P0W4A7gt83sXeU2aGYbzazbzLr7+vqqbNZFulGbiEi4qoLf3Y+7e97dC8BDwLoy9Y4GP3uB75WrF9TZ4u6d7t7Z3t5eTbMmaNatmUVEQlUV/Ga2rGR2A7A3pE6jmeXGp4EPhNW7XHJ6/KKISKiKZ/WY2XZgPdBmZj3AfcB6M+ugOIZ/GLg3qLsc2OruXcBS4HtmNv4+f+7uf3cZPkOoC0M96vGLiExQMfjd/e6Q4m1l6h4FuoLpQ8BbZtW6WdA9+UVEwkX4yl2N8YuIhIls8OeCB65rjF9EZKLIBn86maA+ndRQj4jIJJENfih+wat78ouITBTt4NeN2kREpoh28NenOaMvd0VEJoh08Lco+EVEpoh08C9qSHN6UMEvIlIq2sGve/KLiEwR6eBvqa+jf3iM0Xyh1k0REVkwoh38Dbp6V0RkslgE/2kFv4jIBZEO/vH79egLXhGRiyId/C26UZuIyBTRDv6GOgBOnx+pcUtERBaOaAe/hnpERKaoGPxm9rCZ9ZrZ3pKy+83siJntDl5d06yfNLNnzOzxuWr0TGmMX0Rkqpn0+B8Bbg8pf9DdO4LXjmnW/ySwr5rGzVYyYeSyKd22QUSkRMXgd/edwMlqNm5mK4FfBbZWs/5caGnQ/XpERErNZox/k5ntCYaCWsvU+UPgs0DFS2fNbKOZdZtZd19f3yyaNVFLfR2nB/XlrojIuGqDfzOwBugAjgEPTK5gZh8Cet1910w26O5b3L3T3Tvb29urbNZUi+rTuoBLRKREVcHv7sfdPe/uBeAhYF1ItXcAHzazw8B3gNvM7M+qbmmVFmmoR0RkgqqC38yWlcxuAPZOruPun3f3le6+GrgLeNLdf6OqVs5CS32aMzqrR0TkglSlCma2HVgPtJlZD3AfsN7MOgAHDgP3BnWXA1vdvezpnfOtpaE41OPumFmtmyMiUnMVg9/d7w4p3lam7lFgSui7+z8C/3iJbZsTi+rT5AvOueExctl0LZogIrKgRPrKXSie1QNonF9EJBD54F/UoKt3RURKRT74x+/Xox6/iEhR5INfPX4RkYkiH/wa4xcRmSj6wX/h8Yu6bYOICMQg+LPpJJlUQhdxiYgEIh/8UOz1n9KN2kREgJgEf2tDnb7cFREJxCb4Tw6oxy8iAjEJ/sWNdZzUUI+ICBCj4D+lHr+ICBCT4G9trOP0+VHyBa91U0REai4Wwb+4IY27LuISEYGYBH9rY/Hq3ZMDwzVuiYhI7VUM/uBh6r1mtrek7H4zO2Jmu4PXlHvwm1nWzH5mZv9iZs+Z2e/NdeNnavGF4FePX0RkJj3+R4DbQ8ofdPeO4LUjZPkwcJu7v4XiQ9lvN7Nbq29q9VobxoNfX/CKiFQMfnffCZy81A170blgNh28avLt6pKmYvDr6l0RkdmN8W8ysz3BUFBrWAUzS5rZbqAX+JG7PzWL96uaevwiIhdVG/ybgTUUh3COAQ+EVXL3vLt3ACuBdWZ2U7kNmtlGM+s2s+6+vr4qmxUum07SUJdU8IuIUGXwu/vxINQLwEPAugr1T1N82HrYdwXjdba4e6e7d7a3t1fTrGm1NugiLhERqDL4zWxZyewGYG9InXYzawmm64H3Afureb+5oNs2iIgUpSpVMLPtwHqgzcx6gPuA9WbWQfHL2sPAvUHd5cBWd+8ClgHfNLMkxV8wf+nuj1+ODzETum2DiEhRxeB397tDireVqXsU6Aqm9wA3z6p1c2hxYx2HTpyrXFFEJOJiceUuBLdmPqcev4hIbIJ/cWOagZE8Q6P5WjdFRKSmYhP84/fr0ZO4RCTuYhP8Sxp1EZeICMQo+Mev3tVtG0Qk7mIT/ON36HxdPX4RibnYBP+Fe/Kf0z35RSTe4hP8DXUkTD1+EZHYBH8yYSxuzHBCPX4RibnYBD9AW1Mdff0KfhGJt1gFf3suQ5+u3hWRmItV8Lc1ZTihHr+IxFzMgr+OE+eGca/JEyBFRBaEmAV/huGxAueGx2rdFBGRmold8AOc0Di/iMRYvII/Nx78GucXkfiqGPxm9rCZ9ZrZ3pKy+83siJntDl5dIetdY2Y/NrN9ZvacmX1yrht/qdqDHr9O6RSROJtJj/8Rwh+S/qC7dwSvHSHLx4DPuPsvAbcCv21mN1bf1NlryxVv26Aev4jEWcXgd/edwMlL3bC7H3P3p4PpfmAfsOKSWziHFjfUYYZO6RSRWJvNGP8mM9sTDAW1TlfRzFZTfP7uU9PU2Whm3WbW3dfXN4tmlZdKJljcUKeLuEQk1qoN/s3AGqADOAY8UK6imTUBjwGfcvez5eq5+xZ373T3zvb29iqbVVlbk+7XIyLxVlXwu/txd8+7ewF4CFgXVs/M0hRD/9vu/t3qmzl32nJ1Cn4RibWqgt/MlpXMbgD2htQxYBuwz93/oLrmzT31+EUk7mZyOud24KfAWjPrMbPfAr5kZs+a2R7gPcCng7rLzWz8DJ93AL8J3DbdaZ/zra0pQ1+/btsgIvGVqlTB3e8OKd5Wpu5RoCuY/j+Azap1l0FbU4ah0QIDI3maMhU/vohI5MTqyl2Aq4Krd3vPDtW4JSIitRG74F/anAXg+FmN84tIPMUu+K9eVOzxH1ePX0RiKnbBP97jf03BLyIxFbvgz2XTNNYl1eMXkdiKXfADLF2UVfCLSGzFM/hzWV47o+AXkXiKZfBfvSirs3pEJLZiGfxLm7P09g9RKOjqXRGJn1gG/9XNGUbzzslB3Z5ZROInnsG/KDilU+P8IhJDsQz+qy5cvavgF5H4iWXwX62LuEQkxmIZ/O25DGa6X4+IxFMsgz+dTNDWlOG4xvhFJIZm8iCWh82s18z2lpTdb2ZHKj1gJWzdhWLZoixHz5yvdTNERObdTHr8jwC3h5Q/6O4dwWtHyPLp1q25FS31HDmt4BeR+KkY/O6+EzhZzcZns+7ltrK1niOnzusRjCISO7MZ499kZnuC4ZzWOWvRPFnRUs/wWIET53QRl4jES7XBvxlYA3QAx4AHZtsQM9toZt1m1t3X1zfbzVW0srUBgJ5Tg5f9vUREFpKqgt/dj7t73t0LwEPAutk2xN23uHunu3e2t7fPdnMVrWitB9A4v4jETlXBb2bLSmY3AAvurJ1KLgT/KQW/iMTLTE7n3A78FFhrZj1m9lvAl8zsWTPbA7wH+HRQd7mZ7aiw7oLQnE3TnE3Ro+AXkZhJVarg7neHFG8rU/co0FUyH7bugrGytUFDPSISO7G8cnfcitZ6fbkrIrET7+Bv0bn8IhI/sQ7+la31DIzkOT04WuumiIjMm1gH/6rFxXP5Xz6p4R4RiY9YB/91bY0AHD4xUOOWiIjMn1gH/6olDZjBIQW/iMRIrIM/k0qyoqVePX4RiZVYBz8Uh3sOv67gF5H4UPC3NfJS34BO6RSR2Ih98K9e0kj/8BivD+j2zCISD7EP/uvai2f2vKRxfhGJCQX/EgW/iMRL7IN/ZWs9dckEL/aeq3VTRETmReyDP5VM8Iarmtj/Wn+tmyIiMi9iH/wAb7o6x/7Xzta6GSIi80LBD6y9Osfxs8OcHtSZPSISfRUfxGJmDwMfAnrd/aag7H7gPwDjT0X/b+6+I2Td24GvAElgq7v/rzlq95xae3UOgP2v9XPr9Utmta2RsQL/9+AJfnb4JK+fGyaXTfPmFYt49w3ttDbWzUVzRURmpWLwA48AXwe+Nan8QXf/crmVzCwJ/BHwfqAH+LmZfd/df1FlWy+bX1rWDMD+Y2erDn53Z/vPXuWrTxzgtbNDpJPG4sY6zp4f4/xonoTB21Yv5v03LmX92qtY096Imc3lx5AFxt0ZHMnTPzRG/9AoZ4fGMINsKkk2naC+LklzNk1DXVLHgsyrmTx6caeZra5i2+uAg+5+CMDMvgPcCSy44L8ql6GlIc3zx6v7gvfM4Cibtj/NPx04wdtWt/LFj9zEO9/YRjadpFBw9hw5wxP7jvP3zx3niz/Yxxd/sI/WhjRvvbaVG5cvYk17I9e3NXFVc7EdmVRywvYLBWd4rMDQaJ6hsTxDo8H06KTpsTyjYw4GBpgZCQMzMKz404zJETM5c0prTF023bo2pcwm1LUy5WWmmbqhS95emc8S0uzyHEbyBQZH8gwMjxV/jowxMDwWhPrFYB+f7h8a49zwGPlC5SvCs+kESxoztDXVsaQpQ0t9mub64jOhm4PppkyKdDJBOmnUJROkUwmSCcPdyRcgX3AK7uQLTt6dQsEZf+vxq9IduHiBelDmpXPzb75/3S303691qQS3vWnpZX+fmfT4y9lkZvcA3cBn3P3UpOUrgFdL5nuAXy63MTPbCGwEWLVq1SyadenMjBuXNfPskTOXvO5rZ4a45+GnOHxikC985CZ+45dXTQikRMLouKaFjmta+MwH1vLK64P89NAJug+fYtfLp3hify+T7xaRufCfGvLujIwVZvsR5TJIJoxcNlV8ZdLksilWtNTTnM0F5ekJP5uyxf9uw8Ev7POjec6eH+X1gRFOnBvm9XMj9PYPcaC3nzODo/QPj005NiTa2poydP/uwg3+zcAXKHYUvgA8AHx8Up2w361lD2N33wJsAejs7Jz3w/3mVS184yeHOD+Sp74uWXkF4Mz5Ue55+CmOnh7ikY+9jV95Q1vFdVYtaWDVklX827cVf7kNjeZ55eQgh/oGOHFumDPnRzk9OIJ78ZeGWfEuotl0gvp0kmy6OF0cLkiSSSfIppMXlqWTxd3uHrwo9vzcPejxTdy1k4PFp11Wft2LPUcPXV52vTL1J7bDQ8rCW11+G1O3F1anXI8wnTQaMyka61I0ZlI01CXJpBKXdYimUHD6h8c4e36UgZExxvLOSL7A6FiB0bwzWiiQMCNpRiIBSTOSCSORCMrMpv7FZhf/ChpfduHnPPe/Jx9Pl/39roBfoqnE/PwbVBX87n58fNrMHgIeD6nWA1xTMr8SOFrN+82HW1a1ki84zx45w7rrFlesP5Yv8Ik/28VLJwb45sfWzSj0w2TTSW5YmuOGpbmq1pfoSiSMRfVpFtWna90UiZiqTuc0s2UlsxuAvSHVfg680cyuM7M64C7g+9W833zouKYFgGdemTxiFe5rTx7k/734Ov9jw5urDn0RkVqYyemc24H1QJuZ9QD3AevNrIPiX8qHgXuDusspnrbZ5e5jZrYJ+CHF0zkfdvfnLsunmANLmjJcu6SBp2cQ/D8/fJKvPXmAf33zCv5N5zUV64uILCQzOavn7pDibWXqHgW6SuZ3AFPO71+oOq9dzJP7j5MvOMkyY23DY3n+62N7WNFaz+9/5KZ5bqGIyOzpyt0S77qhjVODo9Oe3bP1n17iUN8Av3/nTTRlZnNSlIhIbSj4S7zzDW2Ywc4X+kKX95wa5GtPHuCOm67mPWuvmufWiYjMDQV/iSVNGd68YhFP7u8NXf6VfzhAweF3P3TjPLdMRGTuKPgn6XrzMna/enrKg1kO9p7jsad7+M1br2VFS32NWiciMnsK/kk23LyChMFju3omlD/4Dy+QTSf5xPo1NWqZiMjcUPBPsrQ5y21vuoo//eeXOTs0CsCzPWf4wZ5jfPwd19HWlKlxC0VEZkfBH+JT77uBM+dH+Z879jE0mue//81eWhvSbHz39bVumojIrOl8xBA3rVjEf3z3Gr7xkxf562eOcn40zx//+i00Z3XpvIhc+RT8ZXz2g2tZvaSB7pdP8a/espx339Be6yaJiMwJBX8ZiYRx17pV3LVufm8RLSJyuWmMX0QkZhT8IiIxo+AXEYkZBb+ISMwo+EVEYkbBLyISMwp+EZGYUfCLiMSMuXut2zCFmfUBL1e5ehtwYg6bcyXTvphI+2Mi7Y+LorAvrnX3Gd1iYEEG/2yYWbe7d9a6HQuB9sVE2h8TaX9cFLd9oaEeEZGYUfCLiMRMFIN/S60bsIBoX0yk/TGR9sdFsdoXkRvjFxGR6UWxxy8iItOITPCb2e1m9ryZHTSzz9W6PXPFzK4xsx+b2T4ze87MPhmULzazH5nZgeBna1BuZvbVYD/sMbNbSrb10aD+ATP7aEn5W83s2WCdr5qZzf8nvTRmljSzZ8zs8WD+OjN7Kvhsf2FmdUF5Jpg/GCxfXbKNzwflz5vZB0vKr6hjycxazOxRM9sfHCdvj+vxYWafDv6f7DWz7WaWjfOxUZa7X/EvIAm8CFwP1AH/AtxY63bN0WdbBtwSTOeAF4AbgS8BnwvKPwf872C6C/hbwIBbgaeC8sXAoeBnazDdGiz7GfD2YJ2/Be6o9eeewX75L8CfA48H838J3BVMfwP4RDD9n4BvBNN3AX8RTN8YHCcZ4Lrg+EleiccS8E3g3wfTdUBLHI8PYAXwElBfckz8uzgfG+VeUenxrwMOuvshdx8BvgPcWeM2zQl3P+buTwfT/cA+igf4nRT/wxP8/EgwfSfwLS/6Z6DFzJYBHwR+5O4n3f0U8CPg9mBZs7v/1ItH/bdKtrUgmdlK4FeBrcG8AbcBjwZVJu+P8f30KPDeoP6dwHfcfdjdXwIOUjyOrqhjycyagXcB2wDcfcTdTxPf4yMF1JtZCmgAjhHTY2M6UQn+FcCrJfM9QVmkBH+K3gw8BSx192NQ/OUAXBVUK7cvpivvCSlfyP4Q+CxQCOaXAKfdfSyYL/0MFz53sPxMUP9S99NCdT3QB/xJMPS11cwaieHx4e5HgC8Dr1AM/DPALuJ7bJQVleAPG3OM1OlKZtYEPAZ8yt3PTlc1pMyrKF+QzOxDQK+77yotDqnqFZZFYn9Q7OHeAmx295uBAYpDO+VEdn8E32PcSXF4ZjnQCNwRUjUux0ZZUQn+HuCakvmVwNEatWXOmVmaYuh/292/GxQfD/4MJ/jZG5SX2xfTla8MKV+o3gF82MwOU/xT+zaKfwG0BH/ew8TPcOFzB8sXASe59P20UPUAPe7+VDD/KMVfBHE8Pt4HvOTufe4+CnwX+BXie2yUFZXg/znwxuDb+zqKX9R8v8ZtmhPBmOM2YJ+7/0HJou8D42defBT4m5Lye4KzN24FzgR/6v8Q+ICZtQY9ow8APwyW9ZvZrcF73VOyrQXH3T/v7ivdfTXFf+cn3f3XgR8DvxZUm7w/xvfTrwX1PSi/Kziz4zrgjRS/xLyijiV3fw141czWBkXvBX5BPI+PV4BbzawhaOv4vojlsTGtWn+7PFcvimcrvEDxW/ffqXV75vBzvZPin5N7gN3Bq4viWOQTwIHg5+KgvgF/FOyHZ4HOkm19nOIXVQeBj5WUdwJ7g3W+TnBh30J/Aeu5eFbP9RT/cx4E/grIBOXZYP5gsPz6kvV/J/jMz1NypsqVdiwBHUB3cIz8NcWzcmJ5fAC/B+wP2vunFM/Mie2xUe6lK3dFRGImKkM9IiIyQwp+EZGYUfCLiMSMgl9EJGYU/CIiMaPgFxGJGQW/iEjMKPhFRGLm/wPxYbmIjix9AgAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "peaks, info = nk.rsp_findpeaks(rsp)\n", "rate = nk.rsp_rate(peaks)[\"RSP_Rate\"]\n", "rate.plot()" ] }, { "cell_type": "code", "execution_count": 26, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "15.0139897440302" ] }, "execution_count": 26, "metadata": {}, "output_type": "execute_result" } ], "source": [ "rate.mean()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Distord the signal" ] }, { "cell_type": "code", "execution_count": 43, "metadata": {}, "outputs": [], "source": [ "data = []\n", "for noise_amplitude in np.linspace(0.01, 1, 200):\n", " distorted = nk.signal_distord(rsp, noise_amplitude=noise_amplitude)\n", " \n", " raw = nk.rsp_rate(nk.rsp_findpeaks(distorted)[0])[\"RSP_Rate\"]\n", " \n", " # Khodadad 2018\n", " cleaned = nk.rsp_clean(distorted, sampling_rate=1000, defaults=\"khodadad2018\")\n", " khodadad = nk.rsp_rate(nk.rsp_findpeaks(cleaned)[0])[\"RSP_Rate\"]\n", " \n", " # Biosppy\n", " cleaned = nk.rsp_clean(distorted, sampling_rate=1000, defaults=\"biosppy\")\n", " biosppy = nk.rsp_rate(nk.rsp_findpeaks(cleaned)[0])[\"RSP_Rate\"]\n", " \n", " # Results\n", " rez = pd.DataFrame({\"Raw_Mean\" : [(rate - raw).mean()],\n", " \"Raw_SD\" : [(rate - raw).std()],\n", " \"Khodadad_Mean\" : [(rate - khodadad).mean()],\n", " \"Khodadad_SD\" : [(rate - khodadad).std()],\n", " \"BioSPPy_Mean\" : [(rate - biosppy).mean()],\n", " \"BioSPPy_SD\" : [(rate - biosppy).std()],\n", " \"Noise_Amplitude\": [noise_amplitude]})\n", " data += [rez]\n", "data = pd.concat(data)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Results" ] }, { "cell_type": "code", "execution_count": 45, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 45, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZUAAAELCAYAAAARNxsIAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJzsvXt4VMX9+P+a7CVZEiKoUSkIsQ1WLVq5qj8/ooUElFqhasVovREFWrGC1Wqr4gUF/NoKCK2gDaJtCahgQYugKNrWG9ciCmqohgKirNwJSTbZvH9/zJ7NXs5uNtnNlXk9z3l2d86cOZOdzbzPzPumRASDwWAwGFJBWkt3wGAwGAztByNUDAaDwZAyjFAxGAwGQ8owQsVgMBgMKcMIFYPBYDCkDCNUDAaDwZAyjFAxGAwGQ8poUaGilHIopTYopV4NfD5FKfWhUqpUKbVQKeUOlKcHPm8NnM8NaeO3gfLPlFJDW+YvMRgMBgO0/ErldmBLyOfHgGki0hPYBxQFyouAfSKSB0wL1EMpdQZwNfAD4GLgT0opRzP13WAwGAwRqJbyqFdKdQOeAx4F7gB+AniBk0SkRil1HvCgiAxVSq0IvH9fKeUEvgZygHsARGRKoM1gvVj3Pf744yU3N7cJ/zJDoqxbt+5bEclJVXtmbFsHZlzbJ4mOq7M5OhOD6cBvgI6Bz8cB+0WkJvB5B9A18L4rsB0gIHAOBOp3BT4IaTP0miBKqdHAaIDu3buzdu3a1P4lhkahlNqWyvZyc3PN2LYCzLi2TxId1xbZ/lJKXQrsFpF1ocU2VaWec/GuqSsQeVpE+olIv5ycRjxAeb2wZo1+NRgMBkNMWkqncj5wmVKqDFgADEKvXDoFtrcAugFfBd7vAE4GCJw/BtgbWm5zTWooKYEePaCgQL+WlKS0eYPBYGhPtIhQEZHfikg3EclFK9rfEpFrgVXAlYFqNwBLAu+XBj4TOP+WaGXQUuDqgHXYKUBPYHXKOur1QlERVFTAgQP6tajIrFgMBoMhBi1t/RXJ3cAdSqmtaJ1JcaC8GDguUH4HdQr6T4AXgM3AcuBWEfGnrDdlZeB2h5e5XLrcYDAYDFG0pKIeABF5G3g78P4LYIBNnUrgZzGufxRtQZZ6cnPB5wsvq67W5QaDwWCIorWtVFoXOTlQXAweD2Rn69fiYl1uMATwlntZs3MN3nKzLWowtPhKpdVTWAj5+XrLKzfXCBRDGCWbSihaWoTb4cbn91E8vJjCXoUt3S2DocUwK5VEyMmB/v2NQDGE4S33UrS0iIqaCg5UHaCipoKiJUUxVyxmRWM4GjBCxWBoJGX7y3A7wg05XA4XZfvLouqWbCqhx/QeFPylgB7Te1DysTFNN7RPzPaXISZer5dnnnmGsrIyampqguVz585twV61HnI75eLzhxtyVPurye2UG1YWuqKpqKkAoGhJEfmn5JOT2fDVr7fcS9n+MnI75TbqeoOhKTFCxRCT4cOHc8EFF5Cfn4/DYeJ0RpKTmUPx8GKKlhThcrio9ldTPLw4aqK3VjSWQIG6FU1DhYLR4RhaO0aoGGJy5MgRHnvssZbuRqumsFch+afkU7a/jCx3Fod9h/GWe8OERaIrGotYKxGz4jG0BYxOxRCTSy+9lGXLlrV0N1o9OZk5bN27lb5P97XVmVgrGo/TQ3Z6Nh6nx3ZFA/F1Lw3R4dRHIjqe5cuX8/3vf5+8vDymTp0adb6qqoqRI0eSl5fHOeecQ1mdU3C2UmqdUmpT4HWQdUIp9XYg99F/AscJDe68oVXTYqHvW4p+/fqJiXiaGB07dqS8vJz09HRcLhciglKKgwcPpqR9pdQ6EemXksZoubH1lnvpMb1H2PaWx+lh2/htUSuNeCuD+tpJ9D6p6K/f7+fUU0/ljTfeoFu3bvTv35+SkhLOOOOM4DV/+tOf+Oijj5g9ezYLFizg5ZdfZuHChSilNgMFIvKVUqoXsEJEuoIWKsCdIpLwQJn/2dZBov+vZvvLEJNDhw61dBfaBInqTHIyc+JO/vW1E0+Hs3PnTrZt2xZmUDFw4MAG3WfDrg109nQmt1MuWz/aSl5eHt/97ncBuPrqq1myZEmYUFmyZAkPPvggAFdeeSXjxo0j8JBaISJWYNdPgAylVLqIVNXzVRraAUaoGOKyb98+SktLqaysDJbFmqzaM/FWGQ3VmcSivna85V7yOuexbvQ6DvsOB3U4t91xG68sfoUzzjgjaFChlIo5Tnb3OeI7wmULLsPtcFNTW8Noz2hOPrkuAHi3bt348MMPw67ZuXNnsI7T6eSYY45hz549kbe7AtgQIVCeVUr5gUXAI3K0bZe0c4xQMcTkz3/+MzNmzGDHjh2cffbZfPDBB5x33nm89dZbLd21ZiXU4qqqpop7B97LmL5jgsIlUSuw+ojXTqTVV1GfIorXF+N2uDn4l4M8ufRJzul+TlDoWY6WoUIwVDCG3udI1RFqpIYafw1Vfj33P7XmKX6WHh5uT6nw9EV2siC0jlLqB+jU30NCqlwrIjuVUh3RQuU64HmbdsIS6xnaDkaoGGIyY8YM1qxZw7nnnsuqVav49NNPeeCBB1q6Wykn3irEzuLq/lX3M/lfk8PMeUOtwGK1E3kusixyJRJaHtmHWatnAejPneC2V26jY8eO1NTWhAkcy+wYIcoUedv4bWzYtYHLFlxGjb8mrL/OTk62btka/Lxjxw6+853vhNXp1q0b27dvp1u3btTU1HDgwAGOPfZYIJgu/GXgehH5r3WNiOwMvB5SSs1HB5CNEioi8jTwNGidSn1jaGg9GKFiiElGRgYZGRmAtvQ57bTT+Oyzz1q4V6kllt+HNeF/ue9L0lS0kaQVkiXUnDdUZxIqMFZ+sTLsHtMunob3sJfJ/55su/Kw+tG/a39A60CcaXH+VV3AbDj03UPggFlLAwJnmBZAo/4+CqVUmFC66e83sWGM1qG4He7gCsWi5qQadi3fxZdffknXrl1ZsGAB8+fPD573lnv54cAfMvvPsznvvPN46aWXGDRokLVScQD/AH4rIu9a1wQS7HUSkW+VUi7gUmBlwoNlaBMYoWKISbdu3di/fz8jRoygoKCAzp07Rz2ttmVi+X0crDzIhBUTqJXaqMk2lFgOjJHbZbXU4vP7gvcY++rYYF3blQdaEByXcRy9u/TmnbJ3OOSLYzTx/cARA0datONqlb+K3nN6M+OSGdTU1kSdf/LHT3Jy/5MZOnQofr+fUaNG8YMf/IBf3/Nrvs7+mpflZVy1Lg5/eJgV3VfQ/aTuLFiwwLr8BKATcL9S6v5A2RCgHFgRECgOtEB5JnbPDW0RY1JsSIh33nmHAwcOcPHFF+OOTFzWSFrapHjNzjUU/KWAA1UHgmVZ7iyq/dVxhYmFnTnvFu8Wes/pndD19ZHpyqSquooaoif9NNJwOpwg4Kv12VxdR4YjI7hSsfsbpg2dxoQVE3CkOaj2VzPjkhmM6Tsmqm7JphJGLRlFpb8yqo3Q76Glx9XQNBiTYkNK+Pe//01paSk33XQTXq+XnTt3csopp7R0txqFt9zLhl0bADj5mJPZ+PVGjlQfCatTVV2lJ2sb0kjD5XCR7kynqrqK6394Pd8e+TY4mZZsKuGmJTelRKAAlFeXxzxXSy1+vx8/ftiDfub3QlD+KMj+TXZQ2X+w8iC3vXYb1bXVYe24HC76dOnDtvHb6vWhKVpaFCVQrDYaE3LG0D4xHvWGmDz00EM89thjTJkyBYDq6mp+/vOfp6Tt5cuXA/RSSm1VSt0TeV4pla6UWhg4/6FSKjeZ+5VsKqHrH7oy9G9DGfq3oZzxpzO45dVboibZaqm2faIHPZG/ef2bXNj9Qqpqq5izbg5n/OkMbl56My98/ALX//16W4HiVE7SHelkujJj9m/U2aPwOD1x60TiJ5A5++9Af/R/843AD8HxQwcvXvki28ZvA4EJKyZEeeMD+Gp87KvYB0D/rv1jCoay/WVIrf2uRmPMpw3tFyNUDDF5+eWXWbp0KZmZeqL7zne+kxKHSL/fz6233grwOXAGUKiUOiOiWhGwT0TygGlo09QG4y33Mu29afx88c+plur6L4hDhjODFz5+gVdKXwkrL95QzMhFI211Ew7lIC0tDWeaE5/fx9DvDo2qk+XOYuQPRjJv+Dx+d8HvyHBmNKxjNcB3A+87AT8CyqCzpzNAUG8UuvLJcmfhUi5qqeWql66KCtUSmfsly51FZW30KgXgdxf8zqxSDEGS3v5SSqWjHZxyQ9sTkYeTbdvQsrjdbpRSQd+D8vLY2zENYfXq1eTl5fHFF1/4RMSnlFoADAc2h1QbDjwYeP8SMEsppRriKFeyqYRrFl+Tkj4DSK3w1LqnGnSNI82Bz+8LOhuu+GIFrjRX2AqpqrqKS+dfGhR6DuXA7XCT4cygwleBKLEVWMF2nEAtcCzwIZANtYdrye2Uy4ZdG2yt1yqrK1FKhfXNsmaLtFYrHl5MXuc8MhwZUdtf7jS3rf7F0Hi83tiJZuOdS+V9kiEVK5Ul6AmgBm3dYR2GNs5VV13FmDFj2L9/P8888wz5+fnccsstSbcb6okdYAfQNaJaV2A7gIjUAAeA4xK9h7fcy40v35hcRwFFnTNfjdRAA+xaHDhwpbmiyv21fjKcGWSnZ5PhyABF2CrKL37SSOOZS5/hr5f/laeGPUW6Iz2sjQxHBjMvmYnH6cHzEw9pNWk4f+zEtdtF2qY0Hpz+ICu/WMnwBcNtdTM1UmOrX9mwa4NtNsssd1aU8yPAk5c8aVYpKaSkBHr0gIIC/VpSkti5VN4nWVKhqO8mIhenoB1DK+POO+/kjTfeIDs7m88++4yHH36YgoKCpNuNsdiILIyewWym9Fie12X7y3QLSRo3SkgDfvHX6THqIcORwfRLpnP7a7dHnaulltraWu664C4GfGcAl79wedQEXyu1XLP4Gvyi7xe6erGU74W9ChnYYyCrd65mQNcBeMSDd4I3qN/oMb2HrWI9FtV+3Qe7uGCHfYeDXvj1WYkZGofXC0VFUFGhD9Cf8/Pr3tudy8lp2Kpjyxa46SaoqrJvK1lSIVTeU0qdKSKbUtCWoZVRUFCQEkESiuWJHVoEfBVRbQdwMrAj4DR3DLA3sq1Ynte5nXKTFiiNwamcPHDRA8EwLocqD3HXyrui6vlqfUz+12TWjV6HvzZaUEWaCfvFj6pV/KHgDww/bXhY6Ba1Q1GxqIJOaZ3Y+/VeNm7cyKQnJuE+3R3T6ADApVw4Hc6wkDC9u/SOGX+sf9f+caMGGJKjrAzc7rqJHsDl0uVgf27DBli9GiZP1ud9PiguhsIYedtKSuoESijWfVqLUPk/4Eal1JdAFYHnQxE5KwVtG1qAjh072m51pCr0ff/+/SktLQVwK6XcwNVApPJjKXAD8D5wJfBWQ/QpOZk5zLt8HtcsSp1OpT7SHeksvXopQ/J0qKuSTSXcv+p+XMplayTgcrhYtHmRrVCxo0ZquOXVW9i4eyMTB04MblPxCvBz2Feyjzlr5zCm3xg+WfMJvlPj+6/MHDaTy0+/PEpIRMYfm3bxtGDOlvoiLRsaT26uFgqhVFfrcog+V1EBw4eDFeu1vlWHtRKKFCiR90mWVOhULgF6oj1mf4IOvfCTFLRraCEGDx7MGWecwX333cfHH3/MwYMHOXjwIIcOHUpJLhWn08msWbMATgW2AC+IyCdKqYeVUpcFqhUDxymltgJ3AFFmx/VR2KuQ3Xfu5ubeN5PWiJ96uiOdcQPG4U5LzNmzyl/FR998BGidzvUvX0+lvzKm1VlldSWP/vvRBlulzVo9i1VfrgoP3XKMfrnttdvY4t2C2+UOJgazM1POcmdxSqdTghkry/aXBS29CnsVsm38NlZet1I7Ri6fEDeZlyE15OToVYbHA9nZ+rW4WJdHnsvIAKXqBEoooasb0MJkzRq9qrHzW3a56u6TClLmUR/I4Ba0hRSR/6Wk4RRjvHMT48CBAyxevJgFCxZQWVnJyJEjufrqq4MBA1NBc3pee8u9zFk3h4fffjhqElco0lRaUH8BehvrqR8/xfDThvP5ns+5cN6FYefjcc/597B+13pe/+L1uPUs66rG8Otzf82sNbO0X8xC4DxgGXAzONY46K/68/5r7wcdPkcsHBGuJwlsfSFQ4a/A4/QAhAXJbGxSMONRnxx2+hGrLCsLDh+GffvgqqvgwIHo6z0e2LZNX1tSolcnbrdeofj9elUSSkYG/O9/9QuVRMc16ZWKUuoypVQp8CXwDlAGvJZsu4aW5ZhjjuGmm27itddeY+zYsUycOJF58+a1dLcaTU5mDmP6jrGNgwVECQxrq+nE35/Ij577UXC14kpz4U5z41Sxd46nvju1XoECJCRQlK29AsxZPwe/P9DnS4E1wCHgCfB/5ec/ff6Dt9xLTmYOQ/KGMG3oNNId6XR0d9RhW9ICASb9WmBY8c+KlhThLffiLfeyrHRZ1PfV2PTFRzPWSsHrjV9mkZMD/fvXTfJz5sDJJ8PgwdC3L2zdCr17R2+HQfjqJlTxf+CAXtXU1kZf43aHr2ySJRXbX5OAc4HPReQUYDDwbrwLlFInK6VWKaW2KKU+UUrdHig/Vin1hlKqNPDaOVCulFJPBryrP1JK9Qlp64ZA/VKl1A0p+HviE+/X0I547733uO222+jTpw/vvvsuL7/8MnfccUdLdyspyvaX2QoViaPRF4Tq2urg5FtdW42/1q/Ni1NIB1cHHCTet8O+w3UxwTLRnmJ3Ab/R790d3cHJv2RTSdCj3uf3Mf688cGVSSQuh4s5a+fQY3oPxr4ylsO+w2Hnjfd8w7Az3bXKBg/WwmLOnNjXz5kDY8fqVcahQ1pAFBXpc5HbYZMm6RVKfn7s7S6PB9LDrdOD+pRUTW1Jb38ppdaKSD+l1Eagt4jUKqVWi8iAONd0AbqIyPpAsp51wAh0kIm9IjI1ELqjs4jcrZQaBtwGDAPOAWaIyDlKqWOBtUA/tK3POqCviOyLde+kltKha8n6zCzaMLm5uXTq1Imrr76aQYMG4XSGP5X36dMnxpUNo7m3SbzlXrpP694gM9tU4kxzkuHMiJqoQetvHhn0CBNXTYxrsRXGstinHGkOdq3cBRC1hZXhyMBf67fV5cQLPgkw+8ezGdMvvhmx2f7SeL1aeIRabFm6kIqIr/eee8B6ZrO2vkALnUjFeseO8OabejUTuVUWud1VWxu+oklPh0cegYkTtS6lulpPYwcPwu236+tqauyntuYMKLlfKZUF/Av4m1JqN9iEVQ1BRHYBuwLvDymltqCd3YYDFwWqPQe8DdwdKH8+YP3zgVKqU0AwXQS8ISJ7AZRSbwAXA6nXJsYzIk+lO2orIDc3F6UUK1as4PXXXw/zK1FKtdnMjzmZOcwdMZcbFt8QnFCdypnyVYcdrjQXrxa+SmdPZ9bvWs/tr91OVW3dbFHlr2LiqonBiMEuhwtfjS8YNt+Wtegg8z8AOuoit8ONiFDURz/OLitdFpWLxelw6u2+iOdJj9PD7y74Hb9/7/e2QiXTlUmfLql5oGiPRE7wdibCDgfYPcdPnQp/+IMWOB6PFgS/+12dcAjF56sTOpYS37p/5BTlcun2RPT2V1qaFijTpkGfPrqdxYv1agjq7pXM1JYKoTIcqADGA9ei7VASDtESCBTYGx1g4sSAwEFEdgWU/xDiXR3A8sCOVZ564hmRtzOh8vbbbydU74033ki5D0tTY2VotKIV9+7Sm4f/+XAwn0kyuNPc2kHSRqHvTHPSu0tvcjJz2Lp3q61rZ5pK45ROpwQzMgJ8uf9LJqyYQJpKi/aMvxP4JHCkgfNMJ/N/O5+Bpw1k5Rcr6TG9B840Z1Qulmp/tXai9NWtVDo4OzA1fyr5381n8r8m2/595dXlLN6y2Pio2GC3iZGfH633sFOUW1jl1jWPPqqFTCQzZujXNWvClfl2U5THA888A9dfrz9b526/XW+PWe8jcTobP7UlrVMRkXK0k9pFIvIc8GcgIZOWwApnETBeROLZqsbyrk7Y61optVYptdbb2A1DOyPyqiptjnGUcvfdd7d0FxqFpbwekjeEnMwcZl4yk3/f9G/bkCoNIU2lxczQOG3otLDUwHbRjMuryxm+YDi3LruVEQtHcNVLVzFhxQSmXTyNxVctjtaDdEBHJ74RHD91MPg7g/nFj3/B4089HvRhCRUoWe4sPE4Pjwx6JGr1c6TmCPe+dS8/fOqHcQ0Ipr479ag2LY6ldA9ViFt6j8WL9VaShdsN06fr1UIiuN16teLx6GkmPR1mz9Y6FDudTCw/l06d9LZbKFVVWtk/Z469mXHoaqihpML66xZ0wD9L3dQVHYy7vutcaIHyNxFZHCj+JrCtZelddgfKLe9qC8sDO1Z5GCLytIj0E5F+OY1dVYQailuaLqW0OYZd4JyjQKHfnhK8nd/9fCZeOLFR13Z0d8Tj9HDvwHtthUoHZ4fgtlHZ/jLbEPQWlf5KXtz8YljsrQnLJ9C7S2+KhxfbKvP5Cvzv+3n95dfxfsfL9P9Nt92+qqqp4qen/ZT73rovWJbhqJttDvkOUS3V9ZpOh1qJHU3Eipc1Z060jsTh0CuA0FVJWhqccgp06JDY/aqrYcwYrXx/6y3Yvh0uv7xOgB06pIXD2LF1guX22/X01LFjnSVYLEuxqiqtX7E7N2NG4zdgUmH9dStwPnAQQERK0Tu9MVHaXbsY2CIiT4ScsryoCbwuCSm/PmAFdi5wILBNtgIYopTqHLAUGxIoaxoKC+E3v6nbeKysrHss2bKlTog0ZbS2VoSd131bZkzfMVGrgfS09KhgjqFkubOYeclMto3fxpi+Y2y94wUJWkzldsptsG+KI81B2f4y8k/Jx+UIWU2tQj/KvQ/0ALlFYDhUH2e/v1JdW838j+dT5a8KrpT84m9QDheLo820ONZq5N13tdVVJNXV0SsAhwP277efxCOJdHy0TIzLyvTWVCS/+AV07ap1M1VVcOSIFjBnn62vmTZN3z+Sqqrw1RTov2tMEiHdUiFUqkQk+DUF4jTV9wh7PnAdMEgp9Z/AMQyYChQE/F4KAp9B27l8AWxF57T+JUBAQT8JbaW/BnjYUto3CXPmwEMPRZfX1urHgYIC6N4dbrwx+tfXjlcs7YWczJygF3p2ejYep4dnf/osz454Fo/TQ5Y7eqvTX+tnWM9hwfAlc0fMDfNhcTvcTLt4Ght2beD1rdp3pXh4MS6V+FbbYd9h1u9aT9n+MtKdIQLuHaAS+AZ4Ey1g/hRyJIBCxc4wWQrMBGagzXBCqPZX08XThZEjR5KXl8c555xDWYizg1LqtwEXgM+UUkNDyi8OlNkmZ0s1sTYMIsvrq2dnnisCP/qRvZCYMCF6si4v17qNwsI6U2BXjJ/BvHn2hqW5ufZhVkTCV0V+vxYwZ5wBF14I48fbGwhYdUOZPz+56SoVivp3lFK/AzxKqQL0hP9KvAtE5N/Y60NA+7lE1hf0isiurbnA3Ab1uDF4vfYaLagbZbvRhnar0M9NVbCgVoSlyI+Mh2WVrd+1PmidZQVhDFNai36Kd6OV9teedS23/eO2oLWZ2+HmyUuexOlwUl2TeHiWCSsmsG70uvBVzvjk/96Y+e1r0Y9y1wHZ6Ee57wMnaCux4uHFLF2wlM6dO7N161YWLFjA3XffzcKFC0FH1rgabZf2HWClUurUQMt/RD8w7gDWKKWWishmmoBYyvM5c3QARqdTl197ra5rWVrde69+Ul+8GH71K13P74+elO1CpIBu50c/gtJSWLQo/FxVFcydC48/rq+3W+VYrFmjdSnbt+sVTqdO+tn1jju0wEiUyK25+kh2ukqFn0oaOkvfELSgWAH8uSHB/5qTRtu8//a3DRvJUELjJrQBFi9eHPf85ZdfnpL7tFV/Bm+51zZSr11YEzvSHek405xxc9BHkp2ezcrrVrL+6/WMfXVs2DmncqJQYX4nDuXAkeYgjTQq51SSMSY6wRbobJaVNTaz43a0Qf91gc+BlcqY8WOY9KNJ5GTmMHToUB588EHOO+88ampqOOmkk/B6vaSlpe0E/igiUwCUUiuoS7j2oIgMDZT/FsCqF4tExnXLFh2td8AAOP10/fnss8NXEQ6H1mvEsr4KxeGIfoJXSusr3G4tEPz+6DqhdePNgE6nvofdc2haWp3aNlIguN3w5JNw222J/R2NIdZ01Wx+KiJSi36OeSbZtlotc+Y0TKC4XPpXE+pd1EYECsArr+iF5u7du3nvvfcYNGgQAKtWreKiiy5KmVBpq8SK1Gsp4esTKg7laLBeJdSTvaO7Y5hVV4ZL51gJTfXiSnOxfsx6DvsOc/3C63n+pueDqywrH4rlbGnLQfQKxSIbHF85ggIFwpOtOZ1OjjnmGPbs2QPgJrapf2T5OXa3j5Unx8LrhVWrdMiSZcu0bsPirLNg06boST2eEIjErp6ILr/wQli+PH5b9T1S19REb4+BFjQOR+zVhc+nt9ZmztSvaWl6972oSE9TyQoah0PrX5KZrhotVJRSH8U7325C38fb9oqkQwf9a7LW2ZYheO/eTde/JuDZZ58F4NJLL2Xz5s106dIFgF27dlm55Q02ZLmz7J/6I6iVWn7R7xc8ufrJeuuGBnq0JvPI1MLV/mrcDneYmbLb6eaw7zC5nXKDhgL9u/YPhrrPcmdx2Hc46GzpSHPYevqHkndsHht2bQj629htRsQx3hDsdbj2Gdti5MkBvVX185/bx7EC+CjuzJQc1dXwStzN/eRwOOx9U0JRSpsIr1unt8ZATzETJ+opZ/58eO656OsGD9bC1+HQAq22NloI+f1aWGVnNz5YSDKK+lr0s9FfgKvQ4e5Dj/aB5VGUCH6/FvOFhbByJYwYoUOJtlELsLKysqBAATjxxBP5/PPPW7BHrZeSTSX0fbpvMMS+x+nB4/QwbsC4MKV8Gjoa8rP/eTahcPwKFTZ52xkTzLhkRpSgqaquYtHmRfSY3oPSvaVB3xLL+bLv030p+EsB45eP5/ZzbmfRzxYx+8ezyXBk0MHZQa9SQj3HDsJnvs8Y+rehdJvWjZKPS8KSrdXU1HDgwAErirWPJFwA4uH16iRTsQRKWyfWllgoR44coMNAAAAgAElEQVTAL38JP/wh/OQncOWV2l9l8WIYMgRiuY9Zq7l77tHC6Ne/tq+XrG1RUjoVpdRpQCFaiGwG5gOvB3KKt0oavO9uF8AHYOhQvf6ONP3wePQjRN++0a6t69bpuNWJ5PxsBYwbN47S0lIKCwtRSrFgwQLy8vKYOXNmStpvqzqVSOx0KemOdDaM2cDpOacHw8+/+vmrzFxT/3fXwdmBIzVHwsoyHBn8b8L/gquVSL1OycclYeFnwpgNjNWCbt3odfR9um/UFp2lfLcMEtZsX8O4YePwFHk4knFEb25fQdBZwOP0MNEzkbLPy5g9ezYLFixg8eLFvPDCCyilPkE/cA5AK+rfROdcUsDnaGOcnWiLzWtE5JN430fouK5ZAwMHxlaSt1bq07Gkiscf19tz8b4jK9T9hg16GrMjO1s/F/fvX1fWLKHvReRTEXlARPqgLb6eByYk02arI9TpMStLP0q4XPDBB/p8pPG3iNYY2tkfWmbHbWTlMmvWLMaMGcPGjRv5z3/+w+jRo1MmUNoTdg6N6c704HZSTmYOvbv05un1T9fbVoYjg6n5U+no7hhWXumvZM66unC2OZk59O/an5zMHGbNmkW/Tv10fhQ7fqpfXA4Xq3eutnW+tBwaAfp37c8vz/0lf/vz3+j4QkfUH5W24zoBeAv4VEcP+L/h/8eePXvIy8vjiSeeYGqd3rESeAH9oLkcuFVE/IGHzXFoY55gcrZ6v5QQcnObZ3JOJQ4HZDbcFSiK4cPrr3PXXfDOO/G/o8pKeCLgHRjLpDmZTJBJKeqVUl3RpoM/BfahBcrLybTZKiksrNORjBhR54NiR2Ul5OVFr2Csx4ZURGxrRi6//PKjXjFfH3YOjZEh4i3BExmexaEcuB3uMBPl/FPy+c0bv4m6z+R/TWZM3zFRRgJff/01g/5vEDXZNXAWkEe4wf6JdX0a0HVATCMBy6HRar/w8kLyh+aHR3bWNhvUSi3fP+n7vPjii7ZticijwKM25cuIG185Pjk58OyzcE3zZYlOGqdTb1AkQ3q61oksWVJ/3fvu05738Z5bp06FP/5Rv3e59ErK59PPzpCcbVGjVypKqXfQqxMXOmT9DcA/0HnHU5cesLWQkwOdOyemX1m1KjzZQXp63WhZROb8bIV88MEH9O/fn6ysLNxuNw6Hg+zs7PovPMqw03NE+q/kdsqN0nsA/HHYH4Ope7eN30Zhr0JyMnO4d+C9UXVjebE/8sgjrPt4HaqPgv8ATwIrgRA3YKtPp+ecTvHw4rDwLBZ2uVIsh85QvZDb4Y72z2lGCgth9+7YWzctgdsdnafEoj4dSSJhW6qqQBvW1Y/DEe0fY8ehQ3pF4nBo44PNm/UqZ9u2JDN6iEijDnSGxy8Dxxchx5fAF41tt6mPvn37SqPZvVvE4xHRq8vYR0aGyIoVIps3i6xerV8jr/N4dHutmL59+0ppaamcffbZUlNTI3PnzpXf/e53KWsfWCutZWxTwO7Du2X1jtWy+7D9uM7fNF88j3gka3KWpE9Kl9lrZ8dty/OIR3iQ4OF5xBOzbRGRKS9NEcd5DlHHK6EforoocfyfQya9Mynqut2Hd8ukdyaJ5xGPZE/JloxJGbb1QuuvKF0hK0pXxO2DSPOO6+bNIvPm6dfZs0XS00U6dBBxOETOPVfE6az/3zXySEtrWH2nU2T+fJGf/az+ukqJuN0iHTvqKWDUKN3n9PT6r3W5ossuv7zhf1/kkZ2tpykRPSWtXm0/NSU6rikb+Jg3gB809T1S9QNNiPnz9a8hK0v/OpSyH6nMTF1v/vzw67Kzw8tbMdZ3deaZZwbLzjvvvJS1396ESiLUJ3hCsYRQ9pRs8Tzikfmb7H8zM2bMkD59+siQIUPklim3iPtBt2Q+minuh9xyfNfj6+3PpLe1cDlmyjFx75MoLTmukZPi7t36+W7FCpHHH9eTt93kHDl511cn8igsTEwwgBZ66em6P4k8o8Y6srJEJk1qnOC0e761pqhjjrGfohId16Q96utDKbVetCK/VZASC6E5c7TvisOh7fviYZlaWEmjQ7P4tHIGDhzIypUrufnmmznppJPo0qUL8+bNY+PGjSlpv71YfzUlsbz3Q5k4cSJFRUV0OL5DtBXa3nS2T94e81o7yzWP08O28dsavb3VmsfV69Wq0eHD41uQWf7LPl/iDpMeT8NColje+YcO1V/XjlhZJOvD8ssO9c3Oz482co30rG8W668EaV+hbL1e7R1khQKtj8rKurjUoeFG2wB/+ctf8Pv9zJo1i8zMTLZv386iRDZrDSkj1MorksrKSqZPn87evXtZvnw5W7/dGm2F1iU9TA/jLfeyZueaYNh6O8u19hyBOCdH+3LMnasnzcg8IxYeD/z97/D++zrJVSwrqVDsPOTj4XLVH7E4lp7G6YSRI+0jFteH2w2PPab/rr//XQsUO3e8xqp9UxFQsj6adinU3NilVwOtbbPcVCN/XY88AldcoYMStSF69OgBgMfj4YEHHmjh3hgiueGGG3C5XFxwwQW89tprrNu4Dl/X2FZoJZtKKFpapK3Qaqq4d+C9XHH6FVHWYL4aH/sq9uEt97bbDI+WQWdZmZ7YBw8OV6hXV2sPgJwcvbHwq1/VHwJl8mQdIjBR4VJTo59Pp0+3XzU5nfD883VBzyOvtfOaDyUjwz7T5JEj2vS4qqrOfmjaNPsEX40xK26OlUr7wi69Wno6vPwy7NgBdpOvlWatDfimAJx55pmcddZZMQ9D62Dz5s389a9/ZcyYMbz00kus/WBtTCs0K+Oklfyr0l/J/avup+/TfSnqUxS8xqVc1FLLVS9d1e4zPFobB+efr82ULWPN0FwmVr1QY06PRzsYhjJuHNx5J3z1lfZYtxJl2Rl+Ql3Mrqee0q92K6H0dJ3UK9RNLlEcDli/Xgseu7YtAWrls58wQQuWWN9BQ2iOlUrDIue1dqxf2A031D0CiGh7v5wcHTN78uToR4uqqjbjm/Lqq68C8MeAIft11+lQtX/729/okGjaOkOT4wqZLZyBfZBYoftjBbusqKmgeH0x60avY/uB7YxYOIKKmorg6qVoSRH5p+S32xWLRejKxU7laXc+MjIy6PIpU3R4+rIyLQj62GiUa2v1s2m87a/aWn2v/v31vZct0+FZInfd09O1kAqdcjIztW9MqIvdZ5/psP52OhyXS/dz27bk1b5JCxWl1F+AfwL/EpFPI8+LyLnJ3qNVEKpkz8/Xa1NLqPh84QKjuFgHKIo0UG8jeVWsba93332Xd0PCv06dOpXzzz+fiRMbl3bXkFo2btxIdnY2lrFNRUVF8LNSioMH64J3xcs46XK4OOw7TGdP5yjBE+kQ2Z6xsiwmev7002PvaIfWvfdeuP/++PfOyNDTiWUU4HZHr5aGDbOPeWYXgDJ068rSI/XuHTsumFW/vu8gEVKx/fUs0AWYqZT6r1JqkVIqwbC+bYTI9MBz5tjnCrW0WoWF+tEgUsuWTOyDFqC8vJx///vfwc/vvfce5eWJ5wAxNC1+v5+DBw9y6NAhDh06RE1NTfDzf3f9N6iQtyzIpg2dFtfpMZHIAIaGM2aM/RZYKErpcP0rVuhjx45oB8ScHG1gELqd5XJpoZXI1lXoNp5loODxJLfVZUcq8qm8FfCu7w/8CBiLjhQ0I9m2WwWhyamt9eWjj0Y/Hhw+rDcxc3PrVjTPPquvbaN5VYqLixk1ahQHAiFpOnXqxNy5TZ9k05AYlZWVzJ49m61bt3LWWWcxatQonE5nmEL+iO8IKO0FXyu1TL9kOt5yL5P/Ndk2e2Xx8GKKlhTFzGxpBccEgiHwDfGxJvOiIr1TXllZNyWEhkWJt/KxCN3OeustnbDr97/XmyXTpmkdDERn27A2WvLz67a4rKySdvWTIRWZH98EMoH30fnh/i0iu1PQtyahwTbva9boFUporK/sbL25GZm4yzJuj8xf2oZ8U+w4ePAgIsIxxxyT0nZbsz9DW2DkyJFh1l89evTgvsn3xc0+6Xa42TFhB0BM/5dYvjElm0q44eUbwlIjzxsxj8Je4Y/UZlzt2bJFT96hu+Lp6VpANNQw1C54ut30U1hon1Y5Xnksmi3zI/AR0BfoBRwA9iul3heRBrrktFLsrL2qq3US6j/+MVzrVV2tD2uki4r0Y0Fo/Og2xj/+8Q8++eQTKkNsHo1OpXWwefNmNm3aBEBRUREDBgyoN/ukz+9jw64NDMkbEnOVYZfZ0lvuZdSSUWGh9X1+31GjyE8Fhw/rbadIodKYYJN2ng2R089NN0H37tEbLUVFOtWyXXkq7IiS1qmIyAQRGYiOVLwHrWPZn2y7rYbI0Pfp6Xqd2bt3/QbpbSBoZDzGjh3LwoULmTlzJiLCiy++yLZt21q6W4YAdtZf8RTyyVC2vwxHmiOqPE2ltVtHyVQT6/k0UTWr16s3Trxe+7YiqarSz76RuFz22TkcDm1h1tjkXBZJCxWl1Dil1EJ0fNQRwFzgkmTbbTV4vTqU/cMP14X0nDBBZ7AJNV7PyIgeJesXE/praEO89957PP/883Tu3JkHHniA999/P5jpz9DyWNZf2dnZdOzYkY8++ojvdfkeMllgik5vHIlLuejdpf4N9EjP+9xOufhro+OV1EqtUeQniJ2/S6Jq1khboUSmHwhfuYSWDRgQLZQOH4bbbktBuqdkg70BdwHnAM5UBpFrqqNBQQetCGsdO0ZHYXO7dRS20Ah2kUEjH39cZOxYHbU4VpS2eGFBW5gBAwaIiMg555wjO3fulMrKSsnLy0tZ+xyFASWbCytw5ew1syVjUoZkPpopGZMyEgoWOf+j+bYBJudvmi+uh1zBqMnuSW7b9sy4xqeh//J2wdGtIJCR049dUMuMDF0eGcs2NDZurCCToSQ6rikbeHReuO7Wkap2U30k/ANNJMz9ihXR11mxuK+5xv4aKyx+ImFBW5iHH35Y9u3bJy+99JKceOKJctJJJ8l9992XsvbN5JMc5eXl4vP5gp8//fRTeeKJJ2Tx4sVh9RoSGbm+kPuJhMA345paVq/WU0SscPWhbN4cLVg8nrosHJGCYvduPV1FPjfbtd9sQgWdn74UKEfnUqkFPkm23aY6Ev6B2o1kfULFEhKZmfGvy8zUwsXtjh79VrhiERGprKyU/fv3p7RNM/kkxwUXXCCff/65iIiUlpZK586dZdy4cTJo0CC5++67G9Xm6h2rpeP9HYXvIhyL8F2k48SOsnpH9Aw2b948ycvLk7y8PJk3b56IaEGH1ql+CnwCTJW6ueJGwIveKv8PcLOYca2XeCsVOxqaZSPR9hP9f02F9dcjwLnAShHprZT6EZBM3rDWQX2aMJcr3Ljb64VRo+LH07aI5UDYSjzuFy9eHPd8MumF9+7dy8iRIynTBgw9lVKdRWRfaB2l1NnAU0A24AceFZGFjb5pO2Xfvn307NkTgOeee47CwkJmzpyJz+ejb9++oTnjEya3Uy4VqyrgFOAC4F9QsaqC3N/khtXbu3cvDz30EGvXrkUpRd++fbnssstI1w6/34jIaUopN/CmUuoSEXktcOlCERnX+L/66CPUz8WKGRZPF1NfyJl47afCpS4VQqVaRPYopdKUUmkiskop9VgK2m1ZIr9pS9vldutYCnPnhudImTQpMYESj1bicf/KK68AsHv3bt577z0GDdKJyVetWsVFF12UlFCZOnUqgwcP5p577kEpdQi4B4gMHnEEuF5ESpVS3wHWKaVWiEj7sSpMASrEAfett97irrvuAsDtdpOW1jgbnJzMHHJ25LBv5D7c6W58/Xx0Xtg5ymR4xYoVFBQUcOyxOnN4QUEBy5cvp1A7OhwCEBGfUmo90K1RnTGEIRL+Go+GhltpqCCKRyqEyn6lVBY6/tfflFK7gQZmFmilRH7TEP6tW95DaWmxVx8N4fbWEd3m2WefBeDSSy9l8+bNdOnSBYBdu3Zx6623JtX2kiVLePvtt62Pe9AWg2FCRUQ+D3n/VeA3lUN7MlVPAWeddRZ33nknXbt2ZevWrQwZMgSA/fuT+5qO7DvC/+7/X9AB8tQ5p0bV2blzJyeffHLwc7du3di5c2dYHaVUJ/T2eGh0jSuUUgOBz4EJImLMCevBCuoR+szaFLFpUxH3C1IT+2s4+slyArAc+C/6h9RsKKUuVkp9ppTaqpS6J6WNhybWCn0fGr4lnkBpSBadmTNTYM+XOsrKyoICBeDEE0/k888/j3NF/XzzzTehbVajDTxiopQaALjRvytDCM888wzHH388ZWVlvP7668EI0ps3b+bOO++Me21+fj69evWKOpYsWQLETw4GWDqSMEJXTkopJ1ACPCkiXwSKXwFyReQsYCUQMyOIUmq0UmqtUmqtt42Z4qeaVCbQag5SEfvLmlFrgeeUUg7gauBvybadCIH7/REoAHYAa5RSS0Vkc4Ma2rJF50TZvVvHlq6u1okWhg+3F9+xknWFkpUFixbBqlXRIV3ssIRTKwmRf9FFFzF06FAKCwtRSrFgwQJ+ZOdNFUF+fj5ff/11VPmjjz7aoPsrpboAfwFuEBGb+Kx68gFGA3Tv3r1B7bd1PB4P99yjn6F8Ph8ff/wxALln5OLq4YqbZGvlypUx2z3xxBPZtWsXXbp0YdeuXZxwQrTc79atW+iKkx07dnDRRReFVnkaKBWR6VaBiOwJOf8MEHObXESeDrRBv3792leivwaSrNNks5OINt/uQCtRfwvMAoag0waPA7YBSxrbbiP6cR6wIuTzb4Hfxqpva0kyblxsS620NHvziURMji0TikmT4teLPGLZC7YAixcvlvHjx8v48eOjTFUbw6mnnipfffWViIgAG4HPJPbvaz3wM7vzdsfRaiX09ttvS/fu3WXgwIFyWp/TRHVSkjk6M8zHpCHceeedMmXKFBERmTJlitx1111Rdfbs2SO5ubmyd+9e2bt3r+Tm5sqePXtERAT4ClgEpEn4mHYJef9T4AMx45oQDbXoagpoapNiYAkwDxgDvAC8AbwDnN3YNhvZjyuBP4d8vg6YFat+1A908+b6J/n0dHv7vciRvuYabSacmamvmT1bX+d0NkyopKfrfrUCvv76a1m6dKm88sor8s033yTdXuiEhV5Z/j+JHlM38CYwPvJcvONonXz69Okjn376aZ2PyTiELtE+Jony7bffyqBBgyQvL08GDRoke/bskd2Hd8u8f8yTa2+4NlivuLhYvve978n3vvc9mTt3roiIbN++XdApxLcQYToMTEGbGW8EVgGniRnXhGlpP+nmECqbQt47gH1Ax8a2l0Q/fmYjVGZG1BkNrAXWdu/ePfybmjev/kne7Y69cogc6dmztVDo2FELmssua5hAsRwkW4Ez5MKFC6V79+5y/fXXy3XXXSe5ubny4osvJtVm6IQFHASOFT1G/axxBH6O1rf8J+So92HlaJ18zjzzTBHRPibHTDlGOy2eoIVK9pRsWx+ThhDLwz4WiU4+iR5H67i2NhId12R0KsFwpSLiV0p9KSI2iSqbnB3AySGfu6GX30Ek3v7sgAH138Hn07lS7KINh5pMeL06LlhVVV0o0qVLE/wzQrDMPG66SYcTbWhc7BTx6KOPsmbNmuCeutfrJT8/nyuvvLLRbR533HG8+eabACilPheRvQAisha4OfD+r8Bfk+z+UUO/fv0oKiriJ1f+hIrSCi2Cv6PPJZtkKzS3vRX52EQmNsQjGeuvHyqlDgaOQ8BZ1nul1MF6r04da9BOdKcEnK2uBhKfyU8/HcYl4Is1YUL9ASHtzDSSoapKO1iGWoM1Y3DK2traMCXtcccdR61dPlNDi/LUU0/xgx/8gOefeZ6zys7CdaKLjld0xOP0RCXZagjeci/LSpfhTAt/9rRSDBsMdjR6pSIi0XGwWwARqVFKjQNWoLfh5orIJw1qZOZMnXTLsv7atg3+8Y+6HPSQmLd7IvGoG0pVlfbUdzj0aunJJxPPqpMkF198cdD6C2DhwoUMGzasye5naBzp6enccccd3HHHHUDsJFsNwcoe6UxzcsgXvgFhUgwb4pEK58cWR0SWAcuSaiQ0l6eVVi1UqCRiw2cX76CoSJc5HPrztdfqlYfLpQXDhRfqpNTxqKyEkSPrPltmzKNGNanp8eOPP86iRYt49913ERFGjx7NT3/60ya5l6HhXHXVVbzwwguceeaZYT4iFh999FGj2g3d8gqlo7sjNbU1Sa1+DO2fdiFUUk4ywXDs4h1MnBj+eerUus9lZfCvf2nfmIZSWQlz5sB99zX82gS54ooruOKKK5qsfUPjmTFDO6q/+uqrKW3XLntkljuLmZfMZFjPYUagGOJihEoskgmGExnvoL7PkoRv16RJcMUVTaLMX7x4MXfffTe7d+8OWnYopTh4sDlVZoZYWJEJevToESz79ttvOe6442xXLolilz3SX+s3AsWQEKkI09J+ycnRAmXDBnj99dQqxy2F++LF9acljofPB716we9/n7q+BfjNb37D0qVLOXDgAAcPHuTQoUNGoLQiPvjgg2CAzw0bNgRDrZx44oksX7680e3mZOZQPLwYj9NDdnp20gp/w9GFWanEY84cuPVWHZUYtIJ83rzkleNz5ujgkU5nagJR1tbCXXfBe+/psDJ5ebqvSYYbPfHEEzm9hcyZDfUzbtw4Jk+ezIEDBxg0aBCvvfYa5557Lp9++imFhYVcfPHFjW67sFch+afkJ63wNxx9GKESizlzYOzY8DKfL/m4XKHtWr4sqeLll/UBWmA5HDBjBlx+eYO28ax8Kv369WPkyJGMGDHCypMBJJdPxZA6ampqgpGJJ06cyLnnngvAaaedlpL2czJzjDAxNBgjVOzwemOHoU9La3wirXjtWrhcWseSzJYY6OtrarQAu+026NChzgzZ0hVlZcHhw1HCxsqnAtChQwdef/314GellBEqrYTQnCkejyfsXDI6FYMhGYxQsaOsTD/p260kamsbHx7Uco60a7djRy0ErEn/iSfg8cfrtt6SoboaDhzQ76+9VgsupXQ/LGfNJ5+EMWMAmDRpEt262edVChU4hpZl48aNZGdnIyJUVFSQnZ0N6NBLlckmjDMYGolR1Nuxfr29rsPlSi7PZm6u/Qrk8cfhzTe102VhoW7/jjtS651vIaJXLJZg8/n0MXasXkV5vQwePNhK9xvGs88+y/jx41PfJ0Oj8Pv9QQOKmpoaDh48GPxcHepjlSDeci9rdq7BW3505y8xJIdZqURixe+K5Fe/0v4gyTgahvq/WM6QM2YEVwhhJJKvJdU8+STMnMm0n/yEgvPOY9k999CzSxfo1Ikp//wn85cs4Z133mm+/hiaDcuD3u1w4/P7KB5eTGGvpovWYGi/GKESid1knpUFP/95ajzXE/V/iRXyJT099Qr+UEQYtnQp6cAl48fzd+DP6ABr/xw/ns6HDmlTaJ8Ptm5NmaWZoeUwQSMNqcQIlUjsJnO/P7Vp1hJJBm3n1T9tGvTpo4Xcpk1w/fVNJmAGo5PlXAT8f+jkJhnTp8P06dGV09L0yuvXv4ZLLzXCpo1h50FvBY00QsXQUIxQiSSZEC2pJt6q5vTTtbArKtKTus8HP/mJdqZMko7oNJ4CVKEFygmBzwqdBCWM2lp9TJ0anjZZKW3wMHky1JMz3dBy2HnQm6CRhsZiFPV2FBZqpfnKlbBunX7qboZQ87bk5Og8LnZCzernqlWwcycsWgSzZ+vJPBSltOBJkENowXEI8AHlIZ8b5E8vooXyXXfBNdc05EpDM2I86A2pxKxUYpGTo4VKUVGzhZpvFJFbaWPGaGfHVavgm2/0Suf44/Vq55134J57UmOm3FBKSvTWnVmxtEqMB70hVZiVSiy8Xi1QKiq0j0dFhf7cUiuWhpCTA1ddpZ0eTz+9brVz552wa5cOtb9woX7dvVsfv/pV9Aon1dx7b9v4/o5ScjJz6N+1P44qBwUFBfTs2ZOCggL27dtnW/+5556jZ8+e9OzZk+eeey701PeVUp8ppf4TOE4AUEqlK6UWKqW2KqU+VErlNvkfZWh2jFCJhV0WRytRV1smJweGDNFCZ8iQupXOjBl6ZbNiBTzzDDz6KNx8s1bApwq3u+1/f0cBU6dOZfDgwZSWljJ48GCmhurJAuzdu5eHHnqIDz/8kNWrV/PQQw9FCp9rReTswLE7UFYE7BORPGAa8FiT/zGGZsdsf8XCzgoskURdbRlL4IQyebKO0rx/v9bbLFsGb7+tVzUNdbBLtRWdoUlYsmQJb7/9NgA33HADF110EY89Fj7/r1ixgoKCAo499lgACgoKWL58eTBLaAyGAw8G3r8EzFJKKZFkcj8YWhtGqMSiNVmBtSSRgmbCBL2FFRo7zOfTUQhqarS117HHakG0bRssWaJ9a/z+o/P7a4N88803wVwtXbp0Yffu3VF1du7cycknnxz83K1bN3bu3Bla5VmllB9YBDwSEBxdge0QTAN+ADgO+LaJ/hRDC2CESjySSdTVnrHzszn//PDP116rXy0BZL6/VkV+fj5ff/11VPmjjz6a0PV2i4uQIJZfiMgPlVId0ULlOuB5tEV6VFN27SulRgOjAbp3755QnwytA3W0rTyVUl5gm82p42k7T0xtqa8Qu789RCRlkibG2Lal76q19LUX8BlQDbiA7wMfR9Q5FsgJ1APogbY630vIuCqlbgT6icg4pdQK4EEReV8p5QS+BnLq2/5qB+MKbau/Sf2/HnUrlVhfilJqrYj0a+7+NIa21Fdovv7ajW1b+q5aS1+VUo8De0RkqlLqHuBYEflNRJ1jgV1AQaBoPdAP7crUKVDHBVwKrAzUWQrcALwPXAm8lYg+pa2PK7St/ibb16NOqBgMhnqZCryglCoC/gf8DEAp1Q8YKyI3i8hepdRX6LBwAA8HyjKBFQGB4kALlGcCdYqBvyiltqJXNFc3359kaC6MULMKwM0AACAASURBVDEYDGGIyB50+LfI8rXAzSFFeyKfaEWkHOgbo91KAgLK0H4xfip1PN3SHWgAbamv0LL9bUvfVVvqK5hxbQhtqb9J9fWoU9QbDAaDoekwKxWDwWAwpAwjVAwGg8GQMo46oaKUujgQ7G5rwFwy8nyrCXqXQF9vVEp5QwL33WzXTnOglJqrlNqtlIr0Z7DOK6XUk4G/5SOlVJ8U39+MaxNgxjVxzLgGEJGj5kCbOP4X+C7gBjYCZ0TU+SUwO/D+amBhK+7rjcCslv5eA30ZCPQBPo5xfhjwGtqr+lzgQzOuZlzNuLa/cT3aVioDgK0i8oWI+IAF6CB3oQwHrDjeLwGDlWrqmPC2JNLXVoOI/BPtexCL4cDzovkA6KSU6pKi25txbSLMuCaMGdcAR5tQCQa0C7AjUGZbR0RqACvoXXOTSF8BrggsT19SSp1sc761kOjf01Rtm3FtGsy4RvQjwFE7rkebUEkkoF3CQe+amET68QqQKyJnoT2Xn4u+pNXQlN+rGdeWw4yrxoxrgKNNqOwAQp8OugFfxaoTCHp3DPGXiU1FvX0VkT0iUhX4+AwxPJlbCYl8903ZthnXpsGMa0Q/Ahy143q0CZU1QE+l1ClKKTdasbc0oo4V9A4aEPSuCai3rxF7nJcBW5qxfw1lKXB9wKrkXOCAiOxKUdtmXFsOM64aM64WLW2F0AJWD8OAz9GWGvcGyh4GLgu8zwBeBLYCq4HvtuK+TgE+QVuarAJOa8G+lqCj1lajn3KKgLHoAISgl9N/DPwtm9Dh0M24mnE149rOxtWEaTEYDAZDyjjatr8MBoPB0IQYoWIwGAyGlGGEisFgMBhSxlGXpOv444+X3Nzclu6GAVi3bt23ksIc9WZsWwdmXNsniY7rUSdUcnNzWbt2bUt3wwAopbalsj0ztq0DM67tk0TH1Wx/JYLXC2vW6FeDwWAwxMQIlfooKYEePaCgQL+WlLR0j9o8o0aN4oQTTgD4gd15pdRFSqkDISHCJzZvD49evOVe1uxcg7fcPEAZGocRKvHweqGoCCoq4MAB/VpUZFYsSXLjjTeyfPny+qr9S0TODhwPN0e/jnZKNpXQY3oPCv5SQI/pPSj52DxAGRqOESrxKCsDtzu8zOXS5YZGM3DgQI499tiW7oYhBG+5l6KlRVTUVHCg6gAVNRUULSkyKxZDgzFCJR65ueDzhZdVV+tyQ1NznlJqo1LqNaWU7TaZIXWU7S/D7Qh/gHI5XJTtL2uZDhnaLEaoxCMnB4qLweOB7Gz9Wlysyw1NyXqgh4j8EJgJ/D1WRaXUaKXUWqXUWq/Zlmw0uZ1y8fnDH6Cq/dXkdsptmQ4Z2ixGqNRHYSFs2wYrV+rXwsKW7lG7R0QOisjhwPtlgEspdXyMuk+LSD8R6ZdjhH2jycnMoXh4MR6nh+z0bDxOD8XDi8nJNN+poWEcdX4qjSInx6xOmhGl1EnANyIiSqkB6IefPS3crXZPYa9C8k/Jp2x/Gbmdco1AMTQKI1QMzU5hYSFvv/02QLpSagfwAOACEJHZ6LwYv1BK1QAVwNViwmk3CzmZOUaYtHO8Xm1rlJurn5UjPyeLESqGZqck4OujlFovIv0iz4vILGBWc/fLYGjvlJRorwinU9sgXXutLnO79efi4uR3+I1OxWAwGI4CQt3uDh2CqiqYOzf1bnhNKlSUUnOVUruVUh+HlB2rlHpDKVUaeO0cKFdKqSeVUluVUh8ppfqEXHNDoH6pUuqGkPK+SqlNgWueVEqppvx7DAaDoa1SVqZXKPFIhRteU69U5gEXR5TdA7wpIj2BNwOfAS4BegaO0cBToIUQes/9HGAA8IAliAJ1RodcF3kvg8FgMGDvdhdJKtzwmlSoiMg/gb0RxcOB5wLvnwNGhJQ/L5oPgE5KqS7AUOANEdkrIvuAN4CLA+eyReT9gBL3+ZC2DAZDA0lF3K9A+J1egd2DeyLPK6XuUEptDuxGvKmU6pFElw0xsIuBm5MDM2bEv27atOSV9S2hUzlRRHYBBF5PCJR3BbaH1NsRKItXvsOm3GAwNJBUxP3y+/3ceuutAJ8DZwCFSqkzIqptAPqJyFnAS8D/S7LrhgjixcAdMwZmz9bbXJFkZUGfPtHlDaU1Kert9CHSiPLoho3XdZOwZ88eZs6cya233sqtt97KrFmz2LPHuJO0NA1dcaQq7tfq1avJy8sD8ImID1iA3oEIIiKrRORI4OMHQLcG3cQAxM7GkUgM3DFjYONGSE8Pv9bvT00EqpYQKt8Etq4IvO4OlO8ATg6p1w34qp7ybjblURiv69SzZcsWevXqxbp16zj11FPp2bMna9as4cwzz+TTTz9t6e4dtTRmxZGquF87d+7k5JND/1Xr3T0oAl5r0E0McVciicbAPf10ePbZpolA1RJ+KkuBG4CpgdclIeXjlFIL0Er5AyKySym1ApgcopwfAvxWRPYqpQ4ppc4FPgSuR8eJMjQD999/PzNmzOCqq64KK1+0aBH33nsvixYtaqGeHb2ErjgqaioAKFpSRP4p+XEdGlMV9yuGf2qs3YOfA/2AC2OcH402wqF79+4N6kd7JnQlUqGHmKIiyM/XAqEhMXALC/V1qXR8hKY3KS4B3ge+r5TaoZQqQguTAqVUKVAQ+AywDPgC2Ao8A/wSQET2ApOANYHj4UAZwC+APweu+S/mqafZ2LRpU5RAAbjiiiv4+OOPba4wNBXWdteGXRtsVxwbdm2Iux2WSNyvRLbUunXrxvbt28OKsNk9UErlA/cCl4lIlV1bZnfBnvpWIqExcLOy9BZXPOV7Tg7075/aKFRNulIRkVi+mYNt6gpwa4x25gJzbcrXAr2S6aOhcWRmZjbqnCG1lGwqoWhpEW6Hm6qaKmqpDTtf4atgxMIRuB1ufH4fxcOLKewV/W+Zf0o+fx+pg0H37tI7TKCE3iNeG/3796e0tBTArZRyA1cD14TWUUr1BuYAF4vI7qhGDHFJZCVSWAgHD8Ltt2sBNGGC3uKK5ymfylAtJkyLoVHs3r2bJ554IqpcRDDGEM2D3XaXUzlJd6Tjdrip9ldTS22922F2QuOSbpdQWlrKvop9jHpjFJXOynq31JxOJ7NmzeLHP/7xqcAWYK6IfKKUehhYKyJLgceBLODFgK/y/0Tksib+qtoN1kqkqEivUKqr63QhlmDIytKCpKpKHxC+RRaJFbolVaFajFAxNIpbbrmFQ4cO2Z67+eabm7k3RyeWgt2a7AFqpAa30sJhwnkTeGrNU2H6EksBbwmEKMFUA9ddfx2ZX2TSvUd3DlQeoPJ/lXAacCngBEeag2WlyxjWc1iUYBk2bBjAx6Ex3URkYsj7/Kb5No4e7HQhoYKhshLSIhQbDoeuHylU6tPRNAYjVAyN4oEHHmjpLhz12CnYAY7UaIvd6e9PJzJyUaQCPkow/RNUreJXJb/iD+v/QJpKg8PAP/Q5BsFh32HGvTYO/z/8MbfCDE3L/9/emYdHVZ6N+35myUIiihopKoJ+aEtVati0UhFEEXFBsUpTNyQWKKCAdatabAWt/VlFBD5BBMSlcUVQVMQIYhUV1LAIVOWToFJaxlaWBJKZzLy/P945me3M5Mxkkpkk576uXJmc9Z28M+d5nz28G4eZYIimqgo++0z7T4xzKivhhx+0IAo/L54Asko25anYtCDmzZtn2M9RSjFq1CgOPfRQevToQUVFRYZH1zYId7AXuGP9WDmuHO48686EDvgYwbQVnJc4eeizhzhYd5BqXzXkAhfqfQZV3iq7j32WYOa8j85BAW0S83giQ5KHDYsVRIYAShVbqNikxIwZM+ga9A6WlZWxYcMGvv76ax5++GFuuummzA6uDVFySgk7Ju1g8ZWLyXflR+zz+X2M6TWGHZN2UH5NOTsm7YjRKqIjv8Qh3H3u3bgcUUaMXEzTjV0Ol93HPsPEq+lVWBj5t9sNFRWRyZE1NRAIxJ5rCKBUsIWKTUq4XC7cwVoPy5Yt49prr+WII47g3HPPpbq6OsOja1sUFRQxuNvguGHBRQVF9DmmT9xcFUMwlV9TTvei7pzb6VwO7jsIB4j8MREqXr/X7mOfYcLDiNu3h7w8LRT8/sjjfD792ywkOS8vdluq1Yptn0qypLtNWgvF4XCwa9cuOnTowDvvvMNdd91Vv+9gPMOuTZMS3Q4YYN3OdZZaAxvCZ/d/d3Pm6WfGSVmMZcaQGXanyAwQHulVVaUd6zt2wNy5cP/98NhjUFenBUhentZk7rwTOneO1WrMvq6NqVZsC5VkSHfsXQvm3nvvpXfv3vj9fi655BJOPvlkAFavXs0JJ5yQ4dG1XQzhYDW3JBxPtYfqcdWoukiJUphTiD/gp+TUEp7d+Cwuh4u6QB0zLpjBmF5jmvLt2JhgPIZAC4T8oNXz3nth2jQdRmwIirw8+O1v4dFH4a9/1QKntFQ/upxOLZCiaShhsiGkrbX+7t27t/rkk0+SP9Hj0d6tcLGen6+XB21UY6mrq2P//v106NChflt1dTVKKQqDBt23336b8847z/R8EfnUrJ1wqqQ8t60MT7WHLo90iQg1znfls2PSDkBHfBXmFFLlrYrQYtbtXMegRwex37EfguaQvG/zOOvAWRz6o0NZdugy3DluvH5vQoFiz2vTYfYYSkRhodY6asPqFuTnw6efwtq1cOONugtkOO3agVKxa2ar82r7VKxi1jbNiL1ro7hcrgiBAjqbvjDMQ3j77bc397BaNOnoaWJWINLpcDJl1RS6PNKFsxeezU//96ec/eTZEUUnux7Wlaq/VYFhHtkFNWU1nPbj03hl1SvULK1hv3c/tf5aJi+fbEd9ZQCzSK9E+HzmPpSqKhg6VJvIojlwoHGthS0JFRFxicjFInJr8OciEWlbprPPPosV6Q3F3sWrT92GMNOER40axVFHHQVwstk5iVpLt2bS0dMEzPNXqrxVzPl0jk5y9OtlrpHweP2S69nq2UpRQRHH5h9L/uHa2e/63MVFIy7iihuuoOCKAtgZul4qVYxtGo+V7o0Gubm6KVe04DD8JeEOfrPKSqk66xsUKiJyNLAZ+B1wNLqU9a3A5uC+1o/Ho8MpzJg0yVxoJKpP3YaITr4DGDlypNEhMB6mraVbM1Z6mljVYsq/LqfOb7IEjUOtv5Yej/Vg7qdzOTT30FAkWFV3rh52Ndt/2E5tILLuYypVjG0aT7ggMHJR4mku77yje6eEC47wEvceD3Trpk1hixeHfDMGqTrrrWgq9wOPKaUGKKUmK6UmKaXOBmYDf07+li2QykrzYG7Qgd4PPxyplVjplNOG6d+/P4cffniiQ+K1lm61NNTTxEyLMRMynmoPI5eOxKd8Sd2/TtUxdtlYik4uYvz145l33zy+/fe3XP3Z1Yx4eQS1e2rBqX0zZkmUNumjIQNHSYl2pIMWFErFWubz8iKFjWEwUEoXm5w2LbTm7dULtm/XBShzc+GQQxrXX8WKCesMpdTI6I1KqUdF5Ivkb9kCWb060tMVzQMPwEMPaQ+XEbsXPcuGLtnGnPpdU4tLjNdCelcahpSVJOppYlY48rrF1+FyumKiuyp2VZiWbrHKByd/wDXOa1j4/kLqfl0XWnZWAYN0Xsqm327iyHZHWg5XtrGGxxMKCTYLME1UMDIaEa1lGOvbmprQvrFjQ68Nh3/4Nr8fZs1KPbDVilBJFGdwIMG+1oHHA3fc0fBxPp/WSgD+9Kf4hsxWwuLFixPuHz58uKXj4pBUq2haQTMnI7O9dGkpTocTn9/H9CHTKSooYtrqaRGRXAA+5cNX54uoHHxax9P4ZFfjoqScDifPBJ6h7vSoz29QT/QrPw+8/wBPT3ya9uPaWw5XtklMWRmMGhV6+EcXdywvT1wwMjdXayG5ufrRY2gZ06ZZjxQzqKvTQmv48KbTVA4VkeEm2wVon/wtWxgVFbGpqQ1hFlLRmMDvLOS1114DdAn8NWvWcM455wCwatUqBgwYUC9UUiReC+kYlFKPA4+DDj1tzE3TgafaU598mOwKvuSUEvbV7GPi8onkOHOYvFz78e77+30NnquUonhuMU5xmu7Pd+Xj9XtxiINcVy5VXpMEBSCgAro3iz++Zv7UxqfAB3tr9SLKSndJm/iYaRMG0aVV4gmI2lptCvN6tXO+pERf976GPzqmOBypG1asCJXVwMVx9r2X/C3bIIWF0LNnq8rGX7hwIQAXXXQRW7ZsoVMnvZTdtWsX48eb9lpLBtPW0o29aFOTSsJhOJ5qD5Pfmkytv7b+oT7xTS1gavyRTxyHOAiokJ8ven84LlwsGbGE4k7FALzx1RtMeHNCjGBxO9zMuGAGk96cZHnMEFtO3yY5jDBhM4Hh9cIXX4R8ItHk5YWEkVEdydAyKiu1UDITVg0RCDRhRr1S6vrULt1KKC4OGTgT4Xbr4/x+PSPhx/v9OvT47LNbXTZ+ZWVlvUAB6NixI19++WXCc0pKSnj33XcBckXkO+AewA2glJqDbi09FN0m+gCQ9Z/BVPvDh2PWH8XtdJv6SMIFSjzynDqDccGlCxjcbXD99pq6mhiBkuPIYf3Y9az/1/qEWooZdiRY44gXJux260fJ738fXzD4/aFSLQYOh9Zutm+PzYJwufR142XTg35EpeqkB2shxReLSJewv6eIyAYReVVEjk/tti2IoiJ48slQTJ4xIwYulzaGhm8bPTpU3S0/X5u+Jk9uldFgAwYM4Pzzz+fJJ59k0aJFXHjhhQwcODDhOWVlZezatQvgM6XUsUqp+UqpOUGBQjDqa7xS6n+UUqcG20ZnNQ1Fb1nBzFl/0HcQfyBJ82uQ0b1G883kbyK0pa2erUxcPjHm2Kt7XM2R7Y5k1NJRKBP3VY4jh1xnbD11OxKs8ZgVhLzjDv1o8XpDGogZubmhQpEG1dW6pP2NN8Yef8stOoR41qzYIpI5Obqcy3ffNW69ayWk+D7AAyAiFwFXA6PQJoo5qd+6hRAdzO1yRfpYnE7429/0UqK6Wv+eP18fX16uy7j07Gme1toKsvFnzZrFmDFj2LBhA+vXr2f06NHMnDkz08NqdhJFb1nFLL/Er/zUKes5J+HM/XRuxN9lm8oonltsqok8vfFpVm1fhdMRXBx9TESIzi+O+wWBsLB653AnUwdONS2nb5M8JSX6UVFeDt98o81XVjLnDad6tICoqYkVNgAzZ+oQ4hUrYrUfr7dxHR8NrPhUlFLKiPIaDsxXSn0KfCoi4xp3+ywnuoDknXfGGj+Vio3rM+ogdO0aigGM1m9bUTTY8OHDG+uYb/GER2+5nW58fl9SK3jDfJZsfkkicpw59b4O4/rxTFu+gI+rF18d0lKq0OEPnYBiWKlWRsTk+Y/yU9SuyNZQ0ojxMI/3yAgnN1evbf1+mD1bm8lycxNnPkBI6/nb38yvGc8klgxWNBURkUIRcQCDgHfC9uXFOaflY5bAeN99sTNtNvO1tfD667rO9KBBemlQWhppEmuM0TKL+Oijj+jTpw+FhYXk5OTgdDpp3771BwWaYfQlefGXL7JkxBLOPd56O/bKPZWxjbEaSV2grl5TMjPPReNTvpBWNAi4CegJrAceBcqB/4aOn7h8ol3/K42EF+EwHhnRGgjoh//o0VpL8fm038TrbVigWCEd61wrQuUR9MfqE2CrYd8WkWJacTKaaeW2nBytrRjCITdXayXR+Hw6V6W2Vs/4wYOxJrFW4KQHmDBhAmVlZZx44okcPHiQJ554ghvNjLlthPKvy7n0+Uu58qUrk6rfFa/ffDzaudpxVuezzPe529X7OkBXHy7MKUzq+jmOHO475z7yDsuDQvSTogZ4AVgRPCaoCdk0HrM17Pz5Or5n6tTI9ei0aTAnjuPB5Yrt+GiVadP0Y6+xrl4r0V8LROQt4ChgQ9iufxEWlSMiJyulNjduOFmEWUiGz6eL6YwZo//7r7+uhUc0ZiVdXC6tW/bp0xSjzSjdunXD7/fjdDq5/vrrOfPMMzM9pIxgFgE28pWROHEy8PiBFBUUsdWzlbU719L3mL50L+pef25RQREzLpjB2GVjY67rxImfSGd9nb+Odf9cF3NsjiOHGefP4LhDj2P7nu10eaRLfYjzr0/9NfMr5lt6L7JWWPj8Qmpra6EYGAw4gQAwU/9tR30lR7yMAo8H3njDvAhHVRXcfXfokfPZZ7qcipm/BLT2cvnlcMop8Ic/NBxObNQPu+oqmDIlPcGplvRtpdROImqUgknewNNoZbl1YIRklJbq2fX5Yk1WDzxg/Xpeb6vxoYTTrl07vF4vp512GrfddhudOnVqs+2EzUKCvQEvI14eQY4zh7O7nM3bX79dv29C3wnMvGBmfcJk/+P6c0e/O5j+0XRcDhc+v4+xfcZyYbcL2bh7I7e9fVu9z8OrvOQQa87yKz+/WfYb8px59bkrxnie2fhMfR+VhvDt9zFl9hT+7fo3t5bfGtrhAH6tXwYIUL693HbUWyBefz9ju8sVG/7r9cIPP2ihYzx2zj67YTPXokVao5k0STvmE30dAwFYtUqb3MKTK41M/lQs9Glr0iUiFUqp4rRcrAlJuuFPvOXFunUwcGDiGQtnzhy93Ghl7Nixg6OOOgqfz8f06dPZu3cv48aNo1u3bg2e29qaOXmqPRzz0DFJOdsfPPdBprw7BRQc9B8k35WPz+9DRHCIg1p/LfmufAKBQEyl4GTJc+ahUInzUHxoQ/d/wd3JjbOXk1+f9uv6jo/VvsjPu9H8K9xh39rmNRXCHxtg3t/v00+17yQ66fGQQ7SGIaI1CSNLvmdP/fA3qkE1hFG6JZHD/5BDtOCZODHyuu3ba0t9uGHF6rym0zOY8RIZTUJRkbm47trVvByLGdddp2MEWyFduugUpvz8fO65554MjybziEMgibSSO1feiS8QEkL1Wk7Ytym67lc4Oc4c3A53zMPejBp/DYLgdrhxipMafw35rnxq62oJEDTZLkFrI13A94UP3799lDnLqBhTwdqda7nxzRvZ7w0tqe1s+lisBI263fqhHd0ZorBQ56jce6/WSAyBMHYsPPig9V4qoM+PNqlFU1cHffumNzjV7vzYGOL1WIlm8eJW11Pl1FNPpUePHnF/2iIVuypwO0wCNxIQLlBSYeW1K7mxj/XACIXCF/Ax6YxJbBm3hSUjlpDrCktq9ACXA72BK4EdWnBUeasYeuJQ6gKRCynbrxKJ1aDR6mqtHRyIKslbWxsSKNHceafeF6+plhmJ1r1GEGr37pHJl40NTk2nUEm93nZLwWh0MHeuFhKPPRaqe2A0IZgwQf8OD8EwIsBaSRY9wLJly3jttdcYMmQIQ4YM4dlnn+XZZ59l6NCh/PKXv8z08Jqdsk1lDHtumCWNwcDRyK/fFT+9gsPzD2fG2hkx+/JceeQ6c+OGET/04UMAdMjvEHlM+JCCeZCG4DBycfJdujOkWTZ9sPnaKcGunTHlvUWkv4h8JiJ1ItLqPigVFbEVhJ1ObawwHtqg80vieR7i+Ux8Pu20nz7dvKmW06lDkK0KnIkTQ8748OTLxganWvapiMjT6AKSf1dK/SP1W2aWlO2ziTxqeXmwdKmuE2a0VHvjDV0nIfxYM0NlC6dfv3588MEHDW4zo7XY3j3VHro80iWhmSoaBw7TQpFWOKvzWazbtY5cZy41dTU4xBFx7wJ3AbOHzqbvMX3pObdn3HvkOnOZccEMJi+fHDr/TwSrsAWp06ZNpzgREfbt2xe3ErPf7+ekk07i66+/3oTWddYBJUqpLcYxItIVXd38FuBVpdRLDb3fluJTiS5fH46Rb3L22fDWW/GvYaXMYH5+SACUlmph4vNpv8vw4bpnoJUYIuM6VjUSq9/XZJZKC9H5tTNF5P9E5GURiS0iZBERqRSRTSKyXkSM3JfDReRtEfkq+LtDcHvcnuUicl3w+K9E5LpUxxODoZVs3aprGhg6bbRAAf1J6NAhJFAqKvSMtfKeKgDV1dW8//779X+vWbOmzUV/WUksjCZAICWBAvD3b/9OTV0Ne2v3UuuvjRFmARVg6IlDWf+v9QnrhtX6a5m8fDLTz59er33k3ZvH1LensuW7Laz9ei279+3mQNUB9u/fz759+wAd/tznmD4xfpS1a9caARpepZQXeA7dxbMepVSlUmoj0HBFzBZEovL1oLfX1CQWKGDNZ2JUeDI6QPp8+hE0ebLWYGbEKq5cbFJnvqkqRVl21CulVorIaqAPMBAYC5wMmLwFywxUSn0f9vcdwDtKqQeCqvMdwO1E9iw/Hd2z/HQRORxd4bY32rX5qYi8qpT6oRFjCmkloAVJQ/UPDGFRVqb1XCOI3OnUs52XZx6S3AqYP38+o0aNYm8wdOSwww5jwYIFGR5V85Js4mJTke/S9hAj6dFK2Re3003PTj3ZMWkHlXsq+VHej3j56ZeZ/cJsevToQfEo6wGdO3fupHPn8DY4fIf+vrZ6zMxe6eCee7TWEf74MR43Hk9sB8jx42PbPxUWaqt8eXlksEBTrXEtCxUReQcoAD4E/g70UUrtTvN4hgEDgq8XAe+ihUp9z3LgIxExepYPAN5WSv03OMa3gSFA6h7xcE+bQTyBUlCgtRGjYfSoUZFZSX6/Xg68+GLINNbK6NWrFxs2bGDfvn0opTj00EMzPaRmJ7pro5U8EDOc4tTRXE43tXW1CJKUNhNQASrGVNC9qDvrdq6zVPYl3F9SVFDEiBEjcLvdnHXWWbzxxhts3ryZGWZLXxPimNJTigptSR09E5m9GkNBAfz735G51C4XXHstfP+9ToyMjioz6yfo9+vHT0Npd+kimZDijUAv4BRgL7BHRD5USiXZrLIeBawQEQXMDXbw62gkVSqldonIUcFj4/Usj7c9dRJ1zAE904GA1j2ffTakd77zjvmnyuGAXbv0rLZSXn/9dTZv3kxN2PufMmVKBkfUvHiqPXTr0I1PR39qGnZrldlDZzO8o9WjqQAAHshJREFU+3Aq91SyunJ1ZNJhkHxXPqU9S5n36byYfJPwjo6f7fqswTGYOdq3bNnCpk2bACgtLaVv376Wx3/sscfy7bffRmwiTsfOhsi2jp7xiGf2cjqTbxgbTXV1bDmWujodJzR3rr6vFXPZnXdq4VFSohMam7pPYDLmr8kAIlKILs+yEPgRENtkwRr9lFL/DAqOt0UkkfM/Xs9yS73Mk1r1xOuYA9qMtXixLhTZq1ek3vnii+bnHDigHfa//W2racwVztixYzlw4ACrVq3ihhtu4KWXXkrqQdTSie72OH3I9JiwWyvc1PcmxvQOJcee/e7ZMcfc0e8Obv75zYB21l+75NoIwWJoHUYHyWjcDjf57ny8dV7u6n8XY3qNifGLuMNq2bkaSnKIok+fPnz11VcAOSKSA/yK+vz71onZGrSgABYsgJEjk+8Pnwzz5+vclSlTtPbh9cb2B8zLi8y5jpd2l04sWwFFZIKIPI8uLnkpsADt60gJpdQ/g793A68AfYF/B81aBH8b5rV4Pcst9TJXSj2ulOqtlOpd1NB/NLxjjhGykZ+vfxYsgMGDQ3qnVYyQ4lGjWk1IscGaNWt46qmn6NChA/fccw8ffvhh9Gq11bLVs5Xrl17PwbqD7K3dy8G6gxGO7wK3xdhOoGNhR174/AVWbFtBxa6KGMd/YU4hw7sPp/zrcro80oUbXruBukAdboc7JrzXLHCgMKeQZSXLKL+mnG8mf8Pd/e82TVjcsGED7du3p3379hxyyCFs3Lix/nVD1addLhezZs0COAnYCryglNosIveKyCUAItIn2O3zCmCuiLToeoFma9C6Om22MnJKzCoNp4uiosg+LEY/QSPfZMGC5re6J7MUyQceBj5VKsWuQUFEpABwKKX2B18PBu5FN/66Dngg+Htp8BTTnuXBQpf3G1Fiwev8vjFjAyL1RK8Xtm3TaafdgwUAE2kz4YQ3kAb9+g9/0GVHjZlu4X3r84PB8u3ateOf//wnRxxxBNu3b8/wqJqesk1lXL/0+hgTVCAQqHd8P7zmYR5YY60+3F0r76p/7Xa4kahUa3/AT2FOYX3Byvrtys/vev+Om39+c72QMAsc8Af8FHcqbjDz3d9Im83QoUMBPg8PPVVKTQl7vQ69+GvxGF9do7Gr06lNVrW1cNNN+hiR+Pko4YQfZ/Uc0I+lcO2juUxcibCsqSilHlRKfQwcLiLHGT8p3rcj8L6IbADWAq8rpZajhcl5IvIVcF7wb9A9y79G9yyfB4wLjum/wFR0PPw64F7Dad9oioq0MDnvPG2+Ki7WhkxjX3QK6tmx5gpTH4uROFlWFtlAoYVm3F900UXs2bOHW2+9lZ49e9K1a1d+9atfZXpYTUqihle1gVpW71gNYJqUOOq0UQ0mPfoCPlRAxSQZVnmrTJ3v0z+azvcHvmfdznV4qj2WkhTjceDAAXxhwSZffPEF06dP55VXXmnw3LaE8dUdNEgnEV52mRYm0cLAqnAIP87qORMmhNa54RQV6VS4TK1Rk0l+vBitqRyNNkt1QfdXObnphpd+LCdSeTyxVeAgsjCksVQpLNRCJzpKzOWKXychL08vSaKrzCWTjZRl1NbWUlNTYzkCrKUmP67buY5BTw2K6wjPdeby6q9e5cqXrmRvbahKX2FOIS9f8TKXPn9pg4mSBe4CFl+5mA75HeqjszzVHjpP7xzroA/2js9z5eH1e5k/bD4lp5TETVJMRP/+/Zk/fz4nnngi27Zto2/fvlx11VVs2bKFPn368ICFrLqWOq9W2brV/OueCrm5OpYnGd/LZZfp0i9mAqUpaYqCktOAM4BypVSxiAwEWpfXOZzKSvNqbBMn6rRVQ+csKtJJkjk5sZ+yRIV3nM7YbUY2UgsQKosXL064vzW3F24ossopTvbU7DE1QQGWQn2NY8MFQryeK4aQMX6XLi3l3OPPrQ8TToYffviBE088EYBFixZRUlLCzJkz8Xq99OrVy5JQaY0Y60ejn0ljBUpeHtx1l+590qtXcucuXx4ymmQjyQgVn1LqPyLiEBGHUmqViPylyUaWaeL5TXJyYh/8yVQsNvD7Y0uU1tam3ratmXnttdcA2L17N2vWrOGcc84BYNWqVQwYMCChUFm+fDkTJ04EXSPqDqVUxJNKREYCDxLq4TNLKfVE2t9ECmz1bGXi8sSFJA7UHeDqxVeHKv+iqwnPHzaf7Xu2Nxjq68BBgABXvnRlhOYBMKbXGPbX7ueulXeR48zB5/fFlGlpTOXgcF/OypUrufVWHdack5ODoymy+1oAiSo0pUJBgQ4iHTxY/11aCjq+QXP66bBxo15jVlfHhiZn+9ozmU/JnmA48XvAsyIyA2iUwz6rKSoyr3dQVxebhhruYyks1LPeUMjHiBHawxceHuJw6GVLC/CtLFy4kIULFyIibNmyhZdffpmXX36ZzZsTB/P4/X7Gjx/Pm2++CbAZKBGRn5oc+rxS6rTgT1YIlLJNZRTPLU7cjySIT/nwq9DTwIGD0zqeZhrqG45b3OS4dHiyEVFWurS0vhd82aYypqyaQq4zF5/fx7RzpsXeuxGVg3v06MEtt9zC9OnT2bZtG4ODT749e/akdL2WTngudDoECuiwXyNtzePRj45wPv5Y+1VuvRU2bYotHJnt1Z6SESrDgAPAZGA58H+ASUWZVsTw4bq5QW5uqApxdBqqUSPs3HO1P2TlStiwIVYLiebFF3XIyL33hjxzRuu1FlTNuLKykk6dOtX/3bFjR7788su4xxv1oU444QTQOUUx9aGyAU+1hxXbVrBi2wq2erayYtuKuM55K7icLp77/Lm4pq88p64q/Lt+v6v3kRgYmkd4u+L93v3U+muZsmpKRO2uZJzyZsybN48jjzySyspKVqxYQbt27QCdFHnLLbekdM1sx/gKm33l5s5Nf65JaWnoEWLkuURTUwP33w9HHpnesvTNQTLJj0aVwACwSESc6OSmZ5tiYBknvNOOCNx2m3bQh89mvB6hoF9ff31846vRSOHuu2P9Mdmu34YxYMAAzj//fEpKShARnnvuOQYOHBj3+CTqQ10uIv2BL4HJSinT5JemKOdRtqmM6165LqJuVo4jB28g9fpeVd4qHvrwobil8Wv8NRTmFPLIh49EmM0gpHmYtSuOrt2VjFPejPz8fO64Q1es93q9fP7554BObDzzzDNTvm62kugr7PFoh3i6mT9fJywWFSXOTnA4dE2xbAgTToYGNRURaS8ivxeRWSIyOFgxeAI6xPfKph9iBojutGMsGxIdE61hlJToT4S7gaZNLlds6PHBg9mt34Yxa9Ysxo4dy4YNG1i/fj2jR49m5syZcY+3WB/qNaCrUqoHUI6uAxfvetYTWy3gqfYwaumomEKMqQgUt8NNYU7IR9ZQr5UqbxU1/hrTcOKigiLT/JPw2l1mlYNTZfXq1Zx44omMHz+ecePGcdJJJ/Hee++l5drZgtlX+PrrdXQXaC3FLCvgjDMad9/w6sCG5dzMWl5dDZdeqgVfpsOEk8GK+etp4MfAJuAGYAU6G3aYUirrzBZpwUwnja4TbeWY7t11A+hEmHnifD5YtCi+Tp5lnHnmmZxzzjkMGjSIfv36JTzWSn0opdR/lFKG6jYPXXOuWajcU4nTYRKZlyQF7gKWlSxj1gWzOCTnkJh9iUrl5+fks2TEEsqvKWfHpB31TvrG5J8ky80338yKFStYvXo17733Hm+99RaTrXY6bSGYfYVra0MpaWZaisMBH33UuPtG+0RKSnQ2/NSpscKlhVnDAWtC5QSl1Eil1Fx0CHFv4CKl1PqmHVoGMdNJoz8JVo4BbTKbM0f7ZeK1ZAuYtJa49VadWZXlSZEvvPACffv25aWXXuKFF17g9NNP56WX4vddMupDBbPuBW1CfTX8GKNUT5BL0CU/moWuh3Wl1mfdb+LAwfOXP0+eK/JpUO2rZvue7aYteAMqwKNDHiXflR+hyRj4/D6KOxWbah4lp5SwY9KOGIGTbnw+Hz/+8Y/r/z7ppJMikiJbA/FMT7W1MG6cuZZi9lVtCKMjYyKfSFGRtoQvXRr7mIheqybyAWUDVoRK/SdJKeUHtiul0hQHkaWYZcxHfxKsHGMwZgx8+y3Mnh0bypGIFtCG+L777mPdunUsWrSIp556irVr1zJ16tS4xxv1oc4//3zQ/Xhi6kMBN4nI5mDFhZuAkU39PsJxOBv+WrRztSPPmcf/XvS/HN/heKYOiH3PRqSXmXYxpvcYdkzawcprVzLnwjlJaR/pNnWZ0bt3b0pLS3n33Xd59913+c1vfkOvZBMqshzjK2yWMpaK8IhHQYEWFlZa9RYXx947fK3aEopwNJhRLyJ+wDAGC7oG2IHga6WUSlxlLstIKjvXSl2uZGp3zZ0LY8cmPsaMLG5DfOqpp9aXSgdd++pnP/tZxLZ4ZFvmtafawxtfvcGENyck7IniwMFtv7iNrod2ZfJbk3Vb4LoaRISautDytn1ue8qvKafPMX0azG5PJfu9KamtrWX27Nm8//77KKXo378/48aNIze34aLk2TavifB44Oijk08zS4ZkC2UYwQPhfU9KSsyLfDRnEY60ZdQrpRpvYG6pWKkTbbWWtMeja4ilQhYHpg8ZMqQ++gvg+eefN4oKtiiMEvYqoKgJJO62FCDAA+8/gNvhxhfwxS25Ep4v0lB2eyrZ701Jbm4uN998MzfffHOmh9JkeDzaypxugeJ26/ibVJthxYv2Miuzn42Bosk1TGhrpLOCcEVFZFdIq+TmZnVg+oMPPsjLL7/MBx98gFKK0aNHc9lll2V6WEkRnv+RDL5A5HzmOfNQKHJdOjGxqZzoTcmVV17JCy+8wKmnnhpTKRlg48aNGRhV+klnt0YzIdLYEGCztapVN26msYVKPBIFsDcXubm6o2ROjhZwWSpYLr/8ci6//PJMDyNlKvdUWqrH1RAiwmejP6PKW5U1ZqxkMVoHL1u2LMMjaRo8Hli1SjfQstK9oiEMV6qZEEn319XwATVHS+DGYAsVM8ID2A1ds7RUf3JSncHi4tAnwSrdu8PAgdqTGAjAo49qp38W9WBZvHgxt99+O7t370YphVIKEWHfvn0ZHZdVPNUeFm9ZnFL7X7e4cTlduJ3ues2ke1Ezl45NM0Z1hC5dutRv+/777zniiCNMNZeWxNy5MH5849v8GlxxhY69aSohYkZLSIRsmxXiGsJKDkqyFBXp3BMLjs561q/XQqimRi+rxo7VZU07d86acOPbbruNV199lb1797Jv3z7279/fYgRK2aYyjpt+nOVGWtHMHDqzWcJ7m5OPPvqoviBoRUUFp5xyCqeccgodO3Zk+fLlmR5eyhgxMukSKACZUuayPRHS1lTMaCrjZUkJHHGErilWnTi7Oi5GyXmjrMuoUY3ToBpJx44d6d7cjR3SgOFHqfGnZlQvzCmkZ6eeWedgbywTJkzg/vvvZ+/evZxzzjm8+eabnHHGGfzjH/+gpKSEIUOGZHqISePxhDoxppNsdJJnA7ZQMaMpjZdmgeiNoaZGL8Puvjv+MU1gLjP6qfTu3ZsRI0Zw6aWXRoSbZns/lbmfzE3aMR+OP+BPuRJwNlNXV1dfmXjKlCmcEaxJ8pOf/CSTw2oUlZU6Ez7dZKOTPBuwhUo8msp4GS6wnE6oip8PYZmpU7VZzExjaKKAA6OfCuj+9CtWrKj/W0SyWqhs9WzlD+/+wfLxboeba352DWWbyiL8J61JQzEI75mSH5Wo2xJ9Kh4P/PBD6k75ggK9Bpw+XV/r/vuz20meDVhuJ9xayJrWpOGt5CZNanxso9MJDzwAPXpAeO+LkSObJFvqu+++49hjjzXd99prr3HxxQ13RchEkpxZBeJEuB1uNozdQPei7lmXoNgUOJ1OCgoKUEpx8ODB+tL3SilqamoslWrJVPJjuEIO8NBD8Mgjush4Kl+vvDydCV9cHPq6ZFGMTLPTFO2EbdKJEYjep4/2sRgV7FIVLn6/zuRqiHBDsMej82cg8ptjgUGDBvHWW2/RNUr/X7hwIdOmTbMkVJobw49iVaAAzLxgZn1EV2vzn5jhT6cnuxkxFHLQ7sZULMyXXAJvvx2piRjdGQ2s5jq3Zezor2zAqCb3zTc6ZLgpOXAA3n9fNx87+mg4/3z9c+yxWrBZrFQ3ffp0zjvvPL766qv6bX/+8595+OGHWb16dVO+g5Qx+pE0RK4jlxxHDnMumsOY3k08HzaNZutWuO66UAaAFYHidIZqvObm6pqvS5dqJd5KjS6b+NiaSjZRVKT9I089lf52cwZ1dWBWesMIWS4o0MdMm6ZNaWCq/w/t04fcOXO44IILWLJkCU888QTr1q3jvffeo0OHDk0z9kbS9bCu1PgSa4IOHLxa8irFnYpbvVbSGvjrX/X6KFkF6ze/0U1XzRIWbU2kcdhCJdsId+Q7HKmHHqeKcb9oU9pVV2l/zJNPavtAXR2D+vfnyZ/+lAF9+3LmSSfxzrRp5NXVhcxqe/bAYYclbVprKooKiuh7bF/+/s3f4x7jdDhtgdIC8HhgxAidHZ8KN91kC5CmwhYq2YgReVZRoVu/NZXWkgzPhnWNrqvjEEDeeQcF1ALvbNjAURdfjEKXr45If3Q6depxU5v2GmDuJ3MTChSAfHc+lXsqbaGSxZSVwa9/nfr5EyaYB0rapAfbp5KtFBVpL2F4z5a8PG0emzOn4TbFTcx+tODYD3jRvRH2hW2PwO/XprW5c5t1jOF4qj1MeHNCg8eFVxa2yT48Ht3yNxWGDYMtWxpuxmrTOGyhku2UlIS8h998ox36Y8bAzp3w/PO6PGpLYeLEjDUbq9hVEdOBEbQPpTna87YVgqVcThGRbSJyR/R+EckVkeeD+z8Wka7JXL+yUocIW0VE+1x274YlS2wNpTloQU+kNoyZ8beoCK68UmsB110XWagyJ0cXn/R4dJiy15veLP5UcTozVtfim73fmG6/rd9t3Pzzm1t9/klz4Pf7GT9+PMCX6Lbj60TkVaXUlrDDSoEflFLdRORXwF+AEVbv0bUrWEmtc7ngsce0dmL7TZoXW6i0dML9L2aO8TFj9IPc64VnnoF589JbVS8ZlMpIXYvSpaUsWL/AdN+1P7u2TeSfNAdr166lW7dufP31116llFdEngOGAeFCZRjwx+Drl4BZIiLKYhZ2UREsXAhXX22+TsrN1fEtmehUYaOxhUprwPC/xNtnCJh+/XQcZbgA6twZvv1Wm9Z274bt23U1ZREtBPr2hTVr4i8PDzssMoM/Hjk5Galr8dcP/hpXoOQ6cxO2DbZJjp07d9K5c+fwTd8Bp0cddgzwLYBSqk5E9gJHAN9bvY+xjlq1CrZtg6OOguOO0x/lqqq2me2eTdhCpa1hJoCiDc333x8ZwB8dIlxQoL/Nffvqc7du1T6fujr4/nsdllxQAMcfrwMMMhRW7Kn2cHv57XH3K5TtlE8jcZSN6I1mHpGYE0VkNDAa4Ljjjos5wbD+2mQfrUKoiMgQYAbgBJ5QSiXfIGPrVnjlFb1aP3BA+yj69WubRtloH46ZIOrXL/S6e/ekPaDhDl1M5kxEcoGngF7Af4ARSqnKZO7x8IcPEyC+L+m+c+6zzV5p5Nhjj+Xbb7+N2AT8M+qw74DOwHci4gIOBf4bfS2l1OPA46BrfzXJgG2ahBYvVETECcwGzkN/YM2cg4m58UaYNSt2+5NPap/EM8/YBto00hwOXU+1h4fWPBR3f2lxKbeceUtK47cxp0+fPkbZnhwRyQF+BURnlLwKXAd8CPwSWGnVn2LTMmgNIcV9gW1Kqa+VUl7AcA5aY+tWc4FiEAjowPgMhcK2RgyHLuBNMGfDgEXB1y8BgySJ2uuVeypxO81zeeZdNI8nLnki6XHbJMblcjFLf5dOArYCLyilNovIvSJySfCw+cARQQ31ZiAm7NimZdPiNRXCHH9BYpyDCe2za9dau4vd4i1tNIdDt+thXQmoWNOXW9wM+4n1NYdNcgwdOhTg8/AS6UqpKWGva4ArMjA0m2aiNWgqDTr+lFKPK6V6K6V6F0ULhr59rd3FbvGWNtLt0BWRT0TkE0+YNllUUMSCSxfglpC24na4WTR8ke1HsbFpQlqDpmI4/gzMnIPx6d5dFwOKZwJzOHRgvK2lpI3mcuiWnFLCucefS8Uu3TPGLhRpY9P0tAahsg44UUSOB3Zi7hxMzMyZMG6cHf3VTDSnQ7eooIjB3eLk8NjY2KSdFi9Ugvb2CcBb6JDiBUqpzUlfKIWwWJvUMBy6F154oeHQXWA4dIFPlFKvoh26Twcduv9FCx4bG5ssp831qBcRD7DDZNeRJJHVm2Fa0lgh/ni7KKXSpgbGmduW9L9qSWMFe16ToSWNt1Hz2uaESjxE5JPwiJVspiWNFTI73pb0v2pJYwV7XpOhJY23sWNtDdFfNjY2NjZZgi1UbGxsbGzShi1UQjye6QEkQUsaK2R2vC3pf9WSxgr2vCZDSxpvo8Zq+1RsbGxsbNKGranY2NjY2KSNNidURGSIiHzRVD2004mFsY4UEY+IrA/+3JCJcQbHskBEdovI53H2i4g8GnwvG0WkZ5rvb89rE2DPq3XseQ2ilGozP+jkyP8DTgBygA3AT6OOGQfMCb7+FfB8Fo91JDAr0//X4Fj6Az3RxQTN9g8F3kTX9DoD+NieV3te7XltffPa1jQVK2XyG1VyPY00rqR/M6OUeg+T2lxhDAOeUpqPgMNEpFOabm/PaxNhz6tl7HkN0taEilmZ/GPiHaOUqgOMkuvNjZWxAlweVE9fEpHOJvuzBavvp6mubc9r02DPa9Q4grTZeW1rQsVKOXVLJdebASvjeA3oqpTqAZQTWrFlI035f7XnNXPY86qx5zVIWxMqVsrk1x+TqOR6M9DgWJVS/1FK1Qb/nIfu556tNK5FQeOvbc9r02DPa9Q4grTZeW1rQqW+TH5YyfVXo44xSq5DZntoNzjWKBvnJeiKv9nKq8C1waiSM4C9Sqldabq2Pa+Zw55XjT2vBpmOQshA1MNQ4Et0pMZdwW33ApcEX+cBLwLbgLXACVk81j8Dm9GRJquAn2RwrGXALsCHXuWUAmOBscH9AswOvpdNQG97Xu15tee19c2rnVFvY2NjY5M22pr5y8bGxsamCbGFio2NjY1N2rCFio2NjY1N2rCFio2NjY1N2rCFio2NjY1N2rCFio2NjY1N2rCFigkiokTkobC/bxGRPzZwzlgRuTbN45ghIjtFJG3zJCJ/FJFbgq/vFZFzg68niUi7FK5Xla6xNTX2vCZ1PXtekx+HPa/YQiUetcBwETnS6glKqTlKqafSNYDgB/MydFG3/um6bjhKqSlKqfLgn5OApD+kLQx7Xlsn9rxmEbZQMacO3ad5cvQOEekiIu8EK42+IyLHBbeHryhuEpEtwWOeC24rEN0YZ52IVIhIQ2WxBwKfA48BJWH3/6OILBKRFSJSKSLDReT/icgmEVkuIu7gcZUi8hcRWRv86WbyXp4UkV+KyE3A0cAqEVkV3FcVdtwvReTJ4OvjReTD4PuYGnW9W4PbN4rInxp4f5nAnld7Xu15DR3XJPNqC5X4zAauEpFDo7bPQvcZ6AE8Czxqcu4dQHHwmLHBbXeh6xL1QX8AHxSRggT3L0GXUngFuMj48AX5H+BCdM+DZ4BVSqlTgYPB7Qb7lFJ9g2N+JN6NlFKPoovFDVRKDUwwJoAZwGPB9/EvY6OIDAZORPeVOA3oJSJNsmJrJPa8mmPPqz2vaZlXW6jEQSm1D3gKuClq18+BvwVfPw38wuT0jcCzInI1ehUFMBi4Q0TWA++iaxYdZ3Zv0QXphgJLguP4OHi+wZtKKR+6Jo8TWB7cvgnoGnZcWdjvn8d5q8nSL+y6T4dtHxz8qQA+A36C/tBmFfa8xsWeV3te0zKvrnRcpBXzCPofvjDBMWbF0y5E21UvAf4gIiejC7RdrpT6wsJ9h6BLeG8S3cSuHXAAeD24vxZAKRUQEZ8KFXALEDmnKs5rK4Qfn5dgn4EAf1ZKzU3yPpnAnleNPa8ae17TiK2pJEAp9V/gBXQFT4M16LLWAFcB74efI9ph11kptQq4DTgMKATeAm6U4KdORIoT3LoEuEEp1VUp1RU4HhgsyUd7jAj7/WEDx+4HDgn7+98i0l1CDkiDD4h8/wZvAaNEpBBARI4RkaOSHG+zYM+rPa8G9rymf15todIwDwHhUSU3AdeLyEbgGmBi1PFO4BkR2YRWLacrpfYAUwE3sFFEPg/+HUPwg3g+oVUOSqlq9Jfh4iTHnisiHwfHGOPEjOJx4E3D8Ye2My8DVqJLZBtMBMaLyDr06swY4wq0meHD4Ht/icgPfbZhz6s9r2DPa9rn1S5930oRkUp0D4TvMz0Wm/Rhz2vrpDXNq62p2NjY2NikDVtTySAicj7wl6jN25VSl5kdb9MysOe1dWLPqzVsoWJjY2NjkzZs85eNjY2NTdqwhYqNjY2NTdqwhYqNjY2NTdqwhYqNjY2NTdqwhYqNjY2NTdr4/z7XKNamzMyBAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "fig, axes = plt.subplots(nrows=2, ncols=3)\n", "\n", "data.plot.scatter(x=\"Noise_Amplitude\", y=\"Raw_Mean\", color='r', ax=axes[0,0])\n", "data.plot.scatter(x=\"Noise_Amplitude\", y=\"Raw_SD\", color='r', ax=axes[1,0])\n", "data.plot.scatter(x=\"Noise_Amplitude\", y=\"Khodadad_Mean\", color='g', ax=axes[0,1])\n", "data.plot.scatter(x=\"Noise_Amplitude\", y=\"Khodadad_SD\", color='g', ax=axes[1,1]) \n", "data.plot.scatter(x=\"Noise_Amplitude\", y=\"BioSPPy_Mean\", color='b', ax=axes[0,2])\n", "data.plot.scatter(x=\"Noise_Amplitude\", y=\"BioSPPy_SD\", color='b', ax=axes[1,2]) " ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.7.4" } }, "nbformat": 4, "nbformat_minor": 2 }