{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Creating and accessing a fit\n", "This example shows you, how you can easily calculate and visualize a fit to your data" ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": true }, "outputs": [], "source": [ "import numpy as np\n", "import xarray as xr\n", "import psyplot.project as psy\n", "%matplotlib inline\n", "%config InlineBackend.close_figures = False" ] }, { "cell_type": "markdown", "metadata": { "collapsed": true }, "source": [ "First we start with some example data to make a linear regression from the equation\n", "``y(x) = 4 * x + 30``" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "\n", "Dimensions: (experiment: 50)\n", "Dimensions without coordinates: experiment\n", "Data variables:\n", " x (experiment) float64 0.0 2.041 4.082 6.122 8.163 10.2 12.24 ...\n", " y (experiment) float64 56.74 74.26 2.13 57.84 69.51 34.99 89.89 ..." ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" } ], "source": [ "x = np.linspace(0, 100)\n", "y = x * 4 + 30 + 50* np.random.normal(size=x.size)\n", "ds = xr.Dataset({'x': xr.Variable(('experiment', ), x),\n", " 'y': xr.Variable(('experiment', ), y)})\n", "ds" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We can show this input data using the ``lineplot`` plot method from the psy-simple plugin:" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX0AAAECCAYAAAASDQdFAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAF0hJREFUeJzt3X9wXOV56PGvJBtiY5k0tuLiGOOQts+M3fArnTpOqQm9\ngYSqNaQ/ZpoJ6YRAcIKnt5dMyA9zMQkhJJSU3rSkNHH5lUw7mUKblNYlkHtDiwoe1RQY6ti83CSQ\n4gn1LyAYgy822vvHrmAtr6TV7p7ds3u+nxkG6Wh39epI+/g9z/u8z+krlUpIkoqhv9MDkCS1j0Ff\nkgrEoC9JBWLQl6QCMehLUoEY9CWpQGbV86CIeAh4rvLpE8DXgC8DB4HvppSuiog+4M+Bk4EDwEUp\npR+1fsiSpEZNG/Qj4miglFL6tapjDwPvTSk9GRGbIuIU4M3A0Smld0TESuB64LysBi5Jmrl6Zvon\nA8dExN3AAPBZ4KiU0pOVr98NvAs4DvgOQEppNCJ+qfXDlSQ1o56c/ovAdSmldwMfBW6pHBu3DzgW\nGAR+WnX8UES4ZiBJOVJPUH4c+CuAlNL/pRzY31D19UHgWeD5ysevvnZKaaxF45QktUA96Z0PAW8F\n1kXEYmAusD8i3gw8Cbwb+AxwPPAbwB0R8XbgP6Z74d2799n4R1JX2nDTKDt27z/i+JKheVx14S9n\n+r2Hhgb7Gn1uPUH/JuCWiBgBxoALKv//a8pXCveklLZExIPAWRFxf+V5FzQ6KEnKu5/sebHm8af3\nHvkPQZ5MG/RTSgeB82t8adWEx5Uo5/wlqectXji35kz/uAXHADC6bSebNj/JT/a8yOKFcxletYyV\nyxe1eZRHcqFVkhowvGrZJMdPYHTbTr565/fZsXs/Y6USO3bv56t3fp/RbTvbO8gaDPqS1ICVyxex\nds0KlgzNY6C/jyVD81i7ZgUrly9i0+Ynaz5n0+Yft3WMtdS1I1eSdKSVyxfVTNnkOd/vTF+SWmzx\nwrk1j4/n+zvJoC9JLTZVvr/TTO9IUouNp3w2bf4xT+/dz3ELjmF41Qm5qN7p6+SN0VuxOSuvZVET\nvfzyy7z//b/D7bffOelj7rzzWwwPr2FgYKCNI5PUbZrZnNXV6Z08l0VNVP7Hderf0ze+cQtjY3au\nkJSdrk7vTFUW1ehs/667/pFNm+6kVCpx4YVreeKJH3Hfffdy4MABjj329VxzzXVcfPEHuf76G5g3\nbx7Dw/+NG27YyM///C/woQ+dz9e+diuzZpVP60svvcRVV/1P9u3bx5vetOTV7/HIIw9xyy0bKZVK\nvPTSi1x55ed55JGH2Lt3L1deuZ6rr76W6667hl27drF37x5OP301F130kYZ+Hkmq1tVBP6uyqMHB\n+XzhC1+iVCrx6KOP8OUv3wjAxz72Bzz22DZWr34no6MPMDT0RhYvfhNbtowye/Zsli494dWAD/Dt\nb/8tJ574c3z4wx9l27atPPTQvwPwxBM/YsOGz7FgwUK+8Y1buPfe/80HPnABt912M1dd9QV27drJ\nihVv5ZOfPJeXX36Z3/qtXzfoS2qJrg76022DbtTSpeUV9r6+PgYGZnHlleuZM2cOe/bs4tChQ6xe\nfSZf//rN/OzPHsfFF1/C7bd/k7GxV3jnO3/tsNd56qkf8453/CoAy5f/IrNmlXP1Q0ND/MmfXMfc\nuXPZvXsXJ510SuUZJUqlEvPnz2f79u/z8MMPMmfOMRw8eLCpn0eSxnV1Tj+rsqj+/vJp+eEPf8DI\nyD/z2c9ew6WXXsbY2BilUokTT3wLTz/9E7Zv/z6rVp3OSy+9yP3338fb3/4rh73OsmUnsnXrowA8\n/vhjHDr0CgBf/OLVXH75Z1i//koWLhxifDG9v7+fsbFX+Kd/+gcGB+dzxRWf4/d+7/0cOHCgqZ9H\nksZ19Uw/67KoJUuWMGfOXC655CJKpRILFgyxZ89uAE455TT+67+ernz8Np588gle97rXHfb88877\nba6++krWrfswS5eewNFHHwXAe94zzCWXXMicOXN5wxve8OprnnTSKVx22f/gYx/7JJ/5zOVs3foo\ns2fP5vjjT2DPnj0sXLiwJT+XpOLq+pJNSSqawpZsSpJmxqAvSQXS1Tl9SepGnewkYNCXpDYa7yQw\nbryTANCWwG/Ql9QS3dIHq9Oy6CQwEwZ9SU3r9Oy1m3T6BisGfUlN6/TstV55uBrJqpNAvazekdS0\nTs9e65GXrrydvsGKM31JTWv37LWRGXterkY6fYMVg76kpg2vWnZYTv+1462fvTa6fpCnq5HJbqje\nDqZ3JDVt5fJFrF2zgiVD8xjo72PJ0DzWrlmRSWCbasY+lTzfrLydnOlLaol2zV4bnbG382okzwz6\nkrpKo+sHnc6l54VBX1JXaWbG3slcel4Y9CV1FWfszbGfviR1mWb66TvTl6Qp5GEX78Sx7Ni9/9A/\n/PG5DcVvg74kTSJPPYUmjGWg0dcx6EvSJDqxi3eyK4vJxjJTBn1JmkS7d/FOdWUx2VhmyqAvSdSe\nYbe7p9BUVxaTjWWmbMMgqfAm68AZS3+m5uOz2sU71ZXFZN05Z8qgL6nwJpthp/98rm09hWDq/kDV\n/Y2AQ41+j7rSOxHxRuBB4F3AK8CtwBiwNaW0rvKYDcAwcBC4NKW0pdFBSVI7TTXDbucu3ul2G4+P\nZWhocHaj32PamX5EzAL+Ahg/K9cD61NKZwD9EXFuRJwKrE4prQTeB3yl0QFJUrvlpQNnO7qV1jPT\n/xJwI/BpoA84LaU0UvnaXcDZQALuAUgpPRURAxGxIKW0t2UjlaSM5KkDZ9ZXFlPO9CPig8CulNJ3\nKQf8ic/ZBxwLDAI/rTr+QuW4JOVeO+8H0GnTzfQvAMYi4izgZODrwFDV1weBZ4HngfkTjj/XwnFK\nUqaK0oGz7oZrEfE94CPAdcAfp5Tui4gbge8BPwSupZzqOR74+5TSqdO9pg3XJGnm2t1w7ePAxoiY\nDWwH7kgplSJiBNhMOQ20rtEBSZKyY2tlST0lT10xs2JrZUkiX10x88oduZJ6xlS9a1Rm0JfUM9rd\nFbMbGfQl9Yy87KzNM4O+pJ4xWSfKTuyszSsXciX1jPHF2k2bf8zTe/dz3IJjGF51gou4VSzZlKQu\n00zJpukdSSoQg74kFYhBX5IKxIVcqUc12o6gCG0MisyFXKkHTWxHMG66HvGNPk/t5UKupMM02o7A\nNga9z6Av9aBG2xHYxqD3GfSlHtRoOwLbGPQ+g77UgxptR2Abg95n9Y7UgxptR2Abg95n9Y4kdRmr\ndyRJdTHoS1KBGPQlqUBcyJWUS7aDyIZBX1LmZhrAJ7aD2LF7/6ufG/ibY3pHUqbGA/iO3fsZK5Ve\nDeCj23ZO+hzbQWTHoC8pU40EcNtBZMegLylTjQRw20Fkx6AvKVONBHDbQWTHoC8pU40E8JXLF7F2\nzQqWDM1joL+PJUPz7OnfIrZhkJS5cvWO/XxapZk2DAZ9Seoy9t6RJNXFoC9JBWLQl6QCsQ2DlBP2\nmlE7GPSlHLDXjNrF9I6UA/aaUbtMO9OPiH5gIxDAGPAR4P8Bt1Y+35pSWld57AZgGDgIXJpS2pLN\nsKXeYq8ZtUs9M/3fBEoppdOBK4BrgOuB9SmlM4D+iDg3Ik4FVqeUVgLvA76S1aClXmOvGbXLtEE/\npfT3wMWVT08AngVOSymNVI7dBZwFnA7cU3nOU8BARCxo+YilHmSvGbVLXQu5KaWxiLgVOA/4XcpB\nftw+4FhgENhbdfyFyvHqY5JqGF+stVWBslZ39U5K6YMR8UZgCzCn6kuDlGf/zwPzJxx/rhWDlIpg\n5fJFbQvylocW17TpnYg4PyI+Vfn0APAK8GBEnFE5dg4wAjwAnB0RfRGxFOhLKT2TxaAlNa6RO1mp\nd9Qz0/874JaI+JfK4/878BjwlxExG9gO3JFSKkXECLAZ6APWZTRmSU2YqjzU2X7vs8umVDAXXXsv\nYzXe9wP9fWz8xJkdGJFmyi6bkupmeWixGfSlgrE8tNjsvSMVjOWhxWZOX5K6jDl9SVJdDPqSVCAG\nfUkqEBdylSm3+7eG51Gt4kKuMjPxblDj1q5ZYcCaAc+jJnIhV7nk3aBaw/OoVjK9o8wU9W5QrU7F\nFPU8KhvO9JWZIm73z6KDZRHPo7Jj0FdmirjdP4tUTJ7O4+i2nWy4aZSLrr2XDTeN2o65C5neUWaK\nuN0/i1RMXs7jxAXl8auY6jEq/wz6ylQ77waVB4sXzmXH7iMDfLOpmDycR/vw9wbTO1IL5SkV02ou\nKPcGZ/pSC+UlFZOFrK5i1F4GfanF8pCKycLwqmU1N4n1wlVMkRj0VXi2OKhPL1/FFIltGFRojbY4\n8B8KdZJtGKQGNVJXn8UGLKldDPoqtEYqUuyFo25mTr9HmX6oTyMVKZYuqps50+9Bph/q10hdvb1w\n1M0M+j3I9EP9Vi5fxNo1K1gyNI+B/j6WDM2bdhG3lzdgqfeZ3ulBph9mZqZ19ZYuqpsZ9HuQOyez\n16sbsNT7TO/0INMPkibjTL8HmX6QNBl35EpSl2lmR64zfXUV9x9IzTHoq2t45yapeS7kqmu4/0Bq\nnkFfXcP9B1LzTO8odybL27v/QGqeM33lylR9g9x/IDXPoK9cmSpv30ifHEmHmzK9ExGzgJuBZcBR\nwOeBbcCtwBiwNaW0rvLYDcAwcBC4NKW0JbNRq2dNl7e3/cGRLGPVTEyX0z8f2JNS+v2I+Bngkcp/\n61NKIxFxY0ScC/wnsDqltDIijgf+FvjlTEeunpRV3r5XA6NlrJqp6dI7fwNcUfXYQ8BpKaWRyrG7\ngLOA04F7AFJKTwEDEbGg9cNVr8sib9/L9xewjFUzNeVMP6X0IkBEDAK3A5cDX6p6yD7gWGAQ2Ft1\n/IXK8epjakCvzlAnk0XfoOnWCbqZZayaqWlLNivpmr8DbkgpfTMi/qjqy4PAs8DzwPwJx59r5UCL\nqKiX7q3O2/dyYLSMVTM1ZXonIhYBdwOfSCndVjn8cESsrnx8DjACPACcHRF9EbEU6EspPZPVoIvC\nS/fW6OXbG1rGqpmabqb/aeD1wBWV6pwS8IfAn0XEbGA7cEdKqRQRI8BmoA9Yl+GYC6OXZ6jtNLxq\n2WFXTK8d7/7AaBttzZStlXNsw02jNS/dlwzN46oLW18c1cvrB+WfzcCYN738N5elZlorG/RzbGJO\nf1wWG5Ka+V6NvnF9wxdbO/++e4399HtUOy/dG61waXSxuaiL1HpNL1dV5ZlBP+fatQO10fWDRt+4\nvuHlmlVn2HtHQOMVLo2+cX3Dq5erqvLMoC+g8dK/Rt+4vuFluWlnGPQF0HAHy0bfuL7hZdfUzrB6\nR01rtBzSMkqpMZZsSl3CMlW1giWbUhewTFV5YE5fahN7KSkPDPpSm1imqjww6EttYpmq8sCgL7WJ\nZarKAxdyVRerTppnG2TlgSWbmpbdEKV8sWSzC3TzTLnIzdG6+fcm1WLQb6HJAkS312cXteqk239v\nUi0u5LbIeIDYsXs/Y6XSqwFi/B+CWrqlPruoVSfd/nuTajHot8hUAaLbZ8pFrTrp9t+bVIvpnRaZ\nKkAsXji35r1uu2WmXNSqk27/vUm1GPRbZKoAMbzqhJrVL900U27XHbzyZHjVsq7/vUkTGfRbZKoA\nUdSZcru1utLG35t6kXX6LWR/+M5xL4GKxDr9nChiCiQviryXQJoJq3fUE6y0kepj0FdPKOpeAmmm\nDPrqCUXdSyDNlDl99QQrbaT6WL0jSV2mmeod0zuSVCAGfUkqEIO+JBWIQV+SCsSgL0kFYslmF/NW\nfpJmyqDfpbyVn6RGmN7pUt7KT1Ij6prpR8RK4IsppTMj4i3ArcAYsDWltK7ymA3AMHAQuDSltCWb\nIQtsMCapMdPO9CPiMmAjcHTl0PXA+pTSGUB/RJwbEacCq1NKK4H3AV/JasAqs8GYpEbUk975AfDe\nqs/fllIaqXx8F3AWcDpwD0BK6SlgICIWtHKgvWx020423DTKRdfey4abRhndtnPa59hgTFIjpk3v\npJS+FRHVkaS658M+4FhgENhbdfyFyvHqY6qh0QVZG4xJakQj1TtjVR8PAs8CzwPzJxx/rolxFUYz\nd3zyTl2SZqqR6p2HImJ15eNzgBHgAeDsiOiLiKVAX0rpmVYNspe5ICupnRqZ6X8c2BgRs4HtwB0p\npVJEjACbKad/1rVwjB3Rro1PixfOZcfuIwO8C7KSsmA//Rom5tnHrV2zgpXLF7X0H4TpvpckTdRM\nP3135NYw3canVu6EdUFWUjsZ9GuYKs/ezMLrZFyQldQutmGoYaqNTy68SupmBv0aptr45E5YSd3M\n9E4N0+XZay28uhNWUjfIXfVON/SIL4/RhVdJndFM9U6ugr7li5I0vZ4p2cyiMqYbrhwkqV1yFfRb\nXRnj3aUk6XC5qt5pdWWMd5eSpMPlKui3uke8NfWSdLhcpXda3ZLAZmaSdLiOVu+c+/E7S61YXJ1s\nsdZqIEm9qGurd8ZKpaYXV+tZrLWmXpLKcpPeabQsc7oyT5uZSdJrcrOQ2+jiqou1klS/3AT9RhdX\nbYAmSfXLTdBvtCyz1WWektTLOprTH+jva3px1cXambEthVRsuWq4pmxZwir1hmZKNnOT3lH2bEsh\nyaBfIFY6ScpNnX4nFC2/bVsKSYWd6Y/nt3fs3n/YzuDRbTs7PbTMWOkkqbAz/Sxu2JJ3VjpJKmzQ\nL2p+27YUUrEVNr3jTl5JRVTYoG9+W1IRFTa9Y35bUhG5I1eSuow7ciVJdTHoS1KBGPQlqUAM+pJU\nIAZ9SSoQg74kFUhL6/Qjog/4c+Bk4ABwUUrpR638HpKkxrV6c9Z5wNEppXdExErg+sqxphWtDbIk\nZaHV6Z3Tge8ApJRGgV9qxYsWsQ2yJGWh1UF/PvDTqs8PRUTT38Pb/ElSa7Q66D8PDFa/fkpprNkX\nLWobZElqtVbn9O8HfgO4IyLeDvzHVA+ut3/EWKn0KPDWicdfGSs9OjQ0eHIjA5WkImppw7Wq6p2T\nKocuSCk93rJvIElqSke7bEqS2svNWZJUIAZ9SSoQg74kFYhBX5IKpO33yLU/D0TELOBmYBlwFPB5\nYBtwKzAGbE0prevU+NotIt4IPAi8C3iFgp4HgIj4FLAGmE35fXIfBTsflffHbZTfH4eAD1PAv4tK\nK5svppTOjIi3UOPnj4gNwDBwELg0pbRlutftxEz/1f48wKcp9+cpmvOBPSml1cA5wA2Uz8P6lNIZ\nQH9EnNvJAbZL5Q3+F8D4DrxCngeAiDgDWFV5b7wTWEoxz8evAwMppV8BPgdcQ8HOQ0RcBmwEjq4c\nOuLnj4hTgdUppZXA+4Cv1PPanQj6mfTn6TJ/A1xR+bif8mzmtJTSSOXYXZRnvUXwJeBG4CdAH8U9\nDwDvBrZGxLeBO4F/pJjn43FgViUrcCzlWWzRzsMPgPdWff62CT//WZRj6T0AKaWngIGIWDDdC3ci\n6GfSn6ebpJReTCntj4hB4HbgcsoBb9w+yn/sPS0iPgjsSil9l9d+/uq/hUKchyoLgbcBvwN8FPgr\nink+XgDeDDwGfBX4Uwr2/kgpfYvyZHBcrZ9/kMNj6QvUcV46EWwz6c/TbSLieOB7wG0ppW9SztWN\nGwSe68jA2usC4KyIuJfyGs/XgaGqrxflPIzbC9ydUjpU2cl+gMPfxEU5H5cC30kpBa/9XRxV9fWi\nnIdqE+PDs5Rj6fwJx6c9L50I+vdTztlRT3+eXhQRi4C7gU+klG6rHH44IlZXPj4HGKn55B6SUjoj\npXRmSulM4BHgA8BdRTsPVf4VeA9ARCwGjgH+TyXXD8U5H8/w2gz2OcoFJw8X8DxUe6jG++IB4OyI\n6IuIpUBfSumZ6V6o7dU7wLcoz+7ur3x+QQfG0GmfBl4PXFFZfS8Bfwj8WUTMBrYDd3RwfJ30cWBj\nEc9DSmlTRPxqRPwb5cv5jwJPAn9ZsPPxv4CbI+I+ylVMnwL+neKdh2pHvC9SSqWIGAE2U/57qaui\nyd47klQghVpAlaSiM+hLUoEY9CWpQAz6klQgBn1JKhCDviQViEFfkgrEoC9JBfL/AUfHpopE0cMy\nAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "raw = psy.plot.lineplot(\n", " ds, name='y', coord='x', linewidth=0, marker='o', legendlabels='raw data', \n", " legend='upper left')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The visualization of the fit is straight forward using the ``linreg`` plot method:" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX0AAAECCAYAAAASDQdFAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xd8W9X9+P+XlmVJtpxhx44ThzBPCJTZNoRC2BQIm7Y/\n+im07BVWIIOET6HfflpGEnZomGF1AmGFUWgpLSkjhAJtGTnMkOHtJLa1pXvv7w/ZxrElW5bt2Jbe\nz8ejD+wr6d4jNX77+Nz3eb9tlmUhhBAiP9iHegBCCCG2Hwn6QgiRRyToCyFEHpGgL4QQeUSCvhBC\n5BEJ+kIIkUecmTxJKTUOeBc4EvABK4FP2x5eprV+Qil1PXAcEAdma63XDMJ4hRBC9EOvQV8p5QTu\nAUJth/YDbtFa39bpOfsCB2utpymlqoAVwHcHYbxCCCH6IZPlnSXAMqC67fv9gZlKqX8ope5XShUB\nBwGvAGitNwAOpdTYwRiwEEKI7PUY9JVSZwH1Wuu/ALa2/60G5mqtDwG+BK4HioHmTi8NACWDMWAh\nhBDZ622mfzZwlFLqNWAf4BHgJa31+22PPwPsC7QA/k6vKwa2DvBYhRBC9JMt09o7Sqm/ARcBjwKX\naa3XKKUuBSYCfwIWA0cBVcCzWut9eztnQ0OrFP4RIgdc9+BqNjYEux2fWFbEL8/N39t7lmURqq3B\n3dqK25GcY5/95AbMFJHPYYPlP6hKdyIK/voK3jtuwd7UxGlQssKyWrIZU0bZO11cBNytlIoCtcAF\nWuuAUup14C2SS0CzshmMEGJkqm4MpTxe09T9F0G+iAQCmHU1FFsWNsc3iyqVfhcbm+Pdnl/pd6U8\nj33DBnxLbsT1zmqsAjehCy6G+5ZlPa6Mg77W+vBO334vxeO/BH6Z9UiEECNWZak35Ux//FjfEIwm\nvdUf1/HCW+uobgxRWepl5vTJTJtaPqDXME2TcPvs3ukAm22bx0+Y4mfZ6qZurzt+in/bA7EYhY89\njOfRh7DFYsQOOJDQnPkYlRO2T9AXQoh0Zk6fzL3PfZTi+A5DMJrUVn9ct80YNzYEO74fqMAfaW3F\nrKvBb7OB05HyOQdM8gLw/NoWqlviVPpdHD/F33EcwLnmHXxLbsKx/mvM0jICV15N/PAjk79A+lkO\nX4K+EKLf2oPmC299TU1TkPFjfcycvsOAz6LbZTNjf+GtdWmOf93vcZqmSaimmsJAIDm778UBk7zb\nBPl2tqZGvHfehvuVP2PZ7UR+9GNCF1wEvqJ+ja8zCfpCiAExbWr5oAX5zrKdsQ/WfYf2tfsSSDu7\n75Vh4H72KTzLlmIPBEjsvgfB+Qsw1O79GlsqEvSFECNKtjP2gb7vYFkWwbo63M1b8WUb7AGHXotv\n0Q04P/4I0+cjOGc+0ZNPA0f25+yJBH0hxIiS7Yx9IO87xCIR4tWb8JsGtmwDfjCA9757cD/5J2ym\nSfSo7xO64iqssaXZnS9DEvSFECNKtjP2gbrvEKivp2DLZvwpMnMyYlm4XnsV321LsDc2YFRNIjjn\nGhLfndb3c2VBgn4XhmFw5ZWXsHHjBi68cBbHHDOTFSse57TTfjTUQxNC0L8Ze3/uO8QiEWLVmyg2\nEtiznN3bqzfhXXIzBW+9geVyET73AsJnngVud1bny4YE/S4aGhoIh8M8/fSLHcceffRBCfpCDBPb\nO1MIINDQQMGWzZQ47GDPog1JPE7h7x/Ds/wBbLEo8e9MIzhnPuak7Z/SOqyDvu8X/4t75TMDes7o\nCScT/MWv0j5+yy03snHjehYvvoFdd1W0tDTT0tLCrbfezFVXzR/QsQghsrO9MoXi0SjR6k0UJ+LY\nHdn1nHK+/x6+RTfgWPcV5pixBK+4jthR389uaWgADOugPxSuvvoarr9+IaWlZdhsNn7603NYseJx\nCfhC5JlAQwOuzU18Uh1lZaeNVCd02UiVjm3LFrxLb8f94vNYNhuRU39I+KJZWMXFWY/p7fUhVq5t\nITb76c1kGb+HddAP/uJXPc7Ktx+pCydEvojHYkQ3baQ4Eeed6ug2JRM2Nsc7vk8b+E0T98pn8Nx9\nF/bWFhK7KYLzFmLssWe/xvX2+tA3Y7HZss7nHNZBf7jo565nIcQIEdzchLOhgRKnA+x2Vq5NXcjy\n+bUtKYO+4/PP8N58A64P/4Pl9RG88mqip/0InP0PtenG0lcS9FOwdVlr23HHnfi//7uOn/9c6skJ\nkYsS8TiR6k0URSM4OmXmVLd0r4aZ8ngohOfB+yj80++xGQaxw44geOUcrHHj+jyW9iWcrstJ6cbS\nVxL0u6ioGM899yzf5tgdd2Rf0U4IMbyFtm7BXl+fzMzpsgs2kzLIrtf/jvfWRTjq6jAqJxCaM5/4\n9G6FiDOyzRIO2y4npRtLX0nQF0LkJcMwCFdvxBcO43Q4Us6weyqDbK+pwXvrIgr++TqW00n4Z+cQ\nPuscKPRkPaaelpPSjaWvJOgLIfJOqHkrtvo6SuzJ2X26GfbF08Zy8bSx25RBPmEXL4e8/gSeB+/D\nFokQ33d/gvMWYE7esd/j6mk5aVqVhzBj+MvnYTY1hRPZXkOCvhAibyTicSI1m/BFIjg7LeX0NMP+\n1dEVHTdtnf9+H+/Pb8T55ReYo0YRnLuA2LEzByznPt0STkWxk5YiPwcfuiszDrNRVlacus1WBjIK\n+kqpccC7wJGAATwMmMCHWutZbc+5DpgJxIHZWus12Q5KCCEGWnBzE47GxpRr973dsLU1b8Vz950U\nrnwWgMhJpxC++DKskpIBHWO6JZxjD96ZooqKAblGr0FfKeUE7gHaS9vdCizUWq9SSi1TSp0ErAdm\naK2nKaWqgBVA/nZDFkIMG9tk5qQpV5z2hm2xk4IXnsN71+3Ym5tJ7LwLoXkLSey196CMtf0vipVr\nm6luSTB+rJfjD9xxQHcfZzLTXwIsAxaQbHq+n9Z6VdtjLwFHAxp4BUBrvUEp5VBKjdVa9/+ugxBC\nZKmn2X1nqWbYVU0b+L+XHqDok39jeTyELruSyI9OB2fWKyu9Mk2TqTv42HP/HfH4/b2/IAs9Bn2l\n1FlAvdb6L0qphW2HOxegaAVKgGKg8ycWaDs+6EF/ezQ6HgixWIyf/OQHPPHEc2mf89xzTzNz5olp\nZyNCiMx0zO5jURwZ1Mzp3Le2sSnAOe+v4Og3V2A3DGKHHEZo9hzM8oFZXknFsiwCpok1thTvmLHd\n9goNpN5m+mcDplLqKGBv4FGgrNPjxcAWoAXwdzm+dQDHmdL2aHQ8UCzLIvmHUnqPPfYQxx57vAR9\nIfphm9l9HypiHjDJy8Eb/oX3nkU4aqoxKsbTetU84gfPGMTRQjCRIFEyGu+4cdizqeDZRz0Gfa31\nIe1fK6X+BlwELFZKzdBavw4cC/wN+AK4WSm1BKgCbFrrzYM37KTBaHT80kvP88ILz2FZFueeeyFf\nffUlr7/+GpFIhJKSUdxww2IuuOAsbr11KUVFRcyceQRLl97PrrvuxjnnnMF99z2Ms23LdTgc5pe/\n/F9aW1uZMGFixzU++OA9HnrofizLIhwOcf31v+aDD96jqamJ669fyK9+dTOLF99AfX09TU2NHHTQ\nDM4776Ks3o8Q+SKTtft0bPV1+G5bQsHf/4blcBA+82eEzz4fPNnn3PcmkjCIFRXhHldOoWvwloy6\nyiZlcw5wv1LKBXwCPKm1tpRSq4C3SE5nZw3gGNMarEbHxcV+brxxCZZl8Z//fNCxI/eqqy5j7dqP\nmTHjUFavfpOysnFUVk5gzZrVuFwuJk3aoSPgAzzzzAp22mkXzj//Yj7++EPee+9fAHz11Zdcd93/\nMXZsKY899hCvvfZXzjzzbB55ZDm//OWN1NfXscce32L+/JOIxWKceupxEvSF6EFydt9AicPRt96y\niQTuJ/6E94F7sIVCxPfeh9DcBRg77zJoY40bCcKFXlwTyvEVFg7addLJOOhrrQ/v9O2hKR7/JbBd\ni9MMdKPjdpPaGhvYbDYcDifXX78Qj8dDY2M9iUSCGTMO49FHl1NRMZ4LLriEJ574I6ZpcOihh29z\nng0bvubAAw8GYOrUPXG21fQoKyvjttsW4/V6aWioZ6+99ml7hYVlWfj9fj755CPef/9dPB4f8fjA\n1NwQItck4nEimza2rd33bXbv+Oi/+G6+Aednn2L6SwguvI7YzBOya5KSAcMwCDpd2MdPxNeP8sr9\nNaI3Zw1ko+PO2tfVvvjic1at+jv33fcw0WiEc889E8uy2GmnnampqWbLls1cdNGlPProct5443Vu\nu+0325xn8uSd+PDD/3DQQTP49NO1JBIGADfd9CueeOI5PB4Pv/71L9rW+5PXNU2DF19cSXGxn7lz\nF7Jx4wZWrny6X+9HiFwU3LIZR0N9t9l9uoJl7WwtLXjuWYr7maewWRbRmScSuvRyrFGjB2WclmUR\nAKyycfhGjxmUa/TFiA76g902beLEiXg8Xi655Dwsy2Ls2DIaGxsA2Gef/aitrWn7en/WrfuKwi5/\nqp188mn86lfXM2vW+UyatANudwEAxxwzk0suORePx8uYMWM6zrnXXvswd+6VXHXVfH7xi2v58MP/\n4HK5qKragcbGRkpLSwfkfQkxkhmGQXjTBoqi3Wf3PRUsO6DKQ8ErL+G94zbsWzaT2GlnQnMXkNhn\n30Eba8AwMEaNwVdWNqgZOX1hs4awWHxDQ6tUqhdCZKy9Zk5RmiWYa1+pTbnJav94Hde++QCuf63B\ncrsJn3M+kR+fAYN0AzWcMIgXF1NYXjEo2XhlZcVZ/wYZ0TN9IUR+6FoRM52u5RQK4lF++M4KTnv3\nKVxGgtj3DiZ09TzM8ZWDMs5owiDq8VBQNR5fQcGgXKO/JOgLIYa1rhUxe9K5nMK+697nolfvpbK5\nli3+UpwL5hM/5LBBaUhuGAZBVwGOCePxFRUN+PkHkgR9IcSwZCQShGs29Tq77+yEKX7+9OpnnPuP\n5czQ/8Sw2Xl6/xPxXj6L7+w28PfELMsiYLNBeQW+klEDfv7BIEFfCDHshLZuwdZQn9HsvoNhcMjq\nlRz12N24wkH0+N1YcdKl7HPYvnwnXRPzfggYBsbosfhKS4fNTdpMyI1cIcSwsU3NnD7kyzvWfoJv\n0Q04P/kYs7iY8MWXET3plEHJuQ8nDOJ+P4XjyoesZIrcyBVCjHjBpkYcTY3JvPsMg7Ut0Irn3mW4\nn3oCm2kSPeY4QpddiTVm7ICPb6h30g4UCfpCiCGV1a5ay6Lg1b/gveMW7I2NGJN2IDhvAYn9vzPg\n4zNNk4DDiaOyatjfpM2EBH0hxJDJpmaOfcMGfLfchGv121gFBYTOv4jIGT+DAU6RtCyLgGVhjS3F\nNwh/OQwVCfpCiG1sjx4V8ViMaEe9+wxn97EYhb99BM8jy7HFYsSmTSc0Zz7mxKoBHRtA0DBI+Efh\nKy8fUTdpMyFBX4gclU3w3h49KgINDbg2N1HizHx273z3HXyLb8Kx/mvM0lICV84hfviRA55z317u\nuLC8gkJnbobH3HxXQuS5bIP3YPSoaBePxYhu2khxIo7dmVmwt21uwnvX7bj//CKW3U7kR6cTuuBi\n8A3s2nqu3KTNhAR9IXJQtsF7sHpUBDc34exLNyvTxP3s03iW3YW9tZXE7lMJzluIMWX3fo2j+2Vy\n6yZtJiToC5GDsg3eA92jIhGPE6nZ1FYRM7M0TMenGt/iG3F++F9Mn4/g1fOJnnJa35qj9CJXb9Jm\nQoK+EDko2+A9kD0qAo2NODf3Ie8+GMT7wL24H/9DMuf+yKMJXXEVVmlZ76/tg1y+SZuJXoO+UsoO\n3A8owCTZJ9cNrAQ+bXvaMq31E0qp64HjgDgwW2u9ZlBGLYToUbbBeyB6VERDIRJ1tcm1+0xm55aF\n6x+v4bt1MfaGeowJEwnOvYbEtOkZXzMT4Xiio9xxrt6kzUQm7/wEwNJaH6SUOgS4gWTAv0VrfVv7\nk5RS+wIHa62nKaWqgBXAdwdj0EKInvUneE+bWp7VTVvTNAnX1VLQ0oLfmdns3l69Ce8tiyh4859Y\nLhfhcy8gfOZZ4Hb3+frpfHOTtirnb9JmIqPaO0opu9baVEr9jGR/3DDJmb+T5Gx/NnA24NFaL2p7\nzb+Ao7XWTanPKrV3hMgV4ZYWrPpaiiCzJZN4nMI//BbP8vuxRaPEv/1dgnOvwZzUv1annSVv0jqw\nl47D4/cP2HmHg0GvvdMW8B8GTgZ+AEwA7tdav6+UWgBcD2wBOgf4AFDS5ZgQIoe0lz/2hkO4HJkt\nmTjffw/v4htxfvUl5ugxBBf8nNjRxwxozn3ANDHHjMU7Zmxertv3JOOFLa31WUqpccA7wHStdU3b\nQ88Ad7X9t/Ov02Jg60ANVAgxvCRLKLSlYWYQ8G1bt+BdeifuF57DstmInPoDwhfOwhrAWXh7BUxP\neQX2QaiwmQsyuZF7BjBRa30TECF5M/cppdTlbTdqjwDeBd4AFiullgBVgE1rvXnwhi6EGArxaDRZ\nQiEey6yEgmlS8PxzeO++E3tLM4ldFcH5CzD2+NaAjSmaMIh6vRRUVQzbNoXDRSYz/aeAh5RS/2h7\n/uXARuBupVQUqAUu0FoHlFKvA28BNmDWII1ZCDFE+lpCwfHF53gX34jr3x9geb0EL7+K6A//Pxig\n7JmONoUTK/H5sttLkG+kiYoQolfxaDRZQsFIZLZsEg7jWX4/hX/4LTbDIHbYEQSvvBprXOZZQW+v\nD7FybQvVLXEq/S5OmOLngLYOWJZlEQAoLcM7anR2b2oEkyYqQohBE9zchLOhITm7zyDgu1a9jvfW\nm3HU1mKMryR09TxWVe3Pyg9aqG7Z0C2Ap/L2+hDLVn+TA7KxOd7x/bcmuEmMGoOvrExu0mZBgr4Q\nIiUjkSBcvZGiSARHBgXS7LU1eG9bQsHrf8dyOgn/7BzCZ53D2/VW2gCeLvCvXNuS8viznwbYb8ae\nFA5Rm8JcIEFfCNFNuLkZ6msza0yeiFP4pz/geeBebJEI8X32IzhvAeaOOwGwcm1typc9v7YlbdCv\nbomnPF7XnOHNY5GWBH0hRAfTNAlVb8IbCmSUd+/89wd4F92A88svMEeNIjjnGmLHHb9Nzn26AJ7u\nOECl38XG5u6PZ1v4TXxDgr4QAoBIIIBZW02JzdZr3r2teSue39xF4XPPJF974smEL7kMq2RUt+em\nC+CVflfKc1uWxZFTinl4dfeM72wKv4ltSdAXIs9ZlkWotgZ3Swvu3tbuLYuCF5/Hu/R27Fu3kthp\nZ0LzFpLYe5+0Lzlhin+bNf12x0/pvikrYJoYo8Zw8M674i6v71fhN5GapGwKkccigQBGXQ3FltVr\nJox93Vf4Ft2I6/1/YRUWEj73AiKn/w84U8/YO3t7fYjnO6VfHt8leyeYSGCUjMIzrlx20magPymb\nEvSFyEOmaRKurcHd2tr77D4SxvPwcgp/9yi2RILYjEMIXTkXc/z4fo8jnDCIFxfjLhuH09X7Lw+R\nJHn6QoiMhZubsepr8dvt0EvAd735T7y3LMJRvQmjvJzQVfOIzzi032OQsglDR4K+EHmivSKmLxzG\n2Uvao62+Ht/tSyh47VUsh4PwT35K+JzzwZt+Q1VGY2grm+CsmoCvn+cS2ZGgL0Qe2LYiZg8BP5HA\n/eTjeO9fhi0UIv6tvQnNW4Cxy679un5H2YTyCnwpMnzE9iNBX4gcFo/FkhUxY9FeNzU5PvoQ3803\n4PxMY/pLknXujz8xs/62PQglDOIl+duTdriRoC9Ejsq0IqattRXPPUtxP70Cm2URPe54QpdeiTW6\nf4XM2tsUFlSNp0jW7YcNCfpCDBOrP67jhbfWUd0YorLUy8zpk7PKS4+Fw8RrqpMVMXu6UWtZFLzy\nZ7x33oZ9cxPG5B0Jzl1AYr/9s38TfNOm0FFZha+oqF/nEgNPgr4Qw8Dqj+u497mPOr7f2BDs+D7T\nwG9ZFqG6Wgqam3ttTG5f/zW+xTfhevcdrAI3oYtmEfmfM6GfaZMBw8AYU4pvrLQpHK4k6AsxDLzw\n1ro0x7/OKOhvs8mqp9l9NIrnsYcpfPQhbPE4senfI3T1PMwJE7MbePv1Ewax4mIKyyukINowJ0Ff\niGGgujGU8nhNU7DH15mmSbimGncggM/p6LG5uPOdt/EtvgnHxg2YZeMIzJ5D/NDD+9WQPGEYhAo9\nuCaU4ysszPo8YvvJpEeuHbgfUCT7414ERIGH277/UGs9q+251wEzgTgwu62HrhCiF5WlXjY2dA/w\nPVWVDDVvxdZQj99m63GTla2xAe+dt+H+y8tYdjuR039C6LwLoR/tBdvX7e0VlfgGsLG5GHyZzPRP\nACyt9UFKqUOAG0j2wF2otV6llFqmlDoJWA/M0FpPU0pVASuA7w7ayIXIITOnT95mTf+b492rSibi\ncSLVm/BFIz1vsjIM3E+vwHPPUuzBIIk99iQ4dwGGmpL1OC3LImBZWGNL8Y0Zm/V5xNDpNehrrZ9V\nSq1s+3YHYAtwpNZ6Vduxl4CjAQ280vaaDUoph1JqrNa6e3k9IcQ22tfte6oqaVkWwcZGXFs297rJ\nyrH2E3yLbsD5yceYxcUE5y4gevKp/cq5DyUMEiUlUhRthMtoTV9rbSqlHgZOBn4IHNXp4VagBCgG\nOgf4QNtxCfpCZGDa1PK0N22joRCJ2mqKDQO7o4eAGwzgvW8Z7icfx2aaRI85jtBlV2J1mZX31HS8\n27UTBlGfD/ekCtxSFG3Ey/hGrtb6LKXUOGAN4On0UDHJ2X8L4O9yfOtADFKIfNVR6761Fa/Dnn6m\nblm4/vZXfLcvwd7YiDFpB4JzryHx7e4rrD01He8c+BOGQajAjXO81MnJJZncyD0DmKi1vgmIAAbw\nrlLqEK31P4Bjgb8BXwA3K6WWAFWATWvdvfWNECIjkUAAsz0Ns4fZvX3jBry3LKLg7TexCgoInX8R\nkTN+Bml2waZrOt7es9ayLFptNmxSJycnZTLTfwp4SCn1j7bnXw6sBR5QSrmAT4AntdaWUmoV8BbJ\nG72zBmnMQuS0brXu06VUxmIU/u5RPA8vxxaLEp92AMGrr8Gsqurx/D31rA2YJsbosbK5KodJExUh\nhpFIa2tydt9LwHX+aw2+RTfiWP815tixhK64mtiRR2eUc3/tK7Wpe9aOdvPL86fLTdoRQJqoCDHC\nGYZBpLaGwmCAgp6Ko21uwnvX7bj//CKWzUbkB/8f4Qsvxioqzvha6XrWnnDwLhLw84AEfSGGWKh5\nK7b6umQnq3QB3zRxP/cMnt/cib21lYSaQnD+tRi7T+3z9Q6Y5MU0TVZ+GqC2NUGlNB3PK7K8I8QQ\nScTjRGo24Yv0vMnK8dmnyZz7D/+L5fURuugSoqf+sOdmKGl0NDMpLcM7qn+lk8XQkeUdIUaYZCer\nBkocPdS6D4XwPHAPhY//EZthED3yaEKXX4VVVpbdNQ2DxKgx+MrK5CZtHpOgL8R2FI9Gk52s4rH0\n1SgtC9c/XsN72xIc9XUYEyYSmjOf+AEHZnXNcDxB3O+nsLyCQqmAmfck6AuxHViWRbChgYKtW3os\noWCvqU7m3L+xCsvpJHz2eYR/ejZkUcGyvXOVa0KVVMAUHSToCzHI2kso+E0z/SarRJzCP/wOz4P3\nYYtGie//bYJzrsGcvGOfr5esgOmUzlUiJQn6YlANVAvAkcg0TcJ1tRS0tCQ7WaVZR3d+8D7exTfi\n/PILzNFjCC74X2JHH9vnOvdSAVNkQoK+GDQD0QJwpAq3tGDV1/ZY6962dQveu+/C/fyzyZz7U04j\nfNGlWCnq0/dWIE0qYIpMSdAXg6a/LQBHIiORIFyzCV84nD4N0zQpePF5vEtvx97cTGLX3QjOW4ix\n57dSPr2nAmn7VrqJer0UVFXgS1NrR4jOJOiLQZNtC8CRKrR1C7b6uh7TMB1ffI538Y24/v0BltdL\n6PLZRH54OjjT/yimK5D27NoWvvWdfWXdXvSJBH0xaLJpATgStXeyKopFcTgcqZdiymx4Hrqfwt//\nFpthEDv0cIKz52CN6/0vnnQF0upaExRKwBd9JEFfDJq+tAAcqQKNjTg3NyZn93Z7yqWYd//wZw59\n40E8DbUYFeMJXT2P+EEzMr5Gpd+VskBarv3yFNuHBH0xaDJpAThSxSIR4rU1FMdj2Dst5XReiilt\nbeD81x7kwM/fxrA7CJ95FuFzzoNCT6pTpnXMbkU8sGZLt+ND8cszn7OxcoXU3hGiDyzLItSWhlmY\nIuf+7Cc3YEskOOH95/mft/6IJx7hwwlTuffIi/h/F32vT9fq2Fw1rpz3v2we8l+eXbOx2l144h4S\n+Lczqb0jxHYQbmnBaqjrsZPV95o/54fPLGXHxnW0FBZz7+Hn8+rUw6kalXlmTXJzlQN7xQR8bemb\n06YWDnlgzcdsrFzUY9BXSjmB5cBkoAD4NbARWAl82va0ZVrrJ5RS1wPHAXFgttZ6zWANWojtKRGP\nE6mtxhsO4XI4U26asjU341m2lHnPPgXAK3seycMH/5RWTzJoHz+le+59KgHTxBymm6vyLRsrV/U2\n0z8DaNRa/1QpNQZ4H/h/wC1a69van6SU2hc4WGs9TSlVBawAundkFmKECTY14mhqu1HrSPHjYlkU\nvPRCMud+yxYSO+3Mmp9dxTNMJtQSp8rv4vguG6lSCScM4n4/nvKKYbu5Kl+ysXJdb0H/ceCJtq9t\nJGfx+wNTlFInk5ztzwYOAl4B0FpvUEo5lFJjtdbd2/MIMQLEIhHiNdUUJeLb3KjtzL7uK3yLb8T1\n3r+wCgsJzbqCyOk/Rjld/CrD67Sv2xdMrMDndg/cGxgE+ZCNlQ96DPpa6xCAUqqYZPD/X8ANPKC1\nfl8ptQC4HtgCdA7wAaCkyzEhhj3LsgjW1VHQvDVZLyfVrDsSwfPIcgp/+wi2RILYQTMIXTUPc/z4\njK+Tat1+uMvlbKx80uuN3LblmqeApVrrPyqlSrTWzW0PPwPc1fbfzv9yi4GtAz1YIQbDN2mIQSr8\nLk6cUsz0SamXLFxvv4l3yU04Nm1is7+Uew49jw37HcQJ8RIOyOBaI70o2rSp5RLkR7jebuSWAy8D\ns7TWr7VvpTeuAAAgAElEQVQdflkpdanW+l3gCOBd4A1gsVJqCVAF2LTWmwdx3EIMiK5piNXNce5Z\nvRkbtm3W4W0NDXjvuAX3q3/BtDt4ev+T+P3004kUeKAl0bEhq6e1+1A8QXzUaLzjxg3bdXuR+3qb\n6S8ARgE/V0pdB1gk1/DvUEpFgVrgAq11QCn1OvAWybX/WYM4ZiEGzMp/fpny+PNrW5IB3DBwP/UE\n3nt+gy0UJL7nXtx88AWs9kxM/5ouEoZBqNBDQdUOFElRNDHEZHOWyEvtaZhX/v4LzBT/Ch02eGSP\n1mRD8rWfYBb7Cc+6jOgJJ3P2U5vSvmb5D6o6vrcsi1abDXtZOZ4Rsm4vRgbZnCW6ke3y6XVuSp6q\nro03GuTiNX/Af+sL2CyL6LEzCV16JdaYMUD6WjiVftc31zBMEqPH4CstlSbkYliRoJ+D8rl5SU/i\n0SixmmqKOtXLOWGK/5sCaZbFwfqfnPeP5YwJbsHYYTLBeQtI7Pftbc6zzWs6OX6Kn0jCIFZURGHF\neGlCLoYlCfo5SLbLb6ujKfmWzd3SMNvX4Fe/8QmnPHs3+3z9bwxXAaELZxH5yZngcnU7X/trnu9U\nPvm43YqYOrkYs7wCn7fnjVhCDCUJ+jlItst/Y5um5KnaFkajHPqX33Lsow9hi8WIHXAgoTnzMSd0\nv1Hb2QGTvB3Bv9U0sUrL8I0eMxhvQYgBJUE/B8l2+bam5LU1uAMBvA57yno5zndW41tyE44N6zFL\nywjMnkP8sCMybkgeThjEpS+tGGEk6OegfN8uH2reCg31yabkKaph2poa8d5xK+6/vIxltxP50emE\nLrgYfJl1oeoonVA1XvrSihFHgn4Oytft8u1pmGmbkhsG7mefwrNsKfZAgMTuexCcvwBD7Z7R+Udi\n6QQhupI8fTHiWZZFsLER55YmfOkakuu1yZz7jz/CLCoifPGlRE86NW0D866ChkFi9FhJwRTDguTp\ni7wVCQQw6muTN2pTBfBgAO999+B+8k/YTJPo0ccQunw21tjSjM4fNxKEPT7cFeMpTJHJI8RII0Ff\njCjbFEcrKeC4XX0cPLmo+81Xy8L12qv4bluCvbEBo2oSwTnXkPjutIyu07GbVpZyRI6RoC9GjG7F\n0bbGeGBNDJfdvk3NG/umjXhvWUTBW29guVyEzruQyBk/gwzr1YcSBvHRY/CVlclSjsg5EvTFiLHy\njV6Ko8ViFP7+MTwPPYgtFiX+nWkE516DWTUpo/N3zsqRwmgiV0nQF8OekUgQqammdnM45ePVLXGc\n772Lb9GNOL5ehzlmLMErriN21Pczyrlvz8pxVFbhK8osbVOIkUqCvhjWAo2NODc34k9THM0fauay\ntx7F/+9XsWw2Iqf+gPBFl2IVF2d2fsPAGFNKUWlmN3aFGOkk6IthZ/XHdTz/xlfUbA5RWezkhN1L\nOGCSd5tCZzbL5Kj//pWzVj1KcTRAYjdFcN5CjD32zOgakYRBrLiYwvIKHFIYTeQRCfpiWHnrwxru\nf/6Tju83puhK9f7f/81pz9zFlBpN3OMleOXVRE/7ETh7/+dsGAbBAjeuCVX4CgsH500IMYxJ0BfD\nRmjrFl54/bOUjz2/toUDSuGwZ+/j2D/9HpthEDvsCIJXzsEaNy6j8wcsC2tcOb5Rowdy2EKMKL31\nyHUCy4HJQAHwa+Bj4GHABD7UWs9qe+51wEwgDszWWq8ZtFGLnBKLRIjX1uCLRaltTaR8zsT3/knJ\nb5bjqKvDqJxAaM584tO/l9H5wwmDuN+Pp7wiJwujScMc0Re9zfTPABq11j9VSo0GPmj730Kt9Sql\n1DKl1EnAemCG1nqaUqoKWAF8d1BHLkY80zQJ19VS0NKSrHOf4mZtWUs9F7z2AAd88Q6W00n4rHMJ\n/+wcyGBppiMFc2IFH36xlRdeXJNzgVEa5oi+6m3a8zjw807PTQD7aa1XtR17CTgKOAh4BUBrvQFw\nKKXGDvxwRa4INW8l8sVn+ENBCjvVuT9hSnL3q8NIcOqap/jNw5dxwBfv0DR1X5of+yPhCy/pNeBb\nlkWLZRGtmIBv0g6898VW7n3uIzY2BDEtqyMwrv64blDf4/bQU8McIVLpcaavtQ4BKKWKgSeAa4El\nnZ7SCpQAxUDn/nGBtuPde8qJvBaPRom1LeU4Uiy1HDDJy+i1/2HS3YuYUP81rb4S9AVzqDrjFMwM\ncu7DCYNYySh85eUdu2lzuZOYNMwRfdXrjdy25ZqngKVa6z8qpRZ1ergY2AK0AP4ux7cO5EDzVa6s\n11qWRbCujoLmrd1aFrazNW/Fc/edTF/5LACRk04lcfGlVJWU9Hr+zks5RV3KLeRyYJSGOaKvelze\nUUqVAy8D87TWj7Qdfl8pNaPt62OBVcCbwNFKKZtSahJg01pvHqxB54v29dqRviwRCQQIffk5/kAL\nnlQtCy2Lgheeo+T00yhc+SyJXXal5d7lhK65FquXgG+aJi02G7HKKnyTdsCVor5OZWnqnrW5EBhn\nTp+c5nh+NMwRfdfbTH8BMAr4eVt2jgVcAdyllHIBnwBPaq0tpdQq4C3ABswaxDHnjZG+LJGIx4nW\n1eIOBvE5HSlLIti/+hLfohtxffAelsdD6PLZRH54ekY59wHDwBxbiq+XMsm53EksXxvmiOxJE5Vh\n7LybX8NM8f+Pw27j/nmHDcGIMpNJUxMiYTwPPUjh7x5N5tzPOJTQVXMxyyt6PX/UMIn6fBRWjM94\nN21ymUwCo8gN0kQlR23v9dps7x90fl3FmEKO3sXLoVWe1E1NANcbq/DesghHTTVGxXhCV88jftCM\nlM/tzDAMgq4CnOMn4POmXrJJZ9rUcgnyw1Cu3LMaSSToD2Pbc1ki23zvbjXum8I83BTGw9htatwD\n2Orr8N22hIK//w3L4eCLk37CnXueyrpaB5Wv1HLCFH+310DyL4eAZWGVjcM3ekw/36kYLmSPwdCQ\noD+Mbc/12mzvH6R7XUeNe4BEAvcTf8L7wD3YQiHie+/DW2dexeJqP7RVS97YHO9WYwfadtOWlOAZ\nV56Tu2nz2Ui/ZzVSSdAf5rbXskQ2aY3BLZupbkz9eHVLclet48P/JhuSf/YpZkkJwYXXEZt5An/4\naz3Jih3bav9lkTAMQoUeCiZW4Muw45UYWXI5lXY4k6AvgL7dP+hcKydVjXuAnV0RvItuwP3MU9gs\ni+jxJxGadRlWW7Gz9l8KXVW3xGm1LCivwFcyqp/vSgxnssdgaMjfywLILN/bNE2CNdXYvl6H30jg\ncDg6yiZ0sCwO+eQf3LDsYgqfXoExeUdafnM/wWuv6wj4AJV+V8rrlY9yU7jTLngl4Oc82WMwNGSm\nL4De7x+Etm6BxgaKAVunDVbt6+/Pr22Bdeu4/O/3MeWrf2O53YQuuYzI6T8BV/cA37khyjbHD95F\n1u7zhOwxGBqSpy96FA2FSNTXJWvlpM25j+B59CEKf/sItnic2PcOJnT1PMzxlT2e++31IZ77pJma\n1gSVpT5J1xMiQ/3J05egL1IyDINIbQ0FgcA2VTC7cr39Jt4lN+HYtAljXDmh2XOIH3JYrw3JE4ZB\n0F2Ie3wlroKCgR6+EDlNNmeJARXc3IS9qRG/3Q5pAr6tsQHvHbfi/usrWA4H4R+fQfi8C6GXTVOW\nZRGw2bCVj6cog0JquUY2I4mhJkFfdIhFIsRqqimKx9Iv5RgG7qeexHPv3diDQRJ7fivZkHzX3Xo9\nf6qyx/lENiOJ4UCCvsCyLEJ1tbiamylp62CVimPtJ/hu/jXOtZ9gFhcTnH8t0RNPTlkmuTPTNGl1\nOCnYoYqiPG5GLpuRxHAgQT/PhVtasBrqKLasbbJyOrMFWvHcuwz3iseTOffHzCR02RVYY3pvjhZK\nGMTHjKWorGyghz7iyGYkMRxI0M9T8ViMWG0N3kgYpyN12WMsi4K/voL3jluwNzVhTNqB4LwFJPb/\nTq/n77hRW7UDRXKjFpDNSGJ4kKCfZzov5fh7WMqxb9iAb8mNuN5ZjVVQQOj8i4ic8TPIIIC3miaU\nV1AkG6y2kct1/cXIIUE/j4Sat0JDfbcNVtuIxSh87GE8jz6ELRYjdsCByZz7iVW9nj+aMIgWF+Op\nGC8brFKQzUhiOJA8/TwQj8WI1lTji0aSSzlpON99B9/im3Cs/xqztJTQFVcTO+IosNl4e32IlWtb\nqG6JU+l3bVMG2bIsWm02HOXjKSwq2l5vS4i8Neh5+kqpacBNWuvDlFL7AiuBT9seXqa1fkIpdT1w\nHMnSibO11muyHZQYGO0drFybm3rMyrFtbsJ71+24//wilt1O5EenE7rgYvAlA/jb60PblEzoXAZ5\n70o3sVGj8Y0bl5NpmJJXL3JNr0FfKTUXOBMItB3aD7hFa31bp+fsCxystZ6mlKoCVgDfHYTxDms9\nBYjtHTwigQBGfS1+00y/lGOauJ95Cs89S7G3tpLYfWoy537K7ts8beXalpQvf1a3ss+BUyjK0dLH\nklcvclEmM/3PgVOAx9q+3x/YTSl1MsnZ/mzgIOAVAK31BqWUQyk1VmvdvaJWjuopQADbLXgYiQSR\n2poem5EDOD7VyTr3H32I6fMRnDOf6MmnpfxrIF0Z5LqWOK4cDfggefUiN/V6t01r/TSQ6HRoNTBX\na30I8CVwPVAMNHd6TgDIqz32PQWInh4bKJZlEWhsJP7l5/ijEdzpZvfBIN47bsV/9hk4P/qQ6FHf\np/mPK4ie9qO0yz/pyiDneqqh5NWLXJRN9s4zWuv2AP8McFfbfzsXVi8GtvZzbCNKTwEi3b3ygQoe\nkWAQo66WYiOBPd2NWsvC9fe/4bttCfaGeowJEwnOvYbEtOm9nj9dGeRcTzWUvHqRi7LJq3tZKfXt\ntq+PAN4F3gC+r5SyKaUmATat9eaBGuRIUFmautDY+LG+Hh/rD8MwCG7aiHPjBvyWmTZN0l69iaI5\nV1C8cB62rVsIn3M+zb97PKOAHzcSTN11DOcfP4WJZUU47DYmlhVx4Yl75PwShzT5ELkom5n+xcBS\npVQUqAUu0FoHlFKvA28BNmDWAI5xROht481Ab8oJbm7C3tiA3+FIWwmTeJzC3z+GZ/kD2GJR4t/+\nLsG512BO6v26HdUwKybg8/uZPhGm79lzffxcI3n1IhdJnv4ASmbopA4QPT3WF7FIhHhNNb6eKmEC\nzvffw7foBhzrvsIcPYbQFbOJHX1sr3XuAaKGmdxkVV4xojZZSXqlyBfSRCUPmKZJuK6WgpaWHpua\n2LZswbv0dtwvPo9lsxE95TTCF87C8vvTvqadZVm02u04Kypx91IXf7jpmj3VLh+WoUT+kSYqOS7U\nvBXq63psaoJpUvD8c3jvvhN7SzOJXRXB+Qsx9tgzs2skDOKjx+ArKxuRm6wkvVKIzEjQH8ZikQjx\nutpkf9oellkcn3+Gd/GNuP7zbyyvl+CVVydTMJ29/99rmCaBAveIr4Yp6ZVCZEaC/jDUeSnH73Sk\nb1ISDuNZfj+Ff/gtNsMgdtgRBK+cgzVuXEbXCRgGZtk4ikaPGcDRDw1JrxQiMxL0h5H2WjmOLU09\nZ+UArlX/wHvrIhy1tRiVEwjNmU98+vcyuk7CMAh6PHjGT8CRwV8DI4GULRYiM7nxE58DQs1bobEB\nv2Vh6yErx15bg/e2JRS8/ncsp5PwWecS/tnZUOjJ6DoB08TKwVr3kl4pRGYke2eIRQIBjIb6XlMw\nScQp/NMf8DxwL7ZIhPi++xOctwBz8o4ZXSduJAh5fcnZfU/XEUIMe5K9MwLFIhHi9XUURsL4HOnL\nHgM4//NvvItuwPnF55ijRhGcu4DYsTMzyrkHCFgWVEygKIO0TSFEbpOgv50ZhkGkrhZXSwt+l7PH\nYG9rbsaz7C4Kn30agMiJJxO+5HKsksxq2bXP7r2VE0fUJishxOCRoL+dtN+kdbbfpHX18NFbFgUv\nPo936e3Yt24lsfMuhOYtJLHX3hlfr9U0scnsXgjRhQT97WCbhia9rKfb132Fb9GNuN7/F1ZhIaFZ\nVxA5/cfgTF3euCuZ3QsheiJBfxDFYzFidbW4Q6EeG5oAEAnjeXg5hb97FFsiQWzGoYRmz8GsGJ/x\n9QKWBeWVFGW4/COEyD8S9AeBYRhE6uuS6/bOnvPtAVxv/hPvLYtwVG/CKC8ndPV84gcfkvH1EoZB\n0OuVzBwhRK8k6A+gb9btN+N39FAnp42tvh7f7UsoeO1VLIeD8E9+Svic86EPxc5yNe9eCDE4JOgP\nkHBLC2ZDXdvmql7W0hMJ3E8+jvf+ZdhCIeJ77U1o7gKMXXbN+HoJw+D1ugR/1UFqmqqllLAQIiMS\n9PupPd/eGwnjdPSybg84PvoQ38034PxMY/pLCC74ObHjT0xfXyeFgGmypqWAR96o7jg2mM3WhRC5\nQ4J+ltrz7QtaW5Pr9r2spdtaW/HcsxT30yuwWRbRmScSuvRyrFGjM75mR82cyon89eF3Uz5HSgkL\nIXqSUdBXSk0DbtJaH6aU2hl4GDCBD7XWs9qecx0wE4gDs7XWawZnyEOrr+v2WBYFr7yE947bsG/Z\nTGLHnQjNXUBi3/36dN2ua/dSSlgIkY1eg75Sai5wJhBoO3QrsFBrvUoptUwpdRKwHpihtZ6mlKoC\nVgDfHaxBD5U+rdsD9vVf41t8E65338FyuwldfCmRH58Brsxy7iF9zRwpJSyEyEYmM/3PgVOAx9q+\n319rvart65eAowENvAKgtd6glHIopcZqrZsGesBDob2ZiTcayWjdnmgUz6MPUfjYw9jicWIHfo/Q\n1fMxKyekfPrb60OsXNtCdUucSr+LE6b4OWCSN7mrNk3evZQSFkJko9egr7V+WinVOZJ0jnitQAlQ\nDHQO8IG24yM66CficaL1dRQEAhmt2wM433kb3+KbcGzcgFk2jsDsOcQPPTztL4q314dYtvqbj2lj\nc5xlq5sIuVwcfOCUtLtqpZSwECIb2dzINTt9XQxsAVoAf5fjW/sxriFlmiah+nqcLVt7bWbSztbY\ngPfO23D/5eVkzv2Pf0L43AvB1/Nyy8q1LSmPv6oDHHJQz0tI06aWS5AXQvRJNkH/PaXUDK3168Cx\nwN+AL4CblVJLgCrAprXePIDj3G6Cm5uwb27CD73WyQHAMHA/vQLPPUuxB4Mk9tiT4LyFGLupjK5X\n3RJPeVxuyAohBkM2QX8OcL9SygV8AjyptbaUUquAt0gu/8wawDFuF8miaHUUGQnsdnvadfbOHGs/\nwbfoBpyffIxZXExw3kKiJ53Sp5z7Sr+Ljc3dA7/ckBVCDIa875yViMeJ1tbgDoVwty3jdF1nb3fx\ntLEcMMmLLdCK575luFc8gc00iR5zHKHLrsQaM7ZP144mDN5sgodfr+722IUn7iFLN0KIlKRzVhYs\nyyLY0IBr65Zu+fbp1tmf/6SZGZ/+E+etSyjc0sSm0ZU8ceKl7HbcwRwwJvN6OZZl0Wq345gwkRmq\nCPeo0XJDVgixXeRl0A+3tGA11rfVt+++FJNqnb1iaw3nPXUfReveJ+Zw8dsDf8yKb59Kwuni1ba/\nCrou/6S8tmESGzUaX1kZtraMHrkhK4TYXvIq6EdDIRIN9b3m23deZ3cm4pz67tP8aPWTuI0YH++8\nH7fPOJ+a0dvWuX9+bUuPQd80TVqdLtwTJ1Hkdg/cmxJCiD7Ii6CfiMeJ1tVSEAxmlG9/whQ/y1Y3\n8a31/+WSV+9h4pZNbPaN5uPzFvILx96YdP9lkS4LByCcMIiNGUtRWVm/34sQQvRHTgd90zQJ19fh\nbG7OqJlJu+lFEfZas4yJq17GxMbfv3s8sYtn8e0p46h8pTZltk2lv3tpBdM0CThduHaooqiwsN/v\nRwgh+itng35wcxP2pkaKbTZsGQZ7TBP3c8/g+c2djG5tJTFld4LzFrLX7lM7ntL+V0BXx0/ZtgF5\nyDCJjxlLUWlpv96HEEIMpGEX9Fd/XMcLb62jujGUVWOQcEsLZmM9RYbRp8bgjs8+Tebcf/hfTJ+P\n4NXziZ5yWreloPZ1++c75fAf3ymHP2EYBN2FFE6agLsPhdWEEGJ7GFZ5+qs/rktZRCyTnPVuzUwy\nFQrheeAeCh//IzbDIHrEUYSuuBori/X3gGlilY3D24ca+UII0Vc5k6f/wlvr0hxP3xikt5u0aXfW\nWhauf7yG97YlOOrrMCZMJDj3GhLTpvd53DHDIOz14RlfKY3JhRDD2rAK+n1pDNJeFM3VvDXtTdp0\nFSw99TV87493UfDGKiyXi/DZ5xH+6dmQxc3WVsvCPn4CRcXFfX6tEEJsb8Mq6GfaGCS4ZTO2xgb8\nvdyk7bqz1mEkOPlfz/G9u/5EQTxKfP/vEJx7DeYOk/s81riRIOwrxjO+sk/3DoQQYigNq6DfW2OQ\nrkXRetM5d37qxo+45NV72aFpPVu9JcSu/V9iRx/be0OUFAKWBRUT8Pn9vT9ZCCGGkWEV9NM1Btlv\nl9EEN6ynMBzC53BkXMWy0u+ipbaJn73+CEd/9ComNl7c6xhePe4srv3+rn0eX+fG5LJ2L4QYiYZV\n9k5XHZurtm7F6+rj7yfTZONjK5j88N34I618WTaZ3xxxMbpSdVTL7IuAYWCVV+Bta0wuhBBDJWey\ndzpr31zlt9uhl4DfNUPnx74mpv/2Nvb64H0Sbg9Pfv88fj/1WCpGFXJxirr4PTFMk4DbjWeHHXE4\nh+3HJYQQGRnSmf5Jc56zum7A6uvmqs4ZOu54hNPffpyT//UsTtMgdujhhK68GrO8IqvxBQ0To6wM\n3+gxWb1eCCEGw4id6ZuWxcaGIPc+9xHxWJR9R5vfVMDMcN2+PUPnO1+u4cK/3Ud5SwN1/nGsOOFi\nTr/0xKzGZRgGgbZdtYWyq1YIkUOyDvpKqff4pvn5V8B9wB1AHPiL1vqXfTnfK2+u45Dvj++1AmZX\nsY3VLHjtQQ78/G0SdgePf/c0Hp/2IxIFbk7v05mSQgmDRGkpRWOlZo4QIvdkFfSVUm7A0lof3unY\n+8ApWut1SqkXlFL7aK0/yPScNa2Jvg0iEafw8T9y9yP3UBiL8OGEqSw74iLWl04CoCpF1cuedNS7\nnzwJn9S7F0LkqGxn+nsDPqXUy4AD+H9AgdZ6XdvjLwNHABkH/VSlidNx/vc/eBfdgPPzz4gVl3D7\nYefz6tTDt8m571r1sidS714IkS+yDfohYLHW+kGl1K7AS8CWTo+3Ajv25YSZBGlbczOeZUspfPYp\nACInnET4ksvYrcXNp2mqXvYk2avWIfXuhRB5I9ug/ynwOYDW+jOlVDPQOcWlmG/W+9Oy22BCJkHa\nsih46QW8S2/HvmULiZ12JjRvAYm99wXggFGZ9aftLGKYREtG4Rs3rqNXbT7ob+lqIcTIlm3QPwf4\nFjBLKVUJeIGgUmpHYB3wfeAXvZ3krh/uQJFl9vgc+7qv8C2+Edd7/8IqLCQ063Iip/8POLPLqknO\n7u04J1ZR5O3bL4qRrmvp6vbMKUACvxB5Itug/yDwkFJqFWACZ7f99/eAHXhFa72mXyOLRPA8spzC\n3z6CLZEgdtAMQlfNwxw/vvfXpjtlwiBWUoK3vCKvZvftsildLYTILVkFfa11HDgjxUN9L0afguut\nN/HechOOTZswyssJzZ5L/JDDsj6fZVkEbDYcVZPw5dnsvrO+lK4WQuSmYVVXwFZfj/eOW3D/7a9Y\nDgfh/zmT8LkXQD8CdTRhEPX78VaM7za7z7f17UxLVwshctfwCPqJBO6nnsB77zJsoSDxPfciNH8h\nxi59r4TZWatlYZ8wEV9RUbfH8nF9u7fS1UKI3DfkQd/x8UfJhuR6LWaxn9A11xI94eSMyzCk0t6+\n0Fs5IW39nnxc305XujpX368QorshDfqjl9xE0YrHsVkW0WNnErr0Sqwx/StulmxwUklRLw1O8nV9\ne9rUcgnyQuSxIQ36xU/+CWOHyQTnLSCx37f7da64kSDk9eEZPyGjBieyvi2EyEdD2tx1y6wraH7s\nj/0O+AHTJFpeSdHESRl3tJo5fXKa47K+LYTIXUNaT3/96g+s3jZn9SRhGAS93oxn910ls3dkfVsI\nMbL0p57+iA36AdPEGlcu7QuFEHlnxDZRyYY0JxdCiOyNqKAfNEzMceUUjRo91EMRQogRaUQE/c7t\nC53SvlAIIbI27IN+0DAwSssoGjN2qIcihBAj3rAN+qZp0uoqoHDSZGlOLoQQA2RYBv1QwiA+tpSi\nUmlOLoQQA2lYBf3OzcmLpDm5EEIMuGET9KU5uRBCDL4BDfpKKRvwG2BvIAKcp7X+sqfXSHNyIYTY\nfgZ6pn8y4NZaH6iUmgbc2nYspahhYpSMwldW1mv7wnxreCKEEINhoAuuHQT8GUBrvRrosZJa8U47\nUTRuXEYB/97nPmJjQxDTsjoanqz+uG7ABi6EEPlgoIO+H2ju9H1CKZX2GpmWUeip4YkQQojMDXTQ\nbwGKO59fa519Gc02+drwRAghBtpAr+m/ARwPPKmUOgD4b09PzrRSnGlZ/wG+1fW4YVr/KSsr3jub\ngQohRD4a0NLKnbJ39mo7dLbW+tMBu4AQQoh+GdJ6+kIIIbavIW2XKIQQYvuSoC+EEHlEgr4QQuQR\nCfpCCJFHtnvBtWzq8+QapZQTWA5MBgqAXwMfAw8DJvCh1nrWUI1ve1NKjQPeBY4EDPL0cwBQSl0D\nnAi4SP6cvE6efR5tPx+PkPz5SADnk4f/LtpK2dyktT5MKbUzKd6/Uuo6YCYQB2Zrrdf0dt6hmOl3\n1OcBFpCsz5NvzgAatdYzgGOBpSQ/h4Va60MAu1LqpKEc4PbS9gN+D9C+Ay8vPwcApdQhwPS2n41D\ngUnk5+dxHODQWn8P+D/gBvLsc1BKzQXuB9przHd7/0qpfYEZWutpwI+BuzM591AE/T7V58lRjwM/\nb/vaTnI2s5/WelXbsZdIznrzwRJgGVAN2MjfzwHg+8CHSqlngOeA58nPz+NTwNm2KlBCchabb5/D\n55gDzf4AAAHtSURBVMApnb7fv8v7P4pkLH0FQGu9AXAopXrtKzsUQb9P9XlykdY6pLUOKqWKgSeA\na0kGvHatJP+x5zSl1FlAvdb6L3zz/jv/W8iLz6GTUmB/4AfAxcDvyM/PIwDsCKwF7gXuJM9+PrTW\nT5OcDLZL9f6L2TaWBsjgcxmKYDso9XlGGqVUFfA34BGt9R9JrtW1Kwa2DsnAtq+zgaOUUq+RvMfz\nKNC5i06+fA7tmoCXtdaJtp3sEbb9Ic6Xz2M28GetteKbfxcFnR7Pl8+hs67xYQvJWOrvcrzXz2Uo\ngv4bJNfsyKQ+Ty5SSpUDLwPztNaPtB1+Xyk1o+3rY4FVKV+cQ7TWh2itD9NaHwZ8AJwJvJRvn0Mn\n/wSOAVBKVQI+4NW2tX7In89jM9/MYLeSTDh5Pw8/h87eS/Fz8SZwtFLKppSaBNi01pt7O9FQtEt8\nmuTs7o22788egjEMtQXAKODnbXffLeAK4C6llAv4BHhyCMc3lOYA9+fj56C1fkEpdbBS6h2Sf85f\nDKwDHsizz+N2YLlS6nWSWUzXAP8i/z6Hzrr9XGitLaXUKuAtkv9eMspokto7QgiRR/LqBqoQQuQ7\nCfpCCJFHJOgLIUQekaAvhBB5RIK+EELkEQn6QgiRRyToCyFEHpGgL4QQeeT/B5a7tDb/fS9iAAAA\nAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fit = psy.plot.linreg(ds, ax=raw.plotters[0].ax, name='y', coord='x', \n", " legendlabels='fit', color='red')\n", "fit.share(raw[0], keys='legend')\n", "fit.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The shaded red area displays the 95% confidence interval. To access the data for the fit, just use the ``plot_data`` attribute:" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "\n", "array([[ 32.122999, 36.32817 , 40.533342, ..., 440.024632, 444.229804,\n", " 448.434975],\n", " [ 6.452861, 11.076892, 15.702347, ..., 406.556079, 410.383088,\n", " 414.318081],\n", " [ 55.35959 , 59.295839, 63.059391, ..., 467.569211, 472.145364,\n", " 476.721518]])\n", "Coordinates:\n", " * x (x) float64 0.0 1.01 2.02 3.03 4.04 5.051 6.061 7.071 8.081 ...\n", " * variable (variable) " ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fit.update(legendlabels='%(yname)s = %(intercept).3g + %(slope).3g * %(xname)s, R$^2$=%(rsquared).3g')\n", "fit.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "To improve the fit, we can also force the line to go through a given fix point. For example here, we know, that the\n", "fit crosses the y-line at 30:" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX0AAAECCAYAAAASDQdFAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3Xd4VGX68PHv9MxMkoEUUiCI9VGxl0VdFbGtiii67v5c\ne1fEBlIEC65rQ7CjYG+7rruKvay6a4GXRcS1InCsCCSBFEKS6TPnnPePCTEhk2QySUi7P9flRebM\nOWeeGZM7d55yPxbTNBFCCDEwWHu6AUIIIbYdCfpCCDGASNAXQogBRIK+EEIMIBL0hRBiAJGgL4QQ\nA4g9lZOUUkOAz4CjAS/wBvBdw9PzNU17USk1CzgBiAGTNE1b3g3tFUII0QntBn2llB1YAAQbDu0H\n3K1p2r1NztkXOEzTtFFKqRJgIfCbbmivEEKITkile2cuMB8oa3i8PzBWKfWxUuoxpVQmcCjwHoCm\naesAm1IqtzsaLIQQIn1tBn2l1HlAhaZp7wOWhv+WAVM1TRsN/ATMArKA2iaX+gFfdzRYCCFE+trL\n9M8HjlFKfQjsAzwDvKNp2hcNz78K7AvUAdlNrssCNndxW4UQQnSSJdXaO0qpD4DLgGeBKzVNW66U\nugIYBvwDmAMcA5QAr2matm9796ysrJfCP0L0Azc9sYz1lYEWx4flZ3LLhQN3eE+PxwmVl+INhbDb\nbACc/9I6jCSRz2aBJ08rSXof+1df4LnrDuw//YgxOIc/1GzyLTTNunTalNLsna1cBjyklIoAG4BL\nNE3zK6UWAUtJdAFNTKcxQoi+qawqmPR4eXXLXwQDRXBzDZaKjfhsNmgI+ADF2Q7W18ZanF+c7Whx\nzFJbi/vhB8h4/VUAwuN/T/CyiXDckWm3K+Wgr2la01f5bZLnbwFuSbslQog+qzjPkzTTL8r19kBr\nWrds5UbeWrqGsqogxXkexh48glG7F3Tpa8RjMcLlpXjD4cbsvqlxu2Yzf1l1i+Mn7tqkh9w0cb7z\nFp4H78W6eTPxHXYkOP164nvtTWcrI6eT6QshRDNjDx7BI69/m+T4dj3QmuSWrdzYrI3rKwONj7sq\n8AdqNmGrrMRnszbL7ps6aLgHgDdX11FWF6M428GJu2Y3Hrf+sgbvnDtx/G85ZkYGwYlXEz79T2Bv\n+ZdAOiToCyE6bUvQfGvpL5RXByjK9TL24O26PIveIp2M/a2la1o5/kun2xmPxQiXlZIZjWCztT8T\n/qDhnsYg3ygSwf3c02Q8+xSWWIzoIYcSnDIdo6i4U23bmgR9IUSXGLV7QbcF+abSzdi7a9whsKka\nW1VVIru3plfZxv7Zp3jn3Ilt7S8Yefn4J08ldsSRYLF0qm1JX6vL7yiEEN0o3Yy9q8cdGrP7SBhb\nK1057bFs2oTnwXtx/ettTKuV8B9PJ3jJBPBmpnW/VEjQF0L0Kelm7F057pDI7itbzMxJmWHgevM1\n3PMewFpfR3zX3QhMvx591906fq8OkqAvhOhT0s3Yu2LcIR6LES5d39B3n152b/vpRzyzb8Px9VeY\nHi+BSVOI/P6P6f3ySIMEfSFEn9KZjL0z4w6dzu7DIdxPPk7G889h0XWiY44icM0UzCFD0mpPuiTo\ni21m9eqVhEIhVq5cwZlnntvTzRF91LaeKRSLRok0zsxJLxt3LF2CZ+5sbGWl6IWFBK+dTuzQw7u4\npamRoC+2mdWrV3HCCeNYuvT/EQqFcLvdPd0k0Udtq5lCgeoq7FVV+OzpZfeWqko8992N6z/vY9ps\nhM48h9CFl0APfu9L0O8GhmEwe/atrF37C1arlSlTZjBixPbcffed/PDD9zidTqZPv4GhQ4eldf+a\nmk1ceOHZ3HffwwwfnvxP2m+/XcGCBQ/y4IOPEIvFuP32P1NWVorXm8m1105v8doVFRtZsmQxw4aV\nsNNOuzB48OAOt+uLL/7HTTfNYPvtdwAgEAgwdOgwbrrpL9jtdsaP/z2GYWAYZpsB3zTNNj+r1t7P\nBRecSWZmFgBFRcXMmHFT0s9DiPbEolEipevJisf4tCzCG00WUo1rspCqVbqO65WFuBfMwxoIENtj\nL4LTZ6LvtHOn2vXJ2iBvrK4jOumVTaQZvyXod4MlSxZhsViYP/8Jvvjifzz66EMcd9yJRKNRFix4\nkm+/XcG8efdyxx13d/je8XicOXPuICMjo9Vznn/+Wd59923c7sQ35uuvv4LH4+GRR55i7dpfuPvu\n2dxzz4PNrvnll595/fWXyc72cfnlV6cV9AH23/9Abr75tsbHf/7zDSxZsojRoxNVPD788D+cffZ5\nxONx7Pbk336LFn3U5meV7P3ceWfi+QceWNDu5yFEW/yVlTg2VeOz2/hkfbhZyYT1tbHGx60Fftt3\nGt47b8O+6luMrCwC02YSOfmUtOfwb/HJ2uCvbbFY0h71laCfxMsvv8g333zFrFm3ctttNzNy5B6M\nH38aAB999B8WLvwnliaLJi6//Cp23XX3xseHHXYEv/1tor9uw4ZysrKy+eabrxg16hAARo7cg9Wr\nV7X6+k8++Sg5ObmMH//7Fs899ND9nHLK73nuuadbvX7o0BJuv30uf/lLItNds+ZnDjoo8drDh2/H\nL7/83OKaXXbZlWOOOZ7a2s3suONOKX0WyTStCxKLxaiuriIrK1FT5L33/sXnny/ns88+ZerUGa1+\nll9//WWbn1Wy9/PDD98RDoeZPPkKdN3gkksuZ+TIPZJ+HkIkE4tEiJaXkRWLYrUnYuobq5MXsnxz\ndV3LoB8M4nlsAa5//h2LYRA59jiCV0/GzOma/aRaa0tH9eqg7735BlxvvNql94yMG0/g5lvbPOfU\nU//AZ599yu23/5l4PN4syB1xxFEcccRR7b6O1WrltttuZvHij/jLX2bzwQfvk5n564ILm82GYRhY\nm/z2/89/3ufVV19iw4YNOBx2Pvjgfc455wIOOCBRmvbtt99g8ODBHHjgQTz77FOtvvbo0WPYsKG8\n8fHOO+/Cf//7/zjssCNYseIbqqoqMU2zWbD1+QZxxhlnd+izSObzzz/jqqsuY9OmTVitFk4++VT2\n2+8AAI499jiOPfa4xnNb+yxfe+3lNj+rZO/H5crgjDPO5sQTx7Nu3VqmTLmKv//9ZaxWa4vPQ4it\n+SsrcdZsInurVbVldS2rYSY77lj0EZ6778JWsRF96DAC02YQ/81BabVlSxfO1t1JrbWlo3p10O9J\nZ555LhMmXMATTzzX7PiW7HQLi8XSItPf4vrrb6amZhMXX3wuo0cfSTD469zirQM+wFFHHcNRRx3D\nk08+Sm5uHieffGqz599++w0sFgvLly/j+++/49ZbZzF79j0MHpzT5nsZO/YkfvnlZyZOvJg999wb\npXZrFvDT/SyS2dK9U1dXy6RJV1BUNLTVc1v7LD0eb5ufVbL3s912Ixg2LFGLvKRkONnZPqqrq8jP\n37bT4UTfEg2HiZWXkRWPYU1SM6e9MsjWDeV47pmDc/HHmHY7ofMvInTO+dBG92tbmnXh8Gt3Usww\nKPQ5KdscTeu+TfXqoB+4+dZ2s/LuEIvFeOCBu5k6dSZz5tzBww8/3tj/nEqm/+67b1NRUcHZZ5+H\n0+nEarWy6667sXTpEsaMOZoVK75p1oWSqnnzHm38+sorL2Xq1JntBnyAVatWsv/+v+HKKyezevWq\nDmW9bX0WbcnO9nHjjbdw1VWX8fTTz5OT5E/c1j7LjRs3sGTJ4lY/q2Tv5803X+XHH3/k2munU1VV\nSSgUJDc3r9l1nS1JK/oP0zQJbMnu7TawWpNm2K2VQR63s4eM55/D/fgjWEIhYvvuT2DaDIwR23eq\nXa114bzzc4Rxh++cdH1CR3VuZKGfWrDgQQ499HDGjRvPwQf/lgUL5nXo+tGjj+T77zWuuOISpky5\niquvnsLRR/8Op9PJhAkX8NBD93LllZNbvf6CCy5pkeVvrWmmXldXxw03TGv13JKSEv75z+e57LIL\neOKJBVx55aSU30uyz2LTpmpmzZrZ7rUjRmzPH/5wOvfdNzfl1wM4/PAxLT6rpu8x2fs58cTx+P31\nXH75Rdx88/Vcd91NLf6S6shfN6L/iobDBNf8TFZtDe6GvvstGfb62hiG2XzAdsKoXEp8DmwWKPE5\nuC5nA8fMuhTPg/dhOp34b5hF/UOPdDrgQ+vdSRs2hRi1ewGXnjSSYfmZAPF0XyPl7RK7g2yX2Dfp\nus6CBfOYOPHqnm6KECkzTZNARQXOzb8G+y2uf29D0m6cEp+DW48tBMDir8c9/yFcr7yExTSJjD2J\n4BVXYQ5Kb6ZbMq21Y+ttJ/Pzs9LOYFLq3lFKDQE+A44GdOBpwABWaJo2seGcm4CxQAyYpGna8nQb\nJXo30zSTDvoK0VtFQyFi5WVkGzoWe8vZjm0O2Jomzv+8j+e+uVirq9FHbE9g2kzi++7XpW2MxHWO\n3iOHp5dsbPFcV25G027QV0rZgQXAltJ29wAzNU1brJSar5Q6GVgLHK5p2iilVAmwEBi4uyH3c3a7\nPaWxBCF6mmmaBDZuxFVXm5iZ00oXX2sDtnvGq8icNBvnsqWYThfBSycSPvNscHTNLlaQmKjgt9mx\nDS3icJWJKzevW0tMpJLpzwXmAzNIbHq+n6Zpixueewc4FtCA9wA0TVunlLIppXI1TWs5AiKEENtA\nJBgkvqGMbMPA0s5uVlsP2Nr1GKd89ipnfvoitliU2KiDCEy5DqNhhlhX8RsGRm4e3iYTHbq7xESb\nQV8pdR5QoWna+0qpLSN3TT+9esAHZAFNA7y/4bgEfSHENtWY3dduTszMSWEAv+m+tYNWfskVHzxC\nceVajJxc/NdcS/ToY7t0F6tQXCeWnY27oLDFhIPu1l6mfz5gKKWOAfYGngXymzyfBdQAdUD2Vsc3\nd2E7hRCiXc2y+yR992052BflqOWP4nrzNUyLhfCppxG67ArMrKwua19MjxPK8OAcVojX5eqy+3ZE\nyrN3lFIfAJcBc4C7NU1bpJSaD3wA/AjMJtHVUwK8pmnavu3dU2bvCCG6gmmaBDduwFlbS0YHgz2m\nifPtN/HMuw/r5s3Ed96FwPSZ6CP37LL2JfrtbVjzhuDOzm7/gnZ0++ydrUwBHlNKOYBVwEuapplK\nqcXAUhL9/hPTbZAQQnREJBgkXl5Klml2OLu3rvkZ75w7cHz+P8yMDIJXXkP4j3+CFBYgpsqv6+iD\nc/Hm5fWKtSIyT18I0Sd1KruPRHA/+xQZzz2NJRYjethogpOnYhQWdVn7wnGdaFYWGQWFaW++0ppt\nnen3KstWbuStpWsoqwpSnOdh7MEjtsnmCh0VjUY588zTePHF11s95/XXX2Hs2JO6/BtEiP4mEggQ\n31CWNLtvrWDZFvZPl+Gdeye2dWvRhxQQnDyV2OgxXda2uK4TdGXgGFqCN80aPN2pTwf9ZSs3NqtF\nsb4y0Pi4twX+xF9Ubf9yfu65pzj++BMl6AvRCtM0CW4ox1VfjyfJvPvWCpYBHJwZwvPAvbjefQfT\naiV8+pkEL7oUvG1vqN6RtvktFigoxOsb1CX37A59Oui/tXRNK8d/STvov/POm7z11uuYpsmFF17K\nzz//xKJFHxIOh/H5BnH77XO45JLzuOeeeWRmZjJ27FHMm/cYO++8CxdccBaPPvp0Y0GyUCjELbfc\nQH19fbOdn7788nOeeuoxTNMkFAoya9ZtfPnl51RXJ2ra3HrrbObMuZ2Kigqqq6s49NDDueiiy9J6\nP0L0F2G/H2NjeSK7b2XefbKCZRbToOav/8D30TNY6+uJ77Z7YqBW7dZlbQvoOvFBOXjz83tFv31b\n+nTQL6sKJj1eXh1IejxVWVnZ3HHHXEzT5Ouvv+T+++cDMHnylaxevZLDDz+CZcv+S37+EIqLh7J8\n+TIcDgfDh2/XrALlq68uZIcdduLiiyewcuUKPv/8fwD8/PNP3HTTX8jNzeO5557iww//zdlnn88z\nzzzJLbfcQUXFRkaO3JPp008mGo1y6qknSNAXA5ZhGIQ2bsBVV4ernXn3W5dT2K5yDZf/ZwG7l63G\n9HgJTJ5G5NTT0trvNplIXCfi9eIqKCSjC1fpdqc+HfSL8zysr2wZ4ItyO/fn2pZ9Zy0WCzabnVmz\nZuJ2u6mqqiAej3P44WN49tknKSws4pJLLufFF1/AMHSOOOLIZvdZt+4XDjnkMAB2330P7A19j/n5\n+dx77xw8Hg+VlRXstdc+DVeYmKZJdnY2q1Z9yxdffIbb7SUW65rNE4Toa5pl9ykM1m4pp+CKhTl9\n6T8Y//nr2A2d/408lB3uuAEzP7/de6RC13UCDif2kqF4PX1rG84+XVp57MEjWjneueJEW1bI/fjj\nDyxe/BF//vPtTJo0FcMwME2THXbYkfLyMlat+paDDz6UUCjIkiWLOOig3za7z4gRO7BixdcAfPfd\nauJxHYA777yV66+/mZkzZ5GXl99Y591qtWIYOm+//QZZWdnceONfOP30MwmHw516P0L0NYZhECgr\nxV66nixSL4s9btds9v/pMx565kpO++wVqjNzuXn8DWy86Y4uCfimaVJvGATyh+DdfgdcfSzgQx/P\n9Lf023dXcaJhw4bhdnu4/PKLME2T3Nx8qqoqAdhnn/0aNyPZZ5/9WbPm5xablY8f/3tuvXUWEyde\nzPDh2+FyOQE47rixXH75hbjdHnJychrvudde+zB16jVMnjydm2++nhUrvsbhcFBSsh1VVVXk5TXf\nFESI/ihcX4+xsZxsiwU6MBXTUlnJ0QvmcMKH/0G32lh44KksOvZMfrfXkFY3Me+IUFwn5vPhHlKw\nzUsndCWZpy+E6BUMwyBYVkpGIJDou0+VruN6+UU8Cx7GEgwQ23NvgtNnoqexO10yjaUTCotwOJ1d\ncs/OGtDz9IUQfV+org6zYgO+Dmb3Nm0V3tm3Y1+1EiMrm+B11xMZN77Z5ubp2lI6wVZcgjczs9P3\n6y0k6AshesyW7N4dDODsyIyaQADPo/NxvfQPLIZB5LixBK+8BjOna/Z58Ot6ouRxbv/rUpWgL4To\nEaG6OtiS3aca8E0Tx8cf4r1nDtbKCvSS4QSmzSB+QNfs2dSdpRN6Cwn6Qohmuru0ia7rhMpL8QQD\nOGyphyBreTmeu2fjXLIY0+EgeNGlhM86F7qgRHFvL53QlSToC9FPpRO8u7u0Sai2Fio3NmT3KYaf\neIyMf/wd9+OPYAmHie1/AIGpMzC2G9Hp9pimST1gKSzCm+3r9P36Agn6QvRD6Qbv7ihtAr9m995g\nEHsHuk1sK77BO/s27D98jzFoEIFpM4ked0KX7GIV0HXivajk8bYiQV+Ifijd4N0dpU2CtZuxVGzE\nZ7Wm3HdvqavDvWAerldfxmKahE8aT+jyqzB9nc/GI7pBxOMho7CIjC6sm99XDLx3LMQAkG7w7srS\nJrquEypbjzcUSj27N02c7/0LzwP3Yt1UTXz7HQhOn0l873Y34kupPQGnC3tRQZ8rndCVJOgL0Q+l\nG7zHHjyiWbfQr8c7VtokuLkGS2VFh7J767p1eOfegePTZZhOF8FLJxI+82zoZCEz0zTxAwwpwDto\ncKfu1R+0G/SVUlbgMUABBol9cl3AG8B3DafN1zTtRaXULOAEIAZM0jRtebe0WgjRpnSDd2dLm+jx\neKLvviPZfTRKxt+exf30E1iiUaIHHUJwynSMJuXI0xWM68R8g/AWFAyofvu2pJLpjwNMTdMOVUqN\nBm4nEfDv1jTt3i0nKaX2BQ7TNG2UUqoEWAh0zeRZIUSHdCZ4j9q9IK1B2+DmGqwVFfhsqWf39i8+\nxzv7Nmy/rMHIzSVwzRSiRx3T6YHaqK4TdntwDS/C1UdKHm8rKdXeUUpZNU0zlFLnAkcAIRKZv51E\ntj8JOB9wa5p2V8M1/wOO1TStOvldpfaOEP1BPBYjXF6KNxxOObu3bK7BM+8BXG+9jmmxEDn1D4Qu\nuxwzM6tTbUmUTrBjG1JARj8qnbC1bq+90xDwnwbGA6cBQ4HHNE37Qik1A5gF1ABNA7wf8G11TAjR\nj3Q4uzdNnG+9gWfefVhra4nvrBK7WI3co9Nt8RtGonRCTm6n79WfpTyQq2naeUqpIcCnwMGappU3\nPPUq8GDDv9lNLskCNndVQ4UQvUc8FiNcVkpmNIKtla0Lt2Zd8zPeu+7A8cX/MN1ugldNIvyH06GT\n0yZDcZ1YdjbugsI+XfJ4W2n3E1JKnaWUuq7hYZjEYO7LSqkDG44dBXwGLAF+p5SyKKWGAxZN0zZ1\nR6OFED0nsKma+M8/4YvHsKUSZMNh3I88jO/s03F88T+ih4+m9vkXCf/prE4F/LiuU+dwwnYj8BYV\nS8BPUSqf+MvAU0qpjxvOvwpYDzyklIoAG4BLNE3zK6UWAUsBCzCxm9oshOgBjdl9JJxyMTL7p5/g\nvesObKXr0YcUEJw8ldjoMZ1qh2ma1FssWAuL8WZnt3+BaEY2URFCtCuwqRpbVSXeVAdqq6vw3H8P\nrvffxbRaifzxTwQvuhS8qS/y+mRtkDdW11FWF6M428G4XbPZc2gG8cE5A650wtY6M5ArQV8I0ap4\nLEa4dH1D330KAd8wcL36Mu75D2L1+4nvNpLA9JkscQ9vEcDb2sLwk7VB5i9rOQfk4hN34+A9ijrz\nlvoF2TlLCNHlAtVV2Kqr8NlsKc3Msf3wfaI42opvMLxeAtdOJ3LK7/mkNNIsgK+vjTU+bi3wv7G6\nLunxd5atk6DfSRL0hRDNxKJRIo0zc1LI7kMh3I8/QsY/nsei60SOPpbgVZMx8/OB1gP4m6vrWg36\nZXWxpMc7U/hNJEjQF0I08ldV4aiuwmdPLbt3LF6E557Z2DZsQC8eSnDKdcQOPqTZOa0F8NaOB+M6\nhYNclNVEWjyXTuE30ZwEfSEEsUiESFkpWfEY1hQ2JrdUbMR771ycH32AabMROud8QudfCBnuFucW\nZztYX9sywBdnNy+PEInrRDyJ0gnj9MFdUvhNtCRBX4gBzl9ZiWNTdSK7b2+uezyO66V/4nlsPpZg\nkNje+xCcNhN9hx1bvWTcrtlJB2VP3DUx3dIwDPx2B7ZhxXgbZvd0tvCbaJ3M3hFigIqGw0TLSsnS\n4yktbLKtWpkYqNVWY2T7CE68iuiJJ7X/i4LEbJw3m8zeObFh9o7fNDFycqV0QgfJlE0hRIf4Kytx\n1mzCnUoJhYAfzyMP41r4IhbDIHL8WIJXTsIcnH5t+lBcJ+bz4R5SICtp0yBTNoUQKYmGw8TKyxJ9\n9+0FfNPE8dEHeO+Zg7WqEn34dgSmXkf8gPQrpsf0OKEMD86SIrxOZ9r3EemToC/EAGCaJoGG7D47\nhb57a3kZnrmzcf73/2E6HAQvupTw2edBmoE6UfLYhrVoGN6szpVPFp0jQV+Ifi4aDhPb0nff3syc\neIyMF57H/fgjWCIRYvsfSGDaDIzh6c+a8es6ek4e3tzcAV06obeQoC9EP2WaJoGKCpyba1LK7u1f\nf4Xnrtux//gDxuDBBK67gejvjk97F6twXCealUVGQWHKBdpE95OgL0Q/FA2FiJWXkW3oWNrJ7i21\ntbjnzyPjtZcBCJ98CqEJV2L6fGm9dlzXCboycAwtwZuRkdY9RPeRoC9EL7Fs5UbeWrqGsqogxXke\nxh48osPz0k3TJLBxI87azYnsvq0s3TRxvvcOnvvvxVqzifgOOxKcNpP43vuk1X7TNPFbLFBQiNc3\nKK17iO4nQV+IXmDZyo3NVqCurww0Pk418EeCQeIbysg2jHaze+vaX/DOuRPHZ59iulwEL7+S8J/O\nBHt6m4gHdJ34oBy8+fnSb9/LSdAXohd4a+maVo7/0m7Q35Ldu1LJ7qNRMp57GvezT2GJRoke8luC\n107HKB6aVrsjcZ2I14uroJAMR3q/MMS2JUFfiF6grCqY9Hh7VSU7kt3b/7c8sYvV2l8w8vLwT5pK\nbMxRaQ3U6rpOwOHEXjIUr6f1uvii92k36CulrMBjgCKxP+5lQAR4uuHxCk3TJjacexMwFogBkzRN\nW949zRaifynO87C+smWAb62qpGmaBDduwFlb2252b6mpwfPgvbjeeQvTYiF82v8RunQCZmbH58ub\nponfNDHzh+AdnNPh60XPSyXTHweYmqYdqpQaDdxOYg/cmZqmLVZKzVdKnQysBQ7XNG2UUqoEWAik\nv3RPiAFk7MEjUq4qGQkGiZeXkmWabWf3hoHzzdfxPPQA1rpa4rsoAtOvR999ZFptDMZ14lI6oc9r\nN+hrmvaaUuqNhofbATXA0ZqmLW449g5wLKAB7zVcs04pZVNK5Wqa1rK8nhCimVSqSnYku7f+/BPe\nu27H8eUXmB4PgauvJXLaH8He8R5dKZ3Qv6T0HaBpmqGUehoYD/wBOKbJ0/WAD8gCmgZ4f8NxCfpC\npGDU7gWtDtpGAgHiG8raz+7DYdxPP0HG357FEo8THT2GwOSpmEOa3zfZpuNb72K1pXSCrbgEb2Zm\np9+f6B1S/rWvadp5SqkhwHKg6U4JWSSy/zoge6vjm7uikUIMVKZpEtxQjqu+Ho/N2mZ27/jkv3jm\n3omttBS9sJDg5GnEDhvd4rytNx1Ptmet3zAwcvOk5HE/lMpA7lnAME3T7gTCgA58ppQarWnax8Dx\nwAfAj8BspdRcoASwaJq2qfuaLkT/Fvb7MTaWJ7L7NipiWqoq8dx/D65/v5fYxeqMswldeAm0Mqum\nrT1r9y52EcvOxl1QKP32/VQqmf7LwFNKqY8bzr8KWA08rpRyAKuAlzRNM5VSi4GlJAZ6J3ZTm4Xo\n1wzDINSQ3bva6rvXdVyvLsQ9fx7WQID4yD0SA7U779Lm/Vvbm7a0LgbbjZDSCf2cbKIiRC+yJbvP\nNM02V7bavtMSu1it/BYjM5PQhCuJjD81pV2srn9vQ9I9a4flZ3LLhTLhri+QTVSE6OMMwyBYXkaG\n3992dh8M4n58ARn/fAGLrhM55ncEr56MmZuX8mu1tmetbDo+MEimL0QPC9fXJ/ru21kZ61j0EZ57\n7sK2cSP60KEEp1xH7KBDOv56cZ2llQb/1uoprw7KpuN9kOyRK0QfZBgGwbJS3MEAzjbqzVs3bsBz\nz104F32dkzbHAAAgAElEQVSMabcTPutcQudeAB3se9d1nYDThX1IAS4pndCnSfeOEH1MqK4Os2ID\nPosFWgv48TiuF/+B57H5WEIhYvvsS2DaTIztd+jQa5mmiR9gSAHeQelvZi76Bwn6QmxDqWb3tm9X\n4J19O/bvNQyfj8CkqUTHjktpoLapYFwn5huEt6BASh4LQIK+ENtMqK4Oc2M5Pqu11eze4q/HveBh\nXC+/iMU0iYwdR/CKqzE7mKFH4joRjwfX8CJcUvJYNCFBX4hupus6ofJSPMEADlsrP3KmifM/7+O5\n/26sVVXow7cjMH0m8f0O6NBrJUon2LENLZLSCSIpCfqiW3XFFoB9Wai2Fio3NvTdJ/9xs5aux3P3\nXTiXLsF0OgleMoHwmedAB4ubSekEkQoJ+qLbdMUWgH3VluzeGwxib63vPhYj4/nncD/5OJZohNiB\nowhMvQ6jZHiLU9sqkBaK61I6QaRMgr7oNp3ZArAvC9ZuxlKxsc2+e/tXX+C56w7sP/2IMTiHwMwb\niR57XNJFWa0VSNMNgz13zsE5rBCvy9Vt70f0LxL0RbdJdwvAvkrXdUJl6/GGQq1m95baWtwP3U/G\nG68BED75VEKXX4mZnZ30fGi9QNpbP4U56ChZRSs6RoK+6DYd3QKwLwtursFSWYHPauWT0kjLrpgS\nN85/vY3nwXux1tQQ33EngtOvJ77nXu3eu7UCaRtqwl39NsQAIEFfdJuObAHYV+nxeCK7D4ex22xJ\nu2Je+9fX/Gb5E2Su+BzT5SI48WrCp/8J7KlNpSzOdiQtkNYff3mK7idBX3SbVLYA7MuCm2uwVlTg\ns/3ad9+0K8YRj3Lapwv5w/KFOPQ40UMOJThlOkZRccqvEdd1jh45mKf/W9HiuZ745TnQZ2P1B1J7\nR4gOisdihMtLG7P7ps5/aR2GCXut/ZoJ/1nAsJoyqr05PH7kRVw8449t7nzVlGma+C0WyMvH4xvU\nEGx79pfn1rOxtrj0pJES+Lcxqb0jxDYSqNmErbKyWXbflLIGOP7NRxmz6mN0i5XX9z2Rvx5yBnlD\nfCkH/ICuEx+Ugzc/v7F0Qlv7524rA3U2Vn/TZtBXStmBJ4ERgBO4DVgPvAF813DafE3TXlRKzQJO\nAGLAJE3TlndXo4XY1uKxGOGyUjKjEWzJti40DFxvvsatD9+PM1DP9wU78tDRE/ixYCcATty19dk5\nW0TiOhGvF1dBIRm9sHTCQJuN1V+1l+mfBVRpmnaOUioH+AL4M3C3pmn3bjlJKbUvcJimaaOUUiXA\nQkC24BH9QmBTNbaqSnw2W9KCZ7affsQz+zYcX3+F6fHy7XlX8/B2R1HqNyjJdnBik4VUyei6TsDh\nxF4yFG8vLnk8kGZj9WftBf1/Ai82fG0hkcXvD+yqlBpPItufBBwKvAegado6pZRNKZWraVrL7XmE\n6CMas/tIGFuyeffhEO4nHyfj+eew6DrRMUcRuGYKRUOG8JcU7m+aJn7TxMwfgndwTpe3v6sNhNlY\nA0GbQV/TtCCAUiqLRPC/AXABj2ua9oVSagYwC6gBmgZ4P+Db6pgQfUaz7D5JwHcsXYJn7mxsZaXo\nhUUEp0wn9tvDUr5/KK4T8/lwDynoM6UT+vtsrIGi3YHchu6al4F5mqa9oJTyaZpW2/D0q8CDDf82\n7bTMAjZ3dWOF6A5NpyEW5bo5emcvo4scSbN7S2Ulnvvm4vrg3+hWKwsPPIVFx57F70qGcFAKrxXT\n44QyPDhLivB2sKBab9AbBpRF57Q3kFsAvAtM1DTtw4bD7yqlrtA07TPgKOAzYAkwRyk1FygBLJqm\nberGdgvRJbaehlhaFeSZqiCeUbnN++F1HdfLL+F+5CGsgQCrihQPHz2BNfkjIETjgqzW+u4TJY9t\n2IpLpOSx6FHtZfozgEHAjUqpmwCTRB/+/UqpCLABuETTNL9SahGwlETf/8RubLMQXaa1aYhvrq5r\nDOA2bXViF6tV32JkZfG3cVfwj52OxLRYW72mKb+uY+Tm4cnJld2rRI+TxVliQLto9gcYSb4LbRZ4\n8vgcPI8twPXiC1gMg8ixxxG8ejLnfRBs/ZrTShofh+M60awsMgoKkw8EC5EmWZwlWpDl8m2LRSJE\nykopzrKzvi7e4vnjSz/Dd8Zj2Co2og8rITD1OuK/SfTaF2fHktbCKc5OzK2P6zpBVwaOoSV4MzK6\n940I0UES9Puhgbx5SSr8lZU4NlXjs9sYt5uvWYG0/LpKLvnwMQ768VNMu53QeRcSOvcCaBK8x+2a\n3eyaLcbumkWdaWItKMLr822T9yJER0nQ74dkuXxy0XCYWHkZWfEYVnuiu2VLH/zb39aw/8cv86el\nL5ARDRPbd38C02ZgjNi+xX22XPNmk/LJx+ySye4ji/E0KZ0gRG8kQb8fkuXyLfkrK3Fuqibb3nJV\n7W/rf+TYv9+O/fvvMHw+/NOmEz1hXJu1cg4a7uGg4R4iukHE603029vlx0n0fvJd2g/JcvlfJcvu\nt7DU1+NeMA/XKwuxmCaRsScRvOIqzEGD271vY+mEot5dOkGIrUnQ74dkuXyixEGgshJnzaaW2b1p\n4vz3e3juvxtrdTX6iO0JTJtJfN/9UrqvH2BIAd4UfjkI0dtI0O+HBvpy+Wg4TKyslCw93iK7t65f\nh2fubJzLlmI6XQQvvZzwmedAClUtg3GdmG8Q3oIC6bcXfZbM0xf9RtPs3r1VsCcWI+Nvz+J+6gks\n0QixUQcRmHIdxrCS5DdrIqrrhN0eXIVF2HthyWMx8Mg8fTHgRUMhYuVlZBs6lq0Cvv3LL/DOvg3b\nmp8xcnIJXDOL6NHHtrupSaJ0gh1bcZGUThD9hgR90adsvejshIO2Y49ccNXVkm2zNgvkltrNeOY9\ngOvN1zAtFsKnnkbosisws7LafR2/YWDk5uHNye3GdyPEtifdO6LPaG2P1stG5XDw8CYzk0wT59tv\n4pl3H9bNm4nvtDOB6dej77Fnu68RiuvEsrNxFxT2mZLHYuCR7h0xILS66Gx1fWPQt675Ge+cO3B8\n/j/MjAyCV1xN+P/+BPa2++Ljuk4ww41jaIGUThD9mgR90We0tuisrC4GkQjuZ58i47mnscRiRA89\nnODkaRhFRW3e0zRN6i0WrIXFeLPb38dWiL5Ogr7oE0zTpHCQk7KaSIvnxlSuwHfWI9jWr8PIH4J/\n8lRio8e0O1Ab0A3ig3Pw5uXJFEwxYEjQF73O1oO1x+5fxL5ZEU7axcuCZb8G/UGBzVzw8VOMWf0x\nptVK+P/OIHjxZeBte+VxOK4Tzcwko7CIDCl5LAYYGcgVvUprg7UTGnay+mRtkLdWbmbkf9/mvMXP\n4gkHiO+2O4HpM9HVbm3eWzeMROmEIQW4pHSC6MNkIFf0G+3tZPXbaCnHvHwHjm++wvR4CUyeRuTU\n05JuXr5Fs9IJvkHd0Wwh+oz29si1A08CIwAncBuwEngaMIAVmqZNbDj3JmAsEAMmaZq2vNtaLfqt\n1gZrq6vrcT/0TzL+/lcsuk7kyKMJXjMFMz+/zfsF4zqxQYPxDhnSb/vtZcMc0RHtZfpnAVWapp2j\nlBoMfNnw30xN0xYrpeYrpU4G1gKHa5o2SilVAiwEftOtLRf9TtjvpzDbTtlWu1Id8NNnXPHRo7g3\nV6AXFROcMp3YIYe2ea9IXCfi8eAaXoTL4ei3gVE2zBEd1V7Q/yfwYsPXViAO7Kdp2uKGY+8AxwIa\n8B6ApmnrlFI2pVSupmkttxcSYiuGYRDaUI6rvp6Tm+xKlVNfzcUfPcGh3/8Xw2YjdPa5hC64GDLc\nbd7Lb7NjG1aMt2FAtz8HRtkwR3RUm0Ff07QggFIqi0Twvx6Y2+SUesAHZAFNA7y/4bgEfdGmsN+P\nsbGcLNPEYrcldqUydPzP/p0T//MsnmiITbvsge3GG9B32rnNe/lNM2nphP4cGGXDHNFR7Q7kNnTX\nvAzM0zTtBaXUXU2ezgJqgDoge6vjm7uyoQNVf+2WMAyDYHkZGX4/LrutcU69TVvFsbNvx75qJUZW\nNoHJk2HcePQ2SiKE4joxnw/3kIKkpRP6c2CUDXNER7VZXEQpVQC8C0zTNO2ZhsNfKKUOb/j6eGAx\n8F/gWKWURSk1HLBomrapuxo9UGzpllhfGcAwzcZuiWUrN/Z00zolXF9P+Kcf8IVDiYAPEAjguXcu\n2Recg33VSiLHnUDtCwuJnHxqi+0Nt4jpceocTizb74C3sKjVWjnFecmnZ/aHwDj24BGtHB84G+aI\njmkv058BDAJubJidYwJXAw8qpRzAKuAlTdNMpdRiYClgASZ2Y5sHjP7WLWEYBsGyUtzBAM4tUyxN\nE8fHH+K9Zw7Wygr0kuEEplxH/Dej2ryP32bDWjg0pdIJ/XknsYG+YY7oOFmc1YtdNPtDjCT/f2xW\nC49NG9MDLUpfqK4Os2IDWU2mTVrLy/HcPRvnksWYDgfhs88jdM754HK1eh+/rqPn5OHNze3QFMxE\nN5kERtE/yOKsfmpb99emO37Q1nVbsntP0I/D1vDtFo+R8cLfcT/xCJZwmNh++xOYNhNjuxGtvkY4\nrhPNyiKjoBBbGqUTRu1eIEG+F+qvY1a9mQT9XmxbdkukO62xrev2GuaGig34LBZoCPj2b77GM/s2\n7D/+QCTLx9+Pv5xXtj+MYs3JOEswMXunibiuE3Rl4BhaIiWP+5n+PJW2N5Og34tty/7adMcPWrvu\njUXfc+CRuY3ZvaW+Hvf8B3G9+jIW02TdmLFM3/X/qHcn+uTX18Ya5+cfNNyTKJ1gsUBBoZRO6Kf6\n25hVXyFBv5fbVt0S6U5rbO26jbXRRMA3TZzvv4vn/nuwbqomvv0OBKfN5PaKIuq3WnkLiRo7ew51\nER+Ugzc/v9+WThD9eyptbyZBXwDpjx+0dl1xtgPrunV4774Tx7JPMJ0ugpdNJHzG2eBwUPbSuqT3\nK62LYd9+RzIcbe90Jfo+WWPQM2QTUAGkP9872XX2eIyrv30F31l/xLHsE6IHHULt8/8kfO4F0BDM\ni7OTB/XivEzsEvAHBFlj0DMk0xdAeuMHuq4zMivMJQcO5p3v/JTVxRi9aTUXvzufzLJfMHJzCVwz\nhehRx7TYxWpckxo7TckP/MAhawx6hszTF2kJbq7BUllBZsMqWMvmGjzz7sf11huYFguRU04jNGEi\nZmZW8uvjOsuq4d+raimvDsoPvBAd0Jl5+hL0RYfo8TihsvV4w2HsNltioPatN/DMuw9rbS3xnRWB\n6TPQR+6Z9PqYHieU4cFZWITD6dzGrReif5DFWWKbCG6uwVpRgc9mBZsN65qf8c6+HceXn2O63QSv\nmkT4D6eDveW31ZbSCbbiEryZmT3Q+t5BFiOJniZBX7QrHosRLi/9NbsPh3E/8yQZf30GSzxO9PDR\nBCdNxSgsSnq93zAwcnLx5uZt45b3LrIYSfQGEvRFmwI1m7BVVjZm9/ZPP8F71x3YStejFxQQnDyN\n2OFHJL02FNeJZWfjLihstQLmQCKLkURvIEFfJBWPxQiXlZIZjWCzWbFUV+G5/x5c77+LabMR+tOZ\nhC66DDwtyxbHdZ1ghhvH0AIpndCELEYSvYEEfdFCYFM1tqpKfA2FzVwvv4R7/oNY/X7iu48kMP16\n9F1Ui+tM06TeYsFaWJxSyeOBRhYjid5Agr5o1JjdR8LYbDZsP3yPd/Zt2Fd8g+H1Erh2OpFTfg9J\nqlwGdJ344Fy8eXlSOqEV/bmuv+g7JOgLYKvsPhrF/fgjZPzjeSy6TuSoYwhecy1mXn6L68JxnWhm\nJhkFhWQkmbUjfiWLkURvIPP0B7h4LEa4dH1D370Nx+JFeO6ZjW3DBvTioQSnTCd28G/5ZG2QN1bX\nUVYXozjbwdhdMtljx0HYhxTgStKvL4ToPt0+T18pNQq4U9O0MUqpfYE3gO8anp6vadqLSqlZwAlA\nDJikadrydBslto1AdRX26mp8DQO13nvm4Pz4w8RA7TnnEzr/Qshw88naYLOSCetrYzyyvIZLi4oZ\n1c8DvsyrF/1Nu0FfKTUVOBvwNxzaD7hb07R7m5yzL3CYpmmjlFIlwELgN93Q3l6trQDRm4JHLBol\nUlZKViyK1TRwvfACnsfmYwkGie21N8Hp16PvsGPj+W+srkt6n/4+1VDm1Yv+KJVM/wfgFOC5hsf7\nA7sopcaTyPYnAYcC7wFomrZOKWVTSuVqmtayolY/1VaAAHpN8PBXVeGorsJnt2HTVicGarXVGFnZ\nBGbcSPTEk2CrOfVldS3r3kP/n2oo8+pFf9TuihlN014B4k0OLQOmapo2GvgJmAVkAbVNzvEDvi5s\nZ6/XVoBo67ltJRaJ4P/5J7JqqvFEQnjuuYvsi87Frq0mctxYal9YSPSk8c0CvmEY1FmsFOa4k96z\nv081lHn1oj9KZ7rFq5qmbQnwrwIPNvzbdGJ2FrC5k23rU9oKEK2NlW+r4OGvrMSxKdF37/joA7z3\nzsVaVYk+fDsCU68jfkDLnji/YWDk5uHNyWVcOHNATjWUefWiP0pnbfy7SqkDGr4+CvgMWAL8Till\nUUoNByyapm3qqkb2BcV5yQc0i3K9bT7XnaLhMP6ffiRr8yYyKzeSOeVqsq6fjqV2M8GLLqX2uRda\nBPxQXKfO48W14854c3KBRBfUpSeNZFh+JjarhWH5mVx60sh+38Uhm3yI/iidTH8CME8pFQE2AJdo\nmuZXSi0ClgIWYGIXtrFPaG/hzbbOlP2VlThrNuEzdTL+/jfcTzyKJRIhtv+BBKbNwBje/LW3lE5w\nDivE63K1uN+22qu3N5F59aI/knn6XSgxQyd5gGjrua4UDYeJlZeRGY/hXPENntm3Yf/pR4zBgwle\nNZno745vtotVY+mE/ALcfbx0Qm+aISVEd5JNVASmaRJoyO49AT/u+Q+S8dorAIRPPoXQhCsxfc3H\n1gO6QXxwTr8onbD17KktBkI3lBh4ZBOVAS4aDhMrKyVbj+P6z7t47r8Ha00N8R12JDhtJvG992l2\nfjiuE83KSpROSFJHpy+S6ZVCpEaCfh9mmiaBigqcm2sYVF6Kd86dOJYvw3S5CF5+JeHTzwSHo/F8\n3TAIOF04hpb0u5LHMr1SiNRI0O+joqEQsfIyssMh3H99BvezT2GJRoke8luC107HKB7aeK5pmvgt\nFhhSgNc3qAdb3X1keqUQqZGg38eYpklg40ZcdbXkfPU53tm3Y1v7C0ZeHv5rphA78uhmA7UBXSc+\nKAdvfn6f77dvi5QtFiI1EvT7kEgwSHxDGb7qarzz7sP1zluYFgvh0/6P0KUTMDOzfj03rhPxenEV\nFJLRpIunv5LplUKkRmbv9AGN2X3NJrL/9Raehx7AWldLfBeV2MVq95GN5+q6TsDhxF5QKCWPhein\nZPZOP7Ylux/8w/d4596J48svMD0eAldfS+S0P0LDxiWmaeI3Tcz8IXgH5/Rwq4UQvZUE/V7KNE2C\nGzfgrNhIwXNPk/HXZ7DoOtHDjyA4eSpGQWHjuaG4Tsznwz2kAKs1ncoaQoiBQoJ+LxQJBomXl5Kz\ndAneubOxlZWiFxQQvHY6scNGN54X0+OEMjw4S4rwOp092GIhRF8hQb8XMU2T4IZyMn76ifyH7sP1\n7/cSu1j96SxCF10KDX30hmHgt9mwFZfgzczs4VYLIfoSCfq9RNjvRy9bT94rC/E88hBWv5/4yD0S\nA7U779J4nl/X0XPy8Obm9uspmEKI7iFBv4dtye69n3/G4Dl3YF/5LUZmJoGp1xEZ//vGTU2alk6w\n9ZPSCUKIbU+Cfg8K+/0YP/1AwROPkvHPv2PRdSLH/I7g1ZMxc/OAhpLHrox+WTpBCLHtSdDvAYZh\nECwvw/fOW/jun4tt40b0oUMJTJ1BfNTBwK+lEywFRXh9yXeelFLCQoiOkqC/jYX9fixffc7Qe+fi\nXPQRpt1O6LwLCZ17ATRk8qmUTmhrI3YJ/EKI1kjQ30YMwyC4bi25zzxJ1hOPYAmFiO2zL4FpMzG2\n3wH4tXRCRmERGfa2/9dIKWEhRDpSCvpKqVHAnZqmjVFK7Qg8DRjACk3TJjaccxMwFogBkzRNW949\nTe57wvX12D98n5K77sD+nYaR7SMwaSrRsePAav21dELJULwplk6QUsJCiHS0G/SVUlOBswF/w6F7\ngJmapi1WSs1XSp0MrAUO1zRtlFKqBFgI/Cb5HQcOwzAIaavJv/9uPK+8hMU0iZxwIsErrsEcPLix\ndAJDCvAOGtyhe0spYSFEOlLJ9H8ATgGea3i8v6Zpixu+fgc4FtCA9wA0TVunlLIppXI1Tavu6gb3\nFaHaWtz/+BvD75uLtaoKffh2BKbPJL7fAQAE4zox3yC8BQV8uqqCtxYu69CArJQSFkKko92gr2na\nK0qpppGk6chiPeADsoCmAd7fcHzABX3DMIh++gkFt/8Z1ydLMZ1OghdfRvisc8HpJKrrhN0eXMOL\ncDkcaQ/ISilhIUQ60hnINZp8nQXUAHVA9lbHN3eiXX1SqLKS7AfvZciTj2GJRogdOIrA1OswSoYn\nSidYrNiKi5qVTujMgOyo3QskyAshOiSdoP+5UupwTdMWAccDHwA/ArOVUnOBEsCiadqmLmxnr6br\nOsbbb1B8683Yf/4JY3AOgZk3Ej32OLBY8BsGRm4e3pzcFtfKgKwQYltKJ+hPAR5TSjmAVcBLmqaZ\nSqnFwFIS3T8Tu7CNPSLVhU/hX9Yw+JYbyXzjtcTj8b8nNOEKzOzsRMnj7GzcBYWtljyWAVkhxLYk\nO2clsXU/+xaXnjSyMfDr8Ti2x+aTe+8crJs3E99xJ4LTrye+516J0gkZbhxDCnC2UzohldcSQoim\nZOesLtZeP/vnr37IdnfcyK4/f03U4eKnMy5jyITzMG126gFrYTHe7Oyk99iaDMgKIbYlCfpJtNbP\nXrWxhsorruHIhc/i0OMs3/4AFhx5MRW+Ai5YH2GfvfLx5uV1uOSxDMgKIbYVCfpJJOtn32vt11z5\n4SMUVpdS7c3h0TEX8d+dD4aGAP/emiiHHpXfE80VQoiUSdBPounCJ19wMxd+/BRjVn2MabHyxr4n\n8twhZxByNS+XUF6d/K8DIYToTSToJzFq9wIwDGoeeJjx/3qCrLCf2h0U3HADb6/LIVQba3GNzLYR\nQvQFEvSTWfENR0+6goyvvsD0eAlMmoL++z8SNOGoHAvPLCptcYmUPxBC9AW9Luj36MYgwSCO2/+M\n74lHseg60TFHEbhmCuGcXCJuD67CIkYrBxmDBslsGyFEn9Sr5un35Jx127vvkDV9Mo6yUvTCQoLX\nTidyyKH4bXZsBYVkeKX7RgjRO/SbefrdsTFIe385WDeUkzFtEt5/vY1psxE68xxCF16CPyMDIyc3\naekEIYToq3pV0O/qOjRtVrBUeTifeISsO27FGvAT32NPAtOvxz9iB2I+H+4hBa2WThBCiL6qVwX9\nrq5D09pfDl8t/DfHfDQf1zdfY2RlEZg2E/+J4wh5MnEWFuF1OtN6PSGE6O16VdDv6o1Btv7LwR0N\nceZ/n+fEL97CZhpEjj0O/5XXUJ8/BGvekJRLJwghRF/Vq4J+V9ehafqXw0E/fMIlHzxGvr+aipxi\nMmZdz+b9D0TPycObm9vh0glCCNEX9arZO11t2cqNvPzXj7jkw8c46MdPiVntLDzwVLyXX8See5eQ\nUVCIzWbrziYIIUSX6zezd9KVdIbOLrkc8dE/Oe6vt2IPh1gxbCQvjb+CPY/Ym/1+s1O7JY+FEKI/\n6tFM/+Qpr5udXYCVbG7/LuXfcdOnT+D7UcPw+QhecQ3VY8dB/hA8vkFd0XQhhOgxfTbTN0wz5Y3A\nW9N0ho437OecJX/luK/exYpJZOxJVE24guiIHfDm50u/vRBiwEs76CulPufXzc9/Bh4F7gdiwPua\npt3SkfuluwCrrCoIpsmh3y3h4o+eICdQw9qcYcw/ZgKXTh+Pq6CQTIejw/cVQoj+KK2gr5RyAaam\naUc2OfYFcIqmaWuUUm8ppfbRNO3LVO+Z7gKsvay1jH/pXvb/5QsiNifP/fZMXj5gPAUFPrzDStK6\npxBC9FfpZvp7A16l1LuADfgz4NQ0bU3D8+8CRwEpB/0OL8CKRvE8/AA33z8bWzTC59vtw/yjLmXD\noCIATjxk+47dTwghBoB0g34QmKNp2hNKqZ2Bd4CaJs/XAx2Kuh1ZgOX45L9kTr0Gu7YaPSeHT8+c\nzhP5B1FZF2OYVL0UQohWpRv0vwN+ANA07XulVC2Q0+T5LH7t72+VzWrp0AIsy6ZqvLfchPv55zAt\nFgLjT6X66msZtrPiVimdkJIeLV0thOhx6Qb9C4A9gYlKqWLAAwSUUtsDa4DfATe3d5PHpo1J7dVM\nE9c/nifzzzdgra4muuNOVM24Ef3IY3BnZqb5FgaeNgvQSeAXYkBIN+g/ATyllFoMGMD5Df8+D1iB\n9zRNW94VDbR9/x2Z0ybhXLIYIyODmsuvovayiXgLi5A5OR3THaWrhRB9S1pBX9O0GHBWkqcO7lxz\nmgiH8dx/N54H78USjRI45FA23fhnnPvuj1dKHqelq0tXCyH6nl5ZhsHx8YdkTpuE/eefiOcPoXrq\nDOL/dwYZbneXvs5A69/u6tLVQoi+p1elzJaKCrImXMSgP5yM7Zc1bPrD6ZS/vwjreRfi7IaA/8jr\n37K+MtBsZfCylRu79HV6k7EHj2jluGzqLsRA0TsyfcMg46/P4P3LLKy1mwntuhubbrkDx+gxZHRT\n6YSB2L/d1aWrhRB9T48HfdvKb8maeg2O5cvQvV42Tp2BedVknC5Xt77uQO3fHrV7gQR5IQawHg36\n3r/Mwj3/QSzxOLVHHo3/9jk4d9hxm7y29G8LIQaiHu3T9zx4L7EhBZTOf4LoCy9vs4AP0r8thBiY\nejTTr55wJfGpM3D2wAIr6d8WQgxE/Xq7RCGE6I86s4lKr5qyKYQQontJ0BdCiAFEgr4QQgwgEvSF\nEPtDRUoAAAQWSURBVGIAkaAvhBADiAR9IYQYQCToCyHEACJBXwghBpAuXZGrlLIADwN7A2HgIk3T\nfurK1xBCCJG+ri7DMB5waZp2iFJqFHBPw7FOG2gbngghRHfo6u6dQ4F/AWiatgw4oCtuOhA3PBFC\niO7Q1UE/G6ht8jiulOr0a7S14YkQQojUdXXQrwOymt5f0zSjszcdqBueCCFEV+vqPv0lwInAS0qp\ng4Bv2jo51Upxhml+Dey59XHdML/Oz8/aO52GCiHEQNSlpZWbzN7Zq+HQ+ZqmfddlLyCEEKJTerSe\nvhBCiG1LFmcJIcQAIkFfCCEGEAn6QggxgEjQF0KIAaSrp2y2S+rzgFLKDjwJjACcwG3ASuBpwABW\naJo2safat60ppYYAnwFHAzoD9HMAUEpdB5wEOEj8nCxigH0eDT8fz5D4+YgDFzMAvy8aStncqWna\nGKXUjiR5/0qpm4CxQAyYpGna8vbu2xOZfmN9HmAGifo8A81ZQJWmaf+/nXsH0aOK4gD+28RHIauC\nBkEwYHVqTSEi7poiamxEsLEQTJlKBBUfpBKsLEQjGiJCBEFUiIgSo6jgGgVBYyHoEQtBsBCMAaNY\nCLG487GTzcJuszvs3vOr5sE33Pkz93C/edwF7MdhLYenMnMROyLi3ikbuFmGDv4KZl/gdZkDRMQi\nbh36xh3Yrc887sHOzLwNz+BZneUQEY/hKC4fNl10/hFxExYy8xY8gJfWc+wpiv6GzM+zxbyFQ8Py\nDm00c3NmLg3bTmij3h48h5fxG+b0mwPche8j4l28h/f1mcdPuGS4K3CVNortLYefcd9ofc+K89+n\n1dKPIDN/xc6IuGatA09R9Ddkfp6tJDP/ycy/I2Ieb+NpreDN/KVd7NtaRDyE3zPzY8vnP74Wushh\n5Frswf04iDf0mcc53IgfcQQv6Kx/ZOZxbTA4s9r5z7uwlp6zjlymKLYbMj/PVhMRN+BTHMvMN7V7\ndTPzODtJwzbXAeyLiM+0ZzyvY9dofy85zPyBk5n53/Al+78u7MS95PEIPszMsHxdXDba30sOYyvr\nw59aLb1yxfY1c5mi6J/S7tlZz/w821FEXIeTeDwzjw2bT0fEwrC8H0ur/ngbyczFzNybmXvxHR7E\nid5yGPkCd0NEXI8r8Mlwr59+8jhjeQR7Vnvh5HSHOYx9u0q/+BJ3RsRcROzGXGaeWetAm/72Do5r\no7tTw/qBCdowtSdxNQ4NT9/P42G8GBGX4ge8M2H7pvQojvaYQ2Z+EBG3R8TX2t/5g/gFr3aWx/N4\nLSI+195iegLf6C+HsYv6RWaej4glfKVdL+t6o6nm3imllI509QC1lFJ6V0W/lFI6UkW/lFI6UkW/\nlFI6UkW/lFI6UkW/lFI6UkW/lFI6UkW/lFI68j8HY3NwYKuxcgAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fit.update(fix=30)\n", "fit.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "That works for any other point as well. E.g. we also know, that the line goes through ``y = 4 * 10 + 30 = 70``:" ] }, { "cell_type": "code", "execution_count": 8, "metadata": { "nbsphinx-thumbnail": { "tooltip": "Make a linear regression to your data." }, "scrolled": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX0AAAECCAYAAAASDQdFAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xd8FHX++PHX9pYCJCEQiohlFPSsJ6IIWE8PFNudnh4q\noPQinOjZy1e986eHBRVFBRT17AXFgp6ecohYsCEwitJTSC/bd3Z+f2wSE5JNNpsN2WTfz8fjHpfM\n7sx8dmTf+cxnPp/326DrOkIIIVKDsbMbIIQQYt+RoC+EEClEgr4QQqQQCfpCCJFCJOgLIUQKkaAv\nhBApxBzLmxRFWQ9U1P66FVgMPAgEgQ9UVb1DURQD8ChwBOADrlRV9dfEN1kIIUS8Wg36iqLYAF1V\n1VMabPsGOE9V1W2KoqxUFOVIYH/ApqrqCYqiDAMWAOd2VMOFEEK0XSw9/SMAl6Io7wMm4HbAqqrq\nttrX3wdOA/oC7wGoqrpOUZRjE99cIYQQ7RHLmL4HuFdV1T8A04CltdvqVAOZQDpQ2WB7SFEUeWYg\nhBBJJJag/BPwHICqqj8TCey9GryeDpQDVbU/1x9bVdVwgtophBAiAWIZ3pkIHA7MUBQlD3ACbkVR\n9ge2AX8AbgMGAGOBVxRFOR74obUDFxdXS+IfIUSXdMtT69hV7G6yvX9OGndMOq5Dz52Tk26Id99Y\ngv5TwFJFUVYDYWBC7f8/T+ROYZWqql8qivIVcLqiKGtq95sQb6OEECLZ5Zd4mt1eUNr0D0EyaTXo\nq6oaBP7azEvD93qfTmTMXwghur28bGezPf2+WS4A1m0sYuXabeSXeMjLdjJm+CCGDcndx61sSh60\nCiFEHMYMHxRl+36s21jE4yt+ZFexm7Cus6vYzeMrfmTdxqJ928hmSNAXQog4DBuSy5RzhtI/Jw2T\n0UD/nDSmnDOUYUNyWbl2W7P7rFy7fZ+2sTkxrcgVQgjR1LAhuc0O2STzeL/09IUQIsHysp3Nbq8b\n7+9MEvSFECLBWhrv72wyvCOEEAlWN+Szcu12Ckrd9M1yMWb4fkkxe8fQmYXRZXGWEELETtM0vAW7\nmXr0YRnv63p1PMeQnr5IGps3b8Tr9bJx4wYuvfTyzm6OEEnFXVaKqaSEDKOBNIh7Ra6M6YuksXnz\nJoYOPZzKygq8Xm9nN0eIpOD3eHBv/ZW00hJcpvaHbOnpd4BQKMQ//nEHhYUFBINBLrtsItXVVbzz\nzlsYDAb8fj9btvzEihXv43KlNdrvrrtuo7AwH5PJzLXX3sjAgfE9+CkvL2PSpPE88MCjjY7RXNtG\njBgZdXudPXuKWLNmNf37D+DAAw+mZ8+ebW7TN998zS23XM/++w8GwO12069ff2655f8wm82ce+4F\nhMNhwmEdh8MR9Ti6rvOvf/2TLVt+xmq1ct11N9GvX//614PBIHfffTv5+btxudKYN+9a+vcfQDgc\n5p577mTHju0YjUauuSbSlokTLyUtLZIrsG/fPK6//pY2fzYhEi0cDuMtLMBaXU2G2QTGxPTRJeh3\ngFWr3qVHjx7cfPMdVFVVMmHCpbz66tucddZYABYsuIezzx7XKOADfP75GsJhjUWLlvDll+tYvPgR\n7rzz/7X5/KFQiHvv/Qd2uz2mto0YMTLq9jrbt29lxYrXyMjIZPr0OXEFfYBjjvk9t912V/3vt99+\nE2vWfMqoUZEaPR9//B/Gj7+CUCiE2dz8P89PP/0vgUCAxx5bwo8/buDhh+/nH//4V/3rK1a8jtPp\n5PHHl7Jjx3YWLPh/LFiwkDVrPsVgMLBo0VN8883XLF78KLfffjcADz30WFyfR4iO4C4vw1hSTIbR\nCGZTQo8tQb8Zr732Mj/88B233nond911G0OHHsa5514IwH//+x9effUlDIbfhtSmT5/NIYcMqf/9\nlFNO5+STTwMivdKGwWvz5o1s27aVefOua3LeAQP2Q9M0dF3H7a7BbLY0274lSxbTq1cW5557QbOv\nP/LIg5x33gUsX76syWvR2tZSmwEOPvgQTj/9LCorKzjggANjvl57azhxIBgMUlpaQnp6BgCrVr3H\n+vVf8tVXXzB//vVRr/X333/LsGEnADB06GFs3ryp0Tm2bdvK8cdHXh84cD+2b98KwEknjebEEyN/\nyAoLC0hPT2fLlp/w+XzMmzcTTQszefJ0hg49rNm2C9HRAj4fwcICXAE/JlNig32dpA76rttuwvbW\nGwk9pv/sc3HfdmeL7zn//D/x1VdfcPfdtxMKhRoFsNGjT2X06FNb3L+uh+3xuLn55r8zefL0+teW\nL1/KhAlXNbufw+EgPz+fSy65gKqqSu6554FGr//nPx/wxhuvUFhYiMVi5qOPPuCyyyZy7LG/pXF9\n55236NmzJ7///fE888zSmNvWUpsBMjN7cMkl45ttd0vXa2/r13/F7NlTKSsrw2g0MG7c+Rx9dKTI\n2hlnnMkZZ5xZ/95o1/rNN18jLe23uySTyUQ4HMZYe/t70EEH89ln/+Okk0azYcMPlJQUo+s6BoMB\no9HIXXfdxurV/+X//u8e7HY7l1wynrFjz2Xnzh1cc81s/v3v1+qPJcS+EA6H8RYVYq2qigzldFDA\nhyQP+p3p0ksvZ9q0iTz11PJG2+t6n3UMBkOTnj5AUVEhN954LRdc8GdOPfUMAGpqatixYztHHXVM\ns+d88cXnGTZsOFOmzKC4eA+zZk1l+fIXsVgiPf5TTz2dU089nSVLFpOVlc24cec3OUbdc4Mvv1zH\nzz//xJ133so99yygZ8/f6t4017aWtrfneu2tbninqqqSuXNn0rdvv6jvjXatnU4XHs9vy9kbBnyA\nMWPOYfv2rcyYcRWHH34EinJoo7uFG2+8jfLyMq666nKeeeYF+vUbAMCAAQPJyMiktLSEnJzebfr8\nQsTLU1kBe4o6ZCinOUkd9N233dlqr7wjBINBHnroX8yffwP33vsPHn30yfrhjlh6+mVlpfztb7OY\nN++6+l4swLffrufYY4dF3S8jI6P+PGlp6WiaRjisAc0P8zTn4YcX1/88a9YU5s+/oVHAj9a2aNtj\n0dL1iiYjI5Obb76D2bOnsmzZ8/TqldXkPdGudVFRIWvWrObkk09jw4Yfmgw3bdq0kWOOOY5Zs+ax\nefMmCgsLAHj//XfYs2cP48dfgdVqxWg0smrVe/zyyxb+9rfrKCkpxuv1kJWV3abPL0Q8gn4//oJ8\n0oIBTPvwzlLuYZvx2GMLGTFiJGeffS7Dh5/IY4893Kb9ly9fRnV1NcuWPcmsWVOYPXsqgUCAHTu2\nk5fXuGdbVVXFTTddC8Cf/3wJqrqZGTOu4uqrpzN16gxstqYPYydOnNxsL39vdb3bhueI1rZo22PR\n3PUqKyvl1ltvaHG/QYP2509/upgHHrgvpvPUGTnyZKxWK9OmTeSRR+5n1qx5jT7jgAEDeOml55k6\ndSJPPfUYs2bNBWDUqFP4+WeVmTMnc801s5kz5xrGjh1HTU0106dfyW233cjf/36LDO2IDqXrOu7C\nAvRtW8nUQvs04IOsyBUdRNM0HnvsYWbMmNPZTREiaXgrK9GLi9iw08Pbm6vJrwqSl2Hh7EMyOH5g\n80na9qbrOpNPODbzVV2viqcNST28I7ouXdejPvgVItUEAwEChQU4fV6+2u3nsXVl9a/tqgyyaF0p\nQMyBvz0k6IsOYTabGz1LEN1fspYH7Ey6ruMuLsZSXlY/K+etzc130N/eXNVy0HfX4HhmGbbXX2lX\nmyToCyHara48YJ268oBAygZ+b1UVeskeMsJhDA1m5eRXBZt9f7TthELY3l6BY/EijOVlaL1zoTqu\nXGuABH0hRAK0VB4wmYL+vrgbCQWD+ArzcXo9WExmMDTOjZaXYWFXZdMAn5fRdJae+YvPcT64APOv\nv6A7HHgmT8N78aVwyoi42ydBXwjRbslcHrBOR9+N6LqOu6QEc3kpmSYTmJoPr2cfklE/ht/Q2EMy\n6n82btuKc+H9WD9bg24w4B87Ds+UaejZOdDOyTcS9IUQ7ZaX7WRXcdMA31HlAePpsXfk3YivpgZt\nT2FkKKeV1bR14/Zvb66qn70ztnb2jqG8HMdTi7G98SoGTSN4zLF4Zs9DO1hpV/sakqAvhGi3McMH\nNepF/7Y98eUB4+2xd8TdiBYK4S3Mx+724DKbmgzlRHP8QGfjh7aBAPbnnsG+7CmMNTVoA/fDM/Nq\ngiNOivmYsZKgL4Rot31ZHjDeHnui70bcpSWYSkvJNLUjfYKuY/n4PzgfeQhT/m7CGZm4583Hf94F\nECXhYnt1+aDfVaaJBQIBLr30Ql5+eUXU96xY8TpjxpzTYdn1hOhIw4bk7pPvXrw99kTdjfg9HkKF\nBaRpIYztKGpi2vgjzocWYPnuW3SzGd/Fl+KdcCV6RkbrO7dDlw76XWmaWGTlc8u3acuXL+Wss8ZK\n0BeiBfH22Nt7N6JpGr6iQmw1NThNxriLmhgLC3A89gi2998FIDDqZDwzZhMeMDCu47VVlw76HfFg\n5t1332blyhXous6kSVPYuvVXPv30Y3w+H5mZPbj77nuZPPkKFix4mLS0NMaMOZWHH36Cgw46mIkT\n/8rixcvqk415vV7uuOMmqqurG1V2+vbb9Sxd+gS6ruP1erj11rv49tv1lJZG8tXceec93Hvv3ezZ\ns4fS0hJGjBjJlVdOjevzCNHdtKfHHu/diLusFGNpSSQTZry9e7cbx7NPY3/+WQwBPyHlEDyz5xE6\nuvmsux2lSwf9jpomlp6ewT/+cR+6rvP999/y4IOLAJg3bxabN29k5MjRrFv3GTk5vcnL68eXX67D\nYrEwcOB+jbJLvvHGqwwefCBXXTWNjRs3sH791wBs3fort9zyf2RlZbN8+VI+/vhDxo+fwNNPL+GO\nO/7Bnj1FDB16ONddN45AIMD55/9Rgr4Qtfbl84OEFDXRNKwr38K5+FGMpaWEs3NwT51B4KwxCSuB\n2BZdOuh31DSxupqyBoMBk8nMrbfegMPhoKRkD6FQiJEjT+aZZ5bQp09fJk+ezssvv0A4rDF69CmN\njrNz53ZOOOEkAIYMOQxz7cOenJwc7r//XpxOJ8XFe/jd746s3UNH13UyMjLYtOlHvvnmKxwOF8Fg\nlJV6QqSojn5+kKiiJuYvv8D50ALMW35Gt9vxTpqM99LLoIUa0B2tSwf9jpomVpda95dftrB69X9Z\nvHgZfr+PSZPGo+s6gwcfQEFBPuXlZUydOpNnnlnCmjWfcv/9jzY6zqBBg9mw4XtGjBjJTz9tJhTS\nAPjnP+/k5ZdX4HA4uOuu2+pLCBqNRsJhjXfeeYv09Azmz7+BXbt28tZbr7fr8wghYuepKIfiPfVF\nTT7f4eGtBnPqY8mIady2FefDD2JdszqyuGrM2XgmT0fv3b7iPHVtCcx9vYw443eXDvodfZvXv39/\nHA4n06dfia7rZGXlUFJSDMCRRx5dX5zjyCOPYdu2rU0KkZ977gXceeetzJhxFQMH7ofNZgXgzDPH\nMH36JBwOJ7169ao/5u9+dyTz51/NvHnXcdttN7Jhw/dYLBYGDNiPkpISsrOluIcQHSXg8xEsKowM\n5dR2/D7f4Wm0era1jJiGyorI4qrXXoksrjrqGDyz56Idcmi729eoLQZD3LM9JJ++ECKl6bqOp6gQ\nS2Uljr3m29+4qrDZPDkDMi3ceUaf3zYEg9hfeRH70icxVlej9R+AZ+YcgiNHt3lxVbQ7i73b8ta/\nxsW1aqtL9/SFEKI9PJUVULyHdGiUCbNOqxkxdR3LJx/jfPhBTLt3EU5Pxz3nb/gv+BNY2r64qqU7\ni6hZONtIgr4QIuU0LGpirn1I21wPu6WMmKbNm3A+uADLt+vRTSZ8f74Y78Sr0DN7xN2ulnLt9003\ns7sqFPex60jQF0KkDF3Xce/Zg6WivNGsnGg97NMOTGsS9LOqS7jhq1fI/PQ9AAInjcIzcw7hge3P\nMxStN7+7KsjEPyo8uVJt9zkk6AshUoK3qgq9uIgMXW8ylBOth60W+5k2LIu3N1dRVlLF+O9XcOba\n1zAF/IQOOhjP7LmEjj0uYW2MemeRncYJh/fDZDKzcu12dhXXxN3ljynoK4rSG/gKOA3QgGVAGNig\nquqM2vfcAowBgsBcVVW/jLdRQgiRKKFgEH9hAQ6vu9miJtDy2P3x/e2M/P5DnP9+BGNJCeGsLGqm\n/j2yuCrBKVOi5dqvm4Zetz4hJyc97mxsrQZ9RVHMwGNA3fLXBcANqqquVhRlkaIo44AdwEhVVYcp\nijIAeBVI3J8/IYRoo4ZFTTJaKGoC0XvYo0o3kTHhOsw/qeg2G96JV0UWVzkTX8DcH9IYckAPJvXK\n4v2v8jtstXEsPf37gEXA9UQyhh2tqurq2tfeBc4AVGAVgKqqOxVFMSmKkqWqatM/WUII0cHaUtQE\nmvaw+5bnM+HTpxn+yzoA/GeOwTt1OuHcPtEOEbegFsJrd2Lum4PL6eTEgXDiEf1b3zFOLQZ9RVGu\nAPaoqvqBoig31G5umCyiGsgE0oGGAb6mdrsEfSHEPtOeoiYAH63fzcgPnmPMt+9gCmsEjzwqUrnq\n0CEJb2tI0/BYbZj65OFydUyFsea01tOfAIQVRTkdOAJ4Bshp8Ho6UA5UARl7ba9IYDuFEKJFNSUl\nmMviLGoSDDJq7Ruc+dQTGKur0Pr1o3rGHIKjT0l45aqQpuGxWDH27oMrPT2hx45FzCtyFUX5CJgK\n3Av8S1XVTxVFWQR8BPwC3ENkqGcA8Kaqqke1dkxZkSuEaK9GRU3amrVS17Gs/iSyuGrnDsJpafgm\nXInvwovAak1oOzVNw222YMzOwdHOQik5Oelx/yWKZ8rmNcATiqJYgE3AK6qq6oqirAbWEhn3nxFv\ng4QQIhaapuErLIgUNTGb2pym2KRuxrnwfixffxVZXHXhRXgnXYXeo2dC2xkOh6kxmTDk9sHVjoVb\niSK5d4QQXU5dUZO0ZgJ9a1kxDcXFOB9/FOs7b2HQdQInnIhn1lzCg/ZPaBvD4TA1BgNkZePq2Suh\nx97XPX0hhOgUfo+HUFEhrmCg2aImLWbFzDFgf345jmefxuDzETrgwMjiquOOT2gbdV2nBtCzsnH2\n7IUhwc8E2kuCvhAi6cVa1KS5lbUGPcyeF9+gx+rlGIv3EO6Vhefqv+EfOy6hi6t0XadG1wn3zMKV\nlZV0wb6OBH0hRFJzl5dF6tMaDK3Oytl7Ze3QXT8y6ZMlHFT0C7rVivfyiXjHXwEJnCKp6zpuXUfL\n7IkrJydpg30dCfpCiKRUX9TE74u5Pm3dyto+FQVcsfoZTvx5LQBfHj6Kg26fT7hv34S2sUbT0DJ6\n4Ozdu+0zhzqJBH0hRFIJh8N49xRhqaxsc33a8/obcK9Ywthv3sESDrEx7xCeGjWRU88dzgF9E5c6\nwa1paBmZOHrndplgX0eCvhAiabRW1CSqUBDb669x2lOPY6yspKRHLktPHM/W349m7KGZrda0jbl9\nIY1QRgb23rnYE5xsbV+RKZtCiE4XDATwF+Tj8vvqi5rERNexrFmNc+EDmHZsR3e68E6YhO9PF4PN\nlrD2eUMawfR0bDm9McdRESvRZMqmEKJL0nUdd1ERlsoKMts4lGPa8jPOh+7H8uU6dKMR3/kX4p00\nBb1X4ubEe4MhgmlpWAfk4krwCt3OIkFfCNEpvFVVhPcUkkHbhnIMpSU4Fi/C9tabkcVVw4bjnT0X\nbfABCWubL6QRcDqx9h+IK4F3DMlAgr4QYp8KBYP4Cnbj9HkjRU1i5fNhf+E5HM8sxeD1Ehp8AN6Z\nVxMcfkLC2uYPafgdDiz9cnHZ7Qk7bjKRoC9EN7VuYxEr124jv8RDXraTMcMHxVSMI979WlNX1MRS\nVlo7lBNj+AmHsX7wHo5FD2MqKiLcsyeemVfjP+dcMCcmhO2d0747kwe5QnRD6zYW8fiKH5tsn3LO\n0BYDeLz7tcZXU4NWVEC6rrdp8ZL5u29wPng/5k0/olss+C6+BN9lE9DTEpOSuD6nfXYO9rS0hBxz\nX5AHuUKIRlau3RZl+/YWg3e8+0UTCgbxFRXg8HhwmWIvamLM343jkYewffQhAP7TzsA7bSbhvH5t\nbkOz7erknPadSYK+EN1Qfomn2e0Fpe4O2a85kaImJWSaYp+VY6ipxr5sCfaX/o0hGCQ09DA8s+cR\n+t0RbT5/c+py2hv69MaVkZmQY3Y1EvSF6Ibysp3sKm4aqPtmtZxzJt79GqorapKuhTDGOgUzFMK2\n4nUcTzyGsaICrU8fvNNnEzjtjIRUrkq2nPadqWutHxZCxGTM8EFRtu/XIftBbS969y5MO3eQoYdj\nTk9gWbuGzPEX47r3nxgCQTzTZlL571cJnP6Hdgf8cDhMNVCTnYNr8IE4Uzzgg/T0heiW6sbfV67d\nTkGpm75ZLsYM36/Vcfl493OXlWIsKSbDZIq5Pq3ply04Fj6Add3ayOKqcefhnTwNvVdWTPu3JNlz\n2ncmmb0jhIhbXVGTtFAw5p69oawUxxOPY1vxOoZwmOBxw/DMmot24EHtbk9XyWnfXjJ7RwixT4XD\nYbyFBVirqyOZMGMJ+H4/9hefx/H0UgweN9p+g/DMnktw+IntHsbpajntO5MEfSFEm7jLyyJDOUZj\nbEM5uo71w1U4Hl2IqbCAcGYmnmuuwz/uPDC3P3lZjaahZfbEmZPT5dIcdwYJ+kKImAR8PgIF+aRF\nqU/bHNOGH3A9+C/MG35At1jw/uWv+CZciZ6AufF1Oe3tOb1jbo+QoC+EaEVdfVpLVVXMmTCNBfk4\nFj2M7YP3AQicfCqe6bMI9x8Q83k/3+Hhrc1V5FcFycuwcPYhGRw/0Nktctp3JnmQK4SIqq6oSRrE\nNk7ursHx9FLsLz6PIRAgdOgQPHPm8b+eSrMBPJrPd3hYtK60yfYrTurLCccdmBQ57TuTPMgVQiRU\n0O8nUFSI0+eNrahJKITtrTdxPLEIY3k5Wu9cvNNmEjjjTD7f5WsUwHdVBut/jxb439pc1ez2DzdX\nM/LE1A747SVBXwhRr66oibWyIub6tOZ1a3E+dD/mX39BdzjwTJ6G7y+Xgt0BRA/gb2+uihr086uC\nzW6PJx2EaEyCvhACiBQ10YuLyND1mIqaGLf+inPhA1jXrkE3GPCdPS6yuCo7p9H7ogXwaNv9IY0+\nmVbyKwJNXmtLOgjRPAn6QqS4UDCIvyD/t6GcVsbuDeXlOJ58HNubr2HQNILH/B7PnHloBx3c7Pvz\nMizsqmwa4PMyGg/TNMxpf/bIns2meI4lHYRomQR9IVKUruu4i4uxlJfFNpQTCGB/6QXsy57E6Haj\nDdwPz8yrCY44qcU/FGcfktHsQ9mxh2QADXLa5/bFVZvTftiQyLBPW9NBiNbJ7B0hUlCbiproOpaP\nPsT56EJM+bsJZ2TivXIy/vMuiHlx1ec7PLzdYPbO2EMyOLafDY/Fiimnd5cqYJIM2jN7R4K+ECkk\nFAziK8zH6fXEVJ/W9OMGnA8uwPLDd+hmM74/XYTviivRMzLibkNdTntjdg6OdhwnlcmUTSFEi+rq\n05rLS2uLmrT81TcWFkQWV616D4DA6FPwTJ9NeEDsi6v2Jjntk4MEfSG6OV9NDdqeQtI1rfWiJm43\njmefxv78sxgCfkLKIZHKVUcfE/f5w+EwNQYDZOfg6tkr7uOIxJCgL0Q3pWkavsICbDU1uFrLhKlp\n2FauwPH4IoxlpYSzc3BPnUHgrDGxZdBshuS0T04S9IVIEus2FrFy7TbySzzkZTsZM3xQ3LNV2lLU\nxPzFOpwL78e85Wd0ux3vpMl4L70MHI64zp0qOe27Kgn6QiSBdRuLGs1L31Xsrv+9LYG/UVGTVoZy\njNu24nz4QaxrVqMbDPjHnI1n8nT03r3j+gyS075rkKAvRBJYuXZblO3bYwr6bSlqYqiswPHkYmyv\nvxJZXHXUMXjmzEVTDo2z9ZLTvitpNegrimIEngAUIAxMBfzAstrfN6iqOqP2vbcAY4AgMFdV1S87\nptlCdC/5JZ5mt8eSaybmoiaBAPZXX8K+9EmM1dVo/QfgmXU1wZNGxV25qi6nvaN3rgT7LiKWnv7Z\ngK6q6ghFUUYBdwMG4AZVVVcrirJIUZRxwA5gpKqqwxRFGQC8ChzXYS0XohvJy3ayq7hpgG8p10zA\n5yNYWIAr4G+5iIiuY/nkY5wPP4hp9y7C6em4Z8/Df+GfIc4UxZLTvutq9U+zqqpvApNrf90PKAeO\nVlV1de22d4HTgRHAqtp9dgImRVHaX9ZeiBQwZvigKNub5poJh8O4C/IxbN9GhhZqMeCbNm8iffpk\n0q+fj7GwAN+fL6by5Tfw/+XSuAK+N6RR5XBiGnwArr55UrGqC4ppTF9V1bCiKMuAc4E/EQnydaqB\nTCAdaJhgo6Z2e9OkG0KIRurG7VvLNeOprIA9Ra0O5Rj2FOF87BFs764EIHDSKDwzZhPeb1Bc7fMG\nQwTT0rAOyMVltcZ1DJEcYn6Qq6rqFYqi9Aa+BBrO5Uon0vuvAjL22l6RiEYKkQqGDcmN+tA26PcT\nKCzA6fe1XNTE48Hx3DPYn3sGg99P6CAFz+yrCR3beKQ1WinCvflCGgGnE2v/gbhstnZ9PpEcYnmQ\n+1egv6qq/wR8gAZ8pSjKKFVVPwHOAj4CfgHuURTlPmAAYFBVtazjmi5E9xdzURNNw/ruSpyPP4Kx\npIRwdjbuv11H4I9jm+yzdynC5ipZ+UMafocDS79cXHZ7x3w40Sli6em/BixVFOWT2vfPBjYDTyqK\nYgE2Aa+oqqorirIaWEvkQe+MDmqzECkh1qIm5q+/xPng/Zh/VtFtNrwTr4osrnK2rRTh25urOKaf\ntT6nvSvK/qJrkyybQiSZJkVNojDu2B5ZXLX6EwD8Z47BO3U64dw+LR5/wis7CTfzzTMa4OGZw7C7\npDpVspMsm0J0A7EWNTFUVuJY8gS2V1+KLK468ig8s+ehHTokpvNErWSVnSYBPwVI0BciCdQVNWlx\nKCcYxPbayzieegJjdRVav354ZswhOPqUNi2uilbJSkoRpgYZ3hGiE4WCQfyFBTi87uhFTXQdy+pP\nIourdu5F6L7CAAAdoklEQVQgnJaGb8KV+C68CNo4fbIup/3XZUY+/K5EShF2UVI5S4gupmFRE1dL\ni6vUzTgfWoBl/dfoJhP+8y7EO+kq9B4923S+cDiM22hE75UlOe27ARnTF6ILiaWoiaG4GOfjj2B9\n520Muk7ghBF4Zl1NeND+bTqX5LQXe5OgL8Q+omkavoJ8bG539KImXi/255fjePZpDD4foQMPwjNr\nLqHjhrXpXJLTXkQjQV+IfcBdWoKptJQMU5T0CeEw1vfewfnYIxiL9xDulYVn7jX4x5wTfUFWMySn\nvWiNBH3RoRJZDaoralzUpPn8huZv1uN8aAHmzZvQrTa8V0zC+9fLocH0yVjSJtRoGlpGD5y9e0ua\nYxGVPMgVHWbvalB1ppwztNsHfk3T8BUVYq2uxh5lCqZx506cjz6E9b8fAeA//Q94p88i3Kdvo/ft\nnTahzrRhWRw/0Ck57VOQPMgVSam91aC6KndZKcbSkqiZMA1VVdiXPYn95RcxhEIEDz8iUrlq6OHN\nHi9a2oQVmyoZckiO5LQXbSJBX3SY9lSD6orqipps+KWClT/VNB2KCQWxvf4qjicXY6yqROubh3f6\nLAKnnt7i4qr8qqarZwEKajRcffM66uOIbkqCvugw8VSD6orC4TDeokKsVVVszPfz+Jfl9a/tqgyy\n6PMSen+tctQLizDt2E7Y5cIzYza+P10MMaQrjpo2oZtdR7FvSNAXHWbM8EHNjul3p+X+nopyKN5T\nP5Sz91DMoOKtTPpkKUfu+B7daMR33gV4r5yK3iv2BVJnHpzGkw3+kNTpjOuY6g/muwMJ+qLDxFoN\nqisK+HwEiwoj9WkbPDytG4rp4S5n/JrnOG3DfzCi8/X+R3PgnX9HG3xAzOeoy2l/7ImHYupb2enX\nce8H87uK3fW/d4f/pqlCgr7oUC1Vg+qKdF3HU1SIpbIykglzr9ky+znC/P7jV7nwi1dxBH1szxrI\nU6MmUHLEcdw5uOWUx3WCWqhJTvthQ+ydfh1T9cF8dyNBX4gYeSsr0YuLSIemmTDDYawfvMc9ixbi\nKN1DuTOTp0ZN4IPDTiNsNDHtkIxmj9lQSNPw2OyY+uThSsIUx6n2YL67kqAvRCuCgUCkPm2Uoibm\n776NLK7a+CO61cov4y7lkSHn8mvAQr8MC2Oj1J+tE9I0PBYrxt59cKWnd+RHaZdUeTDf3UnQFyKK\n1oqaGHfvwvHoQmwffQiA/7Qz8E6fRc++edwUw/E1TcNttmDsk4sro/U7gc6WCg/mU4EEfZHympuR\n8rv+DsJR6tMaaqqxL1uC/aV/YwgGCQ09DM/seYR+d0RM56vLaW/I7YMrs0cHfKKO0Z0fzKcSScMg\nUlq0VBGTf9+DEwftNdQSCmF78zUcTz6OsaKC0swclo64jF9+fzJnH5rZ4hAOSE57kTiShkGIOEWb\nkfLuT+7fgr6uY1m7BufCBzBt20rQ4WT5iPGsOGosAYsNqkL1uXGaC/yS014kEwn6IqVFm5FSN9/e\ntOVnnAvvx/LFusjiqnHnc8vB57NJS2uyz9ubqxoFfclpL5KRBP1uSlZOtk4LheiTaSG/ItDkNcVY\ng/Ofd2F76w0M4TDB44bhmTUX7cCDUF/Z2ezx6v5Q1Oe079ELV3a2BHuRVCTod0OycrJ1NSUlmMtK\nGaekN0pbbAkFGLf+LS75+lUsXg/aoP3xzJ5L8PgT6pOiRc2Fk2GJ5LTP7IkzJ0fSHIukJEG/G5KV\nk9H5PR5ChQWkayGMJmP9cMzbmyoZ/OXHTFiznKyKPYR79MA9czb+c84Dc+OvydmHZDSb3/7UI3Kw\nHXiwBHuR1CTod0OycrIpTdPwFRZgq6nBuVf6hBGVWzjjlQWYN/yAbrHgvfQyfJdPRI+yUKr+D8Xm\nKnZXBenTw8aYEwcz/LC+zb5fiGQiQb8bkpWTjUUramIsyI8srvpwFQD+U06LVK7q17/VYx6RZ2OI\nMhhbTm/MFkuHtV2IRJOg3w3JysmIgNdLsLAAVzCAqeFqWncNjqeXYn/xeQyBAKFDh+KZM4/QEUe2\nekxvMEQwLQ3bwD64JNiLLkiCfjeU6isnGxY1aZQ+IRTC9tabOJ54DGN5GVpuLt5pswic/ocm2TL3\n5gtpBJxOrAP2w2W17oNPIUTHkBW5olvxVJRjKCkmba9pkpbPP8Ox8AHMv/6C7nDgHX8Fvr9cCnZH\ni8ery2lv6Z2L1W7vyKYLETNZkStSRrT1B9GKmhi3/opz4QNY165BNxjwnXMu3qumomfntHie5nLa\nC9EdSE9fdBnR8uRcMTKP47PA0eAhraG8HMeTj2F783UMmkbwmN/jmTMP7aCDWzxHSNPwWG2Ycnpj\nT8Kc9kKA9PRFioi2/uA/3xVz8hm1Van8fuwvv4B92VMY3W60gfvhmXU1wRNPql9c1Zy6nPam3L64\n0pqmWBCiu5CgL7qMFvPk6DqWjz7E+ehCTPm7CWdk4p43H/95F4A5+iybrpbTXoj2kqAvkk60cfto\n6w9OqPqF9Cm3YvnhO3SzGe9fLsV3xZXoLQTxrprTXoj2kjF9kVSijdtPOH0welUlyxqkP8ipKuay\n/y1n9OZPAQiMPgXP9NmEBwyIevxwOEyN0QiS0150YTKmL7qNaOP2H3yxg7vO6IsD+OC7Qk74z4uM\n+3oF1lCA0CGHRipXHXV01ONKTnshIloM+oqimIElwCDACtwFbASWAWFgg6qqM2rfewswBggCc1VV\n/bLDWi26rah5g6pCoGmM/HYVZy5+FGNpKeGc3tRMnUHgzD9GXVyVCjntJY22aIvWevp/BUpUVb1M\nUZSewLe1/7tBVdXViqIsUhRlHLADGKmq6jBFUQYArwLHdWjLRbcUbdz+5OINZFwxH/OWn9HtdjxX\nTsF3yXhwNL+4au+c9l9s2sPKN7/odoFR0miLtmot6L8EvFz7sxEIAUerqrq6dtu7wBmACqwCUFV1\np6IoJkVRslRVbZp/VogotFCIU5U0nm4Q9PuX7WLCp8s47tev0A0G/GPOwTNlOnpO9MVVe+e0786B\nUdJoi7ZqMeirquoBUBQlnUjwvxG4r8FbqoFMIB1oGOBrardL0G+nVLl1rytqckpfK85hWXz89U5G\nr3qWM79/D1M4TPCoY/DMmYumHBr1GG5NQ8vIxNE7t1FO++4cGCWNtmirVh/k1g7XvAY8rKrqC4qi\n/L8GL6cD5UAVkLHX9opENjQVdeceap29i5oQCDD6f69y5tInMdbUoPUfQPXMOQRHjo66uMoT0ghl\nZGDvnYu9YTbNWt05MEoabdFWLaYWVBQlF3gfuFZV1adrN3+jKMrI2p/PAlYDnwFnKIpiUBRlIGBQ\nVbWsoxqdKlrqoXZ1mqbh3r0L084dZOhhjAYDlo//Q+YlF+Jc+AAYjLiv/huVz79McNTJzQZ8b0ij\nyuHENPgAXH3zGqdPbiAvu/ncOd0hMI4ZPijK9tRKoy1i11pP/3qgB3Bz7ewcHZgDLFQUxQJsAl5R\nVVVXFGU1sBYwADM6sM0po7v2UPcuamLatBHnQwuwfPsNusmE76JL8E64Ej0zs9n963LaWwfkxpTm\nuDvXF0j1NNqi7WRxVhK75al1zd66989J445JiZ8c1dHPD/weD6GiQtJCQYxGI4Y9RTgXPYLtvZUA\nBE4ahWfmHMIDmw/G9Tnte+disdnadO7IZ5PAmGxS5ZlVorVncZYE/SQWbXXqlHOGJvyL0Z5ztfbF\nDYfDeAsLsFZXYzebwOPB8dwzWJ99BlPAz685+/PmmKtQzjqpvv5sQ5LTvnval/++uxtZkdtN7ctb\n93hnuLT2sLmuqEmGwQAGsL69Aufjj2AsKaHU1ZPlJ0/m40NHEzaa+Kg2xUJd4Jec9t1bd55Vlcwk\n6Ce5YUNy98kXIN7nB9G+uG9/tpXDXF5cfh8mkwnz11/ifOh+zD+p6DYbK0f9hWWHn4PP2nhx1dub\nqzi2nw2PzY6pTx4uyWnfbXXXZ1bJToK+AOKf+hf9i+shIxTEuHsXzocfxLr6EwD8Z47BM20Giz8N\nEG5mcG93VRB/Xn/JaZ8CZLpp52i5GrRIGfFO/Ys2HfIgsw/n/feRecmfsK7+hOCRR1O5ZDnuW+9A\n751LXkbzOe7zstOwS8BPCTLdtHNIT18A8T8/2Hs6pFkL8sfv3uXyL1/G6q5G69cfz8w5Teban31I\nBovWNV2wLV/41CHTTTuHzN4RcdN1HXdREd/9WMAqtZoB6//HpP89Q++yfMJpafgmXoXvwovA0rhX\nHw6HqTEY+LrcxIffl8oXXog2kimbYp/zVlWhFxeRpuuYf1Iji6vWf41uMuE/70K8k65C79Gz0T71\nOe17ZaVsTnuZly4SQaZsin0mFAziL8jH6fNiKSvD+fgjWN95G4OuEzjxpMjiqkH7N9onFXLaxyIV\ncimJ5CdBX8RE13XcxcVYysvICAawP78cx7NPY/D5CB14EJ5ZcwkdN6zpPg1y2qdqsK8j89JFMpCg\nL1rlq6lBKyogQ9OwrXoX52OPYCzeQ7hXFu658wmMORv2SnZWo2loGT1w9u7dKM1xKpN56SIZSNAX\nUYWCQfyFBdi9HjK/+xbnQwswq5vRrTa8l0/EO/4K2GvxVLSc9kLmpYvkIEFfNCtS1KSEHvn5OB95\nEOsnHwPg/8NZeKfOINynb6P3t5bTXnTvbJ+i65CgLxqpK2qSUV6G6+kl2F55EUMoRNnBh/HEqAl8\nlnEAed8bODvg4fiBTrwhjWB6Orac3tgszS+4EhEyL10kA5myKYBIURNfYQG2igoy3nodx5OLMVZV\nouX147uLpnArhzcpZHLFiD6cMOwgzBLshdinZMpmF5DM87PdZaUYi/eQ9flnOBc+gGnHdsIuF54Z\nc/D9+WKe+G8ZVAab7PehWsPIEd074Cfzfzch4iFBP4GiBYhknZ9dV9Qkc9OPpD38IJavvohUrjr/\nT3ivnILeM7K4Kr+qacCH7j/rJFn/uwnRHhL0E6SlAJFs87PriprYt24le8nj2N5eEVlcNfxEPLOu\nJrz/4Pr3BrUQfTKt5FcEmhynu886Sbb/bkIkggT9BGkpQCTT/Gx3eRmmXTvp/eLzOJYvw+D1Ehp8\nAN5ZVxM8/oT694U0rT6n/dkje6XkrJNk+u8mRKJI0E+QlgJEMszPDvh8BPN302vlClyLH8VUVES4\nZy88s+fhH3sOmCP/FEKahsdixZTbtz6n/bAhkXam2qyTZPjvJkSiSdBPkJYCxJjh+3VaTzkcDuMt\nKsS1+hNyH3kQ86aN6FYr3vFX4L18ArgigV3TNNxmC8Y+ubgyMpocZ19V8EomMq9edEcS9BOkpQDR\nWfOzPZUVmL9dT99HF2L7+D8A+E//A95pMwn3zQNq0xybTBhy++DK7NGh7eloiZ5pI/PqRXck8/QT\nKBJ0Oj9ABP1+Qlt+ImvxIpyvvIghGCR02OF4Zs8jdPjvgEiwd5tM6D174erZa5+3MdH2fpBeZ8o5\nQyVIi25H5uknic4eAtF1Hffu3fR8dhkZSxZjrKhA69MX74zZBE49HQyG33LaZ2V3q5z2MtNGiNhI\n0O8mvJWV2N56nQEP3Y9521Z0pwvPtJn4LroEbLZIsA+Hu21Oe5lpI0RsJOh3caFgEG31x+QuuBfb\nF+vQjUZ8516A96op6L2yIr3/cLjb57SXmTZCxEaCfhel6zreTZvIWnAPrrffxBAOEzxuGJ5Zc9EO\nPAiozWmf2RNnTk63T3MsM22EiI0E/S7IX1KM66EF5DyzFIPHgzZofzyz50YWVxkMKZnTXmbaCBEb\nmb3ThYQCAcxPL6HXQwswFRUS7tED75VT8I87H8zmRjntTZLTXohuS2bvdHO6rhP6cBVZd9+O7ccN\n6BYL3r9eju/yCehp6XiCIUIOp+S0F0K0SoJ+kgtu3kTmHTfj+nAVAP5TTsM7fRbhfv3xBkMEbXZs\nA/tIsBdCxESCfpIKV5Rjv/sOsp9fjiEQIHToUDxz5hE64kh8IY2A1YZ1wH64rNbObqoQoguRoJ9s\nQiEMixeR/cB9mCrK0Xrn4p02k8AZZ+IP6/gtViz9cnHZ7Z3dUiFEFyRBP5m8u5LM22/C+usv6A4H\nninT8V18CUGLBa/Njjk7B5fT2dmtFEJ0YRL0k4Bh00acN8zHuWY1usGA75xz8V41lWDPXpGc9tk5\nuFxNFxlJKT8hRFtJ0O9EhuJirHfeSvpL/8agaQSPPQ7P7Ln4Bx8QyWmf07s+p/3epJSfECIeEvQ7\ng8+HddFC0h+6H6O7Bm3gfnhmXo1v+Am4LVaM2TnN5rRvSBKMCSHiEVPQVxRlGPBPVVVPVhTlAGAZ\nEAY2qKo6o/Y9twBjgCAwV1XVLzumyV2YrmN54zVct9+EJX834cxM3POuxTvuPGpsNgxZ2THntJcE\nY0KIeLQa9BVFmQ+MB2pqNy0AblBVdbWiKIsURRkH7ABGqqo6TFGUAcCrwHEd1eiuyLz+Kxw3zMe+\n/mt0sxnvX/6K57IJVPfoAb2y2pzTXhKMCSHiEUtPfwtwHrC89vdjVFVdXfvzu8AZgAqsAlBVdaei\nKCZFUbJUVS1NdIO7GuOunTjuuAXnG68CEBh1Mu4Zs6nqPwC9VxbOnr34YtMeVr62rk0PZCXBmBAi\nHq0GfVVVX1cUpWEkaZjzoRrIBNKBhgG+pnZ7ygZ9Q001jgcX4HjsYYx+P6FDDsU9ay4VRx7VKKd9\nvA9kJcGYECIe8TzIDTf4OR0oB6qAjL22V7SjXV2XpmH/97M4774DU0kx4ewcqq+dQfkZZ6H1ym6S\n0749D2Q7u1KXEKLriSfor1cUZaSqqp8CZwEfAb8A9yiKch8wADCoqlqWwHZ2CZZPPsZ1yw1YNv2I\nbrfjuXIKpRddQqhPXtSc9vJAVgixL8UT9K8BnlAUxQJsAl5RVVVXFGU1sJbI8M+MBLaxU7Rl4ZPp\nJxXXbTdh+/B9dIMB/5izKblqKsHBB+LonYu9hZz28kBWCLEvST79Zuw9zl5nyjlDGTYkt/4PQvWO\nQq785mVGfrESY1gjePQxlM24Gu+xv485p31r5xJCiL1JPv0Ea2mcHeCp177l7G9X8ud1L5Hm97C7\nRx47J06nz/jzsfXOxdWGNMfyQFYIsS9J0G9G1HH2khryn1jOo28/Tt/KIqptaSwePYl3jziTPjk9\nuKNf/7jOJw9khRD7igT9ZjQ3zn5g4Ramf/Y0B237gZDRxJtHjeWF4y+ixpEOQEFp838ohBAimUjQ\nb0bDhU9Z1SVc9r9nOWXTfwH4VjmeRSeMJ79nv0b7yINXIURXIEG/GcOG5GLyuAnfdx+nffoStlCA\n8v0OxH/99WwdPIz897Y02UdWwgohuoKkC/qdniNe07C/+Dxn3H0Hpj1FhLKy2TN9Fv7LJ2HPyGAE\nYLHa5MGrEKJLSqopm509fdHyv08ji6s2fE/YZqf8kvF4512LLVcCuhAieXSbKZsdkSM+ljsH05af\ncd1xM7b33gGg/I9j8Vx/C1blEGxxnVUIIZJTUgX9RKckaC2ZmaGsFOe/7sGx9EkMoRA1Rx1D5Y23\nYh05GmtcZxRCiOSWVEE/0SkJot05vLd6C6NXv4LzX/dgrKjAP2Agpdf8HcvFl2I1xH3XJIQQSS+p\ngn6ic8Q3uXPQdY7/ZR0TPn2atIoCtPR0Cq/5O8bZ87DY7XGdQwghupKkCvqJTknQ8M7hgKJfmPjJ\nUn63awOa0UjpRZcQvOk2TLl9EtZ+IYRIdp06e2fcNSv0REzLjPawdt3GIl5+7lPGf/Ycp/z4MUZ0\n1g3+PZXX38RR405O3AcRQoh9qMvO3gnresyVoqKJ9rDW5PVw8kcvcObyBzD7vGzNGcQbYyZz4OUX\ncvxQ6d0LIVJT0gzvxDstc++HtQY9zMkb/8vwJ5/DVVVKKCuLPddci33yDP4q4/ZCiBSXNEE/3mmZ\nDR/WHrZzA5M+WcKBe37Fb7ZSfvkkaubOx56XlzwfVAghOlHSxMJ4p2XmZTvRftrCFauf5oQtnwPw\n8aGjeH/MJObNH4ddpmAKIUS9pAn68UzLNFSUM3/98wx4dTmWcIiNeYfy5KgJ/Nz3YKaMHdqoALkQ\nQohODvomoyG+aZnBII5lT+K4959kV5Tj7t2XpaOvYOWA48jLTmeKJECLqtMT2gkhOlVSJVxrla5j\nXfUezltvwPLrL4RdLqrHT6Bi0mQcA/eTnn0rOjuhnRAiMbrslM22MP3wPa5brse2ZjW6yYT3/AvZ\nM2UG5kOH4pRZOTHpiIR2QoiuJemDvrGoEPudt+N86XkMuk5g+ImUzJhN8Ljjcfbo2dnN61ISndBO\nCNH1JG/Q93iwP/IgrocfxOj1EBp8AJUz5lB92hk4cvtgMRrbfYpUG99OdEI7IUTXk3xBPxzG8vIL\npN15G+aiQsI9e1E9cw6l512Ipf8AXAkaymkt7XJ3lOiEdkKIriepgr5p7RpcN16LbcMP6FYr3ssm\nUPLXy9AGDsLVKyuh50rF8e1EJ7QTQnQ9SRH0DVt/xXHT33F98B4A/tP/QOXk6XgOOhh7bh/sJlPC\nz5mq49vDhuRKkBcihXVu0K8ox3L3/5H53NMYgkGCh/0O96yrqTjiKMx9+uJyOjvs1DK+LYRIRZ0a\n9HsddySminK0Pn3xTp9F+cmnoWVn48rK7vBzy/i2ECIVdWrQNwQCeKbNpPKCi/D36oW9b16HDOU0\nR8a3hRCpqFNX5O54/yM93KMnpt59sKeldVo7hBCiK+myK3LDBx6MMz1D0icIIcQ+0qlB35WR2Zmn\nF0KIlNP+Za1CCCG6DAn6QgiRQiToCyFECpGgL4QQKSShD3IVRTEAjwJHAD7gSlVVf03kOYQQQsQv\n0bN3zgVsqqqeoCjKMGBB7bZ2S7U0yEII0RESPbwzAngPQFXVdcCxiThoXRrkXcVuwrpenwZ53cai\nRBxeCCFSRqKDfgZQ2eD3kKIo7T5HS2mQhRBCxC7RQb8KSG94fFVVw+09aKqmQRZCiERL9Jj+GmAs\n8IqiKMcDP7T05ljzR4R1/Xvg8L23a2H9+5yc9CPiaagQQqSihCZcazB753e1myaoqvpTwk4ghBCi\nXTo1y6YQQoh9SxZnCSFECpGgL4QQKUSCvhBCpBAJ+kIIkUL2eREVyc8DiqKYgSXAIMAK3AVsBJYB\nYWCDqqozOqt9+5qiKL2Br4DTAI0UvQ4AiqL8HTgHsBD5nnxKil2P2u/H00S+HyHgKlLw30VtKpt/\nqqp6sqIoB9DM51cU5RZgDBAE5qqq+mVrx+2Mnn59fh7geiL5eVLNX4ESVVVHAmcBDxO5DjeoqjoK\nMCqKMq4zG7iv1H7BHwPqVuCl5HUAUBRlFDC89rsxGhhIal6PPwImVVVPBP4PuJsUuw6KoswHngBs\ntZuafH5FUY4CRqqqOgz4C/BILMfujKDfIfl5upiXgJtrfzYS6c0crarq6tpt7xLp9aaC+4BFQD5g\nIHWvA8AfgA2KorwBrADeJjWvx0+AuXZUIJNILzbVrsMW4LwGvx+z1+c/nUgsXQWgqupOwKQoSlZr\nB+6MoN8h+Xm6ElVVPaqquhVFSQdeBm4kEvDqVBP5x96tKYpyBbBHVdUP+O3zN/y3kBLXoYFs4Bjg\nQmAa8BypeT1qgP2BzcDjwEOk2PdDVdXXiXQG6zT3+dNpHEtriOG6dEaw7ZD8PF2NoigDgI+Ap1VV\nfYHIWF2ddKCiUxq2b00ATlcU5WMiz3ieAXIavJ4q16FOKfC+qqqh2pXsPhp/iVPleswF3lNVVeG3\nfxfWBq+nynVoaO/4UE4klmbstb3V69IZQX8NkTE7YsnP0x0pipILvA9cq6rq07Wbv1EUZWTtz2cB\nq5vduRtRVXWUqqonq6p6MvAtMB54N9WuQwP/A84EUBQlD3AB/6kd64fUuR5l/NaDrSAy4eSbFLwO\nDa1v5nvxGXCGoigGRVEGAgZVVctaO9A+n70DvE6kd7em9vcJndCGznY90AO4ufbpuw7MARYqimIB\nNgGvdGL7OtM1wBOpeB1UVV2pKMpJiqJ8QeR2fhqwDXgyxa7HA8ASRVE+JTKL6e/A16TedWioyfdC\nVVVdUZTVwFoi/15imtEkuXeEECKFpNQDVCGESHUS9IUQIoVI0BdCiBQiQV8IIVKIBH0hhEghEvSF\nECKFSNAXQogUIkFfCCFSyP8H2FL5ThOMMIsAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fit.update(fix=[(10, 70)])\n", "fit.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "For more informations, look into the formatoptions of the ``regression`` group" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "yrange\n", " Specify the range for the fit to use for the y-dimension\n", "fit\n", " Choose the linear fitting method\n", "line_xlim\n", " Specify how wide the range for the plot should be\n", "ideal\n", " Draw an ideal line of the fit\n", "xrange\n", " Specify the range for the fit to use for the x-dimension\n", "fix\n", " Force the fit to go through a given point\n", "ci\n", " Draw a confidence interval\n", "nboot\n", " Set the number of bootstrap resamples for the confidence interval\n" ] } ], "source": [ "fit.summaries('regression')" ] }, { "cell_type": "code", "execution_count": 10, "metadata": { "collapsed": true }, "outputs": [], "source": [ "psy.close('all')" ] } ], "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.6" } }, "nbformat": 4, "nbformat_minor": 2 }