{ "metadata": { "name": "end-of-day-data-analysis" }, "nbformat": 3, "nbformat_minor": 0, "worksheets": [ { "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "% for i in `find . -type f`; do grep ^Year $i | cut -d/ -f3; done > months.txt\n", "% wc months.txt\n", " 351 351 1053 months.txt\n", "% find . -type f | wc -l\n", " 354\n", "\n" ] }, { "cell_type": "code", "collapsed": false, "input": [ "pwd" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "pyout", "prompt_number": 1, "text": [ "u'/Users/t/dev/2012-scripps/python'" ] } ], "prompt_number": 1 }, { "cell_type": "code", "collapsed": false, "input": [ "cd ../1-Shell/data" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "/Users/t/dev/2012-scripps/1-Shell/data\n" ] } ], "prompt_number": 2 }, { "cell_type": "code", "collapsed": false, "input": [ "ls months.txt" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "months.txt\r\n" ] } ], "prompt_number": 3 }, { "cell_type": "code", "collapsed": false, "input": [ "months = numpy.loadtxt('months.txt')" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 6 }, { "cell_type": "code", "collapsed": false, "input": [ "print len(months)" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "351\n" ] } ], "prompt_number": 7 }, { "cell_type": "code", "collapsed": false, "input": [ "months[0]" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "pyout", "prompt_number": 8, "text": [ "2.0" ] } ], "prompt_number": 8 }, { "cell_type": "code", "collapsed": false, "input": [ "months[1]" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "pyout", "prompt_number": 9, "text": [ "1.0" ] } ], "prompt_number": 9 }, { "cell_type": "code", "collapsed": false, "input": [ "hist(months, bins=12, range=(1,12))\n", "axis(xmin=1)" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "pyout", "prompt_number": 15, "text": [ "(1, 12.0, 0.0, 45.0)" ] }, { "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAXIAAAD9CAYAAAChtfywAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAFZRJREFUeJzt3X9MVff9x/HXWWVZnD+AjntxxQa+aUm5iHKdm8tSFyy9\nuLWB4Y+Y2lVvFJrFZcmUZq4mM8Ml02vUdMj8wyw0g2ypM0uk1FAynbnGzTjSCV3jzNiMpGgvt3OI\nwtApcL9/yHCoXC6Xe+/hc3k+khOv13s/7/dRzstzP/ee+7FCoVBIAABjfcbuBgAAU0OQA4DhCHIA\nMBxBDgCGI8gBwHAEOQAYLqIgHxoaktvtVmlpqSSpurpaWVlZcrvdcrvdamlpiWuTAIDxzYrkQTU1\nNXK5XOrr65MkWZalqqoqVVVVxbU5AMDEJjwjv3r1qpqbm1VZWan/XjsUCoXEdUQAMD1MeEa+fft2\n7d+/X7du3Rq9z7Is1dbWqqGhQcuWLdPBgweVmpo65nmWZcW+WwCYASZ7ohz2jPzEiRNyOBxyu91j\nBt66dauuXLmi9vZ2LViwQG+88ca4zSR6+/GPf2xLXTs39nlmbOzzzNiiETbIz507p6amJuXk5GjD\nhg06ffq0Nm3aJIfDIcuyZFmWKisr1draGlVxAMDUhQ3yPXv2qKurS1euXNHRo0f1wgsvqKGhQYFA\nYPQxx48fV0FBQdwbBQA8XkSfWpHuT5P8d957x44d+vDDD2VZlnJycnTkyJG4NThZRUVFdreQcOzz\nzMA+YzxWKNpJmYkGtqyo53uiMW9euvr6biSs3ty5abp1qydh9YDpimMvtqLJzqQJ8vuvFhL5kcjE\n7h8wXXHsxVY02ckl+gBgOIIcAAxHkAOA4QhyADAcQQ4AhiPIAcBwBDkAGI4gBwDDEeQAYDiCHAAM\nR5ADgOEIcgAwHEEOAIYjyAHAcBEF+dDQkNxut0pLSyVJPT098ng8ys3NVUlJiXp7e+PaJABgfBEF\neU1NjVwu1+gKQT6fTx6PRx0dHSouLpbP54trkwCA8U0Y5FevXlVzc7MqKytHv+y8qalJXq9XkuT1\netXY2BjfLgEA45pwzc7t27dr//79unXr1uh9wWBQTqdTkuR0OhUMBh/73Orq6tHbRUVFrL8HAA/x\n+/3y+/1TGiPsUm8nTpzQ+++/r8OHD8vv9+vgwYN67733lJaWphs3HqzRl56erp6esWvosdQbMDNw\n7MVWNNkZ9oz83LlzampqUnNzs+7cuaNbt25p48aNcjqd6u7uVmZmpgKBgBwOx5QaBwBEL+LFl8+c\nOaMDBw7ovffe044dO/Tkk0/qhz/8oXw+n3p7ex95w5MzcmBm4NiLrbgvvvzfT628+eabOnnypHJz\nc3X69Gm9+eabkyoKAIidiM/IJz0wZ+TAjMCxF1txPyMHAEw/BDkAGI4gBwDDEeQAYDiCHAAMR5AD\ngOEIcgAwHEEOAIYjyAHAcAQ5ABiOIAcAwxHkAGA4ghwADEeQA4DhCHIAMBxBDgCGCxvkd+7c0fLl\ny1VYWCiXy6WdO3dKkqqrq5WVlSW32y23262WlpaENAsAeNSEKwQNDAxo9uzZGhwc1PPPP68DBw7o\n97//vebOnauqqqrxB2aFIGBG4NiLrbisEDR79mxJ0t27dzU0NKS0tDRJSuq/SAAwyayJHjA8PKyl\nS5fq8uXL2rp1q/Lz8/Xb3/5WtbW1amho0LJly3Tw4EGlpqY+8tzq6urR20VFRSoqKopl7wBgPL/f\nL7/fP6UxIl58+ebNm1q1apV8Pp9cLpcyMjIkSbt27VIgEFBdXd3YgZlaAWYEjr3Yiuviy/Pnz9fL\nL7+sDz74QA6HQ5ZlybIsVVZWqrW1ddLNAgBiI2yQX79+Xb29vZKk27dv6+TJk3K73eru7h59zPHj\nx1VQUBDfLgEA4wo7Rx4IBOT1ejU8PKzh4WFt3LhRxcXF2rRpk9rb22VZlnJycnTkyJFE9QsAeEjE\nc+STHpg5cmBG4NiLrbjOkQMApieCHAAMR5ADgOEIcgAwHEEOAIYjyAHAcAQ5ABiOIAcAwxHkAGC4\nCb/GFuOZNXJFW/zNnZumW7d6ElILgHkI8qgNKlGXJff1JeY/DABmYmoFAAxHkAOA4QhyADAcc+RA\nkpk3L119fTfsbgMJxPeRR18xgfWS+/uXEVvJfSzcr5fMx0PMv4/8zp07Wr58uQoLC+VyubRz505J\nUk9Pjzwej3Jzc1VSUjK6HBwAIPEmPCMfGBjQ7NmzNTg4qOeff14HDhxQU1OTvvCFL2jHjh3at2+f\nbty4IZ/PN3ZgzshjWiuZz0AQW8l9LNyvl8zHQ1xWCJo9e7Yk6e7duxoaGlJaWpqamprk9XolSV6v\nV42NjVG0CwCIhQnf7BweHtbSpUt1+fJlbd26Vfn5+QoGg3I6nZIkp9OpYDD42OcWFRWN3s7OzlZ2\ndnZMmgaAZOH3++X3+6c0RsRvdt68eVOrVq3S3r17tWbNGt248eBd8fT0dPX0jL2E/P7Lu42S/m9K\nDUZut5L35WRyv5REbDG1YrZoplYi/vjh/Pnz9fLLL+vPf/6znE6nuru7lZmZqUAgIIfDMc6zvJKK\nJ9VQ9HYnqA4ATC9h58ivX78++omU27dv6+TJk3K73SorK1N9fb0kqb6+XuXl5fHvFADwWGHPyAOB\ngLxer4aHhzU8PKyNGzequLhYbrdb69evV11dnbKzs3Xs2LFE9QsAeEhcLwiSTilxUyvJPC+Y3HOC\nyc6eKy2T9Vi4Xy+Zj4e4zpEDiM79EE90sGIm4UuzAMBwBDkAGI4gBwDDMUeOGYevecVkmPDzQpBj\nxuHNR0yGCT8vTK0AgOEIcgAwHEEOAIYjyAHAcAQ5ABiOIAcAwxHkAGA4ghwADEeQA4DhCHIAMFzY\nIO/q6tLKlSuVn5+vRYsW6dChQ5Kk6upqZWVlye12y+12q6WlJSHNAgAeFXaFoO7ubnV3d6uwsFD9\n/f360pe+pMbGRh07dkxz585VVVXV+AOzQlBMayXziiiJNhNWmU/2eok8Huz4eYnpCkGZmZnKzMyU\nJM2ZM0d5eXm6du2aJBEsADBNRDxH3tnZqba2Nn31q1+VJNXW1mrJkiWqqKhQb2/vOM9qkFQ9svmn\n1ikAJCW/HuRkdVQjRLT4cn9/v4qKivSjH/1I5eXl+vTTT5WRkSFJ2rVrlwKBgOrq6sYOzNRKTGvx\nCih2mFoxvx5TK2NNeEZ+7949rV27Vq+99prKy8slSQ6HQ5ZlybIsVVZWqrW1Nbp+AQBTFjbIQ6GQ\nKioq5HK5tG3bttH7A4HA6O3jx4+roKAgfh0CAMIK+2bnH//4R/3qV7/S4sWL5Xa7JUl79uzRO++8\no/b2dlmWpZycHB05ciQhzQIAHhXRHHlUAzNHHtNazJHHDnPkptdLkTSYwHrSdJ8jZ81OAIYZ1HRf\nQzPRuEQfAAxHkAOA4QhyADAcQY5HzJuXPnqdQCK2efPS7d5lwGi82YlH9PXdUCLfTOrrm/5vJgHT\nGWfkAGA4ghwADEeQA4DhCHIAMBxBDgCGI8gBwHAEOQAYjiAHAMNxQZARZo189SoAPIogNwJf2wlg\nfGGnVrq6urRy5Url5+dr0aJFOnTokCSpp6dHHo9Hubm5KikpUW9vb0KaBQA8KmyQp6Sk6K233tLF\nixd1/vx5HT58WJcuXZLP55PH41FHR4eKi4vl8/kS1S8A4CFhgzwzM1OFhYWSpDlz5igvL0/Xrl1T\nU1OTvF6vJMnr9aqxsTH+nQIAHiviOfLOzk61tbVp+fLlCgaDcjqdkiSn06lgMDjOsxoknR25XTSy\nAQAe8I9s0YsoyPv7+7V27VrV1NRo7ty5Y/7sv98p/XiblLjFlwHAREUae5K7e9IjTPg58nv37mnt\n2rXauHGjysvLJd0/C+/u7pYkBQIBORyOSRcGAMRG2CAPhUKqqKiQy+XStm3bRu8vKytTfX29JKm+\nvn404AEAiWeFQqFxP6D8hz/8QV//+te1ePHi0emTvXv36itf+YrWr1+vjz/+WNnZ2Tp27JhSU1PH\nDmxZkk4pcVMrlhL/WetE1UvmfbtfL8yPYeyrWcn/90k9s+tN9ngIG+RTaoUgN7SWPfUIcupR70G9\nyR4PfNcKABiOIAcAwxHkAGA4ghwADEeQA4DhCHIAMBxBDgCGI8gBwHCsEIRpgKXsgKkgyDENsJQd\nMBVMrQCA4QhyADAcQQ4AhiPIAcBwBDkAGG7CIN+yZYucTqcKCgpG76uurlZWVpbcbrfcbrdaWlri\n2iQAYHwTBvnmzZsfCWrLslRVVaW2tja1tbXpG9/4RtwaBACEN2GQr1ixQmlpaY/cn8gVXQAA44v6\ngqDa2lo1NDRo2bJlOnjw4CNrdt7XIOnsyO2ikQ0A8IB/ZIteRGt2dnZ2qrS0VB999JEk6dNPP1VG\nRoYkadeuXQoEAqqrqxs7MGt2GlqLetSjnt31ErJmp8PhkGVZsixLlZWVam1tjWYYAEAMRBXkgUBg\n9Pbx48fHfKIFAJBYE86Rb9iwQWfOnNH169e1cOFC7d69W36/X+3t7bIsSzk5OTpy5EgiegUAPEZE\nc+RRDcwcuaG1qEc96tldLyFz5ACA6YMgBwDDEeQAYDiCHAAMR5ADgOEIcgAwHEEOAIYjyAHAcAQ5\nABiOIAcAwxHkAGA4ghwADEeQA4DhCHIAMBxBDgCGI8gBwHBhg3zLli1yOp1jlnLr6emRx+NRbm6u\nSkpK1NvbG/cmAQDjCxvkmzdvVktLy5j7fD6fPB6POjo6VFxcLJ/PF9cGAQDhhQ3yFStWKC0tbcx9\nTU1N8nq9kiSv16vGxsb4dQcAmNCEiy8/LBgMyul0SpKcTqeCwWCYRzdIOjtyu2hkAwA84B/Zojfp\nIP9flmWNLLI8nk1K3OLLAGCiIo09yd096REm/akVp9Op7u5uSVIgEJDD4Zh0UQBA7Ew6yMvKylRf\nXy9Jqq+vV3l5ecybAgBEzgqFQqHx/nDDhg06c+aMrl+/LqfTqZ/85Cf61re+pfXr1+vjjz9Wdna2\njh07ptTU1EcHtixJp5S4qRVL0ri7Yni9ZN436lGPeg/XCxPLj39GuCCfUisEuaG1qEc96tldb7Kx\nzJWdAGA4ghwADEeQA4DhCHIAMBxBDgCGI8gBwHAEOQAYjiAHAMMR5ABgOIIcAAxHkAOA4QhyADAc\nQQ4AhiPIAcBwBDkAGI4gBwDDTWnx5ezsbM2bN09PPPGEUlJS1NraGqu+AAARmlKQW5Ylv9+v9PT0\nWPUDAJikKU+txGmlOABAhKZ8Rv7iiy/qiSee0He+8x29/vrrDz2iQdLZkdtFIxsA4AH/yBa9KS2+\nHAgEtGDBAv3zn/+Ux+NRbW2tVqxYcX9gFl82tBb1qEc9u+sldPHlBQsWSJIyMjK0evVq3uwEABtE\nHeQDAwPq6+uTJP373//W7373OxUUFMSsMQBAZKKeIw8Gg1q9erUkaXBwUN/+9rdVUlISs8YAAJGZ\n0hx52IGZIze0FvWoRz276yV0jhwAYD+CHAAMR5ADgOEIcgAwHEEOAIYjyAHAcAQ5ABiOIAcAwxHk\nAGA4ghwADEeQA4DhCHIAMBxBDgCGI8gBwHAEOQAYLgmD3G93Azbw292ADfx2N2ADv90N2MBvdwNG\niDrIW1pa9Nxzz+nZZ5/Vvn37YtnTFPntbsAGfrsbsIHf7gZs4Le7ARv47W7ACFEF+dDQkL73ve+p\npaVFf/3rX/XOO+/o0qVLse4NABCBqIK8tbVVzzzzjLKzs5WSkqJXXnlF7777bqx7AwBEIKrFl69d\nu6aFCxeO/j4rK0t/+tOfHvPIF6PtK0rWyK+7E1xvOtSK9T4nct+irTeVfTZh/x4n0n02df8e53H7\nnEz7N3VRBfn9hZXDi9OazgCAh0Q1tfLUU0+pq6tr9PddXV3KysqKWVMAgMhFFeTLli3T3//+d3V2\nduru3bv6zW9+o7Kyslj3BgCIQFRTK7NmzdLPf/5zrVq1SkNDQ6qoqFBeXl6sewMARCDqz5F/85vf\n1N/+9jf94x//0M6dO2PZU1S6urq0cuVK5efna9GiRTp06JDdLSXM0NCQ3G63SktL7W4lIXp7e7Vu\n3Trl5eXJ5XLp/PnzdrcUd3v37lV+fr4KCgr06quv6j//+Y/dLcXcli1b5HQ6VVBQMHpfT0+PPB6P\ncnNzVVJSot7eXhs7jL3H7fMPfvAD5eXlacmSJVqzZo1u3rw54ThJc2VnSkqK3nrrLV28eFHnz5/X\n4cOHZ8xn22tqauRyuSJ6EzoZfP/739dLL72kS5cu6S9/+UvSvxrs7OzUL37xC124cEEfffSRhoaG\ndPToUbvbirnNmzerpaVlzH0+n08ej0cdHR0qLi6Wz+ezqbv4eNw+l5SU6OLFi/rwww+Vm5urvXv3\nTjhO0gR5ZmamCgsLJUlz5sxRXl6ePvnkE5u7ir+rV6+qublZlZWVM+KTQjdv3tTZs2e1ZcsWSfen\n+ebPn29zV/E1b948paSkaGBgQIODgxoYGNBTTz1ld1sxt2LFCqWlpY25r6mpSV6vV5Lk9XrV2Nho\nR2tx87h99ng8+sxn7kfz8uXLdfXq1QnHSZog/1+dnZ1qa2vT8uXL7W4l7rZv3679+/eP/sMnuytX\nrigjI0ObN2/W0qVL9frrr2tgYMDutuIqPT1db7zxhp5++ml98YtfVGpqql58MdHXaNgjGAzK6XRK\nkpxOp4LBoM0dJdbbb7+tl156acLHJd3R39/fr3Xr1qmmpkZz5syxu524OnHihBwOh9xu94w4G5ek\nwcFBXbhwQd/97nd14cIFff7zn0+6l9sPu3z5sn72s5+ps7NTn3zyifr7+/XrX//a7rYSzrKsGTN9\nKEk//elP9dnPflavvvrqhI9NqiC/d++e1q5dq9dee03l5eV2txN3586dU1NTk3JycrRhwwadPn1a\nmzZtsrutuMrKylJWVpa+/OUvS5LWrVunCxcu2NxVfH3wwQf62te+pieffFKzZs3SmjVrdO7cObvb\nSgin06nu7m5JUiAQkMPhsLmjxPjlL3+p5ubmiP/DTpogD4VCqqiokMvl0rZt2+xuJyH27Nmjrq4u\nXblyRUePHtULL7yghoYGu9uKq8zMTC1cuFAdHR2SpFOnTik/P9/mruLrueee0/nz53X79m2FQiGd\nOnVKLpfL7rYSoqysTPX19ZKk+vr6GXGC1tLSov379+vdd9/V5z73ucieFEoSZ8+eDVmWFVqyZEmo\nsLAwVFhYGHr//fftbith/H5/qLS01O42EqK9vT20bNmy0OLFi0OrV68O9fb22t1S3O3bty/kcrlC\nixYtCm3atCl09+5du1uKuVdeeSW0YMGCUEpKSigrKyv09ttvh/71r3+FiouLQ88++2zI4/GEbty4\nYXebMfXwPtfV1YWeeeaZ0NNPPz2aY1u3bp1wHCsUmiGTqwCQpJJmagUAZiqCHAAMR5ADgOEIcgAw\nHEEOAIYjyAHAcP8PolP4SSc1uVIAAAAASUVORK5CYII=\n" } ], "prompt_number": 15 }, { "cell_type": "code", "collapsed": false, "input": [ "months == 0" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "pyout", "prompt_number": 13, "text": [ "array([False, False, False, False, False, False, False, False, False,\n", " False, False, False, False, False, False, False, False, False,\n", " False, False, False, False, False, False, False, False, False,\n", " False, False, False, False, False, False, False, False, False,\n", " False, False, False, False, False, False, False, False, False,\n", " False, False, False, False, False, False, False, False, False,\n", " False, False, False, False, False, False, False, False, False,\n", " False, False, False, False, False, False, False, False, False,\n", " False, False, False, False, False, False, False, False, False,\n", " False, False, False, False, False, False, False, False, False,\n", " False, False, False, False, False, False, False, False, False,\n", " False, False, False, False, False, False, False, False, False,\n", " False, False, False, False, False, False, False, False, False,\n", " False, False, False, False, False, False, False, False, False,\n", " False, False, False, False, False, False, False, False, False,\n", " False, False, False, False, False, False, False, False, False,\n", " False, False, False, False, False, False, False, False, False,\n", " False, False, False, False, False, False, False, False, False,\n", " False, False, False, False, False, False, False, False, False,\n", " False, False, False, False, False, False, False, False, False,\n", " False, False, False, False, False, False, False, False, False,\n", " False, False, False, False, False, False, False, False, False,\n", " False, False, False, False, False, False, False, False, False,\n", " False, False, False, False, False, False, False, False, False,\n", " False, False, False, False, False, False, False, False, False,\n", " False, False, False, False, False, False, False, False, False,\n", " False, False, False, False, False, False, False, False, False,\n", " False, False, False, False, False, False, False, False, False,\n", " False, False, False, False, False, False, False, False, False,\n", " False, False, False, False, False, False, False, False, False,\n", " False, False, False, False, False, False, False, False, False,\n", " False, False, False, False, False, False, False, False, False,\n", " False, False, False, False, False, False, False, False, False,\n", " False, False, False, False, False, False, False, False, False,\n", " False, False, False, False, False, False, False, False, False,\n", " False, False, False, False, False, False, False, False, False,\n", " False, False, False, False, False, False, False, False, False,\n", " False, False, False, False, False, False, False, False, False,\n", " False, False, False, False, False, False, False, False, False], dtype=bool)" ] } ], "prompt_number": 13 }, { "cell_type": "code", "collapsed": false, "input": [], "language": "python", "metadata": {}, "outputs": [] } ], "metadata": {} } ] }