{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Plot a fit over a density plot\n", "Use the ``densityreg`` plot method to combine fits and their raw data.\n", "\n", "This example uses artifical data to show you the capabilities of the ``densityreg`` plot method." ] }, { "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": {}, "source": [ "First we define our data which comes from multiple realizations of the underlying equation ``sin(x)``" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "\n", "Dimensions: (experiment: 1500)\n", "Dimensions without coordinates: experiment\n", "Data variables:\n", " y (experiment) float64 0.375 -0.474 0.4116 1.099 1.278 0.2544 ...\n", " x (experiment) float64 -5.01 -4.805 -4.601 -4.396 -4.192 -3.987 ..." ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" } ], "source": [ "all_x = []\n", "all_y = []\n", "for i in range(30):\n", " deviation = np.abs(np.random.normal())\n", " all_x.append(np.linspace(-np.pi - deviation, np.pi + deviation))\n", " all_y.append(np.sin(all_x[-1]) + np.random.normal(scale=0.5, size=all_x[-1].size))\n", "x = np.concatenate(all_x)\n", "y = np.concatenate(all_y)\n", "ds = xr.Dataset({'x': xr.Variable(('experiment', ), x),\n", " 'y': xr.Variable(('experiment', ), y)})\n", "ds" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "This dataset now contains the two variables ``x`` and ``y``. A scatter plot of the data looks like" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "psyplot.project.Project([arr0: psyplot.data.InteractiveList([ arr0: 1-dim DataArray of y, with (experiment)=(1500,), ])])" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXUAAAD9CAYAAABDaefJAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3X1wHOd9H/Dv3uF4EO4AkEegkIgrKTmUngps/MaUMC1X\nEU0pHRsyHSVRG1X1TCXSw6Sa6diNIjd2TdNMnJpjx0nbKK5Mk45HUjoJncqlgjiWLMGuwmEwNSWl\nNiEvw1gkA5AFAUIk8f521z8Oezoc9tn33duX7+cfCYc73OKW+O2zv+f3/B6lXC6DiIjiIdXoAyAi\nIu8wqBMRxQiDOhFRjDCoExHFCIM6EVGMMKgTEcVIU1BvNDY2ydpJIiKbOjtbFTvP50idiChGGNSJ\niGKEQZ2IKEYY1ImIYoRBnYgoRhjUiYhihEGdiChGGNSJiGIksMVHRHEzODSK/lPncWl8Bps6WtC3\n81b09nQ1+rCohp/naGFhAQ8//Cs4fvyE7vdPnHgOfX17kE6nPXk/qxjUiRwYHBrFUyfOVL8eHpuu\nfq0XNHgBCJ7dc2RXZYMh+WLPp5/+Bj70ofsZ1ImioP/UecnjF6oBQwvkI2PTqO2R4XVwIX1WzpFd\ns7OzOHToP2FychLd3UUAwOuvv4pvfOMIyuUyZmdn8LnPfQGvv/4qrl69is997tP4nd85jC996Xdx\n5coVXL06jg984G7s2/drDn8rcwzqRA5cGp/Rffzy1WkAa0eJerTgwlG8P8zOkRPf/vaf4x3v2IqP\nf/zXMTT0Y7z66g9x/vxPceDAb2Pjxg48/fQ3MDDwPXzsY4/gm988hkOH/jOuXBnFtm0/i0996qNY\nWFjAL/3ShxnUicJmU0cLhsfWBodbNuYAyEeJtS5fnfY9RZBkZufIiX/4hwt4//v/OQCgp+efIp1u\nQkdHJ37/97+ElpYWjI1dwTvf+e6VZ5dRLpfR1taGN944g9de+yFuuimHxcVFx+9vBatfiBzo23mr\n5PEtAOSjxFq3bMwZpgjIHbNz5MStt74DP/7x/wUAnD37EywvL+Hw4S/gM585iE9/+nPo6OhcybUD\nqVQKpdIy/vIvn0draxs++9nfxq/+6sOYm5tz/P5WpA8ePOjrG2hmZhaCeSOiABQ787i50ILRiVlM\nzy2iuyOPh+69vTq6/qF6BTdmjEdkD917O37w+iXo9aS+MbOA0+oVtDRnUOzM+/AbxJ/ZOXLi9tsF\n+vtP4Fvf+lNcvnwZ4+NjuOee3fj617+KF174DjKZJiwvL+Ouu+7G2bM/wXPPfQu//Mv/Cs888018\n73vfxeDgKTQ1ZXD33bvQ0tJi6T1zuezn7Ryjol1V/MZ+6pQkspx6SgE2deTRt3MLenu6cODooG6K\noNb+PduYikkwu/3UmVMn8oEWhPtPXcDlq9O4ZWOuGshr9e281fKEKpEVjoK6EKIJwDEAtwJYB+AL\nqqo+7+FxEUVeb0+XaTCuDf7DY1O6z3FTrUHJ43Si9N8AGFdV9W4AHwbwh94dElGy9PZ04dDeHSh2\n6ldluKnWoORxGtT/DMBnV/5fAeBvjQ5RAvhRrUHJ4yj9oqrqDAAIIVoBHAfwGS8PiijKnC4mspqH\nJzLiuPpFCPGPAfxPAH+oquo3zZ7P6pfo4ApH+2QtATQKgO7OHD9Lss1u9YujoC6E6AIwAOAxVVUH\nrLyGQT14ToKzrBSPZXVyVloC1OJnSXbYDepOc+q/BWA9gM8KIQaEEC8LIbIOfxb5QAs0w2PTKJXL\n1eXng0Ojhq/jCkf7rLQEWP18fpbkH6c59U8A+ITHx0Iectqhzo8mSHE2ODRqunio3vDYFAaHRjla\nJ1+w90tMOQ3Omzr0ly6zrG4tu2mXWlbumoicYFCPKafBmWV11tlNu6x9PdMw5D22CYgp2fJzs+DM\nsjprnKRd6g2PTWHvF19GU1rBUqmM7o4cxOYNUC++xcojcowNvWKsUv2yNjizZNEdu2mXTDqFxeWS\n4/djtUyysaEXVen1HuGmDO7ZTbsslZwH9Mr7saEXWcecesKwZNE9Kxtg1OruyGP/nm0oduaRTtka\ndAFg5RHZw5F6wrBk0T3ZNmmFtiwmbsyveVxLe2mjbSs91Gux8ojs4Eg9YViy6Nzg0CgOHB3EiCQg\nv+f2Tks/R1ZhJH8+K4/IOo7UI87upKfTqpiks7KTkZUFX9r5UhQgpShYLsnrB4qdeVYekW0M6hHm\nZNKTJYvOyAL2po48Du3dYVjiqKW26s/XskHlGSteyCkG9Qhz2grAyo48tJrRXISVEsfHnzyJicm1\n+XZNJp3CUqmE7g6OzskdBvUI82rSk3Xr5mSToylFMQ3oy6WyYUAHUK1jZ0AntzhRGmFeTHo67eaY\nNLLJTTeLivSwtJTc4kg9wryY9HSawkma+rmIlKJ4HtCB1fl33j2REwzqEebFpCfr1q2rnYvYd9jS\n3jC23bIxx1W/5AqDesS5nfSU5YpZt25M9rm5NTI2ha89r5+j590TWcGcesxoC2T2HR7AgaODprlx\nttp1xu4CIplMOoXaxgFlALJKx+GxKUvnlJKNXRpjxOn+orJujmRscGgUx/rf8CW3boZ17MkRyMbT\nTjCo+6sSYIawuLz2Yy60ZdGSbeKkmwv1E5da3/ORsWk04h92sbOy6Inij613E8hs8cvEjXlMoFIn\nzUk367RAPjI+vSolMjw27Us+3Q5OZJMMg3oMONlW7akTZ9B/6jxH7RJu9h8NAieySYYTpTFgt7+3\nhguN5NzuP+o3TmSTDEfqMSArr8ukU2htyZguUT/W/wYAcKu7Gk4vlH7Yvb0I9eI1TmSTJemDBw8G\n8kYzMwvBvFECtTRncFodW/P4xz/Sg223FXS/V6tULuO0Ooap2UUc//7f48bMIsoAbsws4rQ6hpsL\nLSh25n06+nD6oXoFN2YWG30YAIDFpRIO7d2Brg0t+Lvha/jB65fwQ/UKWpoziTsvSZTLZT9v5/kM\n6jFQ7Mzj5kILRidmMT23iO6OPB6693b09nSh2JnHK397CbMLy6Y/Z/jKFEo61VCjE7PY9d5uPw49\ntGQXyka4MbOAmwsteOrEmTUX3Ff+9hLa81kG9xizG9SZU4+V8kqVxurA/J47rO3II6u3ZqVF4z3z\ngqr7+MTkPOdFaBXm1GPAqFcIALx0etjVz09ipUXYJkqn55YMv88WAqRhUI8Bo06L9aN2I4rk2Ums\ntAjTRKkVvJsiDdMvMTAyrv8HfWl8ylZw6u7MY/+ebSh25pFOKSiufJ3EEaCsV31YJfFuivQxqMdA\nU0p/FXE6lbIVnIbHpnB84Bz6dm7Bvvt7AJRx5PmhRDaR8qphV1DE5vWNPgQKCaZfIkavjnxJp98L\nUJn4bM9nbS1p1ybeaiW1tUCuuUk3l51OKfgnWzbgzJsTDTgqfS+dHsbW7vZEnR/Sx4ZeESJbul5o\nzZouMPJCUppIhalFgGyeQ09Szk/S2G3oxfRLhEgrMmydcueSMhkXpsqXMoBttxVWzXMokvOdlPND\nxph+iRDZpOf1qQXs37PN99FlUibjwlb5cubNiVUT1o8/eVL3zqw9tw4Hjg4mvsVD0jGoR4jR1nO9\nPV3oP3Xe15awM3OL2Hd4ILIBw2pfG7+2qnNDq0MfHBqVptomJuer30vqPAgxpx56tYFofX6d7h+0\nNorzKxdcaMti4ob8faPAaFcoAGs2wHC7YKtWsTPn+iKRUoBNHfKfk04pWC6t/RNjnj36uElGjNQH\nIi2gF9qyuD61sKZjn/ZfLwP7/j3b0H/qfHWTjVpRWsUoy5Mf//65VRcsbQOM3duLeO3smCcT0KMT\n7tM5pTIMLwx6AR1gnj2JXAV1IUQvgC+qqrrLo+OhGrJA1JLN4Mv/7q41j2ujeq/s3l5Eb08Xjjw/\npPv9KAUMWZ5c7w4EANSL1/Dgrq2eXCD1thgMSnueefakcRzUhRC/CeBjAKa8O5x4s9urXBaI9IKp\nH6mXl04P47WzY7qdG4HK6PDA0cFIBAq7efLLV6dDVQXjFLcyTB43I/VzAB4A8LRHxxJrRk23ZH9g\nRhOj2s/ULhJpn4pTzdIPUQkUfTtvtVXjn1IUafuFqJDl2aOUNiP7HIcCVVWfA2DcOo6qjJtu6ZMt\nVe/buaV6kRgem0apXJbe4ks6CHjO6PcIg96eLt2+Ng/u2qr7/MXlEgKqIfAN8+zJxInSgNhJpWi0\n0VT/qQtrtjI7cHTQ0vtu6sijb+cW9J+6gOEx/zJlUQgUvT1d0hHqsf43pP3k9chGwVGQlPUGSeVF\nUA9oLBhtZqkUGVkgsrpAZmZuEUeeH8Kmjhbs3l70tFSvVpQDRW9PF772vP58hIJK98pL41OojeFR\nDehAMlspJ4kXQT26/7oDJMvpOv0DszrxV7sYRSvVUy9e83zUHtZAYXVyuiml6KawyqjchUT5H3km\nnUKpXOam1QnhKqirqnoBwPs9OpZYM0qlOCG7SGgy6ZRuOkG9eA2H9u7A3sMvW8oZpxUFpeo2efrv\n82jfnaEMFHYmp2WdLoHojMrTioJlnRMV1vND/mBOPUBGOV0nPwsAjg+c063eWDLYb3RwaNTyJKBe\nkKgV5oBhNDldf8zdHqz69IvV/H0Z5ZXFYt4MHCiaGNQjTC+wF1qzeM8dndLceXtunSf11+mUgn33\n94Q6YNiZnDa782kkq3cKmzryng4cKJoY1CNMb8HRxOQ8Xjs7Jn1NbdMnN4JISdhdrFX/mnQKKC2v\nfY7epK4fLRaCVrv7kZPPjuKB/dQjTDbiDmLDDKASAP3a6q6+Dl/Lhxu9l9Xa/ZHxKelxZ9LRLeZS\nL14D4Oyzo/hgUI+wMPT99itgOFmsJXtNJp1atbFEubz2uLVA2Mg+LW5paSUnnx3FB4N6hMk2lS60\nZQM+Eu8DhpPFWrLXlMpldHfo19Frxy0LhGnZNkMhlFIUDA6NOvrsKD4Y1CNM1kbgPbd3otAabGAf\nHpvydLQuu2AZLXIyeo1ZoJN9f7lcRq45GlNPi8slPHXiDNbn1+l+P8oLxMg6BvWIqw3ehbZsddVo\nUHn1Wl6mYYz63jh5jVHAHxwaNWyINj0XrRZH0/P6x1s7kUrxlT548GAgbzQzsxDMGyWElgOeXXi7\nvGN2fhnXp+ZXPaZJp5RAGlSNTsxi13u7Xf+cYmceNxdaMDoxi+m5RXR35PHQvbcbVnAYvaalOYPT\n6tqqoH925z/C8e//PSKyvsiSZcm8wJuXb+C0egUtzRkUO/MBHxU5lctlP2/n+dG4r6Q17Fa+BLUq\n0su8rZOaa9lrZCt649Az3Y6otEom5xjUIyoMlS96wpy31Qv4skZeccee6vHFjacj6sDRQd1l7bJN\nH4JSaRj2VmgXvWiLckbGp6VNvJJA28g6rOeJ3mZ342lOlEaUbFLwwV1bA618yaRT1U0ntEnasC56\nqV2UUy43du/QRtM2sg7jeSJ3GNQjSraTT29Pl3Q3Hz+UymUceWIXDu3dAfXiW7rPCcuil+MD5xp9\nCKEWlvNE7jCnHmF2JgXF5vW+bJBRm0MP86KXwaHRhqaloiAM54ncY1CPKb2Av7W73fNt7bTaZ63W\n22oDraAlrcrFiTCcJ3KPQT3i7HTj0wL9vsMDKHk0Qf7a2TFs7W437G4Yhl2RwlotFCZhOE/kHoN6\nyBkFbTs7+9SyuhWeFROT84aNtBqxiYbeZ+bl7xwXhdYsrk8vcDONmGFQDzGzoG1nZ59aXm8IIRsF\nL5VKDQnoep9Z/qaMq5+rKAhkRW4QwrwFIbnH6pcQM2uhKgumw2NTePzJk9h3eEC3b3hvT5dnZY+F\ntqy0r0q5jMDL5GSf2dTsIgDnXRfjEtCBSsUSA3p8MaiHmFHQPnB0UNqND6ikRYxqkL0qe3zwnq3S\nmnkg2DK5waFR0xSL2Z6rScAJ0Xhj+iXEjPLAdvPD9SkZ7f+feUF11YXwmRdUzM7rlLysCKpMTm9r\nP9JndUKUW+JFE4N6yNT+IRmNxDVW2wLIgqvbtrJmrw9qVMiSRW85nYSnxmNQb7D6IF4boLX/Nwrc\n16cXUOzMmY7c9YJrEIEwqDI5lixaZ6WZl9NJeGo8BvUGqh8NyQJ3S3MGLc1NuoFbK0czSz3oBVe/\nA+Hu7UVXAcDO7b9ZyWJKQax6prtRf9em9zmHeXUwGeNEaQNZHSlfvjptuKtPfR+YQmsWhbbsmp4w\n9WRVK17Rdrd3orb5lpWmU0aTtQCgRGivUb/V3rXJPmduiRddHKk3kNWR8i0bc9JNHrTHzTaUqG85\nu7RcxgafuzkOj01h3+EBS5Ns9aPFGUmuXnb7rz0mu2MJapOQdEoJ7L2cqr1rs5uC46rT8ONIvYGs\njpRX/yGVV2qmrQcOvZazZcjTPV6qHf09++JZ0+PTni87NqPb/96eLhQ7GzuSDHNAr71rGxwalfbk\nBypzNbIuoBRuHKk3kGxlZ6Eti+tTq5dvu6lGCEtlyEunh7G1u33N8do5PrPbf69Xy8aJ7N+SHu3u\nkEE8ehjUG8gspVLLaTWClQU5QdI7XjsTtiMrC6/qe+DUpm4quy9dw+Wr00gpChaXS17+CpGlffZW\nLqJMs0QXg3qDWR0NOalG8HtBjpN+KCPja9v+yipXCq1ZtDRncGl8qlq5UsbquxQAa+5ghsemq6mC\nZ18860sf+SjS/q0YXUSLnXk294o45tQjQpZ/N0pH+J12cbLiXq8fjNHWfIf27sCmDv3fsf/UBcM7\nmMGhUQb0Gtq/Fdm/pWJnHof27mBAjzgG9YgwKmmUCeuCnPp+MEZb8wHGdylG3wvLXEJYzMwtYt/h\nAWllEVMu8aCUA2pwNDY2Gd6ygIio5I7N8+8ao+qGRit25iz3EpH9Hpl0Ctl16WoHxlqFtizempyP\nVXdFp2QrktlPPRo6O1ttLbJgUI8Zo7YDYWSlTM7J3EDupiZMz7rraxN3WrqFwo1BPcFkwU8bkbXn\nrAX5TDoVWMVIJp3CcqmMTR0tEJs3QL34lnSXp2P9b8SikqXQlsX1yXksN/gvQlGAo5/6YGMPgkwF\nFtSFEAqAPwLwLgBzAPapqvpT2fMZ1P0nS1PUjsgGh0bx7ItndVMWmnRKQbkMW/uY+rUzUO1I3su9\nVRtl9/ZiqCZvc81NyGbSuDa1EHh7Xbb2tcZuUHczUfqLALKqqr4fwG8B+IqLn0UesFr2aBTQASCl\nKLaDp1+xtnZS1e9eNV5TFKyZ/FUvvtXow1plem7JdEMVP9jt7UPWuQnqHwDwVwCgquoggJ/z5IjI\nMStlj1YqQtykONIpbxtn1V6QZBVAu7cXUezMe/q+Vhn9vt0dlTukI0/sQt/OLeg/dd504trrz8+J\nIHarMtuqkZxzE9TbAFyv+XpJCMESyQayUvZoVOaYSbs/faVy2bD/SiadQjqlWH6v2gvSmm6UbVkU\nWrMYeHUEM3PGdx9+Mer1on3utaNSmUw6hf17tmHf/T2eH6NdQbTXZWtf/7j5K74BoLX2Z6mqGv1Z\nrAgzq/cGjBeeeNGMyiwNUyqXceSJXXi0705LP6++drq3pwuH9u7Avvt7MHFjvpo6CFOVT6E1u+pz\nt3J39GjfndXVxfv3bPP3AE2kFMX3NIiTxXRkjZs2AScB3A/gW0KI9wH4kTeHRG6YtR2QNbyymh6w\nwuhntOfWVY/z+MA502B8fOAcjg+cWzWRBwDH+odcH6fXtAlpbQLwyPND2NTRgpFx+eehtyxf68/S\nqDUGi8sl37euM/p3SO64CerPAbhPCHFy5etHPDge8plZEzG/OxxOTM5jcGgUvT1deHDXVtP3qw36\n9T1fwuby1WndbpoyRnXiQXabVKDfyNnPrevsNLMje1inTqvUrlr1q8NhfYllmAO1HZXJ2rLlEbbe\nwqs1G48rwLXJ+YZsxZdOKTjyxK7g35hWCbKkkWJIy1kfeWKXbxs+DK+0z9VG7I3e2MIrfTu3SCcA\n9cob9QJ6bZnfxOQ8Jm7M4+Mf2YaCT7tUpQ22+WN+O5rYejem3C7sGBwaRToFlJb9Ob7aVEqcNraQ\ntRHWyhuNyCZUj/W/gbvfvcmXRUvLBnfqzG9HE4N6DDnZJalRPWP6T12oBrv+UxcwMj4V2SZc/acu\nuJoAlI3yF5dLeOn0cKCrUQttWduDAK4ODQcG9RgyGvFpFRn1fVVqA1GQ5YHDY1N4/MnKXPu1qQV0\nd+QgNm/Aa2fHqsdRaVtQxqaOfDU49p+6gOGxtRtuNNLlq9OuJgBlo3yNevEaip25QKpiHrxnq+Xn\nutlqkbzHoB5DRiM+YO0fXaP7jtdXuNQHLS23PzO3WL0oic0bQhfUtRy00709zdJQl69OY9/9PbrP\nUZRKfbkX8yApxV4wdrrVIvmDE6UxZLVHirYkO6ybadSr7VESpqZYGrc5aG3hkWy1rbYZ9O7txepz\nMukUdm8v4uinPuhZ2mpTh72WC1wdGi4cqceQ1YlH7Y/O7LafjGXSqeqKULeM1gv07dyyZos+Ld/+\nN2f+n2cdLO1enGT/fuxWzzAv7w2O1GOovl2A0cgPkPeMKbRmYVDxRiuWSiVPg49RuwdZqmNaskVd\nEJxstViPXRu9w5F6TNXmdWULfLQ/OrPJvcefPBmq3ipG8jdlUEY50F2PmlLej41kefkgUmV2c+Fe\nrA5lXt47DOoJUP9H155fB5SBI88Pof/U+eptruyPx8py/kbYvb24qkoGMO8V74flUnB97IJIlTnJ\nhTudHNYwL+8dBvWE0P7o7JSf1eY40ylvKiu8UGjL4paNOfzv10ew2Og94WB/YtENrxZq5ZqboCiK\n7kWwEStJvcrLE4N64li9zV2Tsml87KyaX1jGmTcnGn0YVX6uvBwcGl3VzbLQmsW22wquf//sujQm\nbuin1MTm9ThwdNDShKVXk5vs2ugdBvWEsXqb2+jadSONnBSsV2i1t/LSqvpgrpmYrPSQd7snrCyg\n55qbVlXXmN3JebXoiF0bvcOgnjBWb3OjUrvut6a0gr19PfjaiTO6NytvTXk/gWylc6VfrRRkF0y9\nCUuvJzfd5uWpgkE9wpzc+lq9zQ2ydj2lYFVrWa0twE3NTYFWsegplyvBSxZD/ah8CeNd0sj42tW7\nnNwMJwb1iHJ662v1NjfIzon186/ahGyjAzpQORaji5sflS9+3SUpCrChNStNvRgpl4G9X3wZG1Za\nAF+bWpB28eTkZmMxqEeUm1tfK7e5vT1d+OPv/ATziz713o0JPypf/LpLyjVnHAV0TRmr+/TI2jJz\ncrOxuKI0ooK49WVAN+dHAJOt0HQriBr+3duLzIs3GIN6RHE39saS7V7kBb02Abu3Fz1/Hz+oF681\n+hASj+mXiAqirrfQmo1Me4Ag7d5exMP33eHre+ilyLZ2t1fnQsKyEKweJ0kbj0E9oryq6zWqoAlr\ne4BGa9RotDbQHzg66HnePdfchLmFZVcXDN4pNh6DeoS5res1q6CpvXCEbUMKP2UzKSwtl6XBLQyj\nUT+qkza0NqNv5xYc6x+Stl8wW/TESdLGY049wYwqaDS9PV04tHdHYlrwbrutgK/+xj048sQuFDv1\nR51hGI2abajhxPDYFJ46ccawn05U949NEo7UE8xqBc3g0CiaUkoommf57cybE3j2xbNruj/WCsto\n1GhDDS9l0im0tmQsza+wVW7jMagnmKweuj2/rtrQaX1+XeImS2Vb5RXasnjwnq2JC1qLyyXL/wbC\nkJpKOqZfEkxWDz1xY766A03SArqRlmwmdAE9bC0FwpCaSjqO1BNMr4JmZm6RgVwijKNQOy0Fcs1N\nhi13vRCW1FSSMagnXH0Fzb7DAw08mnAL4yjUTkuBmfkl/LdP3G2pC6QVtY3YCq1ZPLjLODXFjaWD\nwaBOq9gJEikFWJ+PxwKllFLp4yI2r5fm1MM4ChWbN1g+X+VyJbBqgfTrfzHkqCa92JnHob07bL3G\ny97rZIw5dVrFTt+RTR15PLhrK7KZtH8HFJBNHZVA9fB9d2D/nm0otGWr3yu0Zn1rCeDG4NCo9AIk\nc6z/jWpgd7rIyMnFzUr5LHmDI3VaRS/PPjI2pdtPXKtrjoPhsSkcODpougl3mDiZJF1cLtk+Z5l0\nCqVy2dVuROy9HhwGdVplbd5zC/pPnQ9sw4xGilpKwGiSNK0oWDZYKfTMCyrSCmBl6UFrS8Y0X26G\nG0sHh+kXqtLynlo5oxbkRhIQ0GtFJSUg69RZ7Mxj30d6DF87PbdkKaADlR7qT504g8GhUbuHWCVL\n64VxniLqGNSpSnY7H/91pKtFJSVgFigLrVnd7zvl5mKn1044jPMUccD0C1UlabPpdEpB87q07kbL\nUUkJyDp1nhu5bnsC1Qq3F7uozFVEHYM6VQW52XQjaTXVgH7flCilBGoD5eDQKI4PnPOtxDQqF7uk\nU8oBtV0bG5tM2l185Hi1KCUq9u/ZBsB9T/owCOLcFVqzuDa1wIVDAevsbLXVI9VVTl0I8YAQ4lk3\nP4PCozbvKWu1W2jLIi5deN/OEZdXWspGd9xxfOCc7+8xMTm/agLdzcQp+cfxSF0I8QcAfgHA66qq\n/muz53OkHj2V8sa1o1g/dt0Jk6hN4DXqDsvJylKyz+5I3U1O/SSA5wDsd/EzKMRkE1tRmVAttGXR\nkm2yfQGKUk/wRqbMolIllDSmQV0I8SiAT6Jyb6qs/PcRVVWPCyF+3ufjoxCKyoTqxI15TMD+pGGU\ngpXZqtK0oqCpKYX5xWXP3zulKKt6yVA4mAZ1VVWPATgWwLFQRPixP6aX0oqCUrnsOEMepSoPs7um\ndh83OaltOcDAHh5cfES29fZ0hXrP0mUXAR2IVkmjbFWpJogOmlFZgZsUrFMnRzbEpOWuRlGA7o58\n5EoavbhrMusTY2ZkfEr3cfZPbwxXQV1V1R8A+IFHx0IRMTg0GquAvnt7EQ/fd0ejD8MRN5tPa7+3\n241RtD7tAKpBvH5v26g1S4sypl/INtnkXK45/Dd+u7cX1/QfiWpA1/T2dKHYaW8eIH9TBlu72wGY\np3CsOP79c6uawcku+kzV+C/8f4UUOrLJubmFZRQ7c6GsjLGy3VqU2U3DTM0uVp/vRQrH6r6nUaos\niiqO1MndqgH1AAALXUlEQVQ22ciuskDp1mAPxgJFAb782F2xDeiAfhfE2ruSTFr/T12ryddaJvgt\nSpVFUcWROtkmG9nVTjI63f/SD90d+UYfQoDebnmwtbu9mlqS5c2Hx6bw+JMnXb9rodXaxHmUKoui\niiN1ss2sN3ZvTxcC6hNnSRICiWyDE20C0yhvPjE573jiO5NOYf+ebdWul/UKbVn2Tw8YR+rkiFlv\n7DCtOn3qxBn0nzof65I6o42de3u6XOXNjUbhpXJ51Wcah46XUceROnlqcGgUB44OYmQ8HAFdE/fO\ngmYbO9feXVmlja6//Nhd0uoarVXA26Lf8TLq2E+dXNMWmYyMTYf+TzmunQVlnTMz6RQe7btz1YjZ\nSpfNYme+uum4Xt15vW23FXDmzYk1jzPl4p7dLo0M6uRK1DbWSKcUHHliV+xWO5qdh8qK2berk8zO\nmSxIp1OKrQnwuF5EgxToJhlEZl0Cw+aWjTnTScUo0tIrstLFcnn1qk6zVMzZi9d0H7db0cS69OAx\nqJMrUemtrtFSCnqivtqxt6fLUtDVJk8P7d2BlKQz2+JyyZNjYl168BjUyRUvlph7adttBaRT+oGq\n0JZFb0+X6aRilFk5H7W/p+z5shF/oTVr63hm5hax7/AADhwdjPSdUJQwqJMrYVtBeubNCelo9cF7\nKrXURitio87K+aj9PWXPv/vdm3Qff3DXVuzeXtT9Xu0KVi34c1/T4KUPHjwYyBvNzCwE80YUqGJn\nHjcXWjA6MYvpuUV0d+QxObvQ6MNao9CaxaN9dwIAWpozOK2OrXnOQ/febqvkL4xqz4fsPEzOLOC0\negUtzRn09nStOX8P3Xs7PvS+LbqP9/Z04Z0/s1H6ml3v7caeu27DX//oMm7MLK5579GJWex6b7ff\nH0Os5HLZz9t5PqtfyBO11STpFLC4HK7TrQA4+h8/WP1atql23Gi/56XxKejdwPjV6Gzf4QGUdGKL\nVn1E1gW58TQRgLXldCXvt8N0rQys2k/TbEVsXGi/5+NPntStM5+YnPelz7lsRXEcUlxhx5w6uRaV\nssaoV7c4ZWVTE68/G1muPgl9eBqNI3VyzY+yxtxNTZieXbL9OqM+JXGobnHCykXX689GG/UnIcUV\nNgzq5JofzbumZ5ewf8829J+6gJHxKUtdHzPpFN4yGJEm9dbfykXXj88mKSmusGFQJ9e82DlHz/GB\ncwAABQo2tK4DFOD61ALa8+t0d9oxWzCT1Ft/KxddsXk9DhwdjE3bhCRj9Qt5QquyGBmfQlMqhaVS\nCRvy2WogvmVjDu35dbr9ROzQGkSZVXUAlZF7qVxO/K2/rC9MpR9MHmLzerx0enjN99mMKxzY0Isa\nRhY8aoPD7/3p664Ce6E1iy8/dlf1a6OOgyyfe5tRCafsM2QzrnBgSSM1jNlGDQBwfUo/551Jpyz1\nG5mYnMfewy9XOw6OGKQVkppD12OU3zZrmxC3jpZxx5JG8oyVniqy55TKZWTS1gYktR0HJW1eACQ3\nh26XUduEOHa0jDsGdfKMLDikFKXa1Gl9fp3uc27ZmMOSg42qZQtXFXi7mCbOjGrK49rRMs4Y1Mkz\nsuCwuFyqjvJkNeR9O7egu8O7dEl3xHu4BMloI/E4d7SMK+bUyTP1C05SiqKbJ8/d1ISFhRIWl0vI\npFO4+92bqq+1WxpZaMvqljcy9WKPLOcuK4dsz61jCWRIsfqFfCNr6qTHTqli/esArlz0i53tClkC\n6Q9Wv1AoDA6NIp2y3txLy9G+XWWRg9i8AerFa9VgLTavX/V1bfBmMPGH3nL/mblF3TRabZUTNQ5H\n6uQ5J5tRpxTojsw5+gsfttUNFkfq1HCyiolMOoXWlozuKE+BgkqD3NWeOnEGx/qHsLRcRndnjrnb\nEGBb3XBjUCfPySomFpdL0oZbywZ3jNqGG1qNNMB0ix+sLjKS9frh5HQ4MP1CnjNaul8rpQCbOvK4\nNjWPqdm1W5/JpBUF7fl1uDa1wMoLjzz74lnd/i+7txfx8H13rHk8KTtHhQF7v1DDWc2pFzvz6Nu5\nxZMOj8y9O2d2vvjZNpbdoM7FR+S5+sUsMpevTnu2axJXODpndg742UYLc+rki9rFLLJ0zC0bc7g0\n7s3KRK5wdM5sE41L41MBHQl5gSN18p1RbxFZv5hccxOKnXnDhl21zCovBodGceDoYLUHDRtSvU12\nDjSlMvh5RYijkboQog3AMwDaAGQA/Iaqqn/j5YFRfJjtV6mXz33ftpuhXnwLgIJccxrTc8b7lRpV\nXtTnjFlFs5qVnau4sCg6nKZf/gOA76mq+l+FEHcA+B8Atnt3WBQ3st4i9QG/PbcOC0ulVZUYRgFd\nm2w1CjhW+rwnWe05GB7TT7UwvRUdToP6VwBoBccZALPeHA4lkRbw7a5ETacUSzvzsNOgOe0cGM1/\nUDSYBnUhxKMAPonKcj9t2d8jqqqeFkLcDOBpAP/e16OkRLBbCWM10HAFpHVcWBR9pkFdVdVjAI7V\nPy6E+FkAf4JKPv2vfTg2ShizKox6VgMNA5V1ZvMfFH5OJ0p7APwZgH+pquqPvD0kSirZiLpeJp3C\no313Wg40DFT2GO1nSuHnaEWpEOLbAN4J4DwqKZlrqqo+YPQariglM1Zz6lzhSEnCNgEUafU9RYx6\nqBMlAYM6hZrVToBEVMF+6hRaXiwC4kWByBjbBFBgjBYBWaFdFIbHplEql6sXBS5hJ3obgzoFxu0i\nILcXBaIkYFCnwMgaR1ldBMSVoUTmGNQpMEbdGq1we1EgSgIGdQpM/eYZxc68rZpztxcFoiRgSSNF\nCvfGpKRhnToRUYxwj1IiogRjUCciihEGdSKiGGFQJyKKEQZ1IqIYYVAnIooRBnUiohhhUCciihEG\ndSKiGAlsRSkREfmPI3UiohhhUCciihEGdSKiGGFQJyKKEQZ1IqIYYVAnIoqRpiDeRAjxAIBfUVX1\n4ZWvewH8FwCLAF5UVfVQEMcRJ0KIYQBnV748parqZxp5PFEkhFAA/BGAdwGYA7BPVdWfNvaooksI\n8SqAaytfvqmq6t5GHk9UrcTHL6qquksI8TMA/hhACcCPVVV9zOz1vgd1IcQfAPgFAK/XPPzfATyg\nqup5IUS/EOLdqqq+rv8TqN7KiT6tqupHG30sEfeLALKqqr5/5Q/pKyuPkU1CiCyAsqqqH2z0sUSZ\nEOI3AXwMwNTKQ18B8GlVVV8RQnxVCPFRVVX/l9HPCCL9chLAr2tfCCFaAaxTVfX8ykPfBbA7gOOI\nk+0AikKIl4UQfyGEuKPRBxRRHwDwVwCgquoggJ9r7OFE2rsA5IQQ3xVCfG/lIkn2nQPwQM3X21VV\nfWXl/78D4F6zH+DZSF0I8SiATwIoA1BW/vuIqqrHhRA/X/PUNgA3ar6eBHCbV8cRN5LP9TEAv6uq\n6p8LIe4C8AyAHY07yshqA3C95uslIURKVdVSow4owmYAfElV1aNCiNsBfEcIcQc/S3tUVX1OCLGl\n5qHa/UknAbSb/QzPgrqqqscAHLPw1Buo/DFpWvF2Ho7q6H2uQoibACytfP+kEGJTI44tBm6g8u9P\nw4Du3FlURplQVfXvhBBXAdwCYKShRxV9tf8eLcXKwKtfVFWdBDAvhLhtZaLqXwB4xeRltNrnAHwC\nAIQQ7wJwsbGHE1knAXwYAIQQ7wPwo8YeTqQ9CuD3AGBlkNEK4HJDjygeXhVC3L3y/x+ChVgZSPWL\njl8D8CeoXFReUFX1/zToOKLqiwCeEUL0oVJB9G8beziR9RyA+4QQJ1e+fqSRBxNxRwF8QwjxCiqj\ny0d51+OJxwEcEUJkALwB4FtmL2CXRiKiGOHiIyKiGGFQJyKKEQZ1IqIYYVAnIooRBnUiohhhUCci\nihEGdSKiGGFQJyKKkf8PeUF7Og06omsAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "psy.plot.lineplot(ds, name='y', coord='x', marker='o', linewidth=0)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "However, it is hard to see how many data points there are shown. Therefore this is a good candidate for a\n", "``density`` plot:" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "psyplot.project.Project([ arr1: 1-dim DataArray of y, with (experiment)=(1500,), ])" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAD6CAYAAACiefy7AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAFipJREFUeJzt3X20bHVZwPHvIS5Eeq8VodjKxBKezBdaYkiIgviKVki1\nrOVLdRE101xhWvmSrV40yiI1JVcqagKlZlRqKhpGeJepYSQaPUCCpbDshgGXtxt4pz/2PjF3n+Hs\nu2fvmTlzft/PWnfdMzN7Zn6zz8xznnn2bz+/ldFohCRp89tv0QOQJM2HAV+SCmHAl6RCGPAlqRAG\nfEkqxP7zeqKdO3c5HUiSOjrkkK0rQz2WGb4kFcKAL0mFMOBLUiEM+JJUCAO+JBXCgC9JhZjbtEzB\nIffYs+ghSEtj5y3mo0Nzj0pSIQz4klQIA74kFcKAL0mFmOqgbUTsD5wNHAYcALwmMz8w4LgkSQOb\nNsN/FvDfmfkY4CnAm4YbkiRpFqadlvle4H31zyvAHcMMR1p+o69fu+7tK9/+nXMaibS3qQJ+Zt4K\nEBFbqQL/K4cclCRpeFMftI2I+wEXAu/KzPcMNyRJ0ixMe9D2PsBHgRdm5ieGHZK0XNpKOF23t+Sj\nWZm2hv9y4FuBX4uIVwMj4KTM3D3YyCRJg1oZjeaz8qBLHNpLZ7PomtH3VWrGby+dikscSpI6M+BL\nUiFsjyy1aCvh7Pny5b0ef7/7P6jX/aV9ZYYvSYUw4EtSISzpSA2zLuG0PV6zxNMcT6mzdtSfGb4k\nFcKAL0mFsKQjtWgt4Vz1xX5P8MAHr/t8lng0FDN8SSqEAV+SCmFJZxMZosdLieWBzvutpYQzuuqK\nTg+3plFKxxKPtK/M8CWpEGb4S2wWXRs9IDjhIG3HjP6OvGbd7bfEYevevy3jX/P8/s60j8zwJakQ\nBnxJKoQlnU1sX1oAtB0A3Izlgr6lsK4lnKbm9l1LPM1ldDyIq31lhi9JhTDgS1IhLOkskVl0cWy7\nj6f1d3f9FTs7bX9w43KzxNPVpPeJvyeBGb4kFcOAL0mFsKSjdZVwWn/biVZts3K6lnCamvdfU+Jp\nXHbWjqZlhi9JhegV8CPikRHxiaEGI0manalLOhHxMuDZwM3DDUcztS8LdbT0bWna6LN2ZtFvqM3V\n13X7SDzgvvfc63LnEk9Ld01Ym9lttN+T5qNPhn8VcMpQA5EkzdbUAT8zzwfuHHAskqQZcpaO9tYs\n+2zCxTiGnpXTLOFcedvubgO6bu+LzRJPq5bfmbRqiFk6a9p3S5I2niEy/NEAj6GNqmP2uOiDuLM4\nSDt4Rt+w5v7XTd5uVdd5+jBhrn7zPh7ELUKvgJ+ZXwaOHWgskqQZ8sQrSSqEB23Vy0Y/iDuxG2jH\ng7RdfWZXvxLPGh1LPAc88Ihhn1+bhhm+JBXCgC9JhbCkU7hmOaNppVke2GCzdpZxVk6b1pJQo8Rz\n8BGH7HV50u+02X6hadGzqzQfZviSVAgDviQVwpKO1tUsD6wp8WwwbW0TYPhZOU1Hbz1wr8uDz9pp\naJagDp20Jm5jP7hoSpnM8CWpEAZ8SSqEJR11sqbE07i9rVTQdzbIIhY0aTr8oL1LNm2zdpolnqah\nSz6TSlSejCUww5ekYpjhbyLNbHpiW4GOmtnilkkHBNfR1nqhLeNvy+i79raH7v3t2zQzfmmjMsOX\npEIY8CWpEJZ0lkjXcsc0us5JbzuI23W5vY1wULarzksSNrV0wxzC2vMpbLVQIjN8SSqEAV+SCmFJ\nR510nrUz8Cn9XWflTCpRNWfldNVWwml2r5Q2CjN8SSqEAV+SCmFJR3tplmjaZu2sKfE0bm/rrjnE\nyWF9tZ1o1bWE0/XktOaatE1tJaipSkh2zyySGb4kFWLqDD8iVoCzgCOB24HTMvNLQw1MkjSsPiWd\npwEHZuaxEfFI4Mz6Os1J24lYa3rrTHqMxuW2NW67GvrErK6zcvrOyJmkawmnrazVLIOteb59GJO0\nL/qUdI4DPgKQmZ8GHjHIiCRJM9Enw98G3Dh2+c6I2C8zJyWS2qSaGXQzG21mv60Z/xx0PUjbO6Nv\n+RbTtg/avgEMvUTjJLZa2Bz6ZPg3AVvHH8tgL0kbV5+AvwN4CkBEHANcNsiIJEkz0aekcz7whIjY\nUV/ePsB41MM8umm2HQTtXeJpOcC5L60Tuurb7bJrCWeN5vaNA9PNxx/6wPokbQvXaDlNHfAzcwS8\nYMCxSJJmyBOvJKkQtlYoyKSv5WuOsg9cLuhb4mlqK+EMMe++bVZOWwmndwfQlhJPU9d2GND//Ahn\n7SwnM3xJKoQBX5IKYUlHe1lTrhj4pJ62Ek/X+7eZdJJV7zVoB7amBUbPDqJdu3WqHGb4klQIA74k\nFcKSziY2xIlYzfJAswTTLLF07VMzRIlm1vrOymmbwdL6e2k836w7nE7iiVibgxm+JBXCDF9zNXSG\nfuVtu9e9/fCDDmx9jL5LFDb1nZM+9EHceXBe/nIww5ekQhjwJakQlnQK17oM4sAHBNtKMBtSz9YJ\nbTofXJ/HQdxmO4euHUC1IZnhS1IhDPiSVAhLOgWZNHOirXzQt9XCMpRwWrthtpj1jJSus3Ymjb/r\nQjPanMzwJakQBnxJKoQlHXWy5qSknguOfGZXv5LP0Vv3PrGqeaLVpM6YzROt1pjxrJw2fWftTFow\nxRKOwAxfkophwJekQljSUS/N8sisu1k2SzhtJpVvNvqsnDZz6bXjiVabkhm+JBWiV8CPiFMi4tyh\nBiNJmp2pSzoR8XrgicClww1Hi9bWW6fZt6V5ItaaWTHX7X1x6BOx9mVWTlPXBU0WXcLpbFI5ZsLM\nHZWnT4a/A3jBUAORJM1Wa4YfEacCpwMjqgRvBGzPzPdFxPEzHp9mrO8yiH3n5Xc9CNtm6MVMNoK2\n39E+HcTdh7n62vxaA35mng2cPYexSJJmyFk6klQI5+FrXV0XSGltW9DzIO40B2nXWHDrhA3BefZF\n6hXwM/Mi4KKBxiJJmiFLOpJUCEs62kvXWTvNOe1bGrcf3PaE17VtsLdmCadtVs40XSI3+rz7rrN2\noH/7hSLLXpuQGb4kFcKAL0mFsKSjTrq2Xuhc4mnROguoacJslBLLEyW+Zq1lhi9JhTDgS1IhLOlo\nXX177TQ1Z9H0LfGUMCunzdC/oyHGoI3JDF+SCmGGr166HsQdNVox9O1m2bW3fQn2Jdvu+i3ADH5z\nMMOXpEIY8CWpEJZ01EnnA4SNEsuaJRL7mqLro+UJ90GpzPAlqRAGfEkqhCUdDap1fdWBF94ocRaO\nNC0zfEkqhAFfkgphSUe9TLMYx7iuC3N0LeE4G0W6ixm+JBXCgC9JhbCko0F1PTFr6Fk2lnCku2eG\nL0mFmCrDj4htwDnANqpV7H4pM/9xyIFpc5h1r3YzemnfTZvhvwT4eGaeAGwH3jzYiCRJMzFtDf9M\nYHf98xbgtmGGI0maldaAHxGnAqcDI6pmhyNge2ZeEhGHAu8GXjzTUWrT6FvisYQjTW9lNBpNdceI\neChwHlX9/oK27Xfu3DXdE20ih9yjuR6UDPi6OztvcU4JwCGHbB2sq/i0B22/H3gv8PTMvGyowUiS\nZmfaGv5rgQOBN0TECnBDZp4y3LBUCjN2aX6mLul0ZUnHko7UhSWdypAlHfeoJBXCgC9JhTDgS1Ih\nDPiSVAi7Zc6RB6EkLZIRSJIKYcCXpEIY8CWpEAZ8SSqEAV+SCmHAl6RCGPAlqRBza54mSVosM3xJ\nKoQBX5IKYcCXpEIY8CWpEAZ8SSqEAV+SCjGX9sgRcQrwE5n5zPryI4E3AHcAH8vM35zHODaTiPgK\ncEV98VOZ+cpFjmcZRcQKcBZwJHA7cFpmfmmxo1peEfE54Ib64tWZ+ZxFjmdZ1fHxjMx8bER8L/BO\nYA/whcx8YZ/HnnnAj4jXA08ELh27+i3AKZl5TUR8KCJ+IDMvnfwIaqrfBJdk5smLHsuSexpwYGYe\nW3/IzqyvU0cRcSAwyswTFz2WZRYRLwOeDdxcX3Um8IrMvDgi/jgiTs7Mv5728edR0tkBvGD1QkRs\nBQ7IzGvqqz4KPG4O49hMjgK+KyIujIgPRsQRix7QkjoO+AhAZn4aeMRih7PUjgTuEREfjYiP139A\n1d1VwCljl4/KzIvrnz8MPL7Pgw+W4UfEqcDpwAhYqf/fnpnvi4jjxzbdBtw0dnkX8IChxrHZ3M1+\nfSHw2sx8f0Q8CjgHOHpxo1xa24Abxy7fGRH7ZeaeRQ1oid0KvC4z3x4RhwMfjogj3JfdZOb5EXH/\nsatWxn7eBdyrz+MPFvAz82zg7H3Y9CaqD9qqrdxV91PDpP0aEQcBd9a374iI71zE2DaBm6jef6sM\n9tO7gio7JTOvjIjrgfsCX13oqJbf+Puxd6yc+yydzNwF7I6IB9QHzZ4EXNxyN+3t14FfBIiII4H/\nWOxwltYO4CkAEXEMcNlih7PUTgX+AKBOQLYC1y10RJvD5yLiMfXPJ9EzVi5qEfOfA86j+oNzQWZ+\ndkHjWFZnAOdExFOpZjr97GKHs7TOB54QETvqy9sXOZgl93bgHRFxMVVWeqrflgbxUuCtEbEFuBz4\niz4PZrdMSSqEJ15JUiEM+JJUCAO+JBXCgC9JhTDgS1IhDPiSVAgDviQVwoAvSYUw4EtSIQz4klQI\nA74kFcKAL0mFMOBLUiEM+JJUCAO+JBXCgC9JhTDgS1IhDPiSVIj5rWl7640LX0tx9PVrFz0E9nz5\n8kUPoXLVFxc9AkZXXbHoIQBwR16z6CFw/RU7Fz0EAK6+7uZFD4Erb9u96CEA8JldG2McbxndtDLU\nY5nhS1IhDPiSVAgDviQVwoAvSYUw4EtSIQz4klQIA74kFcKAL0mFMOBLUiEM+JJUCAO+JBXCgC9J\nhTDgS1IhDPiSVAgDviQVwoAvSYVYGY0Wvi6JJGkOzPAlqRAGfEkqhAFfkgphwJekQhjwJakQBnxJ\nKoQBX5IKsf+0d4yIFeAs4EjgduC0zPzS2O3PBZ4H3AG8JjM/FBEHA+cB3wxcC2zPzNsnbTvtuBZh\nmn0xdtsvAvfOzFfUl58NvBS4AXhXZp49txfS05TvifsA5wJbgK8Dz8rMWyLimcBLgDuBd2TmW+b7\navqbcn98G3AFcFm92fnAPwCvB0bACnAMcHJmXjCv19JHx/3w25n5txHxh8APUL3m+wL/k5nHRsRJ\nwKvr6z+XmS+a76vpb8r3xWHAu+pNvgw8r46dbwCOBXbVt52cmas/r9Enw38acGBmHgu8HDhzbMD3\nAX4B+CHgycDvRMQWql/UuZl5PHAp8Px1tl0mnfdFRHxzRLwbeMHYtgcDvwU8BjgBeGZEfPfcXkV/\n07wnfpUqoK++J06r7/I64ETgOOCXIuJec3sVw5lmfzwcOC8zT6z//VFm/ktmPjYzTwTeDLx/WYJ9\nrct+OCMitmTm6Zn5WOCJVMnPaRFxT+D3gKfWj3VN/ZlZNtO8L14HnFV/Ti6iSoager88aez9crfB\nHvoF/OOAjwBk5qeBR4zddjTwycy8MzNvAq6k+mv2//cBPgw84W62fViPcS1C133xMKpvOe8CXjO2\n7fcA/5yZN2bmCPgsVTa3LDrvh/qDfW5E7Afcj+rDDfAvwLcBB9WXl/GU8GneF0cBR0XE30fEeyLi\n0NU7RMS3AL8BvHheL2Ag0+yHVS8GLsjMf6XKZC8DzoyIfwC+lpnXz+MFDGya2Pkg7oqdO4Dj6m8K\nhwN/EhGfjIjtbU/cJ+BvA24cu3xn/aGddNsu4F7A1rHrJ10HcHN9/TLpsi9uBu6VmTdk5sepvqKv\nuhJ4cEQcUn+4HwfcY4bjHlrn/QAQEftTfZBPAC6sb/8icEl9/QfrN/+ymWZ/XA68OjNPAP4a+KOx\nbZ4DvDczvz6zEc/GtO+LLVSljd+vb/sOqvfIy4CTgNMj4oGzG/bMdN0f26i+/Z5cX/ejVHHhW4A3\nAs+i+jbw8xHxkPWeuE/Av4kqWP//Y2XmnrHbto3dtg34n8Z9to5dN77tVu7K8pZFl31xt68vM2+g\n+qr2fqq69iXAfw8+2tmZaj/U2cyDgecD746IhwJPBe4PHAbcJyJ+fMZjn4Vp9scngL+vrzufqo69\n6pnA22Yy0tma9vPxeOCisTLF9cBnM3NnZt5CdWxjfP8si2n2x0uBH42IC4E9VHHhVuCNmXl7Zt5M\nlSwdud4T9wn4O4CnAETEMdx1kAngM1RfOQ6oa6/fB3yhvs9T621OAi6mKltM2naZTLMv1oiIbwIe\nnpmPAZ5eb7tjlgMfWOf9EBFvjogT6m1uBr5B9Qa/Fdhdl7b+i6q8s2ymeV+8DVj94/Z4qj/6RMQ2\n4IDM/Oqcxj6kaT8fj6cq/a66BHhIRHx7/a3wGOBfZz34GZhmfzwBeHl9HGcP8DEggE9GxEr9beg4\n4HPrPfHU3TLHjjSv1tu2UwXzKzPzgxHxHKqMbYXqSPNfRcS9qerW96T6C/WMzLxt0rZTDWpBptkX\nY/f9GSDGZum8muqgzm3AH2TmX87vlfQz5XsigLdQvYn3AC/KzIyI5wOnAruBfweem5l3zvcV9TPl\n/jgMWJ2ZdQvVDI6vRcQjgFdk5o/N9UUMYNrPR0R8AHhlZn5+7LGeDvwy1TGd92Tm77NkpnxfHE11\nwP52qnLnCzPzGxHxUqrk8H+BP83MP1nvuW2PLEmF8MQrSSqEAV+SCmHAl6RCGPAlqRAGfEkqhAFf\nkgphwNfCRcTxEfFnY5d/IiI+HxHfNYPn+kREHLHO7c+vz4UY4rnOi4j9I+J+EfHDQzym1IcBXxvF\nCCAifgr4FeDEzPzKYofUT2Y+oz5Z7ETgUYsejzR1P3xpYCsR8Syq1rCPW22WVjeDemO9zfVUZ98+\nHPhdqrNw30rVTOsiqjMX91D3BI+I1wKPpkpszszM90964og4jqrf/PX1/T9VX/8i4Bn1dX+emW+K\niHfUz3sYcCjws5l5aX3991B193xD3QH0auD7qVpAHxQRn6JqhXt4Zo4i4gyq3jATxyUNzQxfG8Wj\ngedS9cwZXw/hrcDP1z1EPkyV/UPVT/z4zDyHqtnUuXWHyWuBkyLiycADMvPRVBn2q9bpqX8W8JOZ\n+STgaoCIeBDwk1SZ+aOBU8ZKQddk5pOBNwHPq/u0Hwf8GFWPqG/U243qn8+g6nH/N1T9o55Ud0c8\niaojpjQXZvjaKK6lahD1XKpOoU+ur38QcFbVcoctVKtBAeTYfUdU7WMB/pNqrYH7U/WVv5CqJ8n+\n9XWT3Dsz/73+eQfwvcBD6u3/rr7/t9bXA/zz2HMdm5k3R8TpVH+ctgLnrPM630bV4/2bgI8tW38g\nLTczfG0UV2Xm/2bmm4HdEfGq+vp/A366zvB/BVhdHnLP2H1XWLtAyuXAhfX9TgTeC3yJyb5aN3ED\n+MH6/wS+UK8i9FjgndzV1XCv56pXKTqqbmz2w8Dv1Z1PV9c62EMV4MnM1T8opwJvv7udIc2CAV8b\n0XOoSiXHUy0B+e56haPfAT4/YftR8+fM/CBwS32/fwJGdc/wSd0Cf65+jo8B313f//PAhfVKQp+l\nWlno2kn3z8yvAYdGxA7gAuB1mfmNsW0vo+pl/vT68rnAoZl5+b7tDmkYdsuU5iwiXgbszMx3Lnos\nKos1fGmO6tk89wV+ZNFjUXnM8CWpENbwJakQBnxJKoQBX5IKYcCXpEIY8CWpEP8HJPkkKPz7SIwA\nAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "psy.plot.density(ds, name='y', coord='x', cmap='Reds', bins=50, density='kde',\n", " clabel='Kernel density')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The ``densityreg`` plot method combines this plot with a fit through the data" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "nbsphinx-thumbnail": { "tooltip": "Make a fit to your data and show the source as density plot" } }, "outputs": [ { "data": { "text/plain": [ "psyplot.project.Project([ arr0: 1-dim DataArray of y, with (experiment)=(1500,), ])" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAD6CAYAAACiefy7AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XeYVOX1wPHvzGzvfZel11eaoCCiggIWxBKNGklsEcVY\nYoxGY9fEJCYaE41GjRoh+lNJxBaNYBcRUREQkPrSQcouy/ZeZu7vj5nVZXbY2ekze8/neXxkZm45\ne3fm7Dvnvvdci2EYCCGE6PmskQ5ACCFEeEjCF0IIk5CEL4QQJiEJXwghTEISvhBCmERcuHZUVlYr\n04GEEMJH+fnplmBtS0b4QghhEpLwhRDCJCThCyGESUjCF0IIk5CEL4QQJiEJXwghTCJs0zIF5Kc6\nIh2CEDGjrF7Go8EmR1QIIUxCEr4QQpiEJHwhhDAJSfhCCGESfp20VUrFAXOBAUACcL/W+n9BjEsI\nIYKipaWFhISESIcRFfwd4V8CHNRanwicATwevJCEEOLwdu/eyR133NKtZZcuXUJjYwMAhmHw978/\n3OXyzz33LEuWfML//d/c757rznqxwt+EPx+4x/VvC9AanHCEiH1Gxb4u/xOB6ddvALfeepfX5crL\nD9LQUE9mZhY1NTXMnz+P1atXHXb5FSu+AmDy5Cm0tbWxZs3qbq3XXRUV5TzyyJ8D3k4g/Er4WusG\nrXW9UiodeAXwfvSFECIIvv12N1u3aq/LLVjwFieeOAWAjIwMZs68mNTU1MMuv3btGoYOVQAMG6b4\n+uvl3Vqvu3JycklJSWXNmsD/ePjL7wuvlFJ9gdeBx7XWLwcvJCFEtLrj8Xd4fdG6oG7zvKmj+NP1\nM7pcZt26tdTUVNPS0szOnTv49NNFHHPMRD7//DNeeGEul112JVu3bmHQoMGccMJkACorK0lMTOp2\nHJWVFSQnJwOQnJxCeXl5t9Zbt24to0aNPuQ5u93ORx99wL59eygsLGLDhvX85CeXcOqp05kz52nG\njDmq23EFk18jfKVUIfAecKvW+vnghiREbPG1ZCMlH98tXfop8fHxTJp0EieccCJpaekAHH/8JBob\nm5gwYSJTpkzjs88+/W6dlpZmn/bhcBhYrVbXv+3YbF2nx8rKCioqylm7dg0VFeVUVlZ899rWrVuY\nOvVkiov7YBgGU6eeTG5uHgMGDGLt2m98iiuY/B3h3wFkAfcope4FDGCG1tq3IyyEiCl/un6G19F4\nKFxwwUz++c9/8Pjjj3DrrXcf8lpGRgY2m434+Hjs9rbvnrfb7T7tIycnh8bGRgDq6+vJysrucvmk\npGTee28hW7ZoUlJSmD79jO9eU+oIANav/4aZMy+muLj3d6/FxUWuo41fe9Za3wjcGORYhIgJoR6F\nu2/fklMc0v3FguXLl3H77ffQ2NjISy8dWlQwDM93T20frXe1bEnJfoqKegFw5JFj2bRpA8cddwIb\nNqxn/PgJXe6jvPwgVVWV9O8/gIqKcsrLD9K7dx8ANm5cT3Fxb3bs2E5xcW/WrFn1XRnHZrP58JMH\nl1x4JYSIeiUl+1m06EOWLv2UhIQEduzYzpo1q/jyy8/ZtWsny5d/yXvvLUTrTezatRPgkPp9Y2Mj\n8+fPY/funcyfP4+mpiZqa2u5777v55uMG3cMVVVVLFr0IRaLhQkTJnpcr12fPn25/PLZ2Gw2Zs26\n6rtkD7Bs2RcsXryIUaOO5NNPP8E5mRG2bdvKyJGH1vvDyXK4v47BJjcxl26ZPUW46+xmHeEH2i3z\n3/9+kbPOOof09PQgRRS4J598jEmTTuTII8d2e51g3sRc2iML4YW3BO/YtTGg7Vv7Dw9ofeHZ2Wef\ny0cfvc8555wX6VAA5zz8hoYGn5J9sElJRwjRI6WlpTFgwCBKS0siHQrgnId/yy23RzQGGeELIXqs\nMWMiN5qORpLwhXAT6hKOt+25l3hk1o4IFinpCCGESUjCF0IIk5CSjhBeeC3hbF0f2A6GjOxyf1Li\nEcEiI3whhDAJSfhCCGESUtLpQYJxBagZywM+HzcvJRxj62afNtfpMkofSzxCdJeM8IUQwiRkhB/D\nQtHTRU4IejhJ6+OIvlXv7HL5eDWgy/W9jfg77V9+Z6KbZIQvhBAmIQlfCCFMQko6PVh3WgB4OwHY\nE8sFgZbCfC3huHNf3tcSj3uTbTmJK7pLRvhCCGESkvCFEMIkpKQTQ0LRxdHbOnJZv+/KN5f5tHyu\n22P3Eo+vPL1P5PckQEb4QghhGpLwhRDCJKSkI7pkhsv6vV1o5W1Wjq8lHHfu63cq8bg9llk7wl8y\nwhdCCJMIKOErpY5VSi0KVjBCCCFCx++SjlLq18ClQF3wwhEh1Z0bdXjp2+Iu2mfthKLfkDc79vv2\nkRjYK+2Qxz6XeLx014TOI7to+z2J8Aikhr8V+CHwQpBiET5qbbOzbFMpn3yzh7KqRurKD9I7K5ER\nxamcMjyXnFT31CDCqRkraxPz2BSfTWlcClbDIN3RwsiWCka3lJNmtEY6RGEyfid8rfUbSqn+wQxG\ndE9NQwuPv7mGJ95aQ2Vds8dlEuMsnDO2gFum92eE2whShFalNZHXUwezMHUAtdYEj8vEG3ZObdjN\nhXVbGBjm+IR5ySydGPP6Z1u54cnFVNQ2kZWWwAWTBnLqUX0Z3i+HjJp97CpvYum2SuavOMD8FaW8\n/vUBfj6lN3edNZiU7uzAvezTA2/GEexZOe0lHAP4ILkvT+aPosGWSIq9iRMqNzG08QB9W2uxWCxU\n25LZkJzPyrT+LEwdyLsp/bmsdhMz67ZgwwA6l3i88vI7E6JdMBJ+p/bdIviamlu5+fFFzH1vA0kJ\nNq47axSzTz+CYX3zsVqdFVrHrhqGFKZx8og87jlrCK+v2M/t/93Oox/v4YONFbw2I5m+GfI3PhSa\nsfJQ9jg+Te5NgqOVM8pWcXLNVpLcJk0WttYyrOkAP6jcwMrUPryWdxTPZYzgy8RCflv5FTkOz9/Y\nhAiGYHz6jSBsQ3ShrqGZC257kcVfb2do7ywenHUsk0b1Ji01+bDrWCwWzj+mmNOPLOD6lzbwytcH\nmTSvkflnZ3Bs76Tu79zH0WOkT+KG4iTt4Ub07b5uhqd7Hc+OpDz6NpZxWemXFNsbutymFYNj6r9l\nRGMp8/LGsTq9H9flTuYX+xeDl5O+vs7TBw9z9d3XkZO4phDQtEyt9S6t9fHBCkZ0Vl3XxIwb5rD4\n6+1MHdOHF26ZxrSj+nWZ7DtKTYxj7qzR3HfWACqbDc59s5plextDHLV51FjiebT3VHYk5TGqdhc3\n7l3kNdl3lOpoYfaBLzi9Yh0V8Wn8pffJ7I6Tcy4iNOTCqyjW3NLGj257gRUb93LWpCN4aPaxjBxY\nQEKC5xOBh2OxWLj59EE8OiWFuhY4761avt4vST9Q9ZY47sw9nv0JmUyo2sLPSr8gyeL7F14LcHbl\nen50YAX1tiRuzTmBA9bu/UEXwhdS0I1SDoeDK3//CktW7+TkCUO472dTGZJmx2az+b3NWWPSaXHA\nzYsbmPl2LYtn2ijO8O2PR6c4o/wkrsduoD6epPXEjoXfZ09gc0I2Y2p2cPHBldgszmLKV7X+1eFT\najdwVBusKh7PTTnHceveD0lxtML+rtdzL/EkDBnm1/5Fzycj/Cj1wPOf8NrH6zhKFfOHa05mcJ98\n4uMDn1d/9VHp/GpcEvvrDWa+XUNjqz0I0ZrPnIwRfJ1UwJjGEq448BVxluDMXRh9YAMnVW2iLCGD\nuQUTO9XehQiEJPwo9PHyrfxhzkcU5abz11/OYFBxDomJiUHb/n2TMpgxII5VB+zc8GF10LZrFouS\ne/Nq2lCKWmu5s3I5cUGcp2YBLij/hmENJWxMLebtrBHB27gwPSnpRJm9ZdVceu9/sFmtPPTLMxhU\nnEVGenBP4lksFp4/M4fjXzrIv3UrJyVs4eL+nmvPFvfyQJTN2gn3rJx9thQezh9LoqOVS/ct4WBb\n90/QdteK2iaOaljMHnUW7+WMYsi+MkY0dYjJrcSTOyz/kMfuJSro3H7BXaRnV4nwkBF+FDEMg6t+\n/yoVNY3ccsmJHD20gML8nJDsKyXewgtnZpJoNbhtYxpbaqR44I0dCw9mj6fZGs9ZB1YyoK02ZPtK\nsjczZecnWDB4oXAiDVZpkyECJwk/ijzz+jIWrdzOCWMGcOEpI+jby/10XHCNzk/gPtVAdZuV675J\nxe6QSyq6Mi9tGBsTchhVu4spdTtDvr/8hnJGl6yhJi6FF/PGhXx/oueTkk6U2L6nnDueeIeM1ETu\nnT2Vopz0oJyk9ea6QW0sKG1hSUUCT2xr4oahhyZ99/JApxJPlPHWNgH8m5WzIy6DeemKnLYGLipb\ngbWLk7QT0g893+LvrB2AMaXrqMruy5r0/iyv+5ZjGvZ2Wsa9BFXk6Z64bsdBbppiTjLCjwKGYXDt\nA6/T2NzGHbOm0Tc/jeysjLDs22KBp8Y2kmIzeGBbOpurZdaOOzvwcNZY7BYrV1WtIdNoC9u+rRhc\nfuAr4hx2Xs0fR71FxmjCf5Lwo8C8d1fx6aqdTBo7gOkTBtKnKLSlHHf9UuC3qoHaNgu3bkjFMKS0\n09GbqYPRCTlMbPiWk1tKw77/wtZaTqtcT11cMq/mjg37/kXPIcOFCKuoaeDXjy0kKSGOOy6fSl5W\nSlhKOe6uHtjGf/a08lF5IvO/bWJmP8/LdSrxuL3urVQQ6GyQcN/QpJQEnk8/gjR7M7+o/gaAocmH\nlmy2NHZdsnEv8bjrTslnepVmZVo/vsoczOqWPYxtrTzssp5KVHIxlgAZ4Ufc3U+8Q2VNI9ecP5F+\nBank52ZHJA6rBR4b04gVg/s2p1HVLKUdgAcYRqM1nh9Xryc/gjcsicPBxWXLAXg8cyzy2xH+kBF+\nBK3evI/nF6xiUO8cfnLaaIpyA6vbu4+mPbYV6MKYTINZ/ZqYszuZB7ck8adRrZ1Gi/GeTgh2wVvr\nBW8jfm8jel9720P3+9uvScjj7bwiBjVXcH7jLucJDw/cR/yhMri5gvE1O1iRMZDXkwfwo8adYdmv\n6DlkhB8hhmFww0P/xWEY3HbZFDLTEkhLS410WPxueAu58Q6e3Z3Mpmrzzs23Y+GJzNFYDINrqtZ8\n1ycn0s6r+IYkRyvzMkdQZZG5+cI3kvAj5OUP1rB8w16mjR/MhBG9KM6PTCnHXUY83K0aaHJY+K1O\nMe0J3HdS+rMzPpPJ9bsYa4+e9hOZ9iYuqd1EvTWBOekylVL4Rko6EdDU3MrdT75LfJyNX108mcy0\nxG6dqPW13OGPVr2TSwx4Jn4oCw6k8JE1mVPSmw67vLeTuL7ebi/cJ2U9qbfE8Xz6cJIcrVxVu8Hr\n8j7fktCdl26Y7gbUb+ft1IF8mDqAc+u3deueuJ2vp5BWC2YkI/wI+PvLS9lbVstPpo+lT34qhXnR\nMbpvZ7PAgwXOC3z+UNGHNpNdgfvvtGFU2xI5u3YzhUZLpMPpJB6Dq6rXYbdYeSZjdKTDETFEEn6Y\nHayq5y8vfkpmaiKzzxlPXmbKd/ekjSaTUho4JbmKb1pSmVcVnovAosFekngjbTC5bQ1cUrc10uEc\n1qTmEoY3H2RVchEfO7IiHY6IEVLSCbNVeh819c3ceumJZKXGkZsTvR/WP+TvZ9HuTB6u7s35GRtI\njbP6PmsnyJf0+zorx9OcdPdZOR39jcG0WmzMrFlP8mHuXuWthOPevTJU/mBs43zy+Kt1GFP4SkZv\nwit5j4TZKROG8MXca7lg2gjys6P73qXDElv4SfpBdrcl8o/KvEiHE3IbSONNiujXUsUPmjr3rIk2\nR1pqOdVRwmZLBq86CiMdjogBkvDDzGKxMGJgAfE2yMqM/lLJ3XkHSLHY+UdNLw5GXzk7qB5iKIbF\nwk+r18XMB+M2yzbiDAePW4fQZK5TLcIPUtKJAKvVSp/C0PS5D5R7iaZA7+TqzAM8UtWLRyoLuL/w\nwCGvdyrxuG3PW3dNXy8OC4Ud++tYlZDHZ3m5jGwqZXLrwUNe97WE4+vFad46J3VVgupPI5dkHuS5\nmgL+k6O4Oa+iezuV7pmmFCsDmR4lLi4uKi6y6q5f5hwk29rG87UF7GrueW8ZA5ib4ZymOKumczvl\naHdHbpnzW1h1EdX26LhATEQnvz+9SimLUuofSqnPlVIfK6UGBTMwET0ybA5+lV1Cg2HjoYqiSIcT\ndJ8lFaMTsjm24VuObKuJdDg+y49r4+rMA1Q44nn4YM8/1yL8F0hJ51wgUWt9vFLqWOBh13MiTLxd\niNWpt46nbbg99tR7BmB2VgX/qMrnlbpcft5UwvCk7rVdCPaFWb7OyumqHALQhoV/pQ/HajiYVbup\nWyH4WsLxVtbydsldd5pl/zLnIM/V5DO3tpBrsg/SK0EK+qKzQL6fTwLeBdBaLwPGByUiEZUSrQZ3\n5pbQipUHyntHOpygeYNe7IlPZ0r9TgbZ6yMdjt8ybA5uzC6h3rDxYA/8FiaCI5ARfgbQsclIm1LK\nqrU2b8etHm5mRhV/qyxgQUM2K+pLGJ/a2mkE7T4adR/9eh3xh0F7N8wWrDxcOIB4q52f1urvXnc/\nSRvwiN7Ltxhvx8DbN4D2bzVXZVXwVFUB/6nN4/qsAwxJCl4TZWm10DMEMsKvAdI7bkuSfc9ms8C9\nuftxYOGPFbH/gX8rdSAHbSlMr91KL8P/+85GiySrwe05JbRg5U8VvSIdjohCgST8pcAZAEqpicDa\noEQkotqZabUclVjHJ01ZfFybFOlw/FZvieM/acNIdrRyaX30tlDw1UWZlQyOa+TN+hxWNcisa3Go\nQN4RbwCnKqWWuh7PCkI8IgDh6KZZsaWMW4wWLrYew+/3FTCf3Vg79IoPuMTj5QRnd1ondMeraUOo\nsSVyYdU6ji5KBPy/iYmvJZxO3Jd3OzHtvv3DnVgH57ewu3NLmFU6kPsrevNqyi7fYnHxduMaEZv8\nTvhaawO4NoixiBgxwVLN8cZBPrfm8a49lzMs3bzYJ0pUWhN5LXUIWfZGLmrYDpkpkQ4pqH6QXsOY\nyno+bsxicW0pJ3XR3lqYS8+7ikaExW04yyCPWYfSFmM3SXkpXdFkjeO8mk2keJysGtssFvhdvvPb\n3f2VxThM1t5aHJ4U+UzE09fyTumui3JBRyOo4wzHfhZae/FvRy8utZR4XC7QEo87byUcb/Pud5HM\ngpQBFLbWcX7DLo9TZLzNyvFWwgm4A6iXEo879/ha9U4mpzQwOamGJU0ZvFGbxo/d1gn0+ohQzdpp\nbW2loaGRzBjoMxWLZIQv/PZryzbiDQdPW4fQGCODyL8yBLvFysU164jv4V0I/pC/DwsGD1b2pjVG\nRvnvfr6RiprGSIfRY0nCF37rQxM/Nr6l1JLEM0a/SIfj1WoyeMdSyODmCqY3e/5G0pOMTmrm3NQK\ntrYm8/R2W6TD8WrBkvVceOfLzP3fykiH0mNJSUccolO5wksJ5ReWnbxhFPOcdSA/MfZSYOn6Yh9v\nJR5vvJVs3LVfZGUA9+WOhUS43baFQYHehzaIOrXACLCDaMcSz/19YeEnBo9sT+XSftVkJkTnGK+t\nzc7tj7+LxQKnHTsk0uH0WNH52xcxI5tWrjV2UGeJ5xEGRzqcw/osqZh1iXkc3bCPKbbaSIcTNn1T\nYHa/Jg60WHlgs/9TT0Ptmde/YOueCs6dMooRAwsiHU6PJQlfBOynlm8pNhp4w9KHDY7ouxirBSv/\nzBiJzXBwVc26SIcTdneqFnLiHTy7O5lttdE3K6m2vpk/PreY5MR4rjv/WOLior/8FKukpNODBeNC\nLPcZIO4lmPLNZSRicIexhV9Yx/CbliH8ufLLw27PvU+NvyUaX/w3bRAlcalMr93KYEcD4Fs5J9BZ\nOd5msHj9vbjtr7sdTttlxMNdwxq5eX0q92xKZt4xvreRCOWFWL99+l3Kqxu47oLjKM5LITMj3ftK\nwi8ywhdBMd1SxtFGBauTi1icED1fycutSbyUpkizN3e7/XFPdOWANo5Ia+N/pUl8FEXnqzfvKuOZ\n/66gOC+DS2YcSa+8rEiH1KPJCF8EhQW4j82cY0zg2awxTCz9kERL56mA/ozQu7KlsevR6qu9xtFo\njefKyq/Joc3jMoHeotBdoHPSg30SF5wtF/4yspGzlqVzx6ZUlhbUEW8N3rxUf+fl//Kvb9Jmd3DL\npSeSn5VCcnL0lQR7Ehnhi6A5gjrOrN9OaVwqL6UNjXQ4bE7KZ1FKXwa2VDovsjK5k/IdnFPYxMa6\nOP62JfJjvf8uWssnK3dw7Ki+TDmqH73ysyMdUo8nCV8E1axaTYa9idcyFLutketR04qVl/PHYTEM\nrqlcTZylh19l1U1/GdVMms3BIztS2FUXuRO4tfXN3PTI28TZrNx66UkU5qRhs8nJ2lCL/J95EVFe\nb4PYzVYL7dKMVn5e/Q1/ypnAI1lH8XDF0kNOMnorwQTL+9nDKUnIZGrddo6yV3tfoSsBtk7wxueT\n6wGcxC1KhtuGNHKPTuXm9Sm8MqERi6c/hu7tHHztAOrFXU8upKS8jqvPO5bh/XPIyc4M6vaFZzLC\nF0E3pWkf45pKWJ+Uz9tJfcO+/33xGbyXPZyM1gauq+m6D40Z/WJIG2MyWnmvLJEXd4c/BXy5dhdz\n3lxJ/17ZXHHWURQXyInacJGEL4LOAtxYvYZERxtzso6k1Bq+E3F2LLxYMAG7xcY5ZSvIJHi3+esp\nbBZ4akwD8RaDuzelUdIQvj47DU0tzLpvPgYG9155MkW5qSQmRu8FYT2NlHRMxNPMCW/lA19bLbQr\ntDdyTfVaHs0+inszj+ZX+xaFZXSxMHsku5JyGV2zk4mN+yHZezLx2g3Ti1Df39XXWTue4ncv84we\nN4xfNVXz4PImrt2ax+ujHFitof8N3fbYAnbur+KyM47m2JFFFOTlhHyf4nsywhchc0bjLkbU72VH\ncgHvZ/qWRP2xNSmP97KHk9Vax0UHpQGXN3ccl8HRBTY+3N3G31eGvt3Ews828uybKxjUO4efXzCB\nPoWS7MNNEr4IGQtwWdlyUtuaWJA7hs2JvrZK675aWyLPFUwE4Mely8gwPM+5F9+Ls1p4/owM0hPg\nd182sWJ/6NoS79pfyRW/e4WEeBv3X3safQszpZQTAVLSET7pdFGSl9YI6fZmZpd+zmPFU5hTdAJ3\n7nmfTPv3t9z7qjawWTsT0hOxY+HZwuOpjE/l/OoNnGet/a6U497KATpfaNVJiGfleBPorB1PN0w5\nXJlqYFYCf5uaypXv1XPJgloWz7RRmJ7gU7zeNDW3MvOOF6mub+Y3V53CuCOKyM6SG5xEgozwRcgN\nayrjnPI11MUl84+iSTRbgjff2gBeyTuarckFDK/7ltl1OmjbNouZw9O4bkwie+oMfrKghpa24M3P\ndzgc/PS3L7NmSwk/OHEE509R9CnKC9r2hW8k4YuwOKV6M+NrdvBtUi5PF56A3dO9Bf2wMHskSzKH\nUNBcxeWly7DJBVZ+eWBKJtP6xvFViZ3LFlRiD8IdsgzD4NbHFvLWpxsZd0Rv7rp8MgP7FHie9y/C\nQko6IiDu5ZHD9cppr+fX2RLZlFrMnIKJzDqwzOf9TUj/vu77YaZiYc4oslvruGbfYtKw4/6W9lS+\nifZZOd6EoteO1WLhpbOyOP3VCt7e0casJQ6emzUaq5/9dgzD4K4n3+WJV75gcO8cHrlpBsP6FxAX\nJyknkmSEL8LGhsHPSj9nUOMB1qT344lek2mxxvu8HQcWXssdyxt5Y0lva+CavZ9Q6GjyvqLoUnqi\njbfPy2ZkrpXXVx3kgqdWU9/s+3UMbXYH1//5vzwy7zMGFGfzxK1nM2pwL5KSpDFapAWU8JVSP1RK\nvRSsYETPl2jY+cX+TxlZt4ctKUUsGDaD8uTuT8+rtiXxdNEkPs5S5LVUc/2ej+hjrw9hxOaSnRzH\nwvOzOaZ/Gu9vqGTKQ1+xubT7x/fbA7WcdvsbzH1rBap/PnPvOpcJI/tJF8wo4ff3K6XU34DTgNXB\nC0dEmrfeOp2+4LtdiNVpVsz+Qx9uaWwmwbBzdennvNZ2JIuzjmDBsDMYeWADI8o2kNzmeaRut1jZ\nmjOYV4qPpsGWyMCGUq4q+YzxSTbg+zKPp1k57ny9oUmkSzg+89T3xsPMncPJTYnj3RtHce2L65m/\n8iDH/vErbjy5LzedOoCMZM8po76placXrOWhV1ZSXd/CqccO5fc/m8qooX2kKVoUCaSgthR4A7g6\nSLEIE7FhcGH5GuIP7mZpvxNYVziKjflH0L9qN0V1JWQ012A1HNQlpFGWmseOrIE0xSeT4GjlzAMr\nOa1mi6sDpiSTUEiMtzHn8tGcfMRe7nlrBw+9v5vHF+3h/HEFTBqSjWrMps3uoKSino9W7+HtZds5\nWN1EenI891w0jtkzp8hVtFHIYhhdn41XSl0B3IRzBpzF9f9ZWuuVSqmTgKu11hd521FZWW34GnZE\nqfzU6LufqDv3Od+dTgi6jRTdL9kv+d/yQx67n8T11C2zxWLjs/QBfJQ1nKr4VI9xJdubGVeznWlV\nmw+p1w91a50w/qSBhzz2dDOTWB/he5uX362TuN5G/B2OSU1jK399bwf/Xn6AfdUtHhfPTE3gR5MH\ncf3ZoxlUnEN8QT/vMXhRVi+nGAHy89ODNq3J6whfaz0XmBusHQrhLsGwM61mG1NqtlGSkIlOyqPc\nlowdK2n2JgY1VzCouYLEzs2bRRhkJMdz37nD+M0PhrJsexUrd1WxrTGJ+DgrOelJTDyikBNG9iYx\nMbgXbIngkzlSImpYgeKWaopbAuxfL0LCarVw3JBsjhuSHfarj0VwSMIXXfL1Bile2xZ4OInrC/cS\nTndO0nYS4dYJUSHINzQRsSGghK+1XgwsDlIsQgghQkjOigghhElISUccwtdOje4zXtyvm/XaEHm/\ntwUO5V7CcS8hBdo2AaJvVo47b78jTyWqQNsvmLLs1QPJCF8IIUxCEr4QQpiElHSET3xtveBziccL\nr7OA3HmYjWLG8oQZf2bRmYzwhRDCJCThCyGESUhJR3TJ5/ureuE+iybQEo8ZZuV4E+zfUTBiENFJ\nRvhCCGHbSpL2AAANvElEQVQSMsIXAfH1JK57d01P3Sx94WvnSzPozmjb128BMoLvGWSEL4QQJiEJ\nXwghTEJKOsInPp8gdCuxBO1ODofZfndIeUKOgVnJCF8IIUxCEr4QQpiElHREUHWatePepTHIN94w\n4ywcIfwlI3whhDAJSfhCCGESUtIRAfHnZhwd+XpjDl9LODIbRYjvyQhfCCFMQhK+EEKYhJR0RFD5\nemFWsGfZSAlHiMOTEb4QQpiEXyN8pVQG8CKQgfMudjdrrb8MZmCiZwh1r3YZ0QvRff6O8H8FfKi1\nngLMAp4IWkRCCCFCwt8a/sNAs+vf8UBjcMIRQggRKl4TvlLqCuAmwMDZ7NAAZmmtVyqlioAXgBtC\nGqXoMQIt8UgJRwj/WQzD8GtFpdRoYB7O+v373pYvK6v1b0c9SH6q+/2ghCR8cThl9TKnBCA/Pz1o\nXcX9PWk7ApgPXKi1XhusYIQQQoSOvzX8PwKJwKNKKQtQpbX+YfDCEmYhI3Yhwsfvko6vpKQjJR0h\nfCElHadglnTkiAohhElIwhdCCJOQhC+EECYhCV8IIUxCumWGkZyEEkJEkmQgIYQwCUn4QghhEpLw\nhRDCJCThCyGESUjCF0IIk5CEL4QQJiEJXwghTCJszdOEEEJElozwhRDCJCThCyGESUjCF0IIk5CE\nL4QQJiEJXwghTEISvhBCmERY2iMrpX4IXKC1vtj1+FjgUaAV+EBr/btwxNGTKKX2AJtdD7/QWt8V\nyXhikVLKAjwJjAGagNla6+2RjSp2KaW+BqpcD3dora+MZDyxypUfH9BaT1VKDQaeAxzAOq31zwPZ\ndsgTvlLqb8BpwOoOTz8F/FBrvVMptUApNVZrvdrzFoQ715tgpdb6nEjHEuPOBRK11se7PmQPu54T\nPlJKJQKG1npapGOJZUqpXwOXAnWupx4G7tRaL1FK/UMpdY7W+k1/tx+Oks5S4Nr2B0qpdCBBa73T\n9dR7wMlhiKMnGQf0UUp9rJR6Wyk1LNIBxahJwLsAWutlwPjIhhPTxgCpSqn3lFIfuv6ACt9tBX7Y\n4fE4rfUS17/fAU4JZONBG+Erpa4AbgIMwOL6/yyt9StKqZM6LJoB1HR4XAsMDFYcPc1hjuvPgT9q\nrV9TSp0AvAhMiFyUMSsDqO7wuE0pZdVaOyIVUAxrAB7SWs9RSg0F3lFKDZNj6Rut9RtKqf4dnrJ0\n+HctkBnI9oOW8LXWc4G53Vi0BucHrV0639f9hBtPx1UplQy0uV5fqpQqjkRsPUANzvdfO0n2/tuM\nc3SK1nqLUqoc6AXsjWhUsa/j+zHgXBn2WTpa61qgWSk10HXSbDqwxMtq4lC/AW4EUEqNAXZHNpyY\ntRQ4A0ApNRFYG9lwYtoVwF8BXAOQdGB/RCPqGb5WSp3o+vcMAsyVkbqJ+TXAPJx/cN7XWi+PUByx\n6gHgRaXUmThnOl0e2XBi1hvAqUqppa7HsyIZTIybA/xLKbUE56j0Cvm2FBS3AP9USsUDG4FXA9mY\ndMsUQgiTkAuvhBDCJCThCyGESUjCF0IIk5CEL4QQJiEJXwghTEISvhBCmIQkfCGEMAlJ+EIIYRKS\n8IUQwiQk4QshhElIwhdCCJOQhC+EECYhCV8IIUxCEr4QQpiEJHwhhDAJSfhCCGESkvCFEMIkJOEL\nIYRJhO+etg3VEb+XolGxL9Ih4Ni1MdIhOG1dH+kIMLZujnQIALTqnZEOgfLNZZEOAYAd++siHQJb\nGpsjHQIAX9VGRxxPGTWWYG1LRvhCCGESkvCFEMIkJOELIYRJSMIXQgiTkIQvhBAmIQlfCCFMQhK+\nEEKYhCR8IYQwCUn4QghhEpLwhRDCJCThCyGESUjCF0IIk5CEL4QQJiEJXwghTEISvhBCmIQkfCGE\nMAmLYUT8viRCCCHCQEb4QghhEpLwhRDCJCThCyGESUjCF0IIk5CEL4QQJiEJXwghTEISvhBCmESc\nvysqpSzAk8AYoAmYrbXe3uH1q4CfAa3A/VrrBUqpXGAekATsA2ZprZs8LetvXJHgz7Ho8NqNQIHW\n+k7X40uBW4Aq4Hmt9dyw/SAB8vM9UQi8BMQDFcAlWut6pdTFwK+ANuBfWuunwvvTBM7P45ENbAbW\nuhZ7A/gU+BtgABZgInCO1vr9cP0sgfDxOPxBa71QKfUIMBbnz9wLqNRaH6+UmgHc63r+a6319eH9\naQLn5/tiAPC8a5FdwM9cufNR4Hig1vXaOVrr9n93EsgI/1wgUWt9PHAH8HCHgAuBXwDHAacDf1JK\nxeP8Rb2ktT4JWA1c3cWyscTnY6GUSlJKvQBc22HZXOD3wInAFOBipVS/sP0UgfPnPXE7zoTe/p6Y\n7VrlIWAaMAm4WSmVGbafInj8OR5HA/O01tNc//1da71Gaz1Vaz0NeAJ4LVaSvYsvx+EBpVS81vom\nrfVU4DScg5/ZSqk04M/Ama5t7XR9ZmKNP++Lh4AnXZ+TxTgHQ+B8v0zv8H45bLKHwBL+JOBdAK31\nMmB8h9cmAJ9prdu01jXAFpx/zb5bB3gHOPUwyx4ZQFyR4OuxOBLnt5zngfs7LDsIWKW1rtZaG8By\nnKO5WOHzcXB9sF9SSlmBvjg/3ABrgGwg2fU4Fi8J9+d9MQ4Yp5T6RCn1slKqqH0FpVQKcB9wQ7h+\ngCDx5zi0uwF4X2u9AedIdi3wsFLqU6BUa10ejh8gyPzJncP5PncuBSa5vikMBZ5RSn2mlJrlbceB\nJPwMoLrD4zbXh9bTa7VAJpDe4XlPzwHUuZ6PJb4cizogU2tdpbX+EOdX9HZbgJFKqXzXh/tkIDWE\ncQebz8cBQCkVh/ODPAX42PX6emCl6/m3XW/+WOPP8dgI3Ku1ngK8Cfy9wzJXAvO11hUhizg0/H1f\nxOMsbfzF9VoezvfIr4EZwE1KqSGhCztkfD0eGTi//Z7jeu4HOPNCCvAYcAnObwPXKaVGdbXjQBJ+\nDc5k/d22tNaODq9ldHgtA6h0Wye9w3Mdl03n+1FerPDlWBz259NaV+H8qvYazrr2SuBg0KMNHb+O\ng2s0MxK4GnhBKTUaOBPoDwwACpVS54c49lDw53gsAj5xPfcGzjp2u4uBZ0MSaWj5+/k4BVjcoUxR\nDizXWpdpretxntvoeHxihT/H4xbgB0qpjwEHzrzQADymtW7SWtfhHCyN6WrHgST8pcAZAEqpiXx/\nkgngK5xfORJctdcjgHWudc50LTMDWIKzbOFp2Vjiz7HoRCllA47WWp8IXOhadmkoAw8yn4+DUuoJ\npdQU1zJ1gB3nG7wBaHaVtg7gLO/EGn/eF88C7X/cTsH5Rx+lVAaQoLXeG6bYg8nfz8cpOEu/7VYC\no5RSOa5vhROBDaEOPgT8OR6nAne4zuM4gA8ABXymlLK4vg1NAr7uasd+d8vscKa5vd42C2cy36K1\nflspdSXOEZsF55nm/yqlCnDWrdNw/oW6SGvd6GlZv4KKEH+ORYd1fwqoDrN07sV5UqcR+KvW+vXw\n/SSB8fM9oYCncL6JHcD1WmutlLoauAJoBrYBV2mt28L7EwXGz+MxAGifmVWPcwZHqVJqPHCn1vq8\nsP4QQeDv50Mp9T/gLq31Nx22dSFwK85zOi9rrf9CjPHzfTEB5wn7Jpzlzp9rre1KqVtwDg5bgP/T\nWj/T1b6lPbIQQpiEXHglhBAmIQlfCCFMQhK+EEKYhCR8IYQwCUn4QghhEpLwhRDCJCThi4hTSp2k\nlPp3h8cXKKW+UUr1CcG+FimlhnXx+tWuayGCsa95Sqk4pVRfpdRZwdimEIGQhC+ihQGglPoxcBsw\nTWu9J7IhBUZrfZHrYrFpwAmRjkcIv/vhCxFkFqXUJThbw57c3izN1QzqMdcy5Tivvj0aeBDnVbj/\nxNlMazHOKxcduHqCK6X+CEzGObB5WGv9mqcdK6Um4ew3X+5a/wvX89cDF7me+4/W+nGl1L9c+x0A\nFAGXa61Xu54fhLO756OuDqA7gBE4W0AnK6W+wNkKd6jW2lBKPYCzN4zHuIQINhnhi2gxGbgKZ8+c\njvdD+CdwnauHyDs4R//g7Cd+ktb6RZzNpl5ydZjcB8xQSp0ODNRaT8Y5wr67i576TwIztdbTgR0A\nSqnhwEycI/PJwA87lIJ2aq1PBx4Hfubq0z4JOA9njyi7aznD9e8HcPa4fwtn/6jpru6IM3B2xBQi\nLGSEL6LFPpwNoq7C2Sn0dNfzw4EnnS13iMd5NygA3WFdA2f7WIBvcd5roD/OvvIf4+xJEud6zpMC\nrfU217+XAoOBUa7lP3Ktn+V6HmBVh30dr7WuU0rdhPOPUzrwYhc/57M4e7zbgA9irT+QiG0ywhfR\nYqvWukVr/QTQrJS62/X8JuAy1wj/NqD99pCODuta6HyDlI3Ax671pgHzge14ttfVxA3gGNf/NbDO\ndRehqcBzfN/V8JB9ue5SNM7V2Ows4M+uzqft9zpw4EzwaK3b/6BcAcw53MEQIhQk4YtodCXOUslJ\nOG8B+YLrDkd/Ar7xsLzh/m+t9dtAvWu9FYDh6hnuqVvgNa59fAD0c63/DfCx605Cy3HeWWifp/W1\n1qVAkVJqKfA+8JDW2t5h2bU4e5lf6Hr8ElCktd7YvcMhRHBIt0whwkwp9WugTGv9XKRjEeYiNXwh\nwsg1m6cXcHakYxHmIyN8IYQwCanhCyGESUjCF0IIk5CEL4QQJiEJXwghTEISvhBCmMT/A0sBl5Rw\n9Q9VAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "psy.close('all')\n", "psy.plot.densityreg(ds, name='y', coord='x', cmap='Reds', bins=50, density='kde', \n", " clabel='Kernel density',\n", " color='Blues_r', fit=lambda x, a: np.sin(a * x), \n", " legendlabels='$\\sin (%(a)1.2f * %(xname)s$)')" ] }, { "cell_type": "code", "execution_count": 6, "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 }