Estuary+Python
Python codes for data analysis
Basic
(c) 2015-10-24 Teruhisa Okada
You can draw a least square line using scipy.optimize.leastsq.
scipy is included in winpython or anaconda.
{ "cells": [ { "cell_type": "code", "execution_count": 10, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[ 1.07660358] 0.987796142524\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXEAAAEACAYAAABF+UbAAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XlclXX6//HXB2lAE7ec1JTSaBHNtTQ1SywF08mc79hi\n+W21Mg+L61RqaY9yzPxaKmmlky2/lmnKyckBBZcINTO3TOGkDi655K64ggKf3x8QuaDCkcPNObyf\njwePzrnPvVycDpfXue7P/bmNtRYREfFNAU4HICIinlMSFxHxYUriIiI+TElcRMSHKYmLiPgwJXER\nER92wSRujJlhjNltjFl72rLxxhi3MWaNMeZfxpjq3g9TRESKcrFK/H2g21nLkoGm1toWwAbgBW8E\nJiIiF3fBJG6tXQQcPGvZPGttXsHTZUADL8UmIiIXcak98SeAxNIIRERESs7jJG6MGQGctNZ+Worx\niIhICQR6spEx5jGgO3DXBdbRpCwiIh6w1prirlviStwY0w0YBtxrrc26SCD6sZZRo0Y5HkN5+dF7\nofdC78WFf0rqYkMMPwO+A240xmwzxjwBxANVgXnGmNXGmKklPqqIiJSKC7ZTrLV9ilg8w0uxiIhI\nCemKzTIQERHhdAjlht6L3+m9+J3eC88ZT3owxdqxMdZb+xYR8VfGGGwJTmx6NDpFxJ8YU+y/F5FS\nVRqFrpK4CKXzxyRSEqVVPKgnLiLiw5TERUR8mJK4iIgPUxIXEfFhSuIifm7Xrl307NmT+vXrExAQ\nwC+//HLB9bds2ULnzp25/PLLCQ8PZ8GCBWe8Pm3aNK677jqqV69OmzZtWLJkSbHi2LJlCwEBAYSE\nhBASEkLDhg155ZVXzlnvww8/5Ntvvz1n+csvv8zWrVsLn2/YsIF7772XK6+8kiuuuIJu3bqxYcOG\nYsVytjfffJN69epRvXp1nnzySU6ePHnedWfPns1NN91ESEgIt912G263u/C17OxsBg0aRP369alV\nqxYul4ucnByPYiouJXERPxcQEED37t2ZOXNmsdbv06cPN998MwcOHGDMmDH07t2bffv2AfDjjz8y\nZMgQvvjiCzIzM3nyySf585//XKLRPZmZmRw5coSZM2cybtw4EhPzZ7OeNm0as2bNAvJHC02bNo1/\n//vfjB07lsWLFwOQk5PDmDFjWLZsGZmZmfTq1YsNGzawe/du2rZty7333luStwaApKQkxo0bx8KF\nC9m6dSubNm1i1KhRRa67ceNG+vbty7Rp08jMzOSee+6hZ8+e5OXl32LhtddeY9WqVaSlpbFhwwZW\nrVrFq6++WuKYSsSLk7hYEV9Qnj+rY8eOtWFhYTYkJMQ2adLEfvXVVx7v69SpU9YYY7du3Xreddav\nX2+DgoLs0aNHC5fdcccd9p133rHWWvvJJ5/Ytm3bFr529OhRa4yxu3btuujxN2/ebI0xNjc3t3BZ\n27Zt7fjx4wvjGz9+vG3ZsqVt3769ff/996211h47dswOHz7cNm7c2Hbt2tV+/fXXRe5///791hhj\nDxw4cNFYTtenTx87YsSIwucLFy60devWLXLd+Ph426NHj8LneXl5tnLlynbhwoXWWmtvueUW+8UX\nXxS+/umnn9rQ0NAi93W+z13B8mLnWlXiIuXYddddx+LFizl8+DCjRo2ib9++7N69G4DFixdTs2bN\n8/589913JT5eWloa1157LZdffnnhshYtWpCWlgbA7bffzubNm/nhhx/Izc1lxowZtGrVijp16hT7\nGLagav/+++9JS0ujTZs2hcsDAvJT0m//PXv5hcZWp6amUq9ePWrWrAnAp59+esH3Z/v27QCkp6fT\nokWLwv00b96c3bt3c/DgwXOOUXA1ZeHzvLw8rLWsW7funN/vt9e3b9/OkSNHivnulJwu9hG5CPNy\n6VyUYUeV/IKi3r17Fz6+//77GTt2LMuWLaNnz5507NixyERzKY4ePUr16mfe+7xatWrs2LEDgNDQ\nUF599VVuu+02AGrWrFnYDimu2rVrk52dTVZWFuPHj6dTp04AvP/++zRq1IiBAwdyzTXXsH79embN\nmoXb7SYyMpLAwED69u3LZ599xpVXXsmtt95auM/t27cTHR3NG2+8UbjsoYce4qGHHirx71ytWjUA\njhw5UvgPwm+6dOnCc889x7fffkv79u0ZN24cJ0+e5Pjx4wB069aNSZMm0blzZ3Jycpg8eTLGGI4f\nP05ISEiJ3qfiUhIXuQhPkm9p+eijj3jzzTfZsmULkJ9w9u/f77XjVa1alcOHD5+x7NChQ4WJ7euv\nv2bChAm43W6uu+46kpKS+NOf/sTq1aupV69esY7xW/yTJk3ijTfe4KmnnqJatWo8/fTTQP6JTWMM\nzzzzDAC9evUCYOHChQQGBjJy5Mgz9rd3714iIyNxuVw88MADl/w7Z2ZmAhSZdG+88UY+/PBDoqOj\n+fXXX+nbty9NmjShQYP8Ww2PGDGCQ4cO0bJlS4KDg+nXrx8//vhjib6plFhJei8l+aEc9xlFTlde\nP6tbtmyxQUFBdsmSJTYvL89aa23Lli3te++9Z621NjU11VatWvW8P4sXLz5jf8XtiQcHB9sjR44U\nLuvYsaN99913rbXWDhgwwA4aNOiMbVq2bGlnzpx50d+nqJ54x44d7ejRoy+67fkcOHDAtmzZ0r7w\nwgvnvPbxxx+f970JCQmx27Zts9Za+9BDD53RE58/f/55e+JnO3jwoK1atapdv359ka+/++67tkOH\nDkW+dr7PHSXsiSuJS4VXXj+raWlpNjg42K5fv97m5OTYGTNm2MDAwMIkXhInTpywR44cscYYu379\nenvixInzrtuuXTs7dOhQe+LECTtz5kxbo0YNu2/fPmttflK64YYb7KZNm2xeXp5NTk62VapUKUxi\no0aNshEREUXut6gk/p///MdeccUV9tixYyX+nTIzM22bNm1sdHR0ibc93dy5c23dunVtenq6PXDg\ngO3UqVOR/yj8ZsWKFTYnJ8fu2bPH3nffffbhhx8ufG3Hjh12x44dNi8vzy5dutSGhobaefPmFbkf\nJXGRUlKeP6sjRoywtWrVsrVr17aDBw+2ERERHiVxY4w1xtiAgIDC//6mf//+tn///oXPt2zZYiMi\nImzlypVt48aN7YIFCwpfy83NtcOGDbMNGjQoHDHz8ccfF77+xBNP2JEjRxYZw+bNm21AQMAZSdxa\na5s2bWonTZpU4t/pgw8+sMYYe/nllxdZYZfEG2+8YevUqWOrVatmn3jiCXvy5MnC1+6++247duzY\nwucdO3a0ISEhtlatWrZ///72+PHjha+lpqbahg0b2ipVqtjGjRvbTz/99LzHLK0krvnEpcI7e8SB\neK5Vq1YsXLjwnBOCcq6iPne5ebkEVgrEevNGySIi57N69WolcA9N+Pxtaj7XoMTbaXSKiIiD9hzb\nw8PvP8Y32xaRm/w28L8l2l6VuIiIQ6b8MIWmU5uyce1+ciftgJ/6lngfSuIiIg75Iv0Lvnn0Gxqu\nvxuyq3m0D7VTREQc8s2j32CMISjoHx7vQ5W4iIhDfpsLJjY2krCwER7tQ5W4CLrjvTirR487AIiP\nf5GkpJJtq3HiIiJesOvoLv467698s+UbJkRO4L4m9xWrWCgYP65x4iIiTsjJy2Hi9xNp9nYzrgq5\nCrfLzf1N7/fatz21U0RESknq1lRciS7qVq3LoscX0bh2Y68f84JJ3BgzA+gB7LHWNitYVgv4HLgG\n2ALcb6095OU4RUTKrZ1HdjJs3jAWbV3EG1Fv8Jfwv5TZeZaLtVPeB7qdtex5YJ619gZgQcFzEZEK\n51TuKSZ8N4HmbzfnmurX4Ha56d2kd5meKL9gJW6tXWSMaXjW4p5Ap4LHHwIpKJGLSAWTsiUFV6KL\n0GqhfPfkd9xwxQ2OxOFJT7yOtXZ3wePdgBdvWSEiUr7sOLyDofOGsnTbUt6MepNejXs5OkT1kkan\n/Db3bSnFIiJSbp3MPcn4JeNp8U4LwmqGke5K58/hf3b8GgNPKvHdxpi61tpdxph6wJ7zrTh69OjC\nxxEREURERHhwOBER70lISGXy5GSyswMJCsohNjay8OKb3yzYtIDoOdE0qtGIpU8u5forri+146ek\npJCSkuLx9he92KegJz77tNEprwP7rbXjjDHPAzWstef0xHWxj4iUdwkJqcTFJZGRMaZwWVjYCCZN\niqJHjzvYfng7g5MGs3znciZGTaTnjT29XnmX9GKfCyZxY8xn5J/ErE1+//sl4N/AP4GrucAQQyVx\nESnvoqJGkpz86jnLu3Ybzl0jqjP+u/G42rh4vuPzVL6scpnEVNIkfrHRKX3O81KXEkUlIlIOZWcX\nkQKvncfiZtMI2taeZf2WEVYrrOwDKwFdsSkiFVZQUM7vT6r/AlGDod4qGm+7i9mvf+5cYCWgJC4i\nFVZsbCT/3fwcm+pWh/ZvwLIYGv3YiFfeuMfp0IpNSVxEKqxKNx4n64lP+eOvVWi05mFqBucQ88Y9\n54xOKc+UxEWkwtl6aCuDkgbx0+6fmP6Xd+l+fXenQ/KYpqIVkQojKyeLV1Nf5eZpN9O6XmvWDVjn\n0wkcVImLSAWRuDGR2DmxNKvTjBVPr6BhjYZOh1QqlMRFxK9tPriZgUkDce9181b3t+h23dkTs/o2\ntVNExC+dOHWCl1Neps30Ntxa/1bWPrvW7xI4qBIXET/0nw3/IW5uHK3qtmLVM6u4uvrVTofkNUri\nIuI3Nh3cRNzcODbs38DbPd4mMizS6ZC8Tu0UEfF5J06dYNQ3o2g7vS23hd7GT/1/qhAJHFSJi4gP\ns9by9fqvGZg0kDZXtWH1M6sJrR7qdFhlSklcRHzSfw/8l9g5sWw6uInp90yny7UVc14+tVNExKcc\nP3WcFxe+SLu/tyOiYQQ/PftThU3goEpcRHyEtZZZP89iUNIg2jVox4/9f6RBtQZOh+U4JXERKfc2\n7N9A7JxYfsn8hRn3zuDORnc6HVK5oXaKiJRbx04eY/iC4XR4rwNdr+3Kmv5rlMDPokpcRModay3/\ncv+LwcmD84cMPvsTV4Vc5XRY5ZKSuIiUK+v3rSdmTgw7j+zkw14fEtEwwumQyjW1U0SkXDh68ijP\nz3+e22bcxt3X3c3qZ1YrgReDKnERcZS1li/Sv2BI8hAiGkaw9tm11Aup53RYPkNJXEQc497rJmZO\nDHuO7eHT//mU26+53emQfI7aKSJS5o5kH2FY8jDu+OAO7rnhHlY9s0oJ3EPGWuudHRtjvbVvEfFN\n1lo+T/ucoclDuevau+hKL/7f2yvJzg4kKCiH2NhIn7pJsTcYY7DWmuKur3aKiJSJtD1pxMyJ4cCJ\nA3ze+3MOrc0lLi6JjIwxhetkZIwAqPCJvCTUThERrzqcfZghSUOI+DCC/wn/H1Y8vYLbrr6NyZOT\nz0jgABkZY4iPn+dQpL5JSVxEvMJayyc/fUL4lHAOZh0kbUAa0W2jCQzIbwBkZxfdCMjKqlSWYfo8\ntVNEpNSt3b2W6DnRHMk+wpf3fUn70PbnrBMUlFPktsHBud4Oz6+oEheRUpOZlcmguYO466O7eKDp\nAyx/anmRCRwgNjaSsLARZywLCxtOTEzXsgjVb3hciRtjXgD6AnnAWuBxa212aQUmIr7DWsvHP33M\nc/Ofo/v13UkbkMYfL//jBbf57eRlfPyLZGVVIjg4l5iYbjqpWUIeDTE0xjQEFgLh1tpsY8znQKK1\n9sPT1tEQQ5EKYM2uNUTPiebEqRNM6T6FWxvc6nRIPq2kQww9baccBk4BVYwxgUAVYIeH+xIRH3Qo\n6xCxc2KJ/DiSvs36sqzfMiVwB3iUxK21B4AJwC/ATuCQtXZ+aQYmIuVTns3jgx8/IHxKONk52aQP\nSOeZW56hUoBGlTjBo564MSYMGAg0BDKBL4wxD1trPzl9vdGjRxc+joiIICIiwtM4RaQcWP3ralyJ\nLnLycvj6wa9pU7+N0yH5vJSUFFJSUjze3tOe+ANAV2ttv4Ln/wu0s9a6TltHPXERP3HwxEFGLhzJ\nl+4vGXPnGJ5o9QQBRoPbvKGseuI/A+2MMZWNMQboAqR7uC8RKafybB4zVs8gfEo4eTYPt8tNv9b9\nlMDLEY/aKdbaNcaYj4AV5A8xXAVMK83ARMRZK3euxJWY/+U64aEEbr7qZocjkqJoFkMROcOBEwcY\nsWAEX/38FX+762881vIxVd5lqKzaKSLiZ/JsHtNXTid8SjiVAirhdrnV+/YBmjtFRFi+YzmuRBeB\nAYHMfXgureq1cjokKSYlcZEKbP/x/QxfMJyvN3zN2LvG8kiLR1R5+xj93xKpgHLzcnl3xbs0mdqE\noMAg3C63et8+SpW4SAWzbPsyXIkuKl9WmeS+ybSo28LpkOQSKImLVBB7j+3lhQUvkLgxkXFdxtG3\neV/yL/MQX6bvTiJ+Ljcvl6nLp9J0alNC/hCC2+Wm1vZr6NbtRSIiRhMVNZKEhFSnwxQPqRIX8WNL\nty0lek40l192OQseWUCzOs1ISEjVDYr9iC72EfFDe47t4fn5z5OUkcTrXV7noWYPFbZOoqJGkpz8\n6jnbREW9yNy5r5R1qHIWXewjUoHl5OXw1g9v0XRqU2oG18TtcvNw84fP6H3rBsX+Re0UET+x5Jcl\nuBJd1Kxck5RHU2h6ZdMi19MNiv2LKnERH7f76G4enfUoD3z5AM93fJ6Fjyw8bwIH3aDY36gSF/FR\nOXk5TF0+lVdSX+GxFo/hdrkJCQq56Ha6QbF/0YlNkTKWkJDK5MnJZGcHEhSUQ2xsZIkT6KKti4ie\nE03tKrWJvzueJn9s4qVopayV9MSmKnGRMnSpw/t2Hd3FsHnDSNmSwoTICdzX5D5dsFPBqScuUoYm\nT04+I4EDZGSMIT5+3gW3y8nLYeL3E2n2djPqh9TH7XJzf9P7lcBFlbhIWfJkeF/q1lRciS7qVq3L\noscX0bh2Y2+FJz5ISVykDJVkeN/OIzsZNm8Yi7Yu4o2oN/hL+F9Uecs51E4RKUPFGd53KvcUE76b\nQPO3m3NN9Wtwu9z0btJbCVyKpEpcpAxdbHhfypYUXIkuQquF8t2T33HDFTc4Ga74AA0xFCkHdhze\nwdB5Q/lu23dMjJpIr8a9VHlXUJo7RcSHnMw9yfgl42nxTgvCaobhdrn5c/iflcCl2NROEXHIgk0L\niJ4TTaMajVj65FKuv+J6p0MSH6QkLlLGth/ezuCkwfyw4wcmdZtEzxt7qvIWj6mdIlJGTuaeZNzi\ncbR8pyWNazcm3ZXOvY3vVQKXS6JKXKQMzMuYR8ycGK6rdR3L+i0jrFaY0yGJn1ASF/GiXzJ/YXDS\nYFb9uopJ3SZxz433OB2S+Bm1U0S8IDsnm78t+hut321NsyubkTYgTQlcvMLjStwYUwP4O9AUsMAT\n1trvSyswEV81979ziZ0TS/gfw1n+1HIa1WzkdEjixy6lnTIJSLTW9jbGBAKXl1JMIj5py6EtDEoa\nxNrda5l892S6X9/d6ZCkAvConWKMqQ7cbq2dAWCtzbHWZpZqZCI+Iisni1e+fYVbpt3CzfVuZt2A\ndUrgUmY8rcQbAXuNMe8DLYCVQJy19nipRSbiAxI2JBA3N45mdZqx4ukVNKzR0OmQpILxNIkHAq2B\naGvtcmPMROB54KXTVxo9enTh44iICCIiIjw8nEj5svngZgYmDSR9bzpvdX+Lbtd1czok8VEpKSmk\npKR4vL1HE2AZY+oCS621jQqedwSet9b+6bR1NAGW+J0Tp07w+pLXif8hnsHtBzOk/RCCAoOcDkv8\nSJncY9Nau8sYs80Yc4O1dgPQBUjzZF8ivmL2+tnEzY2jdb3WrHpmFVdXv9rpkEQ8n4rWGNOC/CGG\nfwAygMdPP7mpSlz8RcaBDOLmxrHxwEbi744nMizS6ZDEj5W0Etd84lJhJCSkMnlyMtnZgQQF5RAb\nG3nBO8wfP3Wc1xa/xtTlUxnaYSiD2g1S60S8rkzaKSK+JiEhlbi4pDPuNJ+RkX+btLMTubWWr9d/\nzcCkgbS5qg2rn1lNaPXQMo1XpLhUiUuFEBU1kuTkV4tY/iJz575S+Py/B/5L7JxYNh/aTPzd8XS5\ntktZhimiO/uIFCU7u+gvnVlZlYD81snIhSNp9/d2dG7YmTX91yiBi09QO0UqhKCgnKKXB+fwL/e/\nGJw0mPah7VnTfw31q9Uv4+hEPKd2ilQIRfXEQ1s+Q+1HVpMddIy37n6Lzo06OxihSD6NThE5j4SE\nVOLj53HsZC47wr5lf8M0Xur8IrG3xnJZpcucDk8EUBIXOS9rLTPdMxmSPISOV3dkfNfxXBVyldNh\niZxBQwxFivDzvp+JmRPDrqO7+KjXR3Rq2MnpkERKhUaniF87evIoz817jo4zOtLj+h6senqVErj4\nFVXi4pestXyR/gVDkocQ0TCCtc+upV5IPafDEil1SuLid9x73UTPiWbf8X18+j+fcvs1tzsdkojX\nqJ0ifuNI9hGGJQ/jjg/uoNeNvVj59EolcPF7SuLi86y1fLb2M8KnhLP3+F7WPbuOmFtjCAzQF03x\nf/qUi09L25NG9JxoDmUd4p/3/ZMOoR2cDkmkTKkSF590OPswQ5KG0PnDzvQO782Kp1YogUuFpCQu\nPsVayyc/fUL4lHAOZR1i3YB1uNq6qBRQyenQRByhdor4jLW71xI9J5qjJ48y8/6ZtGvQzumQRByn\nSlzKvcysTAbOHchdH93Fg00f5Id+PyiBixRQEpdyy1rLR2s+InxKOMdOHiPdlc6zbZ5V60TkNGqn\nSLm0ZtcaXIkusnOzmfXgLNrWb+t0SCLlkpK4lCuHsg7x0jcv8Xna57zS+RWebPWkKm+RC1A7RcqF\nPJvHBz9+QPiUcE7mniR9QDpP3/y0ErjIRagSF8et/nU1rkQXuTaX2X1mc8tVtzgdkojPUBIXxxw8\ncZCRC0fypftLxtw5hidaPUGA0ZdDkZLQX4yUuTybx4zVMwifEo7F4na56de6nxK4iAdUiUuZWvXr\nKgYkDMAYQ+LDibSu19rpkER8mpK4lIkDJw4wYsEIvvr5K8beNZZHWz6qylukFOivSLwqz+YxfeV0\nwqeEUymgEm6Xm8dbPa4ELlJKLqkSN8ZUAlYA262195ROSOIvlu9YjivRxWWVLiOpbxIt67Z0OiQR\nv3Op5VAckA7YUohF/MT+4/t5ZvYz9PxHT1xtXCx6fBE7Vh4mKmokERGjiYoaSUJCqtNhivgFjytx\nY0wDoDswBhhcahGJz8rNy+Xvq/7OSykv8WDTB3G73NQIrkFCQipxcUlkZIwpXDcjYwQAPXrc4VS4\nIn7hUtopbwLDgGqlFIv4sGXbl+FKdFH5ssok902mRd0Wha9Nnpx8RgIHyMgYQ3z8i0riIpfIoyRu\njPkTsMdau9oYE3G+9UaPHl34OCIigoiI864qPmrvsb28sOAFEjcmMq7LOPo274sx5ox1srOL/phl\nZemSepGUlBRSUlI83t7TSrwD0NMY0x0IBqoZYz6y1j5y+kqnJ3HxL7l5uby78l1Gp4zm4WYP43a5\nqR5cvch1g4JyilweHJzrzRBFfMLZBe7LL79cou09SuLW2uHAcABjTCdg6NkJXPzX0m1LiZ4TTdU/\nVGXBIwtoVqfZBdePjY0kI2PEGS2VsLDhxMR083aoIn6vtC720eiUCmDPsT08P/95kjKSGN91PH1u\n6nNO66Qov/W94+NfJCurEsHBucTEdFM/XKQUGGu9k3+NMdZb+5aylZOXwzsr3uHlb1/mkeaPMCpi\nFNWCdD5bxBuMMVhrL14dFdBl93JBS35ZgivRRc3KNUl5NIWmVzZ1OiQROY2SuBRp99Hd/HX+X1mw\naQH/F/l/PND0gWK1TkSkbGkCCzlDTl4Ok5dN5qa3b+LKKlfidrl58KYHlcBFyilV4lJo0dZFRM+J\npnaV2nz72Lc0+WMTp0MSkYtQEq8AEhJSmTw5mezsQIKCcoiNjTxjZMivR37lr/P/SsqWFCZETuC+\nJvep8hbxEUrifu5C85ZEdmvPWz+8xZhFY+jXuh9ul5uqf6jqVKgi4gENMfRzUVEjSU5+9Zzlbe57\nhBOdV1O3al3i746nce3GDkQnImfTEEM5wznzloTshK7D+CnsP3zc6T3+Ev4XtU5EfJhGp/i5wnlL\nAk5B+wnwbHPIvIaOa/rTu0lvJXARH6dK3M/Fxkay7ujD7Gz5IxwOhfe+I6zGBwya1MPp0ESkFKgn\n7se2H97O0OShLNyYQv21t1NtZxMqB+cRE9NV85aIlFPqiQsnc08y8fuJvL7kdfrf0p8ZQ2dQ5bIq\nToclIl6gJO5n5m+aT8ycGBrVaMTSJ5dy/RXXOx2SiHiRkrif2Ja5jcHJg1mxcwUToybS88aeOmkp\nUgFodIqPy87JZuyisbR8tyVNajchfUA69za+VwlcpIJQJe7DkjOSiZkTww1X3MAP/X4grFaY0yGJ\nSBlTEvdBv2T+wqCkQaz+dTWTuk3inhvvcTokEXGI2ik+JDsnmzGpY2j1biuaX9mctAFpSuAiFZwq\ncR8xZ+McYufG0uSPTVjx1Aoa1WzkdEgiUg4oiZdzWw5tYeDcgazbs47Jd0+m+/XdnQ5JRMoRJXEv\nutg83heSlZPF+CXjmbhsIoPaDeIfvf9BcGCwlyMWEV+jJO4lF5rH+2KJPGFDAnFz42hWpxkrn15J\nwxoNvRmqiPgwzZ3iJeebxzsq6kXmzn2lyG02HdzEwLkDce9zE393PN2u6+btMEWknCnp3CkaneIl\n58zjXSArq9I5y06cOsHolNG0md6Gdg3ase7ZdUrgIlIsaqd4SeE83mcJDs494/ns9bOJmxtH63qt\nWf3Maq6ufnVZhCcifkJJ3EtiYyPJyBhxRk88LGw4MTH5FXbGgQzi5sax8cBG3vnTO0SGRToVqoj4\nMPXEvSghIZX4+HlkZVUiODiXmJiudI68hdcWv8aU5VMY1mEYg9oNIigwyOlQRaScKGlPXEm8jFhr\n+ff6fzMoaRBtrmrDhMgJhFYPdTosESlnyuSmEMaYUOAj4ErAAtOstZM92VdFsHH/RmLnxrLl0Bam\n3zOdLtfWAc5bAAAJPUlEQVR2cTokEfETno5OOQUMstY2BdoBLmNMeOmF5R+OnTzGiAUjaP9ee+5s\neCdr+q9RAheRUuVRJW6t3QXsKnh81BjjBq4C3KUYm8+y1vLVz18xKGkQHUI7sKb/GupXq+90WCLi\nhy55dIoxpiHQClh2qfvyBxv2byBmTgzbD2/ng3s/oHOjzk6HJCJ+7JIu9jHGVAW+BOKstUdLJyTf\ndOzkMV6Y/wId3utAVFgUPz7zoxK4iHidx5W4MeYyYCbwsbV2VlHrjB49uvBxREQEERERnh6u3LLW\n8mX6lwxJHsId19zB2mfXUi+kntNhiYiPSElJISUlxePtPRpiaPJv4PghsN9aO+g86/j9EMOf9/1M\nzJwYdh3dxVt3v0Wnhp2cDklEfFxZzZ1yG9AX6GyMWV3wU2Em+zh68ijPzXuO29+/nR7X92DV06uU\nwEXEEZ6OTllMBZw8y1rLP9P+ydB5Q+ncsDNrn11L3ap1nQ5LRCowzZ1STOl704mZE8O+4/v47C+f\n0fHqjk6HJCJS8arpkjqSfYShyUPp9EEnet3Yi5VPr1QCF5FyQ0n8PKy1fLb2M8KnhLPv+D7WPbuO\nmFtjCAzQlxcRKT+UkYqwbs86ohOjyczO5J/3/ZMOoR2cDklEpEiqxE9zOPswg5MGc+eHd3Jfk/tY\n8dQKJXARKdeUxMlvnXz808eETwknMyuTdQPW4WrrolLAubdSExEpTyp8O+Wn3T8RnRjNsVPHmHn/\nTNo1aOd0SCIixVZhK/HMrEwGzh1Il4+60OemPvzQ7wclcBHxORUuiVtr+WjNR4RPCefYyWOku9J5\nts2zap2IiE/y23ZKQkIqkycnk50dSFBQDrGxkTS4uTquRBfZudnMenAWbeu3dTpMEZFL4pdJPCEh\nlbi4pN/vNB98iO9rdSHgx02M6/YaT7Z6UpW3iPgFv2ynTJ6cnJ/ATR60fB9c4Rw+egutlj7O0zc/\nrQQuIn7DLyvx7OxAqLsaerjA5MJns2HnLeR1Gu10aCIipcqrlXhU1EgSElK9eYhzHDxxkE2N/wN9\n74bVT8B7S2HnLQAEB+eWaSwiIt7m1SSenPwqcXFJZZLI82we7616j/Ap4TRvHkrDhIdgVT+w+b9i\nWNhwYmK6ej0OEZGy5NGdfYq1Y2Ms5O87KupF5s59xSvHAVi5cyWuRBfGGKZ0n0Lreq1JSEglPn4e\nWVmVCA7OJSamKz163OG1GERESkNJ7+xTJj3xrCzvnEjcf3w/IxaOYNbPsxh711gebfkoASa/8u7R\n4w4lbRHxe2UyOqW0e9F5No/pK6fTZGoTLgu4DLfLzeOtHi9M4CIiFYXXK/H8XnTp3X5z+Y7luBJd\n/KHSH0jqm0TLui1Lbd8iIr7Gq0k8KupFYmK6lUpbY9/xfQxfMJzZG2Yzrss4/rf5/2JMsdtGIiJ+\nyasnNktj37l5uUxfNZ2XvnmJPjf14eXOL1MjuEYpRCgiUv6UyxObnlq2fRmuRBeVL6vM/Efm07xO\nc6dDEhEpV8plEt97bC8vLHiBxI2JvN71dR5u9rBaJyIiRShXwzly83KZunwqTac2pVpQNX6O/pm+\nzfsqgYuInEe5qcSXbluKK9FFSFAICx9dyE1X3uR0SCIi5Z7jSXzPsT08N/85kjOSGd91PH1u6qPK\nW0SkmBxrp+Tk5RC/LJ6mU5tyReUrcLvcPNTsISVwEZEScKQSX/LLElyJLmpVrsW3j31Lkz82cSIM\nERGf53ESN8Z0AyYClYC/W2vHXWyb3Ud389f5f2Xh5oX8X9f/4/6m96vyFhG5BB61U4wxlYC3gG5A\nE6CPMSb8fOvn5OUw6ftJ3PT2TdS5vA7pA9J54KYHKkwCT0lJcTqEckPvxe/0XvxO74XnPO2JtwX+\na63dYq09BfwDuLeoFVO3ptL63dbM3jCb1MdSeb3r64QEhXgar0/SB/R3ei9+p/fid3ovPOdpO6U+\nsO2059uBW89eqe+/+pK6NZUJkRPo3aR3ham8RUTKiqdJvFiTooRWCyXdlU7VP1T18DAiInIhHk2A\nZYxpB4y21nYreP4CkHf6yc38O/uIiEhJlWQCLE+TeCCwHrgL2An8APSx1rpLvDMREfGYR+0Ua22O\nMSYaSCJ/iOF7SuAiImXPa/OJi4iI93nlsntjTDdjzM/GmI3GmOe8cQxfYIwJNcZ8Y4xJM8asM8bE\nOh2Tk4wxlYwxq40xs52OxUnGmBrGmC+NMW5jTHrBOaYKyRjzQsHfx1pjzKfGmCCnYyorxpgZxpjd\nxpi1py2rZYyZZ4zZYIxJNsZc9A44pZ7ES3ohkJ87BQyy1jYF2gGuCvxeAMQB6RRzdJMfmwQkWmvD\ngeZAhWxFGmMaAk8Bra21zchvzT7oZExl7H3y8+TpngfmWWtvABYUPL8gb1Tixb4QyN9Za3dZa38s\neHyU/D/Wq5yNyhnGmAZAd+DvQIW9YMAYUx243Vo7A/LPL1lrMx0OyymHyS90qhQMlqgC7HA2pLJj\nrV0EHDxrcU/gw4LHHwK9LrYfbyTxoi4Equ+F4/iUgqqjFbDM2Ugc8yYwDMhzOhCHNQL2GmPeN8as\nMsZMN8ZUcTooJ1hrDwATgF/IH+V2yFo739moHFfHWru74PFuoM7FNvBGEq/oX5XPYYypCnwJxBVU\n5BWKMeZPwB5r7WoqcBVeIBBoDUy11rYGjlGMr8z+yBgTBgwEGpL/DbWqMeZhR4MqRwruNH/RfOqN\nJL4DCD3teSj51XiFZIy5DJgJfGytneV0PA7pAPQ0xmwGPgPuNMZ85HBMTtkObLfWLi94/iX5Sb0i\nugX4zlq731qbA/yL/M9KRbbbGFMXwBhTD9hzsQ28kcRXANcbYxoaY/4APAB87YXjlHsmf7KY94B0\na+1Ep+NxirV2uLU21FrbiPwTVwuttY84HZcTrLW7gG3GmBsKFnUB0hwMyUk/A+2MMZUL/la6kH/i\nuyL7Gni04PGjwEULv1K/KYQuBDrDbUBf4CdjzOqCZS9Ya+c6GFN5UNFbbjHAJwVFTgbwuMPxOMJa\nu6bgG9kK8s+VrAKmORtV2THGfAZ0AmobY7YBLwGvAf80xjwJbAHuv+h+dLGPiIjvcuwemyIicumU\nxEVEfJiSuIiID1MSFxHxYUriIiI+TElcRMSHKYmLiPgwJXERER/2/wHJWT/6+T6MygAAAABJRU5E\nrkJggg==\n", "text/plain": [ "<matplotlib.figure.Figure at 0x107134190>" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# (c) 2015-10-24 Teruhisa Okada\n", "\n", "% matplotlib inline\n", "\n", "import numpy as np\n", "import matplotlib.pyplot as plt\n", "\n", "import random\n", "from scipy.optimize import leastsq\n", "\n", "\n", "def residual(param, x, y):\n", " a = param[0]\n", " return y - (a * x)\n", "\n", "def rsq(infodict, y):\n", " ssErr = (infodict['fvec']**2).sum()\n", " ssTot = ((y-y.mean())**2).sum()\n", " return 1-(ssErr/ssTot )\n", "\n", "def plot(param0, x, y):\n", " a, cov, infodict, mesg, ierr = leastsq(residual, param0, args=(x, y), full_output=True)\n", " r2 = rsq(infodict, y)\n", " print a, r2\n", " \n", " plt.plot(x, y, 'o')\n", " l = np.arange(0,11)\n", " plt.plot(l, a[0] * l, '-', label='a={0:.2f}, R**2={1:.2f}'.format(a[0], r2))\n", " plt.legend()\n", "\n", "\n", "param0 = np.array([1])\n", "x = np.linspace(0,10, 11)\n", "y = np.array([x[i] + random.random() for i in range(11)])\n", "\n", "plot(param0, x, y)" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": true }, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 2", "language": "python", "name": "python2" }, "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 }
Teruhisa Okada
Not found.