{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "## pandasプロット(鋭意制作中)\n" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": false }, "outputs": [], "source": [ "# python用のパッケージ\n", "import pandas as pd\n", "import numpy as np\n", "import matplotlib.pyplot as plt \n", "import seaborn as sns\n", "%matplotlib inline\n", "\n", "# jupyter用のdisplayメソッド\n", "from IPython.display import display, Latex, HTML, Math, JSON" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "collapsed": true }, "outputs": [], "source": [ "# 9章のデータをRからJSON形式に変換\n", "d = pd.read_json('''\n", "[{\"x\":3,\"y\":5},{\"x\":3.2105,\"y\":3},{\"x\":3.4211,\"y\":6},{\"x\":3.6316,\"y\":7},{\"x\":3.8421,\"y\":7},\n", " {\"x\":4.0526,\"y\":5},{\"x\":4.2632,\"y\":9},{\"x\":4.4737,\"y\":9},{\"x\":4.6842,\"y\":7},{\"x\":4.8947,\"y\":10},\n", " {\"x\":5.1053,\"y\":12},{\"x\":5.3158,\"y\":8},{\"x\":5.5263,\"y\":7},{\"x\":5.7368,\"y\":4},{\"x\":5.9474,\"y\":4},\n", " {\"x\":6.1579,\"y\":11},{\"x\":6.3684,\"y\":9},{\"x\":6.5789,\"y\":9},{\"x\":6.7895,\"y\":8},{\"x\":7,\"y\":6}]\n", "''')" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/html": [ "
\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
xy
count20.00000020.000000
mean5.0000007.300000
std1.2454862.408319
min3.0000003.000000
25%3.9999755.750000
50%5.0000007.000000
75%6.0000259.000000
max7.00000012.000000
\n", "
" ], "text/plain": [ " x y\n", "count 20.000000 20.000000\n", "mean 5.000000 7.300000\n", "std 1.245486 2.408319\n", "min 3.000000 3.000000\n", "25% 3.999975 5.750000\n", "50% 5.000000 7.000000\n", "75% 6.000025 9.000000\n", "max 7.000000 12.000000" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "d.describe()" ] }, { "cell_type": "code", "execution_count": 28, "metadata": { "collapsed": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYEAAAESCAYAAAAbq2nJAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAEKdJREFUeJzt3X+sZGV9x/H3XN2F/WW4i1dLsS420S81aiJtGqJWAY2U\nqtGWWi2IQaMR0MaarQn9AciPxMZ0jbVW2dhCq1EiVsQ2NggUsT+wTRPb2l98K2qpBapLckt3vcCl\nu6d/zNzu3XXvZeYyZ56Z87xf/9zZw5k932cOs595zrnzfHtN0yBJqtNc6QIkSeUYApJUMUNAkipm\nCEhSxQwBSaqYISBJFWs9BCLieRFxT0RcctT2syPiUNvHlyStrdUQiIitwIeB24/afhxwKXB/m8eX\nJK2v7ZnAI8A5wANHbf814CPAcsvHlySto9UQyMxDmfno6m0R8RzgBZn5OaDX5vElSet7coFjfhD4\npQLHlSQdpTeJtYMi4gpgH3Az8JXB4x7wQuCrmXnmWs9tmqbp9ZwwSNKIhvqHc5IzgV5m3g88e2VD\nRHx7vQAA6PV67Nu3v/XiSllY2OH4Zpjjm11dHhv0xzeMVkMgIk4D9gC7gMci4lzg5zLzvwe7uISp\nJBXUaghk5teANT/pZ+aPtnl8SdL6/MawJFXMEJCkihkCklQxQ0CSKmYISFLFDAFJqpghIEkVMwQk\nqWKGgCRVzBCQpIoZApJUMUNAkipmCEhSxQwBSaqYISBJFTMEJKlihoAkVcwQkKSKGQKSVDFDQJIq\nZghIUsUMAUmqmCEgSRUzBCSpYoaAJFXMEJCkihkCklQxQ0CSKvbktg8QEc8DbgY+mJkfjYgfAa4D\nNgHLwJsy83tt1yFJ+kGtzgQiYivwYeD2VZuvBq7NzDPoh8PuNmuQJK2t7ctBjwDnAA+s2nYxcNPg\n8T5gZ8s1SJLW0OrloMw8BDwaEau3PQwQEXPAO4Er26xBmqTFxUX27L2BxSWY39Kw++LzmD9hvnRZ\n0pqK3BgeBMAngT/LzC+XqEFqw569N3BfEzy89VTu41T2XHtD6ZKkdbV+Y3gN1wOZmVcPs/PCwo6W\nyynL8c221eM7sDxHb3MPgF6vx4HluZkf/6zXv54uj21YEw+BiDgfeDQzrxr2Ofv27W+xorIWFnY4\nvhl29Pi2bzrIQ01Dr9ejaRq2bz400+Pv8vnr8thg+IBrNQQi4jRgD7ALeCwifh54GvBIRHwZaIB/\nycx3tVmHNCm7Lz6PPdcO7glshd0X/WLpkqR1tX1j+GvAmW0eQ5om8yfMc82ll5QuQxqa3xiWpIoZ\nApJUMUNAkipmCEhSxQwBSaqYISBJFTMEJKlihoAkVcwQkKSKGQKSVDFDQJIqZghIUsVK9ROQNMPs\noNYdzgQkjcwOat1hCEga2eJSv3Ma9H8uLhUuSBtmCEga2fyWhqZpAGiahvmthQvShhkCkka2++Lz\nOLmXbFm6m5N7aQe1GeaNYUkjs4NadzgTkKSKGQKSVDFDQJIqZghIUsUMAUmqmCEgSRUzBCSpYoaA\nJFXMEJCkihkCklQxQ0CSKtb62kER8TzgZuCDmfnRiHgG8En6AfQAcEFmPtZ2HZKkH9TqTCAitgIf\nBm5ftfkq4Hcy82XAN4G3tlmDJGltbc8EHgHOAS5dte0M4B2Dx38C7Ab2tlyHtCG2UVTXtToTyMxD\nmfnoUZu3rbr88z3gpDZrkJ4I2yiq60r0E2hWPe4d9edjWljY0V41U8DxTa8Dy3P0Nh9uo3hgee4H\nxjPL4xtGl8fX5bENq0QIHIiI4wYzhJPp3xxe1759+9uvqpCFhR2Ob4pt33SQh5qGXq9H0zRs33zo\niPHM+vgeT5fH1+WxwfABV+JXRG8Hzh08Phe4pUAN0lBso6iu6600i25DRJwG7AF2AY8B9wHnA38I\nHAfcC7wlMw+u89c0XU9rxze7HN/s6vLYABYWdvSG2a/Vy0GZ+TXgzGP8p1e2eVxJ0nD8xrAkVcwQ\nkKSKGQKSVDFDQJIqZghIUsUMAUmqmCEgSRUzBCSpYoaAJFXMEJCkihkCklQxQ0CSKlain4C0rlFa\nOtr+cTaMep48r5PjTEBTZ5SWjrZ/nA2jnifP6+QYApo6i0v9Vo7Q/7m4NJ59Vc6o58nzOjmGgKbO\n/JaGlWZHTdMwv3U8+6qcUc+T53VyDAFNnVFaOtr+cTaMep48r5PTanvJMbG95AxzfLOty+Pr8thg\n+PaSzgQkqWKGgCRVzBCQpIoZApJUMUNAkipmCEhSxQwBSaqYISBJFTMEJKlihoAkVWzi/QQiYhvw\nCWAnsAm4KjNvnXQdkqQyM4ELgbsz80zg9cBvF6hBkkSZzmIPAs8fPN4J7CtQg1Sc3bPGx9dy4x53\nJhARZ0fEUKvRDSMzPwPsiohvAHcCvzKuv1uaJXbPGh9fy40bZibwbuB3I+JTwHWZee8TOWBEnA/c\nm5nnRMQLgN8DfnK95yws7Hgih5x6jm+2bXR8B5bn6G0+3D3rwPLcVL5W01jT0Tb6Ws7C2Nr2uCGQ\nmT8TEfPAzwIfiwiA64GbMvPgBo75YuBLg7/76xFxckTMZeahtZ7Q8TW/Hd8MeyLj277pIA81Db1e\nj6Zp2L750NS9VrNy/jbyWs7K2DZq2IAb6sZwZi4CNwCfBk6gfwnn7yPi9A3Udg9wOkBE7AL2rxcA\nUlfZPWt8fC037nFnAhHxU8BbgLOAm4C3ZubdEXEK8HnghSMecy9wXUTcCTwJeMeIz5c6Yf6Eea65\n9JLSZXSCr+XGDXNP4P3AtcDFmfnoysbM/PeIuHHUA2bm94E3jPo8SdL4DXNP4CXr/Lf3j7ccSdIk\nuWyEJFXMEJCkihkCklQxQ0CSKmYISFLFDAFJqpghIEkVMwQkqWKGgCRVzBCQpIoZApJUsRLtJTVl\nRm3NN8ut/KxdOpIzAY3cmm+WW/lZu3QkQ0AsLvVb8kH/5+LSePefJtYuHckQEPNbGpqmAaBpGua3\njnf/aWLt0pEMAY3cmm+WW/lZu3Sk3soniynWdL0ZtOObXY5vdnV5bAALCzt6w+znTECSKmYISFLF\nDAFJqpghIEkVMwQkqWKGgCRVzBCQpIoZApJUMUNAkipmCEhSxYr0E4iI84H3Ao8Bl2XmLSXqkKTa\nTXwmEBE7gcuBFwGvBl436RokSX0lZgKvAG7LzCVgCbioQA2dZxcq6dhW3hsHlufYvulg9e+NEvcE\nTgG2RcQXIuIrEXFWgRo6zy5U0rGtvDf+Z/NzfG9QZibQA3bSvwz0LODLwK71nrCwsGMCZZXTxvgO\nLM/R23y4C9WB5blir6Pnb7Z1bXzT9N6YBiVC4LvAXZnZAN+KiP0R8dTMfHCtJ3R8ze9Wxrd900Ee\nahp6vR5N07B986Eir2MFa7Y7vhkzLe+Ntg0bbCUuB90KnBURvYh4KrBtvQDQxtiFSjq2lffGU5b/\nzfcGhTqLRcTbgbcBDXB1Zn5xnd3tLDbDHN9s6/L4ujw2GL6zWJHvCWTmx4GPlzi2JOkwvzEsSRUz\nBCSpYoaAJFXMEJCkihkCklQxQ0CSKmYISFLFDAFJqpghIEkVMwQkqWKGgCRVzBCQpIoVWUBOdbW4\ns9WlNL2cCRRSU4s7W11K08sQKGRxqd/aDvo/F5cKF9SimsYqzRpDoJD5LQ0rDX2apmF+a+GCWlTT\nWKVZYwgUUlOLO1tdStOrSHvJEdlecoY5vtnW5fF1eWwwfHtJZwKSVDFDQJIqZghIUsUMAUmqmCEg\nSRUzBCSpYoaAJFXMEJCkihkCklQxQ0CSKlYsBCLi+Ij4ZkS8uVQNklS7kjOBy4AHCx5fkqpXJAQi\nIoBTgS+WOL4kqa9Ue8k9wDuBCwsdX+q8Udt6Hmv/hYUdE6xYJUx8JhARFwB3Zea9g01DLXcqaTSj\ntvW0DWidSswEXgU8KyJeAzwDeCQivpOZd6z1hK5/GnF8s21ax3dgeY7e5sNtPQ8sz61b67H2h+kd\n3zh0eWzDmngIZOYbVx5HxBXAt9cLAKDrjR8c3wyb5vFt33SQh5qGXq9H0zRs33xo3VqPtT909/03\nzeduHIYNOL8nIHXUqG09bQNaJ9tLFlbDpxHHN7u6PL4ujw1sLylJGoIhIEkVMwQkqWKGgCRVzBCQ\npIoZApJUMUNAkipmCEhSxQwBSaqYISBJFTMEJKlihoAkVaxUZzGNaNQuUZI0DGcCM8KuT5LaYAjM\niMWlfrcn6P9cXCpckKROMARmxPyWhpXeD03TML+1cEGSOsEQmBF2fZLUBm8Mz4j5E+a55tJLSpch\nqWOcCUhSxQwBSaqYISBJFTMEJKlihoAkVcwQkKSKGQKSVDFDQJIqZghIUsUMAUmqWJFlIyLiA8BL\ngCcBv5mZny9RhyTVbuIzgYg4A3huZr4IOAf40KRrkCT1lbgc9BXg9YPHi8DWiOgVqEOSqjfxy0GZ\n2QAPD/74duBPB9skSRNWbCnpiHgt8BbglaVqkKTa9Va6VU1SRJwNXAmcnZkPTbwASRJQIAQi4inA\nXwAvz8wHJ3pwSdIRSlwOegNwInDj4IZwA7w5M/+zQC2SVLUil4MkSdPBbwxLUsUMAUmqmCEgSRUr\n9j2BYXR5jaGI2AL8AfB04Djgmsz8YtGixiwijgf+GbgyMz9Rup5xiYiXAZ8F/gnoAV/PzHeXrWq8\nIuJ84L3AY8BlmXlL4ZLGJiLeClxA/5dSesCPZ+ZTylY1PhGxDfgEsBPYBFyVmbeutf/UhsDqNYYi\nYifwd0BnQgB4DfC3mflbEfFM4DagUyEAXAZ09deA78zMXyhdRBsG77fLgRcCO+h/p6czIZCZ1wHX\nAUTESzm8jE1XXAjcnZm/HhEnAXcAP7bWzlMbAvTXGPqbweP/X2OoK0tMZOaNq/74TOA7pWppQ0QE\ncCrdC7YVXV7v6hXAbZm5BCwBFxWup02XA+eVLmLMHgSeP3i8E9i33s5TGwK1rDEUEX8FnAy8unQt\nY7YHeCf9TyVd9NyIuJn+m+yqzLy9dEFjdAqwLSK+AJxA/3LeHWVLGr+I+AngPzLze6VrGafM/ExE\nXBgR36B//l613v5Tf2N41RpD7ypdSxsy88XAa4FPla5lXCLiAuCuzLx3sKlrn5q/AbwvM19HP+R+\nPyKm9gPVBvToh9vr6L/3ri9bTmveRv++XKcM7ufcm5nPBl4OfGS9/ac6BAZrDP0q8NOZub90PeMU\nEadFxDMAMvMfgCdHxFMLlzUurwJeGxFfpf9G+42IOKtwTWOTmfdn5mcHj78F/Bf92VxXfJd+iDeD\n8e3v0P+bq50B3FW6iBa8GPgSQGZ+HTg5Itb8t35qP70M1hj6AP01hrq4yNxLgV3AeyLi6cC2rqyl\nlJlvXHkcEVcA3+7S5YSIOA84KTP3RMQPAU8D7itc1jjdClw/+O28E+nQ/5srBjdM92fm/5aupQX3\nAKcDn4+IXfTHeWitnac2BOj+GkPX0r+M8OfA8cAlhevR8P4Y+PTgUuUm4KIu/WOSmfdHxB8Bf03/\nfdfFS7EnAZ26F7DKXuC6iLiT/q/Xv2O9nV07SJIqNtX3BCRJ7TIEJKlihoAkVcwQkKSKGQKSVDFD\nQJIqZghIUsUMAUmqmCEgjSgi3hMRewePIyL+ddDIQ5o5hoA0ug8Bz4mIF9FfofHtmfn9wjVJG2II\nSCMa9LV4G3Aj8I+Z+ZeFS5I2zBCQNuZEYD/9lWClmWUISCOKiOPprwL7GuDRiHhT4ZKkDTMEpNFd\nCXwuM+8Bfhl4X0T8cOGapA1xKWlJqpgzAUmqmCEgSRUzBCSpYoaAJFXMEJCkihkCklQxQ0CSKmYI\nSFLF/g/f02tjDBPosgAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "d.plot.scatter(x='x', y='y')\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 27, "metadata": { "collapsed": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYYAAAEDCAYAAAAx/aOOAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAFsNJREFUeJzt3XuQpHV97/H3sBBhmWFxYbitAZZEvsSIlpxKHWtBwLBy\niRDkYpGSGKJEDcdTlUSOVcgJJ4pLgptARCwPORaQiFcSLmuKhaCiXCQkeuR4hS8hLLvCru4ss2zN\nZncJ7Pb5o3u0n3Gmu2fpX/f27PtVNUU/l3l+3x/b3Z95nt9zGarVakiSNGmPfhcgSdq1GAySpAqD\nQZJUYTBIkioMBklShcEgSarYs3QDEbE38EPgI5n5mab5S4GrgJeAuzNzWelaJEnt9WKP4QpgwzTz\nrwPOAU4ATo2IY3pQiySpjaLBEBEBHAPcNWX+YuC5zFybmTVgJXBKyVokSZ0pvcdwDfABYGjK/EOA\nsabp9cChhWuRJHWgWDBExDuBhzNzdWNWczhMDYohwHtzSNIuoOTg81uBxRFxFvAqYFtE/Dgz7wOe\npbqHsAhY126DtVqtNjQ0NVO0u3jiiSd454c+z/wFBxVrY8um9dzyF+/g6KOPLtbGXOmHBsqsvjiL\nBUNm/s7k64j4M2BVIxTIzNURMRIRhwNrgTOBd7Tb5tDQEGNjE6VK7rvR0RH718L4+GbmLziI4Vcu\n6mJV07ezM3V22r9dvR/T8b052EZHR2a1fvHTVRtqABFxEfB8Zq4ALgG+2Fj2hcx8ske1SJJa6Ekw\nZOaV08x7CFjSi/YlSZ3zymdJUoXBIEmqMBgkSRUGgySpwmCQJFX06nRVSZpTtm/fztNPP9XVbR55\n5FHMmzevq9vcGQaDJO2Ep59+ij/6yy937Qr2LZvWc90Hf5tf+ZVXd2V7L4fBIEk7qRdXsDd773t/\nnw9/+CoOO2wRY2PrueyyS7nxxlu63o5jDJI0IE4//a187Wv3AvDQQw/wlrecVqQdg0GSBsTSpafy\nwANfB+Dhhx9k6VKDQZJ2a/vtt4DR0YN5/PEfUavVOPDA0SLtGAySNEBOO+0MrrnmY5x8crmHXjr4\nLEk7acum9T3f1vHHn8jy5VcZDJK0qznyyKO47oO/3fVttvO97/0/jj/+RIaHh7vadjODQZJ2wrx5\n83p+zcGNN/4N3/rWv7Bs2fKi7RgMkjQgLr74fVx88fuKt+PgsySpwmCQJFUYDJKkiqJjDBGxD/C3\nwMHAK4BlmXlX0/JVwBpgB1ADLszMdSVrkiS1Vnrw+SzgW5n5VxFxOPAV4K6m5TXg9MzcWrgOSVKH\nigZDZt7aNHk48OMpqww1fiRJu4ienK4aEd8EFgFnTrP4hohYDDyYmZf3oh5J0sx6MvicmccDZwOf\nm7LoCuADwEnAsRFxbi/qkSTNrPTg83HA+sx8JjO/GxF7RsSBmbkBIDM/27TuSuBY4PZW2xwdHSlZ\nct/Zv5lt3FjuFgDNFi4c3uk6O/m9QejHdHxv7j5KH0o6ETgC+JOIOBjYdzIUImI/4FbgrMx8kfpe\nw9+32+DY2ETBcvtrdHTE/rUwPr65i9W0bmdn6uy0f7t6P6bje3OwzTb0SgfDDcCNEfEAsDfw/oi4\nCHg+M1dExF3AIxGxBXg0M28rXI8kqY3SZyVtAy5ssfx64PqSNUiSZscrnyVJFQaDJKnCYJAkVRgM\nkqQKg0GSVGEwSJIqDAZJUoXBIEmqMBgkSRUGgySpwmCQJFUYDJKkCoNBklRhMEiSKgwGSVKFwSBJ\nqjAYJEkVBoMkqcJgkCRVFH3mc0TsA/wtcDDwCmBZZt7VtHwpcBXwEnB3Zi4rWY8kqb3SewxnAd/K\nzJOBC4Brpyy/DjgHOAE4NSKOKVyPJKmNonsMmXlr0+ThwI8nJyJiMfBcZq5tTK8ETgEeL1mTJKm1\nosEwKSK+CSwCzmyafQgw1jS9HjiqF/VIkmbWk2DIzOMj4vXA54DXN2YPTVltCKi129bo6EiXq9u1\n2L+Zbdw43MVKZrZw4fBO19nJ7w1CP6bje3P3UXrw+ThgfWY+k5nfjYg9I+LAzNwAPAsc2rT6ImBd\nu22OjU0Uqrb/RkdH7F8L4+Obu1hN63Z2ps5O+7er92M6vjcH22xDr/Tg84nApQARcTCwbyMUyMzV\nwEhEHB4Re1I/zHRv4XokSW2UPpR0A3BjRDwA7A28PyIuAp7PzBXAJcAXqR9C+kJmPlm4HklSG6XP\nStoGXNhi+UPAkpI1SJJmxyufJUkVBoMkqcJgkCRVGAySpAqDQZJUYTBIkioMBklShcEgSaowGCRJ\nFQaDJKnCYJAkVRgMkqQKg0GSVGEwSJIqDAZJUoXBIEmqMBgkSRUGgySpwmCQJFUUfeYzQEQsB04A\n5gFXZ+YdTctWAWuAHUANuDAz15WuSZI0s6LBEBEnA6/JzCURsRB4FLijaZUacHpmbi1ZhySpc6UP\nJd0PvL3xeiMwPyKGmpYPNX4kSbuIonsMmVkDJvcG3gOsbMxrdkNELAYezMzLS9YjSWqv+BgDQESc\nDbwLOHXKoiuAe4BxYEVEnJuZt7fa1ujoSJkidxH2b2YbNw53sZKZLVw4vNN1dvJ7g9CP6fje3H30\nYvD5NOBDwGmZOdG8LDM/27TeSuBYoGUwjI1NtFo80EZHR+xfC+Pjm7tYTet2dqbOTvu3q/djOr43\nB9tsQ6/oGENE7AcsB87MzE1Tl0XEPRGxV2PWScAPStYjSWqv9B7DBcABwK2NQecacB/w/cxcERF3\nAY9ExBbg0cy8rXA9kqQ2OgqGiBiaZtC4rcz8NPDpFsuvB66f7XYlSeV0eihpTUQsi4ijilYjSeq7\nTg8l/QZwPnBTRLwI3AzclpkvFKtMktQXHe0xZOZPMvOTmXkycEnjZ21jL2LvkgVKknqr47OSIuJN\nEXEjcDfwTer3P3oe+PtCtUmS+qDTwecngaeB/wP8YWa+2Fj0WES8rVBtkqQ+6HSM4QyAzPw3gIh4\nQ2Y+2lj2phKFSZL6o9NDSRdRv3p50mURcTX87H5IkqQ5otNgeHNmvntyIjMvoD7GIEmaYzoNhl+K\niF+anIiIYWCvFutLkgZUp2MMN1AfaP429Sex/Qbw4VJFSZL6p6NgyMwbI+Ir1AOhBvxJZv64aGWS\npL7o6FBS4yK2NwD7AfsDb4mId7f+LUnSIOr0UNI/AduB1U3zasBNXa9IktRXnQbDXpl5UtFKJEm7\nhE7PSvphRBxQtBJJ0i6h0z2GVwFPRsRjwEuTMzPzxCJVSZL6ptNguLpoFZKkXUant92+HxgGjm28\nfgZ4oGRhkqT+6PR01Y8BFwPvasx6B/CJUkVJkvqn00NJJ2XmGyPi6wCZ+dGIeLiTX4yI5dTvqzQP\nuDoz72hathS4ivq4xd2ZuWxW1UuSuq7Ts5K2Nv5bA4iIedS/6FuKiJOB12TmEuq37v74lFWuA86h\nHhynRsQxHdYjSSqk02B4OCJuBg6LiA8A32j8tHM/8PbG643A/IgYAoiIxcBzmbm2cevulcAps6hd\nklRAp/dK+p8RcT6whfqpq3+dmbd38Hs1fr638R5gZdPzGw4BxppWXw8c1Wnhc8H27dt5+umnfja9\nceMw4+Obu94GDDFvXsdPcS3Wxsvt35o1q9uv9DLVduzY6XY67V8v+iG9HJ0+2vMo4DuNn5/Ny8yn\nZv6tyu+fTX3g+tSm2UNTVhuicaiqldHRkU6aHAhPPPEEf/SXX2b+goOKtfHcM4+xz8gBA9/GZDsH\nvOrXiraxdWKMa760gfkL1hVroxf9AFi4cLirn5e59Nmbzlzv32x0Ovj8NX7+pf0K4CDgB9RvrNdS\nRJxG/elvp2XmRNOiZ4FDm6YXAW0/jWNjE+1WGRjj45uZv+Aghl+5qFgbWzb9dE60MdlOL/Ti/1cv\njI9v7trnZXR0ZE599qbaHfo3G50eSlrcPB0Rv0799NWWImI/YDlwSmZumrLN1RExEhGHA2uBM6mf\nBitJ6qNO9xgqMvOHEfFfOlj1AuAA4NbGoHMNuA/4fmauAC4BvtiY/4XMfHJn6pEkdU+nYwxXTpn1\ny9Sfy9BSZn4a+HSL5Q8BSzqpQZLUG52eqrK96ecl4LvAb5UqSpLUP50eSvrodDMjYg+AzNzRtYok\nSX3VaTBsY/ornSfHDdpeBS1JGgydBsNHgB8B91IPgrOAozNz2j0JSdLg6jQYfjMzr2qa/lJE3McM\nh5gkSYOr02A4ICJ+i/q9jwBOBA4sU5IkqZ86DYb3AtdQv+YA6lc9/7ciFUmS+qrTK5//FXhTRAw1\n3QRPkjQHdfoEt9dHxLeBxxrTfxoR/7VoZZKkvuj0ArdPAu/m5ze5uxW4tkhFkqS+6jQYXszM701O\nZOYT1K+AliTNMZ0Gw0uNJ65NPtrzDH7xeQqSpDmg07OSLgVWABERm4Cngd8rVZQkqX86DYYNmfm6\niBgFtk154I4kaQ7pNBg+R/3q57G2a0qSBlqnwZAR8RngYeA/fzYz86YiVUmS+qbl4HNEvK7xcm/q\nz2J4K/Cmxs8JZUuTJPVDuz2Gj1M/hPQugIi4LzPPKl+WJKlf2p2u6impkrSbabfHMPW+SLMOioh4\nLXAncG1mfmrKslXAGmBHo60LM3PdL25FktQrnQ4+T5rVDfQiYj7wCeCrLbZ3emZunWUdkqRC2gXD\nkohY0zR9UGN6CKhl5uFtfn8bcAZw2QzLh/BwlSTtUtoFQ7ycjWfmDuCFiJabuaFxu40HM/Pyl9Oe\nJOnlaxkMmbm6cPtXAPcA48CKiDg3M29v9QujoyOFS+qdjRuH+12C5rCFC4e7+nmZS5+96cz1/s3G\nbMcYuiozPzv5OiJWAscCLYNhbGzu3I1jfHxzv0vQHDY+vrlrn5fR0ZE59dmbanfo32x0enfVbqiM\nJUTEfhFxT0Ts1Zh1EvVHhkqS+qjoHkNEHEf9WdFHAC9GxHnAl4FVmbkiIu4CHomILcCjmXlbyXok\nSe0VDYbM/A7w5hbLrweuL1mDJGl2enkoSZI0AAwGSVKFwSBJqjAYJEkVBoMkqcJgkCRVGAySpAqD\nQZJUYTBIkioMBklShcEgSaowGCRJFQaDJKnCYJAkVRgMkqQKg0GSVGEwSJIqDAZJUoXBIEmqKPrM\nZ4CIeC1wJ3BtZn5qyrKlwFXAS8DdmbmsdD2SpNaK7jFExHzgE8BXZ1jlOuAc4ATg1Ig4pmQ9kqT2\nSh9K2gacAaybuiAiFgPPZebazKwBK4FTCtcjSWqj6KGkzNwBvBAR0y0+BBhrml4PHFWyntmo1Wqs\nWvXv1Gq1Ym2sWbO62La1e6vt2NHV99fGjcOMj2+uzNu+fTswxLx55f6+7EUbAAsXvr7o9gdN8TGG\nFoammW77LTw6OlKmmikmJia49GOfZ6/9FxdrY8umnzJ/wcHFtq/d19aJMa750gbmL/iFnfWuee6Z\nx9hn5ADmLzhooNvYsmk9t/zFMEcffXSxNgZNP4PhWeDQpulFTHPIaaqxsYliBTXbvHmCfUYOZO9X\nLupJe1K3zV9wEMMF37/1P2wGv41Jvfpu6YfZ/kHdy9NVK3sImbkaGImIwyNiT+BM4N4e1iNJmkbR\nPYaIOA64BjgCeDEizgO+DKzKzBXAJcAXqR9C+kJmPlmyHklSe6UHn78DvLnF8oeAJSVrkCTNjlc+\nS5IqDAZJUoXBIEmqMBgkSRUGgySpwmCQJFUYDJKkCoNBklRhMEiSKgwGSVKFwSBJqjAYJEkVBoMk\nqcJgkCRVGAySpAqDQZJUYTBIkioMBklShcEgSaoo+sxngIi4FngjsAP448z8dtOyVcCaxrIacGFm\nritdkyRpZkWDISJOBH41M5dExDHATcCSplVqwOmZubVkHZKkzpU+lHQKcCdAZj4O7B8Rw03Lhxo/\nkqRdROlgOAQYa5re0JjX7IaIeDAi/rxwLZKkDpQeY5i6NzBE/fDRpCuAe4BxYEVEnJuZt7fa4Ojo\nSHcrnMHee8MeQ+7MSLuLXn23DILSwfAs1T2Ew4CfTE5k5mcnX0fESuBYoGUwjI1NdLnE6W3ePMGO\nWq39ipLmhF59t/TDbEOv9KGke4HzASLiDcCzmfkfjen9IuKeiNirse5JwA8K1yNJaqPoHkNm/nNE\n/N+I+CawHXh/RFwEPJ+ZKyLiLuCRiNgCPJqZt5WsR5LUXvHrGDLz8imzvt+07Hrg+tI1SJI655XP\nkqQKg0GSVGEwSJIqDAZJUoXBIEmqMBgkSRUGgySpwmCQJFUYDJKkCoNBklRhMEiSKgwGSVKFwSBJ\nqjAYJEkVBoMkqcJgkCRVGAySpAqDQZJUUfzRnhFxLfBGYAfwx5n57aZlS4GrgJeAuzNzWel6JEmt\nFd1jiIgTgV/NzCXAHwCfmLLKdcA5wAnAqRFxTMl6JEntlT6UdApwJ0BmPg7sHxHDABGxGHguM9dm\nZg1Y2VhfktRHpYPhEGCsaXpDY950y9YDhxauR5LURukxhqFppmsdLOu7oaE9eGnT0+wYerFYGzs2\nbWDbHvsX2z7A1olxfvF/9eC10at2bGP3a2PLpvVFtz+ISgfDs/x8DwHgMOAnTcua9xAWAevabG9o\ndHSke9W1MDo6wtfv/N89aUuSdiWlDyXdC5wPEBFvAJ7NzP8AyMzVwEhEHB4RewJnNtaXJPXRUK1W\n9uhNRPw5cBKwHXg/cBzwfGauiIgTgOXUDyH9Q2b+ddFiJEltFQ8GSdJg8cpnSVKFwSBJqjAYJEkV\nxe+V1C0RsZz6rTPmAVdn5h19LqmrImJv4IfARzLzM/2up5si4kLgg8CLwBWZeU+fS+qaiNgX+Ayw\nENgLuDIz58TZdRHxWup3Lrg2Mz8VEa8CbqH+B+U64J2ZWe5Cn4Km6dsvAzdR/zf8T+B3M3NgL3CY\n2r+m+adRvy9dy52CgdhjiIiTgdc07rl0BvDx/lZUxBXUrwyfUyJiIfC/gCXUT0l+W38r6rrfBx7P\nzDcDb6d+/6+BFxHzqd/b7KtNs68Ers/Mk4B/B97dj9perhn69lHghsw8mfoX6qV9KK0rZugfEfEK\n4DJgbbttDEQwAPdT/9ABbATmR0T5S217JCICOAa4q9+1FLAU+EpmbsnMn2bmH/a7oC7bABzQeL2Q\n6m1eBtk26n+ENV90ejLwj43X/0j933YQTde3S4DbG6/HqP9bDqrp+gdwOfBJ6ntELQ1EMGRmLTO3\nNibfA6xs3HhvrrgG+AC9uK9E7x0J7BsRKyLi/oj4zX4X1E2Z+SXgiIj4N+AbwP/ob0XdkZk7MvOF\nKbP3bTp0NLD3Npuub5m5NTNrEbEH9eutPt+f6l6+6foXEUcDr8vM2+jge2YggmFSRJwNvAv47/2u\npVsi4p3Aw40rwWHuhcMQ9b++3kb93+7m/pbTXY3xk9WZ+Wrqdwf+ZJ9LKqn5j7Fd6t5m3dAIhVuA\nr2Xm1/tdT5ddS/2Pz44MTDA0Bk0+BJyemRP9rqeL3gqcHRH/TP2ZFX86x/6q/in14Ktl5lPAREQc\n2O+iuuh44J8AMvN7wKLGF8xctLlxnBo6u7fZoLkZyMz8aL8L6aaIOAwI4HON75lDI6Jl8A3EWUkR\nsR/1W2eckpmb+l1PN2Xm70y+jog/A1Zl5n19LKnb7gVubpxVdgD1wxFzaZD9SepPKLwjIo4AJjJz\nR59rKuWrwHnUD7OcB8yls8suBF7IzCv7XUuXDWXmWuDVkzMiYlXjZIkZDUQwABdQ/1K5tTHoXAN+\nLzOf6W9Zaicz10bEPwCPUP93mzOHARv+BrgpIr5B/VTq9/W3nO6IiOOoj30dAbwYEecDFwJ/FxHv\nA1YDf9fHEnfaDH07CNjW+Eu6BvwoMwfyvTpN/84Dzs3M5xurtD0E6L2SJEkVc/VYqCRpJxkMkqQK\ng0GSVGEwSJIqDAZJUoXBIEmqMBgkSRUGgySp4v8D+N7Zzb9LQj4AAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "d.plot.hist(x='x', y='y')\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 38, "metadata": { "collapsed": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYwAAAEDCAYAAAAm3zNHAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl8nPV96PvPaN8lSxrJlix5k/X1vmOMbWzqBUMgQBqy\nNBxKk9CFQ+893JzTmyavmzbtSU5ackNPcnryaqBp2qYpBNIb9sWAMeAV77u+trxLsjZLli1Z22jm\n/jEjJIQtjWw9ejQz3/fr5RfW88wzz9c/ntF3frsnEAhgjDHGDCXO7QCMMcZEBksYxhhjwmIJwxhj\nTFgsYRhjjAmLJQxjjDFhsYRhjDEmLAlO30BEngKWAX7gCVXd3e9cMvA0MFNVl/Y7/hDwZ0A38B1V\nfdPpOI0xxgzO0RqGiKwCylR1OfAo8JMBL/khsHfANbnAXwDLgXuBB5yM0RhjTHicbpJaC7wIoKoV\nQI6IZPQ7/63e8/2sA95W1auqWqeqf+JwjMYYY8LgdMIYDzT0+7kxdAwAVW27xjWTgXQReUlE3heR\nNc6GaIwxJhxO92F4rvHzUGuReIBcgk1RU4D3gEkjH5oxxpjhcDphVNOvRgEUAbVDXFMHbFPVAHBK\nRK6ISL6qNl7vgkAgEPB4BuYmY4wxQxjWL06nE8ZG4LvAMyKyEKi+RjOUh08GvRH4hYg8CeQB6YMl\nCwCPx0NDw5WRizqCeb2ZVhYhVhZ9rCz6WFn08Xozh/V6RxOGqm4XkT0ishXoAR4XkUeAS6r6kog8\nD5QA5SKyCXhaVZ8Tkf8AdhBsvvpTJ2M0xhgTHk+ULG8esG8MQfbtqY+VRR8riz5WFn283sxhNUnZ\nTG9jjDFhsYRhjDEmLJYwjDHGhMUShjHGmLBYwjDGGBMWSxjGGGPCYgnDGGNMWCxhGGNMhPujP/oD\namqqAWhoqOfrX3/Ykfs4voGSMcbEkuc3VbKron5E3/OWGQV8cU3Zdc/fddc9vPvuRh5++Kts2fIB\n69dvGNH797IahjHGRLh16+7kgw/eA2Dbtg9Zt86ZhGE1DGOMGUFfXFM2aG3ACVlZ2Xi9hVRUHCUQ\nCJCf73XkPlbDMMaYKLBhw9386Ed/yx13rHXsHpYwjDEmCqxYsYqamipLGMYYYwZ38OB+VqxYRUZG\nhmP3sD4MY4yJcD//+c/YtWsn3/vek47ex/bDiDK21n8fK4s+VhZ9rCz62H4YxhhjHGEJwxhjTFgs\nYRhjjAmL453eIvIUsAzwA0+o6u5+55KBp4GZqrp0wHUpwBHgr1T1X52O0xhjzOAcrWGIyCqgTFWX\nA48CPxnwkh8Ce69z+XeARgfDM8YYMwxO1zDWAi8CqGqFiOSISIaqtobOfwvIBx7qf5GICDADeM3h\n+IwZcYFAgLN1VzhR1UJ9czsEIC87hfKSHCZPyCTOM6yBKcaMGU4njPHA7n4/N4aOVQKoapuI5F/j\nuh8BjwN/4HB8xowYfyDAR0freHX7WWoa2675msLcNO5bPpllswvxWOIwEcbphDHwE+EBBp34ISIP\nA9tU9WywovGp9zBmzLnY0sHTrxzhRFUL8XEebplRwILp+RTlpQNQ13yVA5WNfHSsnmdePcr2I7V8\n7Z6Z5GQkuxy5MeFzOmFUE6xR9CoCaoe45h5gioh8FpgIdIjIeVXdNNhFXm/mTQUaTaws+oxGWRyq\nbOR//PMuWtu7uW3uBL722dmMDyWK/u5ZVUZ901X+928OsFfr+cGv9vLdR5dROj7L8RjBnov+rCxu\njNMJYyPwXeAZEVkIVKvqwLq6h361CFX9cu/fReQvgdNDJQvAZm6G2CzWPqNRFnuPN/APLx0hEAjw\nyF3CqvlFePz+697XAzz+wGxe3Z7Bbz84xTf/fgt//tAiivI/nWBGkj0Xfaws+gw3cTo6SkpVtwN7\nRGQr8GPgcRF5RETuBxCR54FngXIR2SQiXx7k7YwZUw6dushPf3uY+DgPT3xhPqsXFIfVL+HxePjs\n8sk8cpfQ2t7Nj369n6bLHaMQsTE3x9aSijL27amPk2VxpvYyf/urffT4A/y3Ly+gvCTnht7n9R1n\n+c3mk0wtyuLPH1pEQrwz3+HsuehjZdHH1pIyxmEtrZ38+IWDdHX38Mf3zbrhZAFw962lLJtdyKma\nyzy/qXIEozRm5FnCMGYY/P4AP3v5CC1tXXzhd8pYLAU39X4ej4dHNsygKD+dd/ZUUXG2eYQiNWbk\nWcIwZhhe2nKainOXWDg9nw1LS0bkPZOT4vn6PTPxeOAXbxyjs6tnRN7XmJFmCcOYMJ2saeHV7WfI\nz07ha/fMHNGJd1MmZHHX0lIaLnXw8rbTI/a+xowkSxjGhKHb18MvXq8gEICvfWYm6SmJI36P+1dO\nIS8rmbd3nafhUvuIv78xN8sShjFheHnrGWoa21izqJgZk8Y5co+kxHg+f8c0fD0BXth80pF7GHMz\nLGEYM4Sq+lbe2HGO/OwUHrxjmqP3unVmIdOKsthdUc/JmhZH72XMcFnCMGYQgUCAX719HH8gwH+6\ns5yUJGcXR/B4PB8npZe3nHH0XsYMlyUMYwaxq6IePX+J+dPymDftWgsrjzwpHYeU5HDo1EVO1Vwe\nlXsaEw5LGMZcR2dXD7/eVElCvIcvr5s+qve+b8VkAF7ZaiOmzNhhCcOY63hj51mar3SyYWkphePS\nRvXeMyaNY/rEbA6cvMi5OlvGwowNljCMuYaWti7e+ug8WelJ3HPbpFG/v8fj4e5lwfu+s7tq1O9v\nzLVYwjDmGl7deobO7h7uWzHZ8Y7u65k3LY+CcansOFrL5bYuV2Iwpj9LGMYMUN98lc37qynISWXV\n/CLX4ojzeFi/pARfT4DN+6pdi8OYXpYwjBngtx+epscf4HdXT3VsufFwrZg7ntTkBDbtq8bX43c1\nFmMsYRjTz7m6K+w8WsekwkyWzLi5lWhHQkpSArfPm8Dlti72n2h0OxwT4yxhGNPPS1uCw1g/v3oq\ncSO4uODNuD3ULPbBgRqXIzGxzhKGMSHn6q6w70Qj04qymD0l1+1wPlacn05ZcTZHTjfR2GKLEhr3\nWMIwJuSVbWcAuG/llBFdunwk3D5/AgFgy8ELbodiYpjj4wVF5ClgGeAHnlDV3f3OJQNPAzNVdWm/\n408CK4F44G9U9bdOx2liW1V9K3u0gSkTMpkzhmoXvZbOKOTZd06w5dAF7lsxhbi4sZXQTGxwtIYh\nIquAMlVdDjwK/GTAS34I7B1wzR3ArNA1dwP/08kYjYF+tYsVY692AcFd+ZbNKqTpcidHzza5HY6J\nUU43Sa0FXgRQ1QogR0Qy+p3/Vu/5ft4HvhD6ezOQJiJj7xNsokZ1Yxu7K+qZND6TedPy3A7nupbP\nmQDAjiN1LkdiYpXTCWM80NDv58bQMQBUtW3gBaoaUNXenr0/BF5X1YCjUZqY9uq2MwQILvg3FmsX\nvaYVZ5GfncKe4w10dtu+32b0Od2HMfDT5wHC+uUvIvcDXwXuDOf1Xm/m8CKLYlYWfYYqi5rGVnYd\nq2NqUTbrbxubzVH9/c6SEl549wSn69u4fUHxsK6156KPlcWNcTphVNOvRgEUAbVDXSQiGwg2V21Q\n1bCW6mxosBU9IfhBsLIICqcsnn2zAn8A7rxlIo2NraMU2Y2bNyWXF4CN288wozgr7OvsuehjZdFn\nuInT6SapjcCDACKyEKi+RjOUh341ERHJAp4E7lVV26PSOKaltZMth2rx5qSwWLxuhxOW4vx0Sgoy\nOHTqIq3t3W6HY2KMozUMVd0uIntEZCvQAzwuIo8Al1T1JRF5HigBykVkE8EhtplAHvB8qLM7APy+\nqtoaz2ZEvbOnCl+Pn7uWlhIfFzlTkpbNLuSF906yW+u5Y5jNUsbcDMfnYajqtwccOtTv3Bevc9kz\nzkVkDLR3+ti0t5qstERWzJ3gdjjDcuvMYML46GidJQwzqiLna5UxI2jz/mraO32sXVJCUmK82+EM\nS25WCtOKs9Dzl7h81fbJMKPHEoaJOd0+Pxt3nSc5KZ41iyLzG/ri8gICAdh3vGHoFxszQixhmJiz\n/UgtLa1drJ5fRHpKotvh3JDeTvo9agnDjB5LGCam+AMB3tx5jvg4D3feUuJ2ODfMm5PKpMJMjp1t\npq3DRkuZ0WEJw8SUfccbqW26yrLZheRmpbgdzk1ZLF56/AHbWMmMGksYJmYEAgHe2HkWgLtvneRy\nNDfPmqXMaLOEYWLG8fOXOFVzmQVl+RTlp7sdzk2bkJdOsTedw6ebaO/0uR2OiQGWMEzMeHPnOQA+\nsyzyaxe9Fpd78fX4OXTqotuhmBhgCcPEhAsX2zhw8iLTirMom5jtdjgjZuH0YLPU/krrxzDOs4Rh\nYsLbu4Mry2y4pdTlSEZWaWEG4zKTOXTyIj1+v9vhmChnCcNEvStXu9h26AL52SksLM93O5wR5fF4\nWFCWT1uHj8oqW6vTOMsShol6m/fX0OXzs25JSUQtMhiu+WXBJGjNUsZp0ffpMaafbp+fTXuqSE2O\n5/Z5kbXIYLhmTsohOTGe/ZXW8W2cZQnDRLWdR+toaeti1fwiUpMdX5zZFYkJ8cyekktd01UuXPzU\nrsfGjBhLGCZqBQIBNu46R5zHw7rFkbsMSDjml+UBcMBqGcZBljBM1DpwooGqhjaWzPCSlx3Zy4AM\nZf60fDxYP4ZxliUME7VefP8kAHdG2VDaa8lKT2JqcRaVVS22datxjCUME5WqG9vYU1HP9InZTC3K\ncjucUbGgLB9/IMChk9YsZZxhCcNEpbd3nQdio3bRa4ENrzUOc3zYiIg8BSwD/MATqrq737lk4Glg\npqouDecaY4Zyua2LbYdrGZ+XxsLp0TVRbzBF+enkZ6dw5HQTPX5/VM45Me5y9IkSkVVAmaouBx4F\nfjLgJT8E9g7zGmMGtXlfNb4eP/fdPo24OI/b4Ywaj8fD3Kl5XO30cbL6stvhmCjk9FeQtcCLAKpa\nAeSISEa/89/qPT+Ma4y5Ll+Pn/f2VZOaHM/aCN5R70bNnRYcXmur1xonOJ0wxgP9d3dpDB0DQFWv\nNcto0GuMGcxuraelrYuVc4tIi9D9um/GzNJxJMR7rOPbOMLpPoyB7QEeIODANXi9mcMIK7rFclm8\nf2AfHg88uL4ciM2ymDstn33HG4hLSiAvO/Xj47FYFtdjZXFjnE4Y1XyydlAE1DpwDQ0NV4YdXDTy\nejNjtixOX7iMnm1m3rQ8EgPB7xixWBYzSnLYd7yBzbvOsWp+ERDbz8VAVhZ9hps4nW6S2gg8CCAi\nC4HqazRDefhkrSKca4z5lHf3BPe8WLd4osuRuMv6MYxTHK1hqOp2EdkjIluBHuBxEXkEuKSqL4nI\n80AJUC4im4CnVfU5Ednb/xonYzTR4XJbFx8dq6MwN41ZU3LdDsdVheNSKchJ5cjpJnw9fhLibXit\nGRmOz8NQ1W8POHSo37kvXueabzkalIk67++vxtcTYN3iicR5Ymco7bX0Dq99d28VlVUtzJg0zu2Q\nTJSwrx4m4vUOpU1Jimf5HBtQB9YsZZxhCcNEvL3HG7jU2sWKuROids+L4ZpRmkNiQhwHLWGYEWQJ\nw0S8d0Kd3WtjvLO7v6TEeGaUjqO6oY2myx1uh2OihCUME9HO1l6hsqqFOVNzGZ+b5nY4Y8q8ULOU\n1TLMSLGEYSLae/uqAVi7yGoXA82dGhwtZrO+zUixhGEiVnunj51H68jLSmHu1Dy3wxlzCsalUZib\nxtGzzXT7etwOx0QBSxgmYu04Uktndw+rFhTF1Kq0wzFvah6dXT0cPdXkdigmCoSVMETkByIy3elg\njAlXIBBg8/4a4uM83D5vgtvhjFlzpwWbpXZX1LkciYkG4dYwmoF/F5HNIvKwiKQ4GZQxQzlVc5nz\n9a0smJ5PTkay2+GMWVKSQ1JiHHssYZgREFbCUNUnVfUW4GtAMfCuiPxURGY4Gp0x17F5f7Cz+44F\nxS5HMrYlJsQza1Iu5+taabzU7nY4JsINtw+jGCgDMoErwL+IyGMjHpUxg2jr6OajY/V4c1KYOdmW\nvRhK72gpG15rbla4fRh/ISLHgT8HXgYWquo3gduBP3EwPmM+ZdvhWrp9fu5YUBzz60aFo3eZkIM2\nvNbcpHDXURgPrFfVs70HRGSKqp4WkW86E5oxnxYIBHg/1Nm9Yq51docjPzuV0vGZHDvbTFd3D0mJ\n8W6HZCLUkDUMEYkDZgLnRCQu9CeZYE0DVX3T4RiN+diJqhZqGttYLF6y0pPcDidi3DKzkG6fn4pz\nzW6HYiLYoAlDRH4PqABWE9ybwhf60wacczw6Ywb44EANAKuts3tYlswsBOCANUuZmzBok5SqPgs8\nKyLfVdXvjk5Ixlxbe6eP3VpPfnYKUprjdjgRZebkXNKSEzhYeZHA+gAe6/sxN2DQhCEid6vqG8B5\nEfnawPOq+k+ORWbMALsq6unq9rNy7gTr7B6m+Pg45kzN5aNj9dQ0tlHszXA7JBOBhurDmBf670qC\nI6L6/1npYFzGfMqWQxfwAMvn2iZJN2KejZYyN2moJqm/Df33qyLiUdVAqMO7QFXPj0qExgC1TVep\nrGph5qRx5Genuh1ORJozNQ8PwX6Mu5dNcjscE4HCGlYrIt8CWkXkH4E9wGUR2aiqfxHGtU8BywA/\n8ISq7u53bh3wfYId6W+o6vdEJB34VyAXSAT+WlU3DvPfZaLM1kMXAGzdqJuQlZbElKIsKqtaaOvo\nJj0l0e2QTIQJd6b3Z4G/B74IvKKqywijSUpEVgFlqroceBT4yYCX/Bj4XOi91oeWGvkDoEJVfwf4\nQug1Job5/QG2Ha4lNTmBReVet8OJaPOm5eEPBDhy2lavNcMXbsLoVtUAcDfwYuhYOLN/1va+XlUr\ngBwRyYDgxD/goqrWhN77jdDrG4DezQ1yQz+bGHbkTBPNVzq5dWaBTTq7SfOn5QNwoNL6MczwhZsw\nLonIa8BMVd0uIvcSbGIayng++Qu/MXTsWufqgQmq+jwwSUROAJuB/xZmjCZKbTkYbI5aYc1RN620\nMIPsjCQOnbqI3x9wOxwTYcJdGuQrwHpga+jnTuCRMK4bOPbRAwQGOyciDwFnVfVuEZkH/COwdKgb\neb2ZYYQTG6KpLK5c7WLfiUZKCjO4dV7xsOcPRFNZ3Kzeslg6azxvf3SO5g4fMybluhyVO+y5uDHh\nJowegr/o7xWR3k9sCTDUPIxq+moUAEVAbb9z/b8yFgMXgBXAWwCqelBEikUkTlUHrdE0NFwJ598R\n9bzezKgqi/f2VuHr8bNsZiGNja3DujbayuJm9C+L8uJs3gbe332evLTY6/i256LPcBNnuE1SbwH/\nBVjF8OZhbAQeBBCRhUC1qrYBhBYyzBSRUhFJAO4Nvb6S4KgqRGQScGWoZGGi1/ajdXiAZbNt7sVI\nmTV5HPFxHg6ebHQ7FBNhwq1hJKrq6uG+eai/Y4+IbCVYS3lcRB4BLqnqS8BjwHMEay/PqmqliPwM\n+CcR2UywY/2Ph3tfEx0aL7V/PPdiXKbtqjdSUpMTkNIcjp5ppvlKp5WtCVu4CeOIiOSp6rCHVqjq\ntwccOtTv3BZg+YDXtwFfGu59TPTZcTS4reiyWYUuRxJ95k3L5+iZZg6dusiq+UVuh2MiRLgJYyJQ\nKSLHCE6yA0BVVzkSlYl5gUCAHUfrSIiPY7EUuB1O1Jk/LY/n3j3BgcpGSxgmbOEmjL9xNApjBjhf\n3/rxvhdpKeE+piZchblpFOamceRMk22qZMIWVqe3qr4PZABzQ3+vAj5wMjAT23Yc6W2Oss5upyya\nnk9Xt5+jZ2xTJROecPf0/lvg68BXQ4e+wqeX+TBmRPj9AXYeqyMtOeHjFVbNyFsYWmZl7wlbTMGE\nJ9xhtatV9XeBywCq+t+BxY5FZWKanr9E85VOlszwkpgQ7iNqhmtqURZZ6UkcqGy0Wd8mLOF+Gtv7\n/yAi8YS3lpQxw7bjSHBupzVHOSvO42Hh9HyuXO2msrrF7XBMBAg3YWwTkX8GJojIN4D3Ca7zZMyI\n6vb1sFsbGJeZTLltw+q4hdODzVL7rFnKhCHchPFL4HWCk+9WAD9S1W86FpWJWQdPNtHe6ePWWYW2\nDesomDlpHMlJ8ew73kggYM1SZnCDJgwRSRWR3xJMFg8SXP/pduDzIpI0CvGZGLOrIjg66taZNllv\nNCQmxDFvah71l9qpbmxzOxwzxg1Vw/gOwSRRrqpfVNU7gckE+zR+4HBsJsZ0dfdwoPIiBeNSKS3M\ncDucmLFwenCPjH3HrVnKDG6ohHE78A1V7T+7+yrwnwkud27MiDl0qonO7h6WSMGwlzE3N27etDzi\n4zzsPWGLEZrBDZUwfKraNfCgqnYDl5wJycSq3VoPwC0zbCmQ0ZSWksiM0hzO1l6h6XKH2+GYMWyo\nhDFYL5hvkHPGDEu3r4f9lY14c1KsOcoFH0/is2YpM4ihFulZLiLnrnHcA+Q7EI+JUYdPNdHZ1cOS\nhcPfVc/cvEXlXn618Ti7K+pZt6TE7XDMGDVUwpBRicLEvF2h5qgl1hzlipyMZMpLcjgemmVve2SY\naxk0YYR2xTPGUd0+P/tPNJKfncLk8bbXsluWzChAz19it9az3moZ5hpsoR7juiOnm+jostFRblsi\nXjzArop6t0MxY5QlDOO63l9Q1hzlruyMZKQ0h8qqFpqvdLodjhmDLGEYV3X7/OyvbCQvK5kpE6w5\nym29Q5p3Wy3DXIPjW5mJyFPAMsAPPKGqu/udWwd8n+AQ3TdU9Xuh4w8BfwZ0A99R1TedjtO44+iZ\n4NpRt8+bYM1RY8AiKeDf3j7Orop61t9i/RjmkxytYYjIKqBMVZcDj/LpTZd+DHwOWAncKSIzRCQX\n+AtgOXAv8ICTMRp39X6Ttcl6Y0N2ehIzSsdRWd1ik/jMpzjdJLUWeBFAVSuAHBHJABCRKcBFVa1R\n1QDBBQ7XAuuAt1X1qqrWqeqfOByjcYmvx8++E43kZiUztSjL7XBMyBJrljLX4XTCGA/0nzraGDp2\nrXP1QBHBxQ3TReQlEXlfRNY4HKNxydEzzVzt9LG43EZHjSWLy714PLDzmCUM80lO92EM/C3goW+5\nkWud84f+m0uwKWoK8B4waagbeb3WYdorUsri8LuVAKxfNtmxmCOlLEZDuGXh9cKC6V72HW+gGw9F\n3uhbqsWeixvjdMKopq9GAcEaRG2/cxP6nSsGLgBtwLZQM9UpEbkiIvmqOuhSmg0NV0Yu6gjm9WZG\nRFn4evxsP1TDuMxkctMTHIk5UspiNAy3LJaUBxPGax+e5IHbpzoY2eiz56LPcBOn001SGwluvISI\nLASqVbUNPp5FnikipSKSQLCDeyPwNrBGRDwikg+kD5UsTOSpONtMW4ePxeVe21lvDFpYnk9SYhw7\njtTZTnzmY44mDFXdDuwRka0ER0Q9LiKPiMj9oZc8BjxHcI/wZ1W1UlVrgN8AO4BXgT91MkbjDpus\nN7alJCWwuNxL/aV2TtZcdjscM0Y4Pg9DVb894NChfue2EBw+O/CaZ4BnHA7NuMTX42fv8QayM5Io\nm5jtdjjmOm6bM57tR+rYfriWsmL7/2RsprdxgZ67RFuHjyXlBdYcNYbNnDSO7PQkPjpWh6/H73Y4\nZgywhGFGXV9zlNflSMxg4uPiuHVWIW0dPg6dvOh2OGYMsIRhRlWPP9QclZ7E9Ik5bodjhnDb7OAg\nx+1Haod4pYkFljDMqNJzl2ht72aReImLs+aosa60MIPi/HT2VzbS2t7tdjjGZZYwzKj6eO0osdFR\nkcDj8bBy3gR8PQG2H7ZaRqyzhGFGTY/fz57jDWSlJVJeYs1RkeK2OeOJj/PwwcEam5MR4yxhmFFz\n/HwLV652s0gKrDkqgmSlJbGw3Et1QxunL9gM6VhmCcOMmr7mKBsdFWlWzQuu4vPBgRqXIzFusoRh\nRoXfH2DP8QYyUhMpL7XmqEgza3IueVnJ7DxWR0eXz+1wjEssYZhRcaLqEpfbulgsXuLj7LGLNHFx\nHlbOK6Kzq+fjeTQm9tgn14yKjyfr2eioiLVy7gQ8wIcHLrgdinGJJQzjOL8/wG4NNkfNmGTNUZEq\nLzuF2VNyqaxuobqxze1wjAssYRjH9TZHLSq35qhIt2p+EQCb91a7HIlxg316jeM+6h0dNdOaoyLd\nwvJ8xmUms/XwBdo7rfM71ljCMI7y+wPs6W2OstFRES8+Lo7VC4ro6Ophh60vFXMsYRhHHT9vo6Oi\nzer5RcTHedi0t9pmfscY+wQbR+1S21kv2mRnJLNYvFQ3tqHnLrkdjhlFljCMY6w5KnqtXTwRgE17\nq1yOxIwmSxjGMdYcFb3KirMpKchg7/FGmq90uh2OGSWOf4pF5CkR2SYiW0RkyYBz60Rkp4hsFZH/\nZ8C5FBE5KSK/73SMxhm9k/VuseaoqOPxeFizqBh/IMDmfTbENlY4mjBEZBVQpqrLgUeBnwx4yY+B\nzwErgTtFZEa/c98BGp2Mzzgn2BxVT0ZqImLNUVFp2azxpKcksHl/NV3dPW6HY0aB0zWMtcCLAKpa\nAeSISAaAiEwBLqpqjaoGgNdDryeUOGYArzkcn3GInr/E5avdLLHmqKiVnBTP6gXFXLnazY6jdW6H\nY0aB05/k8UBDv58bQ8euda4emBD6+/8LfAOwTRMi1MdrR1lzVFRbu3gi8XEe3vroHH4bYhv1Ehx+\n/4G/8D1AYLBzIvIwsE1Vz4rItV53TV5v5s3EGVXcLoueHj/7TzSSnZHEykUlxMe7V8NwuyzGEifK\nwuvNZNXCYt7bU8X5i+0smVk44vdwgj0XN8bphFFNX40CoAio7XduQr9zxcAF4DPAVBH5LDAR6BCR\n86q6abAbNTTYTmAQ/CC4XRaHT13kUmsnv7OomKYm9xapGwtlMVY4WRar503gvT1VPP+2Mik/zZF7\njCR7LvoMN3E6nTA2At8FnhGRhUC1qrYBhGoQmSJSCtQA9wJfUdWf9l4sIn8JnB4qWZixZfuRYHv2\nbbPGD/FKEw1KCzOZOWkcx842c67uCqWF9u09WjnaVqCq24E9IrKV4Iiox0XkERG5P/SSx4DngPeB\nZ1W10sndHieSAAAPJklEQVR4jPM6u3vYe6KB/OwUphVnuR2OGSUblpYAsHHXeZcjMU5yuoaBqn57\nwKFD/c5tAZYPcu1fORWXccaBykY6u3pYt3giHo+NWYgVc6bmMSEvjZ1H6/j86mmMy0x2OyTjABvv\naEbUjlBz1LLZ1hwVS+I8HjYsLaXHH2DjrnNuh2McYgnDjJjW9m4OnbpISUEGxfnpbodjRtlts8cz\nLjOZzftqaG3vdjsc4wBLGGbE7K6op8cfYNnsyBhaaUZWYkIcG5aW0tndwzu7rS8jGlnCMCNmx5Fa\nPMCtETIW34y81fOLyEhN5J3dVbYjXxSyhGFGxMWWDo5XtVBekkNuVorb4RiXJCfFs/6WEq52+mxR\nwihkCcOMiG2h7TqtOcqsXVRManI8b+06b4sSRhlLGOam+QMBthysISkhjqXWHBXz0lISWbNoIpfb\nuvjw4AW3wzEjyBKGuWknzl+i4VIHS2YUkJrs+NQeEwHWLykhKSGON3eexdfjdzscM0IsYZib1vst\n8vZ5E4Z4pYkVWelJrJpfxMXLnWw5ZLWMaGEJw9yU9k4fuyvqKchJpbzENkoyfT5z2yQSE+J4ddsZ\nun1Wy4gGljDMTfnoWB1dPj8r5k2wpUDMJ+RkJLNmUTFNlzv54ECN2+GYEWAJw9yULYcu4AFWzLGl\nQMyn3b1sEsmJ8by6/YyNmIoCljDMDatubONk9WVmT821uRfmmrLSkli3ZCItrV28Z/MyIp4lDHPD\n3ttbBQRn9xpzPRuWlpKaHM9r28/S0WWzvyOZJQxzQ9o7fWw7XMu4zGQWTM93OxwzhmWkJnLnLaW0\ntnfz7p4qt8MxN8EShrkhO47U0tHVw+oFRcTH2WNkBrd+SQnpKQm8vuOcrWQbweyTboYtEAiwaW81\n8XEea44yYUlLSeCzK6bQ3unjla1n3A7H3CBLGGbYjp+/RHVjG4vFS3aG7axmwrNmUTHenBQ27a2i\nrvmq2+GYG2AJwwzbO7uD7dBrFk10ORITSRLi4/j86mn0+AP8x+aTbodjboDjC/+IyFPAMsAPPKGq\nu/udWwd8H/ABb6jq90LHnwRWAvHA36jqb52O04Snrukqe483MGl8JtMnZrsdjokwt8woYOOu8+zW\nBiqrWygrtmcokjhawxCRVUCZqi4HHgV+MuAlPwY+RzA53CkiM0TkDmBW6Jq7gf/pZIxmeN7adZ4A\ncPetpTaz2wybx+PhS2vKAPj1phMEAgGXIzLD4XST1FrgRQBVrQByRCQDQESmABdVtUZVA8Drode/\nD3whdH0zkCYi9ptpDLjc1sWWgxfIz05hsXjdDsdEqOkTc1hU7uVk9WV2Hq1zOxwzDE4njPFAQ7+f\nG0PHrnWuHpigqgFVbQ8d+0Pg9VBCMS57d08Vvh4/G5aW2lBac1O+tKaMxIQ4fr2p0rZyjSBO92EM\nrBl4gEAY5xCR+4GvAneGcyOvN/MGQ4w+TpTF1Y5uNu+vJjMtiQfWTCclKTL2vbDnos9YKguvN5Mv\nrC3n39+qYOOeah69f86o398Mn9Of+mr6ahQARUBtv3P9N1AoBi4AiMgG4FvABlW9Es6NGhrCelnU\n83ozHSmLV7ad4crVbh64fQpXWtqJhNJ2qiwi0Vgsi9VzC3ln51le+fAUi8vymFiQMSr3HYtl4Zbh\nJk6n2xU2Ag8CiMhCoFpV2wBU9SyQKSKlIpIA3AtsFJEs4EngXlVtcTg+E4arHT7e2nmO9JQE1i8p\ncTscEyUSE+L5yvpy/IEA/7ZRrQM8AjiaMFR1O7BHRLYSHBH1uIg8EmpuAngMeI5gR/ezqloJfAnI\nA54XkfdEZJOI2IB/F729+zxXO33cdWupbcFqRtS8aXksnJ7P8aoW2/87Ajj+6VfVbw84dKjfuS3A\n8gGvfwZ4xum4THha27vZuOscmWmJrF1seduMvIfWl3PsbDO/3nSCOVNsqfyxzIa6mEG9tOU07Z09\n3LNsUsR0dJvIkpuVwpfWlNHe2cO/vmVNU2OZJQxzXdWNbby3t5qCcamssdqFcdCq+UXMmjyOgycv\nsu1w7dAXGFdYwjDX9etNJ/AHAnxpTRkJ8faoGOd4PB7+4K4ZJCfF8+/vnKDxUvvQF5lRZ78FzDXt\nP9HI4VNNzJo8jgVltkGScV5+TipfWTed9k4fP3v5CL4ev9shmQEsYZhPae/08cuNSnych99bO93W\njDKjZuXcCSybVcjJmsu8+OFpt8MxA1jCMJ/ym/dP0nylk3tum0Sxd3QmUxkDwaaphzcIBTmpvL7j\nLIdPXXQ7JNOPJQzzCRVnm3lvbzVF+encc9tkt8MxMSg1OYE/eWA2CfEefvbyEeqabLOlscIShvnY\nlatdPP3KEeI8Hr76mRkkJtjjYdwxeXwWv79hBm0dPn78m4Nc7bB9wMcC+41ggOA+3b94vYJLrV18\nbtUUphXZxjbGXSvnTeCupaXUNl3lH16yTvCxwBKGAeDVbWfYX9nIrMnjuHvZJLfDMQaAB++Yxrxp\neRw+3cTPXzuG32+T+txkCcOwu6Ke3354mrysZP7os7OJs1FRZoyIi/Pw2P1zKJuYzc6jdfzSFil0\nlSWMGHf8/CX+8dWjJCfF838+OJ+s9CS3QzLmE5KT4nniwXmUFmTw/v4afvmWWk3DJZYwYlhldQt/\n98IBevwBHrt/NiWjtB+BMcOVlpLIN768gJKCDDbvr+EfXj5Ct8/6NEabJYwYdexsM3/3/H66u/38\n8X2zmTfNZnObsS0rLYlvfmUR5SU57K6o56lf7+fy1S63w4opljBi0I4jtTz16/10dfv5o/tmsWRG\ngdshGROWtJQEvvHF+Swu96LnL/HX/7yLUzWX3Q4rZljCiCHdvh5+tfE4T79ylKTEeL7xpQUsnVno\ndljGDEtSYjyPfW4Ov7tqKs2XO/nBv+3h5S2nbdjtKLANDmLEyeoW/uVNpaqhleL8dB57YA5F+elu\nh2XMDYnzeLh3+WSmFGXxT68d48Utp9l7vIHfWzcdKR3ndnhRyxMlQ9QCtql70MAN7psud/Dy1jN8\neKCGAMF9B35v3XSSE+PdC3KUDCyLWBbNZXG1w8dzm06wJbTF68Lp+Txw+9TrDuKI5rIYLq83c1hj\n6B2vYYjIU8AywA88oaq7+51bB3wf8AFvqOr3hrrGhKe6sY1Ne6v48EANvp4ARfnpPHxnuX37MlEn\nLSWBr31mJncsKOb5TSfYd6KRfSeCk1DXLSlhzpRc289lhDiaMERkFVCmqstFZAbwT3xyD+8fA+uB\nC8D7IvIboGCIa8x11DdfZXtFPZs+OsfJUEegNyeFe2+bzG1zxtuHxkS1qUVZfPOhRRw8eZG3PjrH\n0TPNHD3TTEZqIktmFLBwej7lE3PcDjOiOV3DWAu8CKCqFSKSIyIZqtoqIlOAi6paAyAirwHrAO/1\nrnE41ogRCARo6/BR23SVs7VXOHPhMpXVLdQ1B3cp8wBzpuayal4RC6bnW6IwMcPj8TC/LJ/5Zfmc\nrb3ClkMX2FVRz+Z91WzeV018nAeZNI6J+emUFGRQWphJ4bhUkmKgiXYkOJ0wxgP9m5MaQ8cqQ/9t\n6HeuAZgG5A1yTVS5crWLltYuunv8dPv8+HqCf7p9frp7/HR29dDW4aO1vZvW9m7a2ru52NJBQ0sH\n7Z2+T7xXSlI8C6fns2JBMVMKMhiXmezSv8qYsWHS+Ewmjc/ky2vLOH7uEofPNFFxtpmKM00cPd30\niddmpSWSm5VCXlYK6amJpKckkJ6aSFpKAqlJCSQmxJEQH0diQuhPfBwJCXEkxHsoyEmNmU3GnE4Y\nA0vRAwSGODfYNVGjtb2b//q/t+LrGd4/LSkhjvycVMonZuMdl8qkwkymTMhifF4acR6PdegZM0B8\nXBwzJ+cyc3IuAOmZKew7Wsv5+laqGlqpb26n6XIHVQ1tnKkd/mfngZVTuG/llJEOe0xyOmFUE6wd\n9CoCavudm9DvXDFQA3QNcs31eLzezJuLdJR5gd8+eZ8z7x1hZeEkK4s+VhZ9ViwqcTuEiOR04/ZG\n4EEAEVkIVKtqG4CqngUyRaRURBKAe0Ovf/t61xhjjHGP4/MwROR/AKuBHuBxYBFwSVVfEpGVwJME\nm5x+o6p/d61rVPWQo0EaY4wZUrRM3DPGGOMwG29pjDEmLJYwjDHGhMUShjHGmLBExWq1IvII8N/p\nm9z3tqr+wMWQRp2tvxUkIquBF4DDBOfwHFTV/+JuVKNPROYQXDHhKVX9qYhMBH5J8EviBeBhVe12\nM8bRco2y+AWwmOCkYIAfquobrgU4ikTkSWAlEA/8DbCLYTwXUZEwQp5T1f/b7SDcEMaaXbFms6p+\n0e0g3CIiacBPgHf6Hf5r4H+p6v8nIt8Hvgb8zI34RtN1ygLgz1X1dRdCco2I3AHMCv2eyAX2Ae8C\nf6+q/xHOc2FNUtHhE2t2ATkiEssbdMfGOg3X1wHcTfAbY687gFdCf3+F4LptseBaZRGr3ge+EPp7\nM5BOcPrCy6FjQz4X0VTDuENEXgcSgT9T1f1uBzSKBluzKxbNEpEXgVzgr1V14LfLqKaqfqBTRPof\nTu/X1FDPJ1dZiFrXKQuAPxWR/wrUAX+qqk2fujjKqGoAaA/9+CjwGrBhOM9FxCUMEfk6wX9s77pT\nAeBZ4C9V9Q0RuQ34F2C+e1GOuphYfytMJ4DvquoLIjIVeE9Epqmqb6gLo1z/5yGWnw+AfyW4UvZB\nEfkm8FfA/+FyTKNGRO4n2PR0J8HPS68hn4uISxiq+nPg54Oc3y4iBSLiCWXUWDDYml0xJbRc/guh\nv58SkVqC65SddTUw97WKSLKqdhIsj5htolHV9/r9+DLwU7diGW0isgH4FsGaxRURGdZzERV9GCLy\nZyLy5dDf5wD1MZQsYJA1u2KNiHwl1NSAiIwnuCFXtbtRjQnvAJ8P/f3zwJsuxuIqEflNaD8eCPbt\nHHYxnFEjIlkEl2K6V1VbQoeH9VxExdIgIlJM39CweOD/irVhpbb+VlCos//fgRyC/VnfVdW33I1q\ndInIIuBHwCSgm2DCfIhgU20ywdrWV1W1x7UgR8l1yuJ/EfyW3Qa0EiyLxuu+SZQQkT8E/hI4Tl/z\n0yMEW2zCei6iImEYY4xxXlQ0SRljjHGeJQxjjDFhsYRhjDEmLJYwjDHGhMUShjHGmLBYwjDGGBMW\nSxjGGGPCYgnDGGNMWP5/Lt7yOJM2SfUAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "d.plot.kde(x='x', y='y')\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 19, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "7.3 5.8\n" ] } ], "source": [ "# yの値は0以上の離散値で、yの平均と分散がほぼ等しい(オーダー的に)\n", "print d.y.mean(), d.y.var()" ] }, { "cell_type": "code", "execution_count": 21, "metadata": { "collapsed": false }, "outputs": [], "source": [ "# statsmodelsを使ってglmを計算します\n", "import statsmodels.formula.api as smf\n", "import statsmodels.api as sm\n", "from scipy.stats.stats import pearsonr" ] }, { "cell_type": "code", "execution_count": 23, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "Intercept 1.566046\n", "x 0.083343\n", "dtype: float64" ] }, "execution_count": 23, "metadata": {}, "output_type": "execute_result" } ], "source": [ "fit = smf.glm('y ~ x', data=d, family=sm.families.Poisson(link=sm.families.links.log)).fit()\n", "fit.params" ] }, { "cell_type": "code", "execution_count": 29, "metadata": { "collapsed": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAESCAYAAAASQMmzAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAHRJJREFUeJzt3XuYZHV95/F3dU93V1dfq7pruocGZkDlOxrN88gm2X3E\nVQQiEmLUGIyZCYoY5aKsGyduCMtFYAzGMGjQdQeNYGKUFcXLs4+GIOKNoIKB3ayJ/BBBVObWl+rb\nVN/77B/n1Jnq63TXVNU51f15Pc88VV3X7/y6+vf5nd85dX4Jz/MQEREBqIu6ABERiQ+FgoiIhBQK\nIiISUiiIiEhIoSAiIiGFgoiIhCoeCmb2YjN7ysyuXHT7+WY2X+n3FxGRtatoKJhZCrgdeGDR7U3A\n1cCBSr6/iIisT6W3FCaBC4CDi26/BvgYMF3h9xcRkXWoaCg45+adc1PFt5nZGcCvO+fuBRKVfH8R\nEVmfLRG8523AVRG8r4iIHEeiGuc+MrMbgH7gK8B3gusJ4KXA951zr1rpuZ7neYmENihERNappI6z\nmlsKCefcAeAFhRvM7JnVAgEgkUjQ3z9W8eJOVDbbpjrLSHWWl+osn1qoEfw6S1HRUDCzM4F9wHZg\nxszeCPy+c244eIhO0SoiEiMVDQXn3GPAilsCzrnTK/n+IiKyPvpGs4iIhBQKIiISUiiIiEhIoSAi\nIiGFgoiIhBQKIiIRePzxf+Haa/98zY//9re/WcFqjlEoiIhEZK1na5idneXzn/9shavxKRRERKrg\nne+8hAMHngOgv/8IH/vYR5iYyHPzzdfx1re+mU9/+m8B+NGPHuHyyy/lqqsu45pr3sfs7Cy3376P\np5/+Gbfd9lcVrzOKE+KJiETqngef4tEnjpT03Pr6BHNzS0/G8Js7t/Kmc56/4vNe85oL+eY37+fi\ni9/GQw99l/POO5977/08n/vcvczNzXHRRb/HJZf8CWNjo7z//R+gt3cbe/fewCOP/IBdu97CT37y\nb7z3vWufbiqVthRERKrgvPNezXe/+y0AHn74e3R1dXPGGTtpbGykubk5fFxnZ5pbbrmZd7/7nTz+\n+L8wMjK80ktWhLYURGTTedM5z191VL+aUk+I197eQTbbwxNP/Due57F161bq6+uXPO6WW27i1ltv\n59RTt/PhD3+opBpPhLYURESq5PzzL2Dfvr/i7LPPZaVlC/L5o/T09DA2NsZjj/2I2dlZEokEs7Oz\nValRoSAiUiVnnfUKDhz4FWeffe6Kj3nDGy7i8ssv5dZb/5Ldu9/KZz5zF3V1dczMzHL99X9R8Rqr\nssjOCfJq5dzlqrN8VGd5qc7yOZEaH3vsR9x339e45pobylzVUtlsW+wX2RER2bQ+9ak7ePTRH7J3\nb/X3E6yHQkFEpAre/vbLePvbL4u6jOPSPgUREQkpFEREJKRQEBGRkEJBRERCCgUREQkpFEREJKRQ\nEBGRkEJBRERCCgUREQkpFEREJKRQEBGRkEJBRERCFT8hnpm9GPgKcJtz7uNmdgpwJ9AATAN/7Jwr\nbbFUEREpq4puKZhZCrgdeKDo5puB/c65s/HDYk8laxARkbWr9PTRJHABcLDotiuALwXX+4FMhWsQ\nEZE1quj0kXNuHpgys+LbJgDMrA54F3BjJWsQiUoul2PfHXeTy0O62WPPFbtId6ajLktkVZHsaA4C\n4TPAN51z34qiBpFK23fH3TznGROpnTzHTvbtvzvqkkSOK6qV1+4CnHPu5rU8OJttq3A55aE6y6vW\n6xyfriPR6C+Tm0gkGJ+ui/T/VOvtGSe1UGOpqh4KZrYbmHLO3bTW58R9IW+ojQXHQXWW22p1tjbM\nMeJ5JBIJPM+jtXE+sv/TRmjPuKiFGqH04KpoKJjZmcA+YDswY2Z/AGwFJs3sW4AH/Ltz7t2VrEMk\nCnuu2MW+/cE+hRTsufyPoi5J5LgqvaP5MeBVlXwPkbhKd6bZe/WVUZchsi76RrOIiIQUCiIiElIo\niIhISKEgIiIhhYKIiIQUCiIiElIoiIhISKEgIiIhhYKIiIQUCiIiElIoiIhISKEgIiKhqNZTEBE5\nIVrZrjK0pSAiNUkr21WGQkFEalIu769oB/5lLh9xQRuEQkFEalK62cPzPAA8zyOdirigDUKhICI1\nac8Vu+hLOJrzT9CXcFrZrky0o1lEapJWtqsMbSmIiEhIoSAiIiGFgoiIhBQKIiISUiiIiEhIoSAi\nIiGFgoiIhBQKIiISUiiIiEhIoSAiIiGFgoiIhCp+7iMzezHwFeA259zHzexk4DP4gXQQuNg5N1Pp\nOkRE5PgquqVgZingduCBoptvAj7qnHsl8DPg0krWICIia1fpLYVJ4ALg6qLbzgYuC67/b2APcEeF\n6xBZMy3zKJtZRbcUnHPzzrmpRTe3FE0XHQG2VbIGkfXSMo+ymUWxnoJXdD2x6OdlZbNtlaumjFRn\neUVV5/h0HYnGY8s8jk/XrVqL2rO8aqHOWqixVFGEwriZNQVbEH34O5tX1d8/VvmqTlA226Y6yyjK\nOlsb5hjxPBKJBJ7n0do4v2Itas/yqoU6a6FGKD24ojgk9QHgjcH1NwL3RVCDyIq0zKNsZhXdUjCz\nM4F9wHZgxsz+ANgN/J2ZXQY8C/xdJWsQWS8t8yibWUVDwTn3GPCqZe56dSXfV0RESqNvNIuISEih\nICIiIYWCiIiEFAoiIhJSKIiISEihICIiIYWCiIiEFAoiIhJSKIiISEihICIiIYWCiIiEFAoiIhKK\nYj0F2UDKtXTl4te55drL0MdTKq1Sn99aXsJVWwpyQsq1dOXi17nx1rvKXKnIUpX6/NbyEq4KBTkh\nuby/ZCX4l7l8eV5nYGyuXCWKrKhSn99SXycOFApyQtLNHp7nL7PteR7pVHlep7utvlwliqyoUp/f\nUl8nDhQKckLKtXTl4te54c/eVuZKRZaq1Oe3lpdwTRTSLca8WlkkW3WWj+osL9VZPrVQI0A225Yo\n5XnaUhARkZBCQUREQgoFEZENYmJqlmcOjvL9Hx8q+TX07SARkRrieR7D49McGjzKgcE8hwbzHBw6\nysHBPLmxqfBxv/eqF5T0+goFEZEYmp2bp394goODeQ4OHg0u8xwaOsrE1NLv8aTbmvi1HWl6u1rY\n1lX6MbEKBRGRCB2dnPFH+4N5Dg35AXBoKM+R3ARz8wuPDq2vS9CTSfGiHSm2daXYlmlhW3eKnnSK\n5qbydOcKBRGRCpuf9xgYmQg6/WMBcGjwKKP5mSWPb27awo7eNrYFo/7erhQndbXQ3Zmkvq6yu4IV\nCiIiZTIxNbtgtH8wmPM/nMszO7dw1J8AujuTvKS3Pez4t2VS9GZStLc0hqfNqDaFgojIOszNzXM4\n53f2h4YKI37/cuTo9JLHJxvrOTnbGnT8LX7H35WiJ91Mw5b4nc5FoSAisozxiZnwyJ7ijr9/eGLZ\nUX+mPcmvnZZhWyYVBkBvJkVna3Sj/lJUPRTMrAX4eyADNAA3Oefur3YdIiIzs/McGZ4IOny/8z88\n5M/9j08sP9d/el8HXW3JcLqnJ+OP+hsb4jfqL0UUWwqXAE845/67mW0DHgReGEEdIrIJzHseudEp\nDgVTPoeH8hzK+ZcDI5MsPv1bXSJBtjPJ805qp7fLn+Pvzfgj//ZUA1u3ttfEuY9KFUUoDAAvCa5n\ngP4IahCpCRtpRa9KG5+YCUb6+YWXuQlmZueXPL69pZEX9HXQkynu+FNkO5vZUh/NyR7i8Ps+biiY\n2fnA/c65spxO1Tn3eTO7xMx+CnQCF5bjdUU2osKKXolUgrznsW//3ey9+sqoy4rM1PQch3P+Mfxh\nx5/zp3yWm+5paqj35/eDTr8QAD3pFKlk/HapxuH3vZZWeQ/wP8zss8CdzrlnT+QNzWw38Kxz7gIz\n+3Xgb4HfWu052Wzbibxl1ajO8lKdMD5dR6Lx2Ipe49N1Jb9frbRnZ7qFw0NHOdB/lOf6xzkwcJQD\n/eMc6B9nYGRyyePr6hL0ZlLs3JHh5K2tnJRtpS/bQl+2lUx7siI7eSvVluX8fZfquKHgnPsdM0sD\nbwD+p5kB3AV8yTlXypqJZwH/FLz2v5pZn5nVOeeWbt8FamH+robOsa46y6jSdbY2zDHieSQSCTzP\no7VxvqT3i1t7zs97DI1OcrhoxH84N8HAyCSHh/LML7POS6a9iRduT4c7dguXK033zE/PMjAwXvba\nK9mW5fp9Q+nBtabtJ+dczszuBqaBK4E/A643s3c4536wzvd8CvhPwJfNbDswtlogiGxme67Yxb79\nwRxzippa0Wve8xgemwo7/MPBNM/h3PKHdQJ0tjZxel87vekUPZlmetL+dE823UzTBjm6ZzVx+H2v\nZZ/CfwbeBpwDfAm41Dn3hJntAL4MvHSd73kHcKeZfRuoBy5b5/NFNo10ZzrW+xAKZ+z0O35/rr8Q\nAEdW2MGbatrCKVtb2Zr2R/q94WGdKbafko7VFk21xeH3vZYthVuA/cAVzrnwvKzOuZ+b2T3rfUPn\n3FHgD9f7PBGJRqHjP7Ko0z88NMGR4TzTM0s7/mRjPSd1tdCTaQ47/8J0T2tzQ019mWuzWcs+hZev\nct8t5S1HRKJQmOo5kpvgyPAEh4eOBcBKHX9TQz296RRbgw5/a9qf7unJpGhPqeOvVfE7JktEKqLw\nJa4juTyHhyc4EszvHxmeoD83wfQyUz3FHX9hxF8IgY4IT9omlaNQENlAZufmGRyd9Ef8wb/+4QkG\nxyY5OJBndm75qZ5tXS1Bx++P9v3L5kjP1inRUCiI1JjpmTn6RybDOf4jw4UAyDM4MrXs4ZwtzQ2c\nnG1ZMOIvdP5tmuqRIgoFkRjKT86EnX1/0eXh3MSCdXiLtbc0cvpJ7eGIf2unHwBb082cdmpmUx/V\nI2unUBCJQOGInkKHf2R4Yee/3CkbEkC6vYmdp3aGI36/4/e/wFWu5Rhlc9OnSKRCCvP7/QumePxO\nv394+R279XUJujubOW1bO1s7m8kGI37/MhnLRVlkY1EoiJyA8fw0Pz80WtTZT4ad/uDo0tMyAzQ3\n1dPblVrQ6ReuZ9qS1NVpfl+io1AQWUXhHD39wxPBzt1gmmd4goHhCY5Ozi77vM7WRp7f10G2s3nJ\niL9NX96SGFMoyKaXn5wNR/f9I/5ofyDo+AdHJpmbXzrc31JfR7YzyQtP66Iz1UA26PCznc10dyQ3\nxXl6ZGNSKMiGNzs3z9DY1LGOv2iaZ7XRfluqgR29bX6HH/wr7NTtaG30V+iK2dlHRU6UQkFqnud5\njOVn6B+ZYGB4koGRhR3/0Ojyx+4XRvvP6+sg29FMtjMZdv7dnUmSjfrzkM1Hn/oqKtdSe3FYsq/a\nJqdnGRieDDv+/mH/3PuFn6dmll/ao7O1kef1tYfTOsWj/sJofzMYHMpx7Qc/sak+M1KahLfc4RHx\n4tXC5vlaphGu/eDH/aX2ggU0+hKupNPknsjrxHW6Y3Zu3t+hO+LP5+dn5nn2wIg/vz8ywVh+6XH7\n4B/J093hd/hP/vRJxr102C4nJZ7kA1dfUdG649qei+398Cf52eTpJ/zZq7RaaM9aqBEgm20racSj\nLYUqyuUhkTq21F4uH+3rVNP8vEdubIqBkWCEH+zE7R+ZZHBkgqGxqWUP39xSn6CrPcmpPcHcfkeS\n7qJRf0tyS3gkz7uuf4hEKgP47TKcj/2Ap2oGxuYWLPNYC58ZiYZCoYrSzR75oqX20qloX6ecPM9j\n5Og0A8FIf2CkMLc/yeDIJIOjyx/FkwA625p4QV8HXUXz+s/fnqEBj862pjVP8cSxXeKiq62OkUm1\njRyfpo/KZC2blLnh3JKl9krap3ACr1Pqpq/neYzmZxgY8Uf4A+G/Yz8vt8oW+OfkyXYk6QpG990d\nSX/KpzNJV3ty2fV1S6mzXO27HrUylVBfP8tffOATVW2bUtRCe9ZCjVD69JFCoUxq6IOybJ2FI3gW\nd/TFPy93WgaAluSWsJMvdPjZTv+yq8Rj9mu9PeNGdZZPLdQI2qcgx+F5HqNHpxnKz/DUs0Nhpz84\neqzjX251LfA7/W1dLXQHo/1wpB/8rBOxiWwc+mveIOY9j5Hx6QWj/GMd/iRDoytP77Qkt7Ato05f\nRBQKNWN2bp7c2FS403ZwZJKB4HK1HbkArc0N9HW30NWR5JTedloa6+lqT6rTF5El1BvExNTM3IIO\nf3HHPzy+/CGbAB0tjWzvbTs20m9P0hXM53e3J2lqPDanXyvzoSISDYVCFXiex9HJ2VU7/eUWVQGo\nSyRItzUGh2z6nX5X+7HL7g6dY19EykehUAazc/McHsrz01/k/M5+1J/mGRr1O/+h0akVT8PQsKWO\nrvYk23vbws6+u6jT72xrpL5u6SGbIiKVoFBYg/zkbFEH74/uh0aPze+vNrXTktxCT6bZ7/Dbk2SK\n5vK72pNaNF1EYmXTh8Ls3DzD41N+Jz9aGN1PLQiBianlR/n+1I7/bdyTetr8Hbgdxzr/rvYmnWlT\nRGrKhu6xCnP5Q4WR/ejCKZ3jjfKbm7YUdfBJMu1NC0b7Ha3Hpna0A1dENoKaDoWpmbmwww87+7Ep\ncoXR/tjkil/IKh7lZxaN7jPtSTJtSVLJmm4eEZF1i6TXM7PdwPuAGeA659x9Kz32yV/kcE8PMDQ2\ntSQAVloxC/xj83szKTJtx0b5xSP+ztYmLZAuIrJI1UPBzDLA9cBLgTbgRmDFUNjzN99dcltTQz2Z\n9iZ2bGsn0xaM7Is6/XRbk9bIFREpQRRbCucB33DO5YE8cPlqD37L77yQ+dm5YErH7/iLz6FfLZtx\ntTMRqT2FvuqHP/7VI//nnz76W+t9fhShsANoMbOvAp3Ajc65B1d68EXnnhGLHbj77rjbX+0slSDv\neezbf3csV64Skc2t0Fed8uKdv1nK86MIhQSQAV4PnAZ8C9i+2hOy2bYqlLW68em6BStXjU/XLakr\nDnWuheosL9VZXrVQZ5xrLO6rShFFKBwGHnbOecDTZjZmZt3OuYGVnhCHLYXWhjlGilb1am2cX1BX\nrRySqjrLS3WWVy3UGfcai/uqUkRx/oT7gXPMLGFm3UDLaoEQF3uu2EVfwtGcf4K+hGPP5X8UdUki\nIksU+qpf/viBR0t5ftW3FJxzB8zsi8APAA94d7VrKEW6M619CCISe4W+KpttW/dOZojoewrOuU8C\nn4zivUVEZGU6/aaIiIQUCiIiElIoiIhISKEgIiIhhYKIiIQUCiIiElIoiIhISKEgIiIhhYKIiIQU\nCiIiElIoiIhISKEgIiKhSE6Itx67r/prWhvmtPxlBWiJURFZLPZbCqONZ/AcO9m3/+6oS9lwCsv2\nTaR2qo1FBKiBUAB/+ctcPuoqNp5cnnB1JrWxiECNhILneaRTUVex8aSbPTzPA9TGIuKLfSi0Tz+p\n5S8rREuMishisd/R/NmPvi/Wi2TXMi0xKiKLxX5LQUREqkehICIiIYWCiIiEFAoiIhJSKIiISEih\nICIiIYWCiIiEFAoiIhJSKIiISEihICIiochCwcySZvYzM3tLVDWIiMhCUW4pXAcMRPj+IiKySCSh\nYGYG7AS+FsX7i4jI8qI6S+o+4F3AJRG9v4iUoBxLuGoZ2Hir+paCmV0MPOyceza4KVHtGkSkNOVY\nwlXLwMZbFFsKFwKnmdlrgZOBSTP7pXPuwZWekM22Va24E6E6y0t1llc56hyfriPReGwJ1/HpunW/\n7vFeoxbasxZqLFXVQ8E59+bCdTO7AXhmtUAAamKRnWy2TXWWkeosr3LV2dowx4jnkUgk8DyP1sb5\ndb/uaq9RC+1ZCzVC6cGl7ymIyJqVYwlXLQMbb4nCwu0x5tVKKqvO8lGd5aU6y6cWagTIZttK2l+r\nLQUREQkpFEREJKRQEBGRkEJBRERCCgUREQkpFEREJKRQEBGRkEJBRERCCgUREQkpFEREJKRQEBGR\nkEJBRERCUa28VjVa5UlEZO02/JaCVnkSEVm7DR8Kuby/uhP4l7l8xAWJiMTYhg+FdLNHYc0Iz/NI\npyIuSEQkxjZ8KGiVJxGRtdvwO5rTnWn2Xn1l1GWIiNSEDb+lICIia6dQEBGRkEJBRERCCgUREQkp\nFEREJKRQEBGRkEJBRERCCgUREQkpFEREJKRQEBGRUCSnuTCzDwEvB+qBDzrnvhxFHSIislDVtxTM\n7GzgRc65lwEXAB+pdg0iIrK8KKaPvgNcFFzPASkzS0RQh4iILFL16SPnnAdMBD++A/h6cJuIiEQs\nslNnm9nrgLcBr46qBhERWShRWJWsmszsfOBG4Hzn3EjVCxARkWVVPRTMrB34HnCuc26gqm8uIiKr\nimL66A+BLuCeYAezB7zFOferCGoREZEikUwfiYhIPOkbzSIiElIoiIhISKEgIiKhyL6nsJiZNQOf\nBnqAJmCvc+5rRfefB3wAmAX+0Tm3N6Z1PgP8ApjH34m+2zl3MIJSMbMk8G/Ajc65vy+6PRZtWVTP\nSnXGoi3N7JXAF4AfAwngX51z7ym6PxbtuYY6Y9GeQS27gfcBM8B1zrn7iu6LRXsGtaxWZyza08wu\nBS4OakgA/8E51150/27gPcAc8Ann3F2rvV5sQgF4LfCoc+5WMzsV+AbwtaL7/wb4beAg8B0z+6Jz\n7okY1ukBr3HOTSz77Oq6DljusN+4tGXBSnXGqS2/7Zx70wr3xak9V6szFu1pZhngeuClQBv+d5bu\nK3pILNpzDXXGoj2dc3cCdwKY2Ss4dhohzCyF//f1G/gh+6iZfdk5N7zS68UmFJxz9xT9eCrwy8IP\nZnYaMOicOxD8/HXgXKDqH5TV6gwkgn+RMjMDdrIwsGLVlsH7L1tnIBZtGVi2jri1J6u3V1za8zzg\nG865PJAHLi/cEbP2XLHOQFzas9j1wK6in/8j8IhzbhzAzB4CzmL5vzcgRqFQYGb/DPQBv1t0cy/Q\nX/TzEeD0ata12Ap1FuwPPtzfc85dU93KQvuAdwGXLLo9bm25Up0FcWhLgBeZ2VeADHCTc+6B4Pa4\ntedKdRbEoT13AC1m9lWgE3/a8MHgvji15w5WrrMgDu0JgJn9BvAL59yRopsXt2c/sG2114ndjmbn\n3FnA64DPFt28OI0LX3qLzAp1gr+p9l7glcBLzOz3q12bmV0MPOyceza4qbj9YtOWx6kTYtCWgZ8C\n73fOvR4/vD5lZoUBVWzak9XrhPi0ZwI/tF6Pf/6zuxbdt/ixUbXnanVCfNqz4E/w93cWW3d7xiYU\nzOxMMzsZwDn3f4EtZtYd3P0cC9OtD3++seqOUyfOuX9wzg045+aBrwMviaDMC4HXmdn38T8o15rZ\nOcF9sWlLVq8zLm2Jc+6Ac+4LwfWngUP47QYxas/j1Bmb9gQO4w8GvKDOsTj+rbN6nXFqz4KzgYcX\n3bbu9ozT9NErgO3An5pZD9BSODeSc+5ZM2sLduwewJ+y2bXyS0VTZ3Bep3uA1zrnZvBHEF+odoHO\nuTcXrpvZDcAzhc3eOLXlanXGpS2DWnYB25xz+8ysF9iK/8cWq/Zcrc44tSdwP3BXsAJjF/H9W1+x\nzpi1J2a2DRhzzs0uuuuHwCeDeueBl+EfibSi2JzmIjgs8VPAKUASf09/NzDsnPuqmb0c+BD+ps8X\nnXMfjmmdV+FvuueBx51z/yWKOgvM7Hrg5/ibjbFqy2Ir1BmLtjSzVuBz+PPKDfi/8x5i1p5rqDMW\n7RnU+g78rUMP2Ivf6caqPddQZ5za80zgZufchcHPf45/JNoPg2mt/4YfCrc75/7Xaq8Vm1AQEZHo\nxWafgoiIRE+hICIiIYWCiIiEFAoiIhJSKIiISEihICIiIYWCiIiEFAoiIhJSKIisk5n9qZndEVw3\nM/uJmbVEXZdIOSgURNbvI8AZZvYy4GPAO5xzRyOuSaQsFAoi6+Sc8/DPh3MP8P+ccw9FXJJI2SgU\nRErTBYzhnzFXZMNQKIisU3Cm3P3463VPmdkfR1ySSNkoFETW70bgXufcU8B/Bd5vZidFXJNIWejU\n2SIiEtKWgoiIhBQKIiISUiiIiEhIoSAiIiGFgoiIhBQKIiISUiiIiEhIoSAiIqH/D53qAOiYZyUA\nAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# λの予測値nuをyhatにセット\n", "d['yhat'] = fit.predict()\n", "ax = d.plot(x='x', y='yhat')\n", "d.plot.scatter(x='x', y='y', ax=ax)\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 35, "metadata": { "collapsed": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAV8AAAFhCAYAAADXxalcAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3X1wZNld3vHvvf0mdetdo3mf0djr2TOzhdd44wSKwcQY\ngnEcCieQQGGcMhQuXlOEIlSREBvWpIoUFVMUEGKKxA6mgMK82ISQAuMsEGfMJlTtFkPwzJkd787M\nal41krpbaqm77+178ke3eqQZSaO37nu79XyqVKPbc9X3XEn96PS5v3uO55xDRES6y4+7ASIiB5HC\nV0QkBgpfEZEYKHxFRGKg8BURiYHCV0QkBh0PX2PMlxljrhljfuCRx99ljIk6fXwRkSTqaPgaY/LA\nLwCfe+TxHPDjwO1OHl9EJKk63fOtAu8G7jzy+L8Bfgmod/j4IiKJ1NHwtdZG1tra2seMMU8Dz1pr\nfw/wOnl8EZGkSsdwzJ8D/kUMxxURSQyvG3M7GGN+EpgFPgP8RetzD3gr8JfW2q/d7Gudc87z1EEW\nkZ7yxNDqZs/Xs9beBs6uPmCMeW2r4AXwPI/Z2cWON64bpqaG++JcdB7J0i/nAf1zLlNTw0/cp6Ph\na4x5DvgoMA0ExphvAf6JtbbY2kVTqonIgdTR8LXWvgRs2rO11r6xk8cXEUkq3eEmIhIDha+ISAwU\nviIiMVD4iojEQOErIhIDha+ISAwUviIiMVD4iojEQOErIhIDha+ISAwUviIiMVD4iojEQOErIhID\nha+ISAwUviIiMVD4iojEQOErIhIDha+ISAwUviIiMVD4iojEQOErIhIDha+ISAwUviIiMVD4iojE\nQOErIhIDha+ISAwUviIiMVD4iojEQOErIhIDha+ISAwUviIiMVD4iojEQOErIhIDha+ISAwUviIi\nMVD4iojEQOErIhIDha+ISAwUviIiMVD4iojEQOErIhIDha+ISAzSnT6AMebLgM8AP2et/WVjzCng\n40AGqAPfaa293+l2iIgkSUd7vsaYPPALwOfWPPzTwMeste+gGco/2sk2iIgkUad7vlXg3cCPr3ns\n+1uPA8wCb+1wG0S2JXKOi5fuMDNb4eRUgQvPHsP3vLibJX2qo+FrrY2AmjFm7WMrAMYYH/hB4PlO\ntkFkuy5eusMLL98C4OpMEYC3v+V4nE2SPtbxMd+NtIL314H/aa39syftPzU13PlGdUm/nEs/nsdc\npU4m7a/b7pXz7JV2bkc/nctWYglf4BOAtdb+9HZ2np1d7HBzumNqargvzqVfz2OykCUIo3XbvXCe\n/fLzgP44l2q1xle/5zumblz67OxW+3U9fI0x7wNq1tqPdPvYIlu58OwxgHVjviLb5ZxjfqHEcuAY\nP35u4En7dzR8jTHPAR8FpoHAGPOtwGGgaoz5M8ABX7TW/lAn2yGyHb7naYxXdqVaqzFXXMJL5chk\ntldE1ukLbi8BX9vJY4iIxGm+WKJSi8hkBnf0dXGN+YqI9LSVlSrz5Uqrt5vZ8dcrfEVEdiAIAuZL\niwQNn/QOe7trKXxFRLapVC5TrgRkcoOk93h/sMJXROQJoiji3oMFnJclk9t9b3ctha+IyBYqy8vM\nl5bJ5PLs583mCl8RkQ0455hbKFILPDK5/L4/v8JXROQRtXqdBwuL+OkBUpnOTK6k8BURWaNULrO4\nHJLO7s/Y7mYUviIiQKPR4P5cEedlSWefeHfwnil8ReTAW1yqUFqqks4O7utFta0ofEXkwIqiiNm5\nIoFLdXyY4VEKXxE5kJolZCukswOkY1ixROErIgeKc44H80VqIft2w8RuKHxF5MBYWakyV6qQygyQ\n7lAJ2XYpfEWk7znnWCiWWQ5c18d2N6PwFZG+tjrRuZ8eIJ1OzmrUCl8R6VvzxRLLtWhPUz92isJX\nRPrO2onO07uY6LwbFL4ifSpyjouX7qxbENSPoaSqm6Io4sF8iXqDRPZ211L4ivSpi5fu8MLLtwC4\nOlME6OsFQkvlRcqV+r5MdN4NCl+RPjUzW9lyu19UazXmi0s4f/8mOu+GHvj7ICK7cXKqsOV2r2s0\nGtx/sMCD4jJ+ZpBUKhV3k3ZEPV+RPnXh2WMA68Z8+0WpXKZUqZPN5XtiiGEjCl+RPuV7Xt+N8QZB\nwOx8GednyXZgdYluUviKSE9Yu3JwP1D4ikii1et1ZhcWoccuqD2JwldEEsk5x0KpTKXaIJOQ+Rj2\nk8JXRBJndfYxP50jk03mHWp7pfAVkcSIooj7DxaoN7zEzD7WKQpfEUmE8uISS9UqkZ/r2fKxnTgA\npygiSRZFEXdn51hcicgO9Hdvdy31fEUkNqtju+nsIL11f9reKXxFJBYLxRJLtagvKxm2Q+ErIl0V\nhiH350rNiXASOtduNyh8RaRrlirLFBdX+r6SYTsUviLScY1GgwcLJYIopeBtUfiKSEeVy0uUKrWe\nmeS8WxS+ItIRqzOQ9ducDPtF4Ssi+26+WGK52tAQwxYUviKyb6rVGnOlpeaqwX06J8N+UfiKyJ45\n53gwX6QWJn/V4KToePgaY74M+Azwc9baXzbGnAR+neatzXeA91trg063Q0Q6o7K8zHxphXR2gHSm\nv5em308dvfZojMkDvwB8bs3DHwF+0Vr794EvAd/dyTaISGc0ZyCbZ2GxuVy75yl4d6LThR9V4N00\ne7ir3gH8YevzPwS+vsNtkAMsco7P//Vtfutzr/D5v75N5FzcTeoLS5Vlbt2bJ/IHSGeycTenJ3V0\n2MFaGwE1Y8zahwtrhhnuA/2zpKokzsVLd3jh5VsAXJ0pAvTdopLdFEURD+ZL1COfTI8vYBm3OC64\nre16eI9sb2hqarhzremyfjmXXjmPuUqdzJrK/rlKfV3be+U8nqQb51GpLDO7UGXs0ERHjzMxUejo\n8ydFHOG7ZIzJWWtrwAnWD0lsaHZ2sfOt6oKpqeG+OJdeOo/JQpYgjNZtr7a9l85jK50+j3Zvt+E1\nhxgqlY4da2KiwPx8554/SeII388B3wL8ZuvfP46hDXJAXHi2Oao1M1vh5FShvS3bUyyVWVwOdGtw\nB3Q0fI0xzwEfBaaBwBjzrcD7gF8zxnwvcAP4tU62QQ423/M0xrsLleVliuUVvHROtwZ3SKcvuL0E\nfO0G//UNnTyuiOxOGIbMFcsEDV+3BneY7nATEZxzLJTKVKohmayGGLpB4StywDWnfKySzg6S0XwM\nXaPwFTmgoiji3oMFIi+rmt0YKHxFDqCVlSoPiktkcvkDt2pwUih8RQ4Q5xxzC0WqoafebswUviIH\nxLrZx9KaBCduCl+RPte8Q61IveGrZjdBFL4ifay8uES5UiOt8rHEUfiK9KEgCHiwUCYio5slEkrh\nK9JHVm+WWF28UpUMyaXwFekDzjnKi4ssLgf4aS1e2QsUviI9rlReZObuXHNcN6uXdK/QELxIjwqC\ngNv35igvR2Ryea2h1mP0Z1KkB5XKZcqV5jy7fkoju71I4SvSQ4IgYHa+jPOzqtntcQpfkR6xdlUJ\n6X0KX5GEW6pUKC5W8bWqRKKVKnWu3Fjg2kxpW/srfEUSqlavM19cpOHSulEigZxz3Jlb5vKNBa7c\nWODWg50t/KnwFUmYKIqYWyhRDSGTGdSLNEHCRsSrt8vtwC1V6o/t42+z6kQ/V+moyDkuXrqzbvXg\nzX45d7JvP6rX6xTLFWpBRCY3SKaL90lEzvGSneXu/DJHJ/I8Z6a2/Dltd99+sFwNsDeLXL6xwNWZ\nIvUgemyfgWyKp0+N8cyZcZ4+NcYHf+xzT3xeha901MVLd3jh5VsAXJ0pAmy6mvBO9u0ny8srlCsr\nBA3IZAfI5LrfhpfsLC9+8R4A1+8uAvC2c4f3vG+velBa4fKNBS7fWODG3UWce3yfieEc56fHOX9m\nnOmjw6T8nd02ofCVjpqZrWy5vdt9+0F5cYnFSg38DKn0AJkYy3Xvzi9vub3bfXtFFDlu3l/k8vUF\nrtxcYLZYfWwfDzh1ZIjz0+Ocmx7n8Njgnm5sUfhKR52cKrR7savb+7FvL1tZqbJQruD8LKmEXEg7\nOpFv92JXt/dj3ySrBQ1emSlx5cY8V24WWa6Gj+2TSfucPTnKudPNwB0a3L+xIIWvdNSFZ48BrBvH\n3Y99e1EQBMyXFgkaPulMMkJ31XNmCmDdOO5+7Js0paUal282L5Z96VaZRvT4eMJwPsO50+Ocnx7n\nqROjZDo0EbLCVzrK97xtj9vuZN9esrpu2krdkckOJHJSc9/ztj1uu5N94+ac4/bcMpevN3u3tzcp\nBzs6kW8PJ5yYKnTlAqLCV6SDyotLlJaqpLODZLL9WxGQJNspB0v5Hm84NtIK3DHGhwe63k6Fr0gH\n1Ot15oqLRGS0SnAXVNaUg72ySTnYYC6FOdXs3T59apSBmKffVPiK7KNGo8F8sdy+QULzjXXObHGl\n3bu9cW/jcrDJkYF273b66AgpPznvPhS+IvtgpVplqbJCNYjIZLt7g8RB0YgcN+8tcqVVf/ugtHE5\n2Okjw+3x26mxgcTOc6zwFdmlaq1GabFCPYjw/AzpTI5MNu5W9ZdavcHVmSJXbixgbxZZrm1eDnZ+\nehxzen/LwTpJ4SuyQ/V6nYXyEkHDI52J5460flZcqrV7t6/e3rwc7Px0sxzsjcc7Vw7WSQpfkW1a\nWalSWlomiDwymWSWjPUi5xy3H1RaF8tKvH5/acP9jk0+LAc7fqg75WCdpPAVeYLy4hJLyzUiUs2e\nrq6i7VkQRrx6u9S8YHazSHmrcrAzzR7u2FB/vcVQ+IpsolReZHG5jpfKksoMarXZPVpaCbj6epHL\n11vlYOHj5WCFgTRnT45x/sw4Z0/GXw7WSf17ZiK7VC4vUV6u4adzmsR8D5xzzJaqXLneHL+9eW+R\nDarBmBwdaI/ffvn5I5SKK11vaxwUvnLg1Wo1qrUa9SCiFjQUunvQiBw37i5y5WYzcOe2KAc7Nz3G\n+TMTHB57+L3e6bSMvUzhKwdWZXmF2/fmiFyKdDYLHqRVKrZj1XrIKzMlLl9fwL6+wEqt8dg+2bTP\n2ZNjnJse66lysE5S+MqBs1SpsFipMjo+iq+x3F1ZWGyWg125uXk52Eghy7nTYz1dDtZJCl85EMIw\npLxUYbka4vkZUulB0pkM8PhVdnlc1CoHW62/vTO38QTqxyfznGuN3x4/VEjs3WVJoPCVvuScY3Gp\nQq0eUAsiIjyy2QHSWb3d3a515WA3FigvB4/tk/I9njoxwrnpcc6d7r9ysE5S+EpfiKKIlZUqtXpA\nPWxQDyLS2QF8P6dx3B1YWgmwrYtlr8yUCDYoB8vn0pjWcMLZk2Pksip83o2uh68xpgB8EpgAMsBH\nrLWf7XY7pLc1Gg0qy8vU6iH1MKIRQSqdJZXKQCqD8mB7nHPcL660hxNev7e0aTnYM627y04fGU7U\n7GC9Ko6e7weAK9banzDGHANeAM7H0A7pMZXKMiu1OrWgQSNqrvTreTlSGTR14w40y8HK7dV558u1\nx/bxvIezg52fHmdqTKV3+y2O8H0AvLn1+QQwG0MbJOGiKCIIAmr1OrV6g2o9xEtlSaezCttdqNZD\nrr5e5MqN4ublYBmfsyead5eZ02MUBjQ+3klPDF9jzLuAz1prN3o3smPW2t82xnzAGPMKMAa8Zz+e\nV3qLc46lSoWwERGGEZFzNCJH1Ppwnofvp0mn03heikwuOQO3kXO8ZGfXLSCZxEleVsvBLt9Y4LU7\nG5eDjRay7eqENx4fIZ3av3KwXvk+xWU7Pd8fBv6jMeY3gI9ba2/s5YDGmPcBN6y17zbGPAv8Z+Dv\n7eU5JZmiKKJer1MPAqIoInIQRY6gERG0L4ilm7c8eeD79ETN7Ut2lhe/eA+gvYR6EhaUjJzj9myl\nvTrvpuVghwqcOz3GM2cmODaZ71g5WFK/T0nxxPC11v5DY8w48I+B/2SMAfgE8PvW2sffuzzZBeBP\nWs99yRhzwhjjW2sfv6zaMjU1vKMDOOeaL/YoIgwbhGFII3J4HuQHB8jEuMzATs8lqR49jzAMWV6p\nUq0FBGFEPWjg8PFTaXK55K4mMDFR2PHXFJfrpFPeuu3dPM9+qAcN/ubaAy5dm+XStQeUlh6vW06n\nPMz0BG85e4g3P3WI8ZHuLBa52+9TXN/LbtvWmK+1dsEY81s0K9J/APhXwIeNMR+01r64w2NeA74S\n+LQxZhpY3Cp4gyDg6rXb+H6K1ddv5BzOPfzXrfl39S2r53l4+Hi+j9/6cM4RBgu4qIHvN/dpfq3D\n9308D3wPUikf3/OaHz54nkc6nSaTTpNKpfB3ef/51NQws7OLu/raTnLOEYYh1VqNMGw03/63hgHA\na39/V/cdHy9QLK72qhxRBPg+qVTmke9No/WRzBsZJiYKzM9vvJT4VsbyWcKGW7e9m+fZrcXlenux\nyGu3dlYO5sJG19q6m+/Tbn8mvWg7Y75vB74LeCfw+8B3W2uvGGPOAJ8G3rrDY/4K8HFjzJ/TvG7y\nvVvtHAQh9ShN2k+zrgam9Va19c+2eJ5HJvvkIvAG0HA0jxc1A6exEkBUIyKCKGqGu/fweZtcK/Tb\nB2w3FSCIAubnl9pf0wz95v7NbXC0xj3dw/FP2OhYzbf1nue1/witLf/xPVpB+PAPC85B6ziNKGr+\nAWsdw/NTpNIZfD/T/v56qfan62Ryg6QyD1/wB+3i13NmCmDdWGYnbbcc7NDq7GBnxjl9eBg/5nKw\nbn+fes12er4/A3wM+H5rbbsmxVp73RjzqZ0e0FpbAb5tp18Xp9We7145P4vzm+HvALxm0D98oMVv\nBl4qtXWwbTf02jHpPfzXSzU3e2GMNWl8z+v42GUjirh+d7E5HePNzcvBpo8O89y5I5yZKnAoYeVg\n3fg+9bLtjPl+9Rb/9zP72xyRg2u1HOxya7HIan3jcrCnT461bucdIz+QOVBv1fuJbi8WidHCYpXL\nN5qr8756u0zkulsOJvFR+Ip0UeQct2Yr7clq7s5vXA52YqrAudPNwO1kOZjER+Er0mH1sMGXbpWb\n89/eWGBx5fHZwdIpj6eOjzaHE6bHGS0k56YS6QyFr0gHrCsHmykRNDYoBxtIt3u3bzo5Sk7LIh8o\nCl+RfeCc497CmnKw+0sb7jc1NsD56QnOT49z6vBQ7OVgEh+Fr8guNaKI63cW27ODLSw+Xg7me3D6\n6DDPTE9wbnqMQ6PJKgeT+Ch8RXZgpfawHOzq6xuXg+UyKc6eGuWZ6QmePjVGfkAvM3mcfitEnmC+\nXG33bq/fWdywHGxs6GE52BuOqRxMnkzhK/KIZjnYEpevNwP33sLKhvudnCq0A/fohMrBZGcUviK0\nysFmWotF3iyytFk52InRdoXCiMrBZA8UvnJglZZq/NWV+1y+vsC1W8V1M3CtKrTKwc5Nj3P25ChZ\nlYPJPlH4yoGx3XKww+ODnG8tha5yMOkUha/0tUYU8VqrHOzKFuVg00dH2otFTo52Z7JxOdgUvtJ3\ntlsO9uY3TfLGYyOYU2MM5vRSkO7Sb5z0hfXlYGU2WCuSsaEs58807y47c3SYw1PDmopRYqPwlZ4U\nOcfM/aV24N5/QjnYM2cmODI+qHIwSQyFr/SMetDg2q2H5WCVTcrB3nRilPPT45jpcUbyKgeTZFL4\nSqKVl+vYVu/22q3SxuVggxnOtRaLfNMJlYNJb1D49rnIOV6ys+sWMfQ3eeu9k307xTnH3fnldnXC\nzOzGY7Ij+QzD+SxnT47yzredJL3LFaU7LQnfU0kmhW+fe8nO8uIX7wFw/W5z2frNFjXcyb77KWys\nnR1snuLS40vN+x6cOdYsBwvDiL95bR4HXJ0pMTEykNiFGuP6nkryKXz73KPL1Gy2bM1O992rlVrY\nmmx8nquvl6gFG5eDPX1qjGfOjPP0mnKw//6F611r515183sqvUXh2+eOTuTbPa7V7f3YdzfmSg/L\nwW7c3bgcbHw415ys5vQ4Z44Nbzg7WKfbuZ96qa3SXQrfPvecmQJYN+a4H/tuRxQ5Xl9TDjZb3Lgc\n7NThoebtvNPj2yoH2+92dlIvtVW6S+Hb53zP2/YY40723Uy7HOz6AlduLlCpho/tk0n5PHVilPPT\nY7sqB9uPdnZLL7VVukvhK3tWrtS5cnOBy9cX+NLtjcvBhgYzreGEMZ46OUo2rXIwOdgUvrJjzjnu\nzFXawwm3NikHO7I6O9j0OCcPD6nESrbNOUcURURRhOd5+L6Pn9Bywt1S+Mq2hI2I1+6UW5PVlJgv\nVx/bx/c8zhwb5pkzzekYJ0Y0O9hB4Jxb/xE1gIfbHuD7Hp7n4eHwfR/PoxmqrX9Xt/OZiGqmged5\npPw0vu8ROUej0cC5Bq61hJNzzY9G5NrLOq17Lrx22/BW93d4ntf6PMLhNUPeQaMREUUOz0+RSme6\nEvQKX9nUcjXAtmYHe2WTcrCBbLMc7Pz0+nIw6V1hGNJoBPg40im//bHKOfcwTFuBB+B7KVIpn1Qq\n1e6p7nQujfGxYcIgnndIzjnCMKRaq9FohDQi1348cg7nmjfNrAZ72wabS/MzG19dXkOvFFnnQWll\nTTnYIhusFcmh0YF24J45Nkyqz94OJkkQBOAcfivQgPbbcVyEcxG1FY9adU39sOe1hng8nIvwcKRS\nPr7X/D8cOBwp32/3FqH5zsXzHNnBLAO5ob57m/8knueRyWTIZDJ7fq5r//f3HzxpH4XvAfewHGy+\nVQ72+HCCB5xslYOdnx7n3FOHWFjQzQL7qdFo0Ajr+F5zcqB0yieTTjGWz+H7PmHYIGw0K0dSfopU\nKtPuYR45Msrw4OCat+Rr3vJ7HqmULm4mkcL3AKoFDa61Fou0W5SDvelka3aw02MMrykH07SMexdF\nEWG9Sibtk834DOezDA6Ob9rbzG5Rjbf683j0X0k2he8BUarUudKarGazcrDhwQym1bt904lRMumD\n9bazExqNBlEUNq/a07z4lM2kGMqlGZ6cVFAeYArfPtUsB3s4O9itBxuXgx2dyDfrb6fHOTFVUDnY\nHgX1KinPkcumSKdSZAbTZDODpFIpBa2so/DtI2Ej4tXb5XbgliobzQ7m8Ybjw+3x2/FhlYPtVRDU\nSBExkEtz6NAI6bReVvJk+i3pccvVoDU72AJXZ4rUg+ixfQayKczpMc6dbo7fDmT1Y9+t1QtjHo50\n2ieT8hkfy5PL5eJumvQYvQp70IPimnKwexuXg00M5zh/pnl32ZmjKgfbLeccYb3KYC5FOuWRHciS\ny+VVQSB7pvDtAVHkuHl/kcvXm4H7oLRxOdipI63ZwU6Pc1iLRe5Jo9EgCmsUBjMcPTqh76XsO4Vv\nQtXqDV6ZKXLl5gJXbhRZrm1QDpb2OXtylHOnmz3cocG9F4cfNGEQEDWC5l8v17wHIZdNMZTPUsir\nGkE6R+GbIKWlGpdvtsrBbpXbtzeuNZzPcO5082LZUyoH25UwDHGNOhkvy+RojsGB0bibJAeQwjdG\nzjluzy1z+fo8V24Wub1FOdj5M83APX5I5WC7EQQBRCG5rM94IUehMMLUoWFmZxef/MUiHaDw7bLt\nlIOlfI83tBaLPDc9zviwrqRvV3NylDpEjXY1wupturnciIYRJDFiCV9jzPuAHwMC4EPW2j+Oox3d\nUmmVg33pzpf421fnNiwHG8ylMKeaYfv0qVGVg21TUK/hogaZtN/8yKTIjwzty+QoIp3U9Ve4MWYC\n+DDwVmAYeB7ou/CdbZWDXdmqHGwk17rZYYLpo8OkfPXKtuKca95B5kM27ZPNpJkYKpDdauIDkYSK\no3v19cCfWmuXgWXg+2Jow75rRI6b9xa5cmPrcrDTR4Y5Nz3G+ekJpsYGdv02OHKOl+zsuoUZ+3Us\nOAwCfEIGc2kOHx5TjW0fipzj4qU7zFXqTBayXHj2WN/+Pq+KI3zPAAVjzB8AY8Dz1toXYmjHntXq\nDa7OFLlyYwF7c+tysPPT43zFsycIa8G+HPslO8uLX7wH0F6avJ8Walytsx3IphgdGWRwUBUJ/ezi\npTu88PItMmmfIGwOy739LcdjblVnxRG+HjABvBd4A/BnwPRWXzAxUUjM/fLz5SqXrj3g0iuzXL25\nsOHsYKNDOZ590yHecvYQZnqczNrFIgv78xa5uFwnnfLWbU9MFPblubejE8dyzlGvVcnnUgzlhxga\n6vz5TE0Nd/wY3dDr5zFXqbfLJjNpn7lKvefP6UniSLR7wBestQ541RizaIw5ZK3ddOb3+flKbOG7\nthzs8o0F7sxtPIn4scnm7GDPtMrBVocTFtesdTYxUWB+fuNysp0ay2fXBf9YPrtvz/0k+3keAPXa\nCumUx2A2xejIML7ns7ISsbLS2TKwqan+KDXrh/OYLGQJwqjd850sZHv6nLbzhyOORPss8AljzM8C\nk0Bhq+CNQxBGvHq7Odn4lZtFypuUg73x+Ejzhocz44wNdbcc7DkzBbBuzLcXhPU6kQtJ+c3VGrJp\nnymN4x54F549BrBuzLffdT18rbW3jTG/C7wIOOCHut2GjSytBFx9vcjl6wu8MlOkHm5UDpbm3Okx\nzk2Pc/ZkvOVgvuf1xBhvFEWEQbVVnZBifHxQM4DJY3zP4+1vOd4XvfjtiiU9rLW/CvxqHMde5Zxj\ntlTlSmuympv3FtmgGozJ0YH2ZDUqB9ue1YtluWyKfC6jFRtENpCMq1hd0ogcN+4ucuVmM3DnNioH\n85rlYKt3lx0eG4yhpb2nOc9tjVwmxfBAhuEhBa7IVvo+fKv1kFdmSly+voB9fYGVWuOxfbJpn7Mn\nxzh/ZpynT41pdrBtajQaENXJZVIM5TOaBUxkB/oyfItLtfbdZa/e3nh2sJF8pr122RuPa3awnYii\niKC2wuhQjpHhybibI9KT+iJ8I+e4/aDSDtytysHOT49z/swExyfz6qXtUBgEeC5kKDdIThOMi+xJ\nz4bvunKwGwuUlx+/c2y1HGx1/Lbb5WD9IgxqpH3HxPAg+fwo42MH54q0SKf0VPgurQTY1sWyazOl\nDcvB8rk05vQY56fHOXtyjFxW9aO7FQYBKS/k0NgQAyoPE9lXiQ/fmftLfP7SHezrJV6/t7RhOdih\nVjnY+TO8hMTCAAAK3klEQVTjnD48jK9ysF1ZnTUsnfLIpn1GhnPk85pTQaQTEh++P/qLX3jssdVy\nsGdawwlTKgfbk2Zdbp3CYJojR8bxtdKxSMclPnxXZTOtcrDpcczpMQoDKgfbq0YY4rmAkcIAw0Oq\nWhDppsSH73d+49N4kePsqXHSKfXI9sNq6I4ND1LIj8TdHJEDKfHh+00XzmBfm1fw7oNGGIJCVyQR\nEh++snerodscXlDoiiSBwrePBUGNlBcxNqSerkjSKHz7TBRFNIIqg7k042N5Td8oklAK3z4RRRFR\nWGM4n2FEUziKJJ7Ct8c55wjrK4wUcoxMar4FkV6h8O1hYVAjl4YjRyZ0Y4RIj1H49qBGGOIRaM4F\nkR6m8O2CyDlesrMUl+uM5bM8Z6bwdzE8sDquOzoUf8lY5BwXL91hZrbCyakCF549tqtzEjmoFL5d\n8JKd5cUv3iOd8trLve9k8cswCPAIGRrMJmZc9+KlO7zw8i0Ars4UAXj7W47H2SSRnqLw7YK788tb\nbm/GOUcjqDIxkk/c7GIzs5Utt0Vka7pK0wVHJ/Jbbm8kDOqkqXPiyAT5fPJmbTs5VdhyW0S2pp5v\nFzxnpgDWjfluZrW3Oz4ySCH/5JCOy4VnjwGsG/MVke1T+HaB73m87dxhJiYKzM9v/vY8qK9QyKUZ\nP5KMcd2t+J6nMV6RPVD4JkAjDEl5AcenxkiltOyRyEGg8I1ZUF9hNJ9jZESTmYscJArfmERRhGvU\nODo5QiajVTlEDhqFbwzCoEY+6zNxSL1dkYNK4dtFzjmC2gqTY3nyg8krHxOR7lH4dklQW2EgnePk\n0eRXMohI5+kmiw6KooigtkLK1Th+eIxDE2MKXhEB1PPtCOccUVBlKJ+cuRhEJFkUvvssDGoMZjwm\neuBGCRGJj8J3n6xO9zg5WmBwcCDu5ohIwil892h1LoahwQyjGmIQkW1S+O5BUK8ymPU4cnhcy/iI\nyI4ofHehORdDyJGJIbLZbNzNEZEepPDdgUYYggsYKcS/jI+I9DaF7zZEUYQLa4wkYO00EekPCt8n\nCIIahZzmYRCR/aXw3UQjDPFcwJSWZxeRDogtfI0xA8DfAs9baz8ZVzseFQYBPiGjQwMMFTTEICKd\nEWfP90PAgxiPv06j0YCozvjwIIWErRQsIv0nlvA1xhjgHPBHcRx/rdWbJEYKOUaGNa4r8Yic4+Kl\nO+sWJPU3uWHn0X3f+86nu9xa2Q9x9Xw/Cvwg8IGYjg8052HIpeGo5mGQmF28dIcXXr4FwNWZIsCm\nC5Q+uu/w8ABf/saJ7jRU9k3Xw9cY837gC9baG80OME9MvYmJAun0/jXVOUdYX+Hw5GHyXZ6HYWpq\nuKvH6xSdx/6aq9TJpP1125u17dF9r98t8w++YrrjbeyWpPxMOi2Onu97gDcYY74JOAlUjTGvW2tf\n2OwL5ucr+xa+q73dQxNjVJYCKkvBvjzvdkxNDTM7u9i143WKzmP/TRayBGG0bnuztj2675mjI4k5\nj71K0s9kL7bzB6Tr4Wut/fbVz40xPwm8tlXw7pfVGyU065gk0YVnjwGsG/Pd7r5f93dPMze31JV2\nyv7p+zrf1akeh/NZRnWjhCSU73mbjvE+aV/f1/WKXhRr+Fprn+/Uc6+O644UtJqEiCRPX/Z8g9oK\ngzmfI0cmNNWjiCRSX4VvIwzxCTh6aIRMJhN3c0RENtU34RvUVxjN5xgZ0biuiCRfz4dvFEW4Ro2j\nk+rtikjv6OnwbQR1BjIwqSoGEekxPRu+QW2FidFBCvl83E0REdmxngvfKIqgUePY1Oi+3nIsItJN\nPZVeYVBjMONpmEFEel7PhG9QW2ZiNK9hBhHpC4kPX8/z8KM6R3XDhIj0kcSH7+DgAMePHoq7GSIi\n+0pdSRGRGCh8RURioPAVEYmBwldEJAYKXxGRGCh8RURioPAVEYmBwldEJAaJv8niIImc4+KlO+tW\nsPW19pxIX1L4JsjFS3d44eVbAFydKQJse0VbEektGnZIkJnZypbbItI/FL4JcnKqsOW2iPQPDTsk\nyIVnjwGsG/MVkf6k8E0Q3/M0xityQGjYQUQkBgpfEZEYKHxFRGKg8BURiYHCV0QkBgpfEZEYKHxF\nRGKg8BURiYHCV0QkBgpfEZEYKHxFRGKg8BURiYHCV0QkBgpfEZEYKHxFRGKg8BURiUEsk6kbY34W\n+GogBfx7a+2n42iHiEhcut7zNca8A3jGWvtVwLuBn+92G0RE4hbHsMNfAP+09fkCkDfGeDG0Q0Qk\nNl0fdrDWOmCltflB4H+0HhMROTBiW0DTGPPNwHcB3xBXG0RE4uI51/1OpzHmXcDzwLustaWuN0BE\nJGZdD19jzAjweeDrrLUPunpwEZGEiGPY4duASeBTrQttDvjn1tqZGNoiIhKLWIYdREQOOt3hJiIS\nA4WviEgMFL4iIjGIrc53O/phDghjzCDwX4EjQA74d9baP4q1UXtgjBkA/hZ43lr7ybjbsxvGmL8P\n/A7w/wAPuGSt/eF4W7U7xpj3AT8GBMCHrLV/HHOTdswY893A+2lefPeAv2OtHYm3VbtjjCkAnwQm\ngAzwEWvtZzfaN7Hhu3YOCGPMBPAy0HPhC3wT8FfW2v9gjDkN/CnQs+ELfAjohxLBP7fW/rO4G7EX\nrdfFh4G3AsM0a+d7LnyttR8HPg5gjPkaHk4/0Is+AFyx1v6EMeYY8AJwfqMdExu+NOeA+D+tz9tz\nQPTarcjW2k+t2TwNvB5XW/bKGGOAc/T2H49V/TCfyNcDf2qtXQaWge+LuT374cPAd8TdiD14ALy5\n9fkEMLvZjokN336bA8IYcxE4AfyjuNuyBx8FfpDmX/de94wx5jM0XyAfsdZ+Lu4G7cIZoGCM+QNg\njOZQ0AvxNmn3jDFvA25aa+/H3Zbdstb+tjHmA8aYV2j+TN6z2b6Jv+C2Zg6IH4q7LXthrb0AfDPw\nG3G3ZTeMMe8HvmCtvdF6qJd7jq8AP2WtfS/NPyT/xRiT2I7IFjyafzzeS/M18ol4m7Nn30Pz+kjP\nao3B37DWngW+DvilzfZNdPi25oD418A3WmsX427PbhhjnjPGnASw1v41kDbGHIq5WbvxHuCbjTF/\nSfNF8m+NMe+MuU27Yq29ba39ndbnrwJ3ab4r6TX3aP5BdK3zWOzR361V7wC+EHcj9ugC8CcA1tpL\nwAljzIY5m9i/9q05IH6W5hwQvTz5ztcA08CPGGOOAIVenNPCWvvtq58bY34SeK1X3+IaY74DOGat\n/agx5ihwGLgVc7N247PAJ1pVQZP06O8WQOvi1KK1Noy7LXt0DfhK4NPGmGma5xRttGNiw5f+mQPi\nYzTf1v4vYAD4gZjbI/DfgN9sDWllgO/rxRe9tfa2MeZ3gRdpvj56eWjuGNCzY71r/ArwcWPMn9Ms\nkf3ezXbU3A4iIjFI9JiviEi/UviKiMRA4SsiEgOFr4hIDBS+IiIxUPiKiMRA4SsiEgOFr4hIDBS+\ncuAYY37EGPMrrc+NMeZyaxJska5R+MpB9PPA08aYr6I569QHrbWVmNskB4zCVw6c1rzQ3wN8Cvgb\na+3/jrlJcgApfOWgmgQWac44J9J1Cl85cFqLgH6M5vp6NWPMd8bcJDmAFL5yED0P/J619hrwL4Gf\nMsYcj7lNcsBoSkkRkRio5ysiEgOFr4hIDBS+IiIxUPiKiMRA4SsiEgOFr4hIDBS+IiIxUPiKiMTg\n/wMgCjlQVRSeGAAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# seabornの線形回帰機能を使って、入力データのプロット(lmplot)\n", "sns.lmplot(x='x', y='y', data=d)\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": true }, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "SageMath 7.3", "language": "", "name": "sagemath" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 2 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython2", "version": "2.7.10" } }, "nbformat": 4, "nbformat_minor": 0 }