diff --git a/module3/exo2/exercice.ipynb b/module3/exo2/exercice.ipynb index e272d97bf92755160406335bb31c40c1657d0eb1..09fc344269139bf063b7fc1107ba7c3d14292d91 100644 --- a/module3/exo2/exercice.ipynb +++ b/module3/exo2/exercice.ipynb @@ -3254,6 +3254,147 @@ "sorted_data['inc'].plot()" ] }, + { + "cell_type": "code", + "execution_count": 49, + "metadata": {}, + "outputs": [], + "source": [ + "first_september_week = [pd.Period(pd.Timestamp(y, 7, 1), 'W')\n", + " for y in range(1991,\n", + " sorted_data.index[-1].year)]" + ] + }, + { + "cell_type": "code", + "execution_count": 50, + "metadata": {}, + "outputs": [], + "source": [ + "year = []\n", + "yearly_incidence = []\n", + "for week1, week2 in zip(first_september_week[:-1],\n", + " first_september_week[1:]):\n", + " one_year = sorted_data['inc'][week1:week2-1]\n", + " assert abs(len(one_year)-52) < 2\n", + " yearly_incidence.append(one_year.sum())\n", + " year.append(week2.year)\n", + "yearly_incidence = pd.Series(data=yearly_incidence, index=year)" + ] + }, + { + "cell_type": "code", + "execution_count": 51, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 51, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY0AAAD8CAYAAACLrvgBAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAGhRJREFUeJzt3X9wXeV95/H3xxaxSWqDDTILNkbewcOsIRMTawzZdLtLaGyz7WB3Ch0tNGg23jFLSDfpP6m9YcezQGZCZ2dpPQwkTAkY0wAuWwZvWS9RzGbKbr225QAFQ1wpgRjHrqWsHGzaWo3gu3/cR/hI6Mcj6Ur3XOnzmrlzz/3e8zz3PDqa+73neZ5zjiICMzOzHLNqvQFmZlY/nDTMzCybk4aZmWVz0jAzs2xOGmZmls1Jw8zMsjlpmJlZNicNMzPL5qRhZmbZGmq9AdV24YUXRlNTU603w8ysrhw8ePDnEdE42nrTLmk0NTXR3t5e680wM6srkn6as567p8zMLJuThpmZZXPSMDOzbE4aZmaWzUnDzMyyOWmY2Ud0nTrD73x7L12nz9R6U6xknDTM7CO27engwNs9bPt+R603xUpm2p2nYWbjd8Vdu+nt++DD10/sO8IT+44wp2EWh++9oYZbZmXhIw0z+9BLX7uOG1dewtxzKl8Nc8+ZxfqVl/DSH1xX4y2zsnDSMLMPLZo/l3lzGujt+4A5DbPo7fuAeXMaWDRvbq03zUrC3VNmNsDP3+vl1msu45bVS/nu/iN0ezDcChQRtd6Gqmpubg5fe8rMbGwkHYyI5tHWc/dUgacZmpmNzEmjwNMMzcxGlpU0JP2+pEOSXpf0pKS5khZKapPUkZ4XFNbfIqlT0mFJawvxVZJeS+9tk6QUnyPp6RTfJ6mpUKY1fUaHpNbqNf2sK+7aTdPm53li3xEiKtMMmzY/zxV37Z6MjzMzq1ujJg1Ji4H/ADRHxFXAbKAF2AzsiYjlwJ70Gkkr0vtXAuuAByXNTtU9BGwClqfHuhTfCJyMiMuB+4H7Ul0Lga3ANcBqYGsxOVWLpxmameXJ7Z5qAM6V1AB8HDgGrAe2p/e3AxvS8nrgqYjojYi3gE5gtaSLgfkRsTcqo++PDyrTX9czwPXpKGQt0BYRPRFxEmjjbKKpGk8ztJnAY3ZWDaMmjYj4GfBfgCPAceDdiPgecFFEHE/rHAcWpSKLgXcKVRxNscVpeXB8QJmI6APeBS4Yoa4BJG2S1C6pvbu7e7QmDal/muGzX/ost15zGd3v9Y6rHisnf2F6zM6qY9TzNFJ30HpgGfAL4M8k/e5IRYaIxQjx8ZY5G4h4GHgYKlNuR9i2YX37C2dnmt274arxVGElVvzCvPe3PlnrzZlSvjSIVVNO99SvA29FRHdE/BL4c+CfAydSlxPpuSutfxS4tFB+CZXurKNpeXB8QJnUBXYe0DNCXVYC9fDr3ZMcPGZn1ZWTNI4A10r6eBpnuB54E9gF9M9magWeS8u7gJY0I2oZlQHv/akL67Ska1M9tw0q01/XTcCLadzjBWCNpAXpiGdNilkJ1EN3h78wPWZn1TVq91RE7JP0DPBDoA94mUpX0K8AOyVtpJJYbk7rH5K0E3gjrX9nRLyfqrsDeAw4F9idHgCPADskdVI5wmhJdfVIugc4kNa7OyJ6JtRim7B66u7wF2aFLw1i1eLLiNiYdZ06w73/402+d+hvOfPLD5h7zizWXvlP+Ppv/LNSfhnfvqOdxnlzB3xhFsewzCz/MiK+YKGNWb39ep/ukxy6Tp3hy0++zAO3XF3afWDThy8jYuPiKcrlUQ9jSzZ9uHvKrE4NHlvqV8axJSs/X+XWbJrzzDCrBScNszpVb2NLNj14INysjnkqrU01j2mYmZnHNMzMrPqcNMzMLJuThpmZZXPSMDMbRj1cyXmqOWmYmQ3DZ9t/lKfcmpkNUk9Xcp5qPtIws1KqZdeQz7YfnpOGmZVSLbuGfLb98Nw9ZWalUpauIZ9tPzSfEW5mpVJvN/maLnxGuJnVJXcNlZu7p8ysdNw1VF7unjIzM3dPmZlZ9TlpmJlZNicNMzPL5qRhZmbZnDTMzCybk4aZmWUbNWlIukLSK4XHKUlflbRQUpukjvS8oFBmi6ROSYclrS3EV0l6Lb23TZJSfI6kp1N8n6SmQpnW9Bkdklqr23wzMxuLUZNGRByOiJURsRJYBfw98CywGdgTEcuBPek1klYALcCVwDrgQUmzU3UPAZuA5emxLsU3Aicj4nLgfuC+VNdCYCtwDbAa2FpMTmZmNrXG2j11PfDjiPgpsB7YnuLbgQ1peT3wVET0RsRbQCewWtLFwPyI2BuVMwofH1Smv65ngOvTUchaoC0ieiLiJNDG2URjZmZTbKxJowV4Mi1fFBHHAdLzohRfDLxTKHM0xRan5cHxAWUiog94F7hghLpsDHzLSjOrluykIeljwI3An4226hCxGCE+3jLFbdskqV1Se3d39yibVx319EXsW1aaWbWM5UjjBuCHEXEivT6RupxIz10pfhS4tFBuCXAsxZcMER9QRlIDcB7QM0JdA0TEwxHRHBHNjY2NY2jS+NXDF/EVd+2mafPzPLHvCBGV+xI0bX6eK+7aXetNs1HU048Sm1nGkjT+DWe7pgB2Af2zmVqB5wrxljQjahmVAe/9qQvrtKRr03jFbYPK9Nd1E/BiGvd4AVgjaUEaAF+TYjVTT1/EvmVl/aqHHyU2M2VdGl3Sx4HPA7cXwt8EdkraCBwBbgaIiEOSdgJvAH3AnRHxfipzB/AYcC6wOz0AHgF2SOqkcoTRkurqkXQPcCCtd3dE9IyjnVXz0teuG/YGMWXj+xLUn7Lctc5sOFlJIyL+nsrAdDH2/6jMphpq/W8A3xgi3g5cNUT8DCnpDPHed4Dv5GznVKi3L2Lfl6A8uk6d4ctPvswDt1w97P9LPf0osZnJN2Eah3r6Iv72F85eHv/eDR/J1zaFil1O9/7WJ4dcp95+lJRFTkK26vBNmMwm2eAup37DdTndvqOdxnlzB/woKSZ/+6i7nn2NP91/hFtXLx02IdvIcm/C5KRhNgE5v3C7Tp0ZtsvJv4onZqwJ2YbnO/eZTYGcWU7ucpo8niE49TymYTYOY53lVE/jYPXECXnqOWmYDTIZs5w8IWHyOCFPLScNs0E8y6m+OCFPLScNs8RdTmaj8+wps8SznGwm8+wpszFyl9P4+OKKM4uThllBf5fTs1/6LLdecxnd7/XWepNKzxdXnFncPTXJfHkDm658Yt304u6pkvCvsLFxV0f98Il1M5OTxiSpp/tuFNX6S9tJtn54DGhmctKYJOP5FVbrL2yo3Zd2vSbZmc5jQDOPz9OYJOP5FZZzUtlkqfXNf3wfifrkE+tmHieNSZR78letv7Ch9l/a7uowqw9OGpMo91dYrb+woRxf2j7D2qz8nDRKoAxf2FD7L213dZiVn5NGSdT6Cxv8pW1mo/PJfWZm5pP7zMys+pw0zGzKlOFcJJsYJw0zmzI+47/+eSDczCZdGc5FsurwkYaZTTpf3HD6cNIws0lXlnORbOKykoak8yU9I+lHkt6U9BlJCyW1SepIzwsK62+R1CnpsKS1hfgqSa+l97ZJUorPkfR0iu+T1FQo05o+o0NSa/WabmZTyRc3nCYiYtQHsB34d2n5Y8D5wB8Cm1NsM3BfWl4BvArMAZYBPwZmp/f2A58BBOwGbkjxLwHfSsstwNNpeSHwk/S8IC0vGGlbV61aFVY+J979h7j5W38VJ079Q603xWY4/y8ODWiPjHww6pGGpPnArwGPpCTzjxHxC2B9Sib9SWVDWl4PPBURvRHxFtAJrJZ0MTA/IvamDXx8UJn+up4Brk9HIWuBtojoiYiTQBuwbrRttvLxrBkri8n4X5ysqcRlnKKcM3vqnwLdwKOSPgUcBL4CXBQRxwEi4rikRWn9xcD/LZQ/mmK/TMuD4/1l3kl19Ul6F7igGB+ijNUBz5qxspjM/8XJuq1BLW+XMJycpNEAfBr4vYjYJ+mPqXRHDUdDxGKE+HjLnP1AaROwCWDp0qUjbJpNtTJcwdcMJud/cbISUZl/bOUMhB8FjkbEvvT6GSpJ5ETqciI9dxXWv7RQfglwLMWXDBEfUEZSA3Ae0DNCXQNExMMR0RwRzY2NjRlNqn9lPGwdimfNWFlMxv/iZE0lLvMU5VGTRkT8LfCOpCtS6HrgDWAX0D+bqRV4Li3vAlrSjKhlwHJgf+rKOi3p2jRecdugMv113QS8mMY9XgDWSFqQZmetSbEZr57GCDxrxsqi2v+Lk/WjqMw/trKucitpJfAnVGZO/QT4t1QSzk5gKXAEuDkietL6Xwe+CPQBX42I3SneDDwGnEtl9tTvRURImgvsAK6mcoTREhE/SWW+CPzHtCnfiIhHR9rW6X6V28GHrf3KcNhqNhPdvqOdxnlzB9zWoHibgbLVO5zcq9z60uh1puvUmWH7ZcvwK8TM6pMvjT5Nlfmw1cymP1+wsA6V4S5/ZjYzuXvKzMzcPWVmZtXnpGFmZtmcNMzMLJuThpmZZXPSMDOzbE4aZmaWzUnDzMyyOWmYmVk2Jw0zsylUL7c1GI6ThpnZFKqn2xoMxdeesrrVdeoMX37yZR645WpfsNFKr8x34xsLH2lY3ar3X2w2s5T5bnxj4SMNqzvT5RebzSzT5bYGPtKwujNdfrHZzDMdbn3sIw2rO9PlF5vNPMXbtd674aoabsn4OWlYXfKNqMxqwzdhMjMz34TJzMyqz0nDzMyyOWmYmVk2Jw0zM8vmpGFmZtmcNMzMLJuThpmZZctKGpLelvSapFcktafYQkltkjrS84LC+lskdUo6LGltIb4q1dMpaZskpfgcSU+n+D5JTYUyrekzOiS1VqvhZmY2dmM50rguIlYWTv7YDOyJiOXAnvQaSSuAFuBKYB3woKTZqcxDwCZgeXqsS/GNwMmIuBy4H7gv1bUQ2ApcA6wGthaTk5mZTa2JdE+tB7an5e3AhkL8qYjojYi3gE5gtaSLgfkRsTcqp6E/PqhMf13PANeno5C1QFtE9ETESaCNs4nGzMymWG7SCOB7kg5K2pRiF0XEcYD0vCjFFwPvFMoeTbHFaXlwfECZiOgD3gUuGKEuMzOrgdwLFn42Io5JWgS0SfrRCOtqiFiMEB9vmbMfWElkmwCWLl06wqaZmdlEZB1pRMSx9NwFPEtlfOFE6nIiPXel1Y8ClxaKLwGOpfiSIeIDykhqAM4Dekaoa/D2PRwRzRHR3NjYmNMkMzMbh1GThqRPSJrXvwysAV4HdgH9s5lagefS8i6gJc2IWkZlwHt/6sI6LenaNF5x26Ay/XXdBLyYxj1eANZIWpAGwNekmJmZ1UBO99RFwLNpdmwD8N2I+J+SDgA7JW0EjgA3A0TEIUk7gTeAPuDOiHg/1XUH8BhwLrA7PQAeAXZI6qRyhNGS6uqRdA9wIK13d0T0TKC9ZmY2Ab6fhpmZ+X4aZmZWfU4aZmaWzUnDzMyyOWmYmVk2Jw0zM8vmpGFmZtmcNMzMLJuThpmZZXPSMDOzbE4aZmaWzUnDzMyyOWmYmVk2Jw0zM8vmpGFmZtmcNMzMLJuThpmZZXPSMDOzbE4aZmaWzUnDzMyyOWmYmVk2Jw0zM8vmpGFmZtmcNMzMLJuThpmZZXPSMDOzbE4aZmaWzUnDzMyyZScNSbMlvSzpL9LrhZLaJHWk5wWFdbdI6pR0WNLaQnyVpNfSe9skKcXnSHo6xfdJaiqUaU2f0SGptRqNNjOz8RnLkcZXgDcLrzcDeyJiObAnvUbSCqAFuBJYBzwoaXYq8xCwCVieHutSfCNwMiIuB+4H7kt1LQS2AtcAq4GtxeRkZmZTKytpSFoC/AbwJ4XwemB7Wt4ObCjEn4qI3oh4C+gEVku6GJgfEXsjIoDHB5Xpr+sZ4Pp0FLIWaIuInog4CbRxNtGYmdkUyz3S+CPga8AHhdhFEXEcID0vSvHFwDuF9Y6m2OK0PDg+oExE9AHvAheMUNcAkjZJapfU3t3dndkkMzMbq1GThqTfBLoi4mBmnRoiFiPEx1vmbCDi4YhojojmxsbGzM00M7OxyjnS+Cxwo6S3gaeAz0l6AjiRupxIz11p/aPApYXyS4BjKb5kiPiAMpIagPOAnhHqMjOzGhg1aUTElohYEhFNVAa4X4yI3wV2Af2zmVqB59LyLqAlzYhaRmXAe3/qwjot6do0XnHboDL9dd2UPiOAF4A1khakAfA1KWZmZjXQMIGy3wR2StoIHAFuBoiIQ5J2Am8AfcCdEfF+KnMH8BhwLrA7PQAeAXZI6qRyhNGS6uqRdA9wIK13d0T0TGCbzcxsAlT5QT99NDc3R3t7e603w8ysrkg6GBHNo63nM8LNzCybk4aZmWVz0jAzs2xOGmZmls1Jw8zMsjlpmJlZNicNMzPL5qRhZmbZnDTMzCybk4aZmWVz0jAzs2xOGmZmls1Jw8zMsjlpmJlZNicNMzPL5qRhZmbZnDTMzCybk4aZmWVz0jAzs2xOGmZmls1Jw8zMsjlpmJlZNicNMzPL5qRhZmbZnDTMzCybk4aZmWVz0jAzs2yjJg1JcyXtl/SqpEOS/nOKL5TUJqkjPS8olNkiqVPSYUlrC/FVkl5L722TpBSfI+npFN8nqalQpjV9Roek1mo23szMxibnSKMX+FxEfApYCayTdC2wGdgTEcuBPek1klYALcCVwDrgQUmzU10PAZuA5emxLsU3Aicj4nLgfuC+VNdCYCtwDbAa2FpMTmZmNrVGTRpR8V56eU56BLAe2J7i24ENaXk98FRE9EbEW0AnsFrSxcD8iNgbEQE8PqhMf13PANeno5C1QFtE9ETESaCNs4nGzMymWNaYhqTZkl4Buqh8ie8DLoqI4wDpeVFafTHwTqH40RRbnJYHxweUiYg+4F3gghHqGrx9myS1S2rv7u7OaZKZmY1DVtKIiPcjYiWwhMpRw1UjrK6hqhghPt4yxe17OCKaI6K5sbFxhE0zM7OJGNPsqYj4BfADKl1EJ1KXE+m5K612FLi0UGwJcCzFlwwRH1BGUgNwHtAzQl1mZlYDObOnGiWdn5bPBX4d+BGwC+ifzdQKPJeWdwEtaUbUMioD3vtTF9ZpSdem8YrbBpXpr+sm4MU07vECsEbSgjQAvibFzMysBhoy1rkY2J5mQM0CdkbEX0jaC+yUtBE4AtwMEBGHJO0E3gD6gDsj4v1U1x3AY8C5wO70AHgE2CGpk8oRRkuqq0fSPcCBtN7dEdEzkQabmdn4qfKDfvpobm6O9vb2Wm+GmVldkXQwIppHW89nhJuZWTYnDTOzaaDr1Bl+59t76Tp9ZlI/x0nDzGwa2LangwNv97Dt+x2T+jk5A+FmZlZSV9y1m96+Dz58/cS+Izyx7whzGmZx+N4bqv55PtIwM6tjL33tOm5ceQlzz6l8nc89ZxbrV17CS39w3aR8npOGmVkdWzR/LvPmNNDb9wFzGmbR2/cB8+Y0sGje3En5PHdPmZnVuZ+/18ut11zGLauX8t39R+iexMFwn6dhZmY+T8PMzKrPScPMzLI5aZiZWTYnDTMzy+akYWZm2Zw0zMws27SbciupG/hprbdjAi4Efl7rjZgk07Vtblf9ma5tm0i7LouIUe+XPe2SRr2T1J4zV7oeTde2uV31Z7q2bSra5e4pMzPL5qRhZmbZnDTK5+Fab8Akmq5tc7vqz3Rt26S3y2MaZmaWzUcaZmaWzUljCkj6jqQuSa8XYp+StFfSa5L+u6T5Kf4xSY+m+KuS/lWhzA8kHZb0SnosqkFzPiTpUkn/S9Kbkg5J+kqKL5TUJqkjPS8olNkiqTO1Y20hviq1uVPSNkmqRZvStlSzXaXZZ2Ntl6QL0vrvSXpgUF2l2V9pe6rZtnreZ5+XdDDtm4OSPleoqzr7LCL8mOQH8GvAp4HXC7EDwL9My18E7knLdwKPpuVFwEFgVnr9A6C51u0ptOFi4NNpeR7wN8AK4A+BzSm+GbgvLa8AXgXmAMuAHwOz03v7gc8AAnYDN0yTdpVmn42jXZ8AfhX498ADg+oqzf6ahLbV8z67GrgkLV8F/Kza+8xHGlMgIv4S6BkUvgL4y7TcBvx2Wl4B7EnluoBfAKWcTx4RxyPih2n5NPAmsBhYD2xPq20HNqTl9cBTEdEbEW8BncBqSRcD8yNib1T+ux8vlJly1WrX1G716Mbaroj4u4j438CAO/qUbX9B9dpWNuNo18sRcSzFDwFzJc2p5j5z0qid14Eb0/LNwKVp+VVgvaQGScuAVYX3AB5Nh8z/qdZdAkWSmqj8ytkHXBQRx6HyT0/liAkq/+zvFIodTbHFaXlwvOYm2K5+pdtnme0aTmn3F0y4bf2mwz77beDliOilivvMSaN2vgjcKekglcPOf0zx71DZoe3AHwF/BfSl926NiE8C/yI9vjClWzwMSb8C/DfgqxFxaqRVh4jFCPGaqkK7oIT7bAztGraKIWI1319QlbbBNNhnkq4E7gNu7w8Nsdq49pmTRo1ExI8iYk1ErAKepNIPTkT0RcTvR8TKiFgPnA90pPd+lp5PA9+lBF0gks6h8s/8pxHx5yl8Ih0O93dldKX4UQYeNS0BjqX4kiHiNVOldpVun42xXcMp3f6CqrWt7veZpCXAs8BtEfHjFK7aPnPSqJH+GRmSZgF3Ad9Krz8u6RNp+fNAX0S8kbqrLkzxc4DfpNLFVTPpsP0R4M2I+K+Ft3YBrWm5FXiuEG9JfazLgOXA/nR4fVrStanO2wplply12lW2fTaOdg2pbPsLqte2et9nks4Hnge2RMT/6V+5qvtsKmYAzPQHlSOJ48AvqWT8jcBXqMyE+Bvgm5w90bIJOExlwOv7VK48CZXZHgeBv6YywPXHpBk6NWzXr1I5xP1r4JX0+NfABVQG8zvS88JCma9TOao6TGH2BpXB/tfTew/0/z3quV1l22fjbNfbVCZxvJf+d1eUbX9Vs231vs+o/AD9u8K6rwCLqrnPfEa4mZllc/eUmZllc9IwM7NsThpmZpbNScPMzLI5aZiZWTYnDTMzy+akYWZm2Zw0zMws2/8HyEmhkOb3rl0AAAAASUVORK5CYII=\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "yearly_incidence.plot(style='*')" + ] + }, + { + "cell_type": "code", + "execution_count": 52, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "2020 272809\n", + "2002 540966\n", + "2018 545389\n", + "1996 588428\n", + "2017 595842\n", + "2019 596798\n", + "2015 598951\n", + "2000 599943\n", + "2012 616931\n", + "1995 637629\n", + "1998 641356\n", + "2001 644041\n", + "1994 645857\n", + "2011 652346\n", + "2004 664779\n", + "2005 671304\n", + "2014 677491\n", + "2006 682039\n", + "2007 687771\n", + "1993 690753\n", + "1997 693627\n", + "2013 739418\n", + "2016 748281\n", + "2008 774750\n", + "1999 787325\n", + "2003 792368\n", + "2009 795973\n", + "1992 807510\n", + "2010 825410\n", + "dtype: int64" + ] + }, + "execution_count": 52, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "yearly_incidence.sort_values()" + ] + }, + { + "cell_type": "code", + "execution_count": 53, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 53, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAW4AAAEICAYAAAB/Dx7IAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAE6JJREFUeJzt3X+QXXV5x/H3Q/hhYCFaA2sN1qVoESUCzapVp7gL1kFDdYq2imjV6sQZFLFNx0kdEW11TNU4WoudpkrR+mNVlFGIvxh1RRSRjVADIooQlVBDqRJcjGDg6R/nrCwhu/dkd++5+72+XzM7uXvvd+95nntuPvfc7z3n3MhMJEnl2KfXBUiS9o7BLUmFMbglqTAGtyQVxuCWpMIY3JJUGINbkgpjcEtSYQxuSSrMvt240+XLl+fQ0FA37npWd955JwcddFDry22DvZXJ3srUi942b958W2Ye2mRsV4J7aGiIiYmJbtz1rMbHxxkZGWl9uW2wtzLZW5l60VtE/LjpWKdKJKkwBrckFcbglqTCGNySVBiDW5IK0zG4I+KoiLh62s8dEfHaNoqTJD1Qx90BM/N64DiAiFgCbAMu7HJdkqQZ7O1UyUnAjzKz8f6GkqSFFXvznZMRcR7wncz81z3ctgZYAzA4OLhqbGxswYpsanJykoGBgdaX2wZ7K9Ni6G3Lth1dud/BpbB958y3r1yxrCvLbUMv1tvo6OjmzBxuMrZxcEfE/sAtwOMyc/tsY4eHh9MjJxeWvZVpMfQ2tG5TV+537cpdbNgy82zr1vWru7LcNvToyMnGwb03UyXPpNranjW0JUndtTfBfRrwsW4VIklqplFwR8SBwJ8Bn+5uOZKkThqdHTAzfwU8tMu1SJIa8MhJSSqMwS1JhTG4JakwBrckFcbglqTCGNySVBiDW5IKY3BLUmEMbkkqjMEtSYUxuCWpMAa3JBXG4JakwhjcklQYg1uSCmNwS1JhDG5JKozBLUmFMbglqTAGtyQVxuCWpMI0Cu6IeHBEXBAR34+I6yLiyd0uTJK0Z/s2HPce4AuZ+byI2B84sIs1SZJm0TG4I+IQ4ATgpQCZeTdwd3fLkiTNJDJz9gERxwEbge8BxwKbgbMy887dxq0B1gAMDg6uGhsb60rBs5mcnGRgYKD15bbB3so01duWbTt6XcqCG1wK23fOfPvKFcvaK2Y38328O/U2k/n0PDo6ujkzh5uMbRLcw8C3gKdm5hUR8R7gjsw8e6a/GR4ezomJib2peUGMj48zMjLS+nLbYG9lmuptaN2mXpey4Nau3MWGLTO/ad+6fnWL1dzffB/vTr3NZD49R0Tj4G7y4eTNwM2ZeUX9+wXAH8+1OEnS/HQM7sz8GfDTiDiqvuokqmkTSVIPNH0vcCbwkXqPkhuBl3WvJEnSbBoFd2ZeDTSae5EkdZdHTkpSYQxuSSqMwS1JhTG4JakwBrckFcbglqTCGNySVBiDW5IKY3BLUmEMbkkqjMEtSYUxuCWpMAa3JBXG4JakwhjcklQYg1uSCmNwS1JhDG5JKozBLUmFMbglqTCNviw4IrYCvwTuAXZlpl8cLEk90ii4a6OZeVvXKpEkNeJUiSQVJjKz86CIm4BfAAn8e2Zu3MOYNcAagMHBwVVjY2MLXGpnk5OTDAwMtL7cNtjbwtiybUcry5kyuBS272x1ka3p1NvKFcvaK2Y3813Pc11v8+l5dHR0c9Np6KbB/fDMvCUiDgMuAc7MzEtnGj88PJwTExONC14o4+PjjIyMtL7cNtjbwhhat6mV5UxZu3IXG7bszYxkOTr1tnX96harub/5rue5rrf59BwRjYO70VRJZt5S/3srcCHwxDlXJ0mal47BHREHRcTBU5eBZwDXdLswSdKeNXkvMAhcGBFT4z+amV/oalWSpBl1DO7MvBE4toVaJEkNuDugJBXG4JakwhjcklQYg1uSCmNwS1JhDG5JKozBLUmFMbglqTAGtyQVxuCWpMIY3JJUGINbkgpjcEtSYQxuSSqMwS1JhTG4JakwBrckFcbglqTCGNySVBiDW5IKY3BLUmEaB3dELImIqyLi4m4WJEma3d5scZ8FXNetQiRJzTQK7og4HFgNvL+75UiSOonM7Dwo4gLgbcDBwN9n5il7GLMGWAMwODi4amxsbIFL7WxycpKBgYHWl9sGe1sYW7btaGU5UwaXwvadrS6yNZ16W7liWXvF7Ga+63mu620+PY+Ojm7OzOEmY/ftNCAiTgFuzczNETEy07jM3AhsBBgeHs6RkRmHds34+Di9WG4b7G1hvHTdplaWM2Xtyl1s2NLxv1mROvW29fSR9orZzXzX81zXW1s9N5kqeSrw7IjYCowBJ0bEh7talSRpRh2DOzP/ITMPz8wh4AXAVzLzRV2vTJK0R+7HLUmF2atJnMwcB8a7UokkqRG3uCWpMAa3JBXG4JakwhjcklQYg1uSCmNwS1JhDG5JKozBLUmFMbglqTAGtyQVxuCWpMIY3JJUGINbkgpjcEtSYQxuSSqMwS1JhTG4JakwBrckFcbglqTCGNySVJiOwR0RD4qIb0fEf0fEtRHx5jYKkyTtWZNveb8LODEzJyNiP+CyiPh8Zn6ry7VJkvagY3BnZgKT9a/71T/ZzaIkSTNrNMcdEUsi4mrgVuCSzLyiu2VJkmYS1QZ1w8ERDwYuBM7MzGt2u20NsAZgcHBw1djY2ELW2cjk5CQDAwOtL7cN9rYwtmzb0cpypgwuhe07W11kazr1tnLFsvaK2c181/Nc19t8eh4dHd2cmcNNxu5VcANExDnAnZn5zpnGDA8P58TExF7d70IYHx9nZGSk9eW2wd4WxtC6Ta0sZ8ralbvYsKXJR0nl6dTb1vWrW6zm/ua7nue63ubTc0Q0Du4me5UcWm9pExFLgacD359zdZKkeWnykvL7wAcjYglV0H8iMy/ublmSpJk02avku8DxLdQiSWrAIyclqTAGtyQVxuCWpMIY3JJUGINbkgpjcEtSYQxuSSqMwS1JhTG4JakwBrckFcbglqTCGNySVBiDW5IKY3BLUmEMbkkqjMEtSYUxuCWpMAa3JBXG4JakwhjcklQYg1uSCtMxuCPiERHx1Yi4LiKujYiz2ihMkrRn+zYYswtYm5nfiYiDgc0RcUlmfq/LtUmS9qDjFndm/k9mfqe+/EvgOmBFtwuTJO1ZZGbzwRFDwKXAMZl5x263rQHWAAwODq4aGxtbuCobmpycZGBgoPXltqHfetuybcdvLw8uhe07e1hMF9lbmeba28oVy+a8zNHR0c2ZOdxkbOPgjogB4GvAWzPz07ONHR4ezomJiUb3u5DGx8cZGRlpfblt6LfehtZt+u3ltSt3sWFLk1m78thbmeba29b1q+e8zIhoHNyN9iqJiP2ATwEf6RTakqTuarJXSQAfAK7LzHd1vyRJ0myabHE/FXgxcGJEXF3/PKvLdUmSZtBxEiczLwOihVokSQ145KQkFcbglqTCGNySVBiDW5IKY3BLUmEMbkkqjMEtSYUxuCWpMAa3JBXG4JakwhjcklQYg1uSCmNwS1JhDG5JKozBLUmFMbglqTAGtyQVxuCWpMIY3JJUGINbkgpjcEtSYToGd0ScFxG3RsQ1bRQkSZpdky3u84GTu1yHJKmhjsGdmZcCP2+hFklSA5GZnQdFDAEXZ+Yxs4xZA6wBGBwcXDU2NjangrZs2zGnvwMYXArbd87tb1euWDbn5bZhcnKSgYGBXpexYKav5/mst8XO3so0197mkyOjo6ObM3O4ydgFC+7phoeHc2JiosnQBxhat2lOfwewduUuNmzZd05/u3X96jkvtw3j4+OMjIz0uowFM309z2e9LXb2Vqa59jafHImIxsHtXiWSVBiDW5IK02R3wI8BlwNHRcTNEfHy7pclSZpJx0mczDytjUIkSc04VSJJhTG4JakwBrckFcbglqTCGNySVBiDW5IKY3BLUmEMbkkqjMEtSYUxuCWpMAa3JBXG4JakwhjcklQYg1uSCmNwS1JhDG5JKozBLUmFMbglqTAGtyQVxuCWpMI0Cu6IODkiro+IGyJiXbeLkiTNrGNwR8QS4FzgmcBjgdMi4rHdLkyStGdNtrifCNyQmTdm5t3AGPCc7pYlSZpJZObsAyKeB5ycma+of38x8KTMfPVu49YAa+pfjwKuX/hyO1oO3NaD5bbB3spkb2XqRW+PzMxDmwzct8GY2MN1D0j7zNwIbGyy0G6JiInMHO5lDd1ib2WytzIt9t6aTJXcDDxi2u+HA7d0pxxJUidNgvtK4NERcURE7A+8APhsd8uSJM2k41RJZu6KiFcDXwSWAOdl5rVdr2xuejpV02X2ViZ7K9Oi7q3jh5OSpMXFIyclqTAGtyQVxuCWpML8zgd3RJwYEUf0uo6F1q99Tenn/uytTG329jv74WR9vpUx4HbgXuCczPxab6uav37ta0o/92dvZepFb78zW9wRcXhEHDLtqucDn8rME6ge9NMi4sm9qW7u9qaviNjTUbCLWr+uN7A3e5u7vg/uiDg6Ij4HXAb8Y0RMnSDr18CB9eVPUJ2X4E9KCbe59JUFvb3q1/UG9oa9zVtfBndEHDTt1+OAmzNzCPgK8M76+p8Dd0XEwZn5c+AHwCAw1GKpe2WWvr7KA/saqPv6IYu8ryn9ut7A3rC3BdU3wR0RD4mI8yPiSmB9RBxav+I9HvhGvcX5WeD2iFhN9eAeDKys7+KHVGcEu7sX9c9klr5Wcl9fn+GBfT2+vosfsAj7mtKv6w3szd6611vfBDdwArALeBbVGQ1fDxxC1ePDpk0TfBB4IfBt4JdUXxBBZl4OnAjc0W7ZHc3U1xLK7mtKv643sDewt+7IzKJ+qALrlcDXqM7/vby+/hPAa+rLRwDr69ufQDUntaS+bQD43/p+VgDXAa8G/hN4H3CgfdmfvdnbYu6txC3uU4BnA28Gngy8vb7+EuAp9eWfAl8HnpmZV1K9ao4CZOYkcAXwhMzcBryYaj7qZ8AbMvNXLfWxu37ta0o/92dv9tZqb02+SKF1U3tARMQTqN6KfB3YlJl3AX8E3JiZX4mIm4B3RMQzgM3AX0TE8sy8LSJ+CNwZEX8AvBd4UUQcRnU+8f+jentDZk4AE/Zlfw17tLeCeuvX5+Si2+KOiCX1A/00qrccvwaeDrytHnIv8IOIWJqZN1G9PXk81RzTLVT7VALcQ/UWZh/gU1SnaTwdWAVszMx7W2oJgLrejIgR4Dz6pK8pEXFIv/YXEQ+t/30KcD791dthEfHQiBimmrPtp97278csARbHHDdwEPCK+kFZS/VBwGuBV9W3PwT4LnA81YO5HhiqbzuF6oFcXl/eAiyj+nDgc8D+05azTw96OxjYRHUec4C/7Ye+6uUeCLwE+DLVAQh90R/3HVH8l1RzmuNU36PaL8/Jg4CXUk0J7ABW91Fv+wFnABcB/wY8CjirH3qb/tPzLe56P8kvUz04/wE8A3gu1QcB9wJk5i+AzwCvoXqyHQYcXd/FpcDTgLsz82LgA8AFwLlUWxC/mVpW9mZrdClwAHBkRCwHjqR69S66r4jYD7gWeB7wjsx8bn3T8dPqKrK/zMyIWAb8FfDuzBzJzOuptrCmxhTZW0Q8imq64CTgDcA24CdUH8IV/7wEXkWVJe+m+m7cU+vf76lrKrm3+/TyVWPaq9eDp11+HdWDejrw7WnXPxy4pb78KqpDSx9CFYwXAX8wbezyXvc0rZaXAO8AzgZeDvw5cGXpfdX1fBo4fbfrng9cUXp/VFtt/1RfntoCP7X03qjC+YBpv59H9eL7nNJ7q2u5CPjr+vLLgTPr52TxWTL9p+db3ACZeXtEHBIR51NNlSynegCPjIjBeswtwLUR8aTMPBe4Afg41ZcZX5aZP5l2f7e13cPu6h32oZoX+xHV1umJmXkRcESpfe3mPOCciNgQEeMR8UbgW1TfUXoYFN3fbcCfRsTpwOaI+CDVVttj6ndORfaWmfdk9cHc1Nx9UJ0c6SKq9fawelxxvdU+D7wkIj4JvAk4FriGqrdDoeje7tPrV47dXi3PoHrLtpFqzvubwBupnly/B7yf+tWQai7rGOBBva67Q0+fpNrKWUa1v+gbqJ5IZ9e3F9nXtP6+SHWAwiOonvxnAZf3wXp7dN3Hv1Btjb0QeBfVfruvo3pBLrK33fq8Cji1vvxh4J9Lf17W2XEe1Qbgm4BzgOvr52lfrLdFeVrXiDiOasf4y6nmn46h2vXm6sx8WS9r2xsRMUA1TXIAVR+Pofqw6/XAy6jCobi+pouIA7PeXzUijqU6UOEbVIf/FrneACLiAOBLwLWZeUZUZ4P7G6otuB8Dw1QHXRTXG0BE7JOZ90bEe4GfZubbI+Joqufl0RS63gAi4l1UtX8oIh5JNWVyL9UG1PEUvN6mLMr9uKn2jTySaif2D0XEi6j+A13V47r21i6qPS9+Q7WlfQ/w+szcAvxdwX39Vt7/IIPbqfaNPTszP1pyf5l5V0SsB95SX7UP1UEa6zPzqpJ7g+rDtYiYOqPdDfV11wGvq6eHvldibxGxhGrK4wTgQ8CdVDs6nJGZN5W+3qYsmi3u+lP8k6jekj6Warrk3Mz8zax/WJB6B/5TgbHM/Fmv61kI9ZbpyVRHjT2Oahes92Xmrp4WtkAi4q3AH1JtqX2J6iT5v+htVQsnIq4H3piZH586WKXXNc1XRBxJlR93U627C4G3ZHWkY19YTMG9L9X0yF3AhzPz1z0uacHUWwH39sN/ij2JiFdSvRX9r35ab1Mi4ijgx/3U27QjCo+j+uB8Vz89P+uNpEcD38zMnb2uZ6EtmuCWJDWzKHYHlCQ1Z3BLUmEMbkkqjMEtSYUxuCWpMAa3JBXG4Jakwvw/Snk26Fyal0MAAAAASUVORK5CYII=\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "yearly_incidence.hist(xrot=20)" + ] + }, { "cell_type": "code", "execution_count": null,