analyse varicelle

module3/exo2/exercice.ipynb → module3/exo2/analyse_varicelle.ipynb

@@ -1201,20 +1201,151 @@
     "Comme l'incidence de syndrome grippal est très faible en septembre, cette\n",
     "modification ne risque pas de fausser nos conclusions.\n",
     "\n",
-    "Encore un petit détail: les données commencent en septembre 1984, ce qui\n",
-    "rend la première année incomplète. Nous commençons donc l'analyse en 1985."
+    "Encore un petit détail: les données commencent en semaine 49 de 1990, ce qui\n",
+    "rend la première année incomplète. Nous commençons donc l'analyse en 1991."
    ]
   },
   {
    "cell_type": "code",
-   "execution_count": 16,
+   "execution_count": 21,
    "metadata": {},
    "outputs": [],
    "source": [
     "first_september_week = [pd.Period(pd.Timestamp(y, 9, 1), 'W')\n",
-    "                     for y in range(1985,\n",
+    "                     for y in range(1991,\n",
     "                                    sorted_data.index[-1].year)]"
    ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 22,
+   "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": 23,
+   "metadata": {},
+   "outputs": [
+    {
+     "data": {
+      "text/plain": [
+       "<matplotlib.axes._subplots.AxesSubplot at 0x7f2bbdfff908>"
+      ]
+     },
+     "execution_count": 23,
+     "metadata": {},
+     "output_type": "execute_result"
+    },
+    {
+     "data": {
+      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY0AAAD8CAYAAACLrvgBAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAG61JREFUeJzt3X+MXeV95/H3x4wzA9QmYxgT/wDsKg6KISnEV7a7yXZL3NpmU2GjQjIbGkaLJTfApqRaCezACi24ElSrTWtFYbFCwUD44bpFeLt4yWA2qrd1bI8DKRjiziQQ4+Aww44DJpKnjP3dP+5z6zuT8cy58+vec+fzkq7uud85z3Of4+OZ731+nHMVEZiZmWUxrdoNMDOz/HDSMDOzzJw0zMwsMycNMzPLzEnDzMwyc9IwM7PMnDTMzCwzJw0zM8vMScPMzDJrqHYDxtsFF1wQCxYsqHYzzMxy5cCBA+9GRMtI+9Vd0liwYAEdHR3VboaZWa5I+lmW/Tw8ZWZmmTlpmJlZZk4aZmaWmZOGmZll5qRhZmaZOWmY1Znu90/wxQf30H38RLWbYnXIScOszmze1cn+N3vZ/EJntZtidajurtMwm6ouvWsnff2n/vX143sP8/jewzQ2TOPQpqur2DKrJ+5plHG33vJs9+1Xcc0Vc2maXvy1bpo+jTVXzGX3HVdVuWVWT5w0yrhbb3k2e2YTMxob6Os/RWPDNPr6TzGjsYHZM5rOWMYflKxSHp7C3XqrH+9+0McNyy7hy0sv5ol9h+kZIRmUf1DadO2nJqmVlmeKiGq3YVwVCoWo9N5T3e+fYNNzr/O9g7/gxIenaJo+jVWXfYw7v/DJYT+lmeXV4A9KJf6gNHVJOhARhZH28/AUo+vWm+WZ5z9stDw8lVTarTfLM39QstFy0kge/MrpXtmmtZdXsSVmk8MflGw0PKdhZmae0zAzs/HnpGFmZpllShqS/lTSQUmvSnpSUpOkWZLaJXWm5+ay/TdK6pJ0SNKqsvgSSa+kn22WpBRvlPR0iu+VtKCsTFt6j05JbeN36GZmVqkRk4akecCfAIWIuBw4C2gFNgC7ImIRsCu9RtLi9PPLgNXAtyWdlap7AFgPLEqP1Sm+DjgWER8Hvgncn+qaBdwNLAOWAneXJyczM5tcWYenGoCzJTUA5wBvA2uArennW4G1aXsN8FRE9EXEG0AXsFTSHGBmROyJ4uz7o4PKlOraDqxIvZBVQHtE9EbEMaCd04nGzMwm2YhJIyJ+Dvw34DBwFHgvIr4HXBgRR9M+R4HZqcg84K2yKo6k2Ly0PTg+oExE9APvAecPU5eZmVVBluGpZoo9gYXAXOBcSX80XJEhYjFMfLRlytu4XlKHpI6enp5hmmZmZmORZXjq94A3IqInIj4E/hb4N8A7aciJ9Nyd9j8CXFRWfj7F4awjaXtwfECZNAR2HtA7TF0DRMSWiChERKGlpSXDIZmZ2WhkSRqHgeWSzknzDCuA14EdQGk1UxvwbNreAbSmFVELKU5470tDWMclLU/13DioTKmu64AX07zH88BKSc2px7MyxczMxsS3hR+dEW8jEhF7JW0Hfgj0Ay8BW4DfALZJWkcxsVyf9j8oaRvwWtr/1og4maq7GXgEOBvYmR4ADwGPSeqi2MNoTXX1SroX2J/2uyciesd0xGZm+Lbwo+XbiNiU1/3+Cf7Tky/xrS9f6Rv2TQG+LfzQfBsRs4z8jY1Ti28LPza+y61NWf7GxqnJt4UfG/c0bMryJ86pq3Rb+Gdu+Sw3LLuEng/6qt2k3HBPw6Ysf+Kcuvz9OaPnpGFTmr+IyKwyXj1ldckroswq49VTNiHyckGUV0SZTQwPT1lFav2CKK+IMptYHp6yTPJyQVT3+yfY9NzrfO/gLzjx4Smapk9j1WUf484vfDL3w1QecrOJ5OEpG1d5WZ5azyuiPORmtcDDU5ZJnv4Y19uKKA+5WS1x0rDM8vLHuN7W4O++/aozDrmZTTYnDcus3v4Y50WeenlW/5w0zHIgL708q39ePWVmZl49ZWZm489Jw8zMMnPSMDOzzJw0zKooL/fyMitx0jCrIl/lbXnjJbdmVeCrvC2v3NMwq4K83MvLbDAnDbMq8FXellcjJg1Jl0p6uezxvqSvS5olqV1SZ3puLiuzUVKXpEOSVpXFl0h6Jf1ssySleKOkp1N8r6QFZWXa0nt0Smob38M3q57SVd7P3PJZblh2CT0f9FW7SWYjquiKcElnAT8HlgG3Ar0RcZ+kDUBzRNwhaTHwJLAUmAu8AHwiIk5K2gfcBvwAeA7YHBE7Jd0CfDoiviqpFbg2Ir4kaRbQARSAAA4ASyLi2Jna6CvCzfLB3w9SWybqivAVwE8i4mfAGmBrim8F1qbtNcBTEdEXEW8AXcBSSXOAmRGxJ4qZ6tFBZUp1bQdWpF7IKqA9InpTomgHVlfYZjOrQV45lk+Vrp5qpdiLALgwIo4CRMRRSbNTfB7FnkTJkRT7MG0PjpfKvJXq6pf0HnB+eXyIMmaWQ145lm+ZexqSPgJcA/z1SLsOEYth4qMtU9629ZI6JHX09PSM0DwzqyavHMu3SoanrgZ+GBHvpNfvpCEn0nN3ih8BLiorNx94O8XnDxEfUEZSA3Ae0DtMXQNExJaIKEREoaWlpYJDMrPJNpqVY75yvnZUkjT+A6eHpgB2AKXVTG3As2Xx1rQiaiGwCNiXhrKOS1qe5ituHFSmVNd1wItp3uN5YKWk5rQ6a2WKWQb+RbNaVenKMc9/1I5Mq6cknUNxbuE3I+K9FDsf2AZcDBwGro+I3vSzO4GbgH7g6xGxM8ULwCPA2cBO4GsREZKagMeAKyn2MFoj4qepzE3AN1JT/iwiHh6urV49ddpdz7zCd/cd5oalF7Pp2k9VuzlmFRs8/1Hi+Y/xl3X1lL+EqQ6N5hfNyx+tFnW/f+KM34/u/6fjy1/CNIWNZqLR3X+rRb5yvvb4hoV1qJJfNC9/tFrn70evLU4adSrrL9ru2686Y/ffrBY8+JXTIyab1l5exZYYOGnUray/aO7+m1klPKcxwfKw7NU3zjOzrLx6aoJ52auZ5UHW1VMenpognmA2s3rk4akJUumy1zwMY5mZOWlMkEonmH2dhJnlgYenJlCWZa8exjKzPPFE+CiM5y03fJsEywPfZqb++TYiE2g8h5J8nYTlgYdPrcTDUxWYqKEk3ybBapWHT20wD09VwENJllW9DOf4//zU4eGpCeChpKmtkmXR9TKc4//zNpiHpyrkoaSpqzwRnOnq/noczvH/eSvn4SmzEVTypVYezqmfobmpxsNTlhu1fjV8JVf3ezinfobmbGgenrIJk/UTZ5Zhn2qqNBFM1eGcehyas1/n4SmbMCPd4Xc032VeLX/8WActM5oGJILy7ywxD83lne9ya1WT9RNnnr410N8eNzIPzU0NntOwcZd1DsB/ZOqPv9Cr/rmnYeOukmQwVcf/61W1e2ReuTXxMvU0JH1U0nZJP5b0uqTfljRLUrukzvTcXLb/Rkldkg5JWlUWXyLplfSzzZKU4o2Snk7xvZIWlJVpS+/RKalt/A7dJlLWT5wPfqXAprWXs3juTDatvdzzBDYmXrk18TJNhEvaCuyOiO9I+ghwDvANoDci7pO0AWiOiDskLQaeBJYCc4EXgE9ExElJ+4DbgB8AzwGbI2KnpFuAT0fEVyW1AtdGxJckzQI6gAIQwAFgSUQcO1NbPRFuNvXkaVFFrRq36zQkzQR+B3gIICL+JSJ+CawBtqbdtgJr0/Ya4KmI6IuIN4AuYKmkOcDMiNgTxUz16KAypbq2AytSL2QV0B4RvSlRtAOrR2qzmU0tlX5Tpo1eluGp3wR6gIclvSTpO5LOBS6MiKMA6Xl22n8e8FZZ+SMpNi9tD44PKBMR/cB7wPnD1DWApPWSOiR19PT0ZDgkM6snXlQxebIkjQbgM8ADEXEl8CtgwzD7a4hYDBMfbZnTgYgtEVGIiEJLS8swTTOzeuWVW5Mjy+qpI8CRiNibXm+nmDTekTQnIo6moafusv0vKis/H3g7xecPES8vc0RSA3Ae0JvivzuozPczHZmZTSnVXrlViTyv8hqxpxERvwDeknRpCq0AXgN2AKXVTG3As2l7B9CaVkQtBBYB+9IQ1nFJy9N8xY2DypTqug54Mc17PA+slNScVmetTDEzs9zK8yqvrNdpfA34blo59VPgP1JMONskrQMOA9cDRMRBSdsoJpZ+4NaIOJnquRl4BDgb2JkeUJxkf0xSF8UeRmuqq1fSvcD+tN89EdE7ymM1M6uqerg/l+89ZWY2SWr5/ly+NbqZWY2ph1Vevo2Imdkkyvutczw8ZWZmHp6qZ7X+TXdmVr+cNHIoz8v1zCzfPKeRI/WwXM/M8s09jRzxTdnMrNqcNHKkHpbrmeWN5xAHctLIGd+UzWxyeQ5xIC+5NTMbwlT7YicvuTUzG4PRzCFOhaEsJw0zsyGMZg5xKgxlecmtmdkZZL3lx1RaDu85DTOzMarlu9dm5TkNM7NJMpWWw3t4ysxsHOT97rVZeXjKzMw8PGVmZuPPScNyZSqsgzerZU4alitTYR28WS3zRLjlwlRaB29Wy9zTsFzwbeHNaoOThuXCVFoHb1bLMiUNSW9KekXSy5I6UmyWpHZJnem5uWz/jZK6JB2StKosviTV0yVpsySleKOkp1N8r6QFZWXa0nt0SmobrwO3/PFt4c2qL9N1GpLeBAoR8W5Z7M+B3oi4T9IGoDki7pC0GHgSWArMBV4APhERJyXtA24DfgA8B2yOiJ2SbgE+HRFfldQKXBsRX5I0C+gACkAAB4AlEXHsTG31dRpmZpWbjOs01gBb0/ZWYG1Z/KmI6IuIN4AuYKmkOcDMiNgTxUz16KAypbq2AytSL2QV0B4RvSlRtAOrx9BmMzMbg6xJI4DvSTogaX2KXRgRRwHS8+wUnwe8VVb2SIrNS9uD4wPKREQ/8B5w/jB1DSBpvaQOSR09PT0ZD8nMzCqVdcntZyPibUmzgXZJPx5mXw0Ri2Hioy1zOhCxBdgCxeGpYdpmZmZjkKmnERFvp+du4BmK8xXvpCEn0nN32v0IcFFZ8fnA2yk+f4j4gDKSGoDzgN5h6jIzsyoYMWlIOlfSjNI2sBJ4FdgBlFYztQHPpu0dQGtaEbUQWATsS0NYxyUtT/MVNw4qU6rrOuDFNO/xPLBSUnNanbUyxczMrAqyDE9dCDyTVsc2AE9ExP+WtB/YJmkdcBi4HiAiDkraBrwG9AO3RsTJVNfNwCPA2cDO9AB4CHhMUhfFHkZrqqtX0r3A/rTfPRHRO4bjNTOzMfCt0c3MzLdGNzOz8eekYWZmmTlpmJlZZk4aZmaWmZOGmZll5qRhZmaZOWmYmVlmThpmZpaZk4aZmWXmpGFmZpk5aZiZWWZOGmZmlpmThpmZZeakYWZmmTlpmJlZZk4aZmaWmZOGmZll5qRhZmaZOWmYmVlmThpmZpaZk4aZmWXmpGFmZpllThqSzpL0kqS/S69nSWqX1Jmem8v23SipS9IhSavK4kskvZJ+tlmSUrxR0tMpvlfSgrIybek9OiW1jcdBm5nZ6FTS07gNeL3s9QZgV0QsAnal10haDLQClwGrgW9LOiuVeQBYDyxKj9Upvg44FhEfB74J3J/qmgXcDSwDlgJ3lycnMzObXJmShqT5wBeA75SF1wBb0/ZWYG1Z/KmI6IuIN4AuYKmkOcDMiNgTEQE8OqhMqa7twIrUC1kFtEdEb0QcA9o5nWjMzGySZe1p/AVwO3CqLHZhRBwFSM+zU3we8FbZfkdSbF7aHhwfUCYi+oH3gPOHqcvMzKpgxKQh6Q+A7og4kLFODRGLYeKjLVPexvWSOiR19PT0ZGymmZlVKktP47PANZLeBJ4CPi/pceCdNOREeu5O+x8BLiorPx94O8XnDxEfUEZSA3Ae0DtMXQNExJaIKEREoaWlJcMhmZnZaIyYNCJiY0TMj4gFFCe4X4yIPwJ2AKXVTG3As2l7B9CaVkQtpDjhvS8NYR2XtDzNV9w4qEypruvSewTwPLBSUnOaAF+ZYmZmVgUNYyh7H7BN0jrgMHA9QEQclLQNeA3oB26NiJOpzM3AI8DZwM70AHgIeExSF8UeRmuqq1fSvcD+tN89EdE7hjabmdkYqPiBvn4UCoXo6OiodjPMzHJF0oGIKIy0n68INzOrA93vn+CLD+6h+/iJCX0fJw0zszqweVcn+9/sZfMLnRP6PmOZ0zAzsyq79K6d9PWfvoTu8b2HeXzvYRobpnFo09Xj/n7uaZiZ5dju26/imivm0jS9+Oe8afo01lwxl913XDUh7+ekYWaWY7NnNjGjsYG+/lM0Nkyjr/8UMxobmD2jaULez8NTZmY59+4Hfdyw7BK+vPRinth3mJ4JnAz3klszM/OSWzMzG39OGmZmlpmThpmZZeakYWZmmTlpmJlZZk4aZmaWmZOGmZll5qRhZmaZOWmYmVlmThpmZpaZk4aZmWXmpGFmZpk5aZiZWWZOGmZmlpmThpmZZeakYWZmmY2YNCQ1Sdon6UeSDkr6ryk+S1K7pM703FxWZqOkLkmHJK0qiy+R9Er62WZJSvFGSU+n+F5JC8rKtKX36JTUNp4Hb2ZmlcnS0+gDPh8RvwVcAayWtBzYAOyKiEXArvQaSYuBVuAyYDXwbUlnpboeANYDi9JjdYqvA45FxMeBbwL3p7pmAXcDy4ClwN3lycnMzCbXiEkjij5IL6enRwBrgK0pvhVYm7bXAE9FRF9EvAF0AUslzQFmRsSeKH7H7KODypTq2g6sSL2QVUB7RPRGxDGgndOJxszMJlmmOQ1JZ0l6Geim+Ed8L3BhRBwFSM+z0+7zgLfKih9JsXlpe3B8QJmI6AfeA84fpi4zM6uCTEkjIk5GxBXAfIq9hsuH2V1DVTFMfLRlTr+htF5Sh6SOnp6eYZpmZmZjUdHqqYj4JfB9ikNE76QhJ9Jzd9rtCHBRWbH5wNspPn+I+IAykhqA84DeYeoa3K4tEVGIiEJLS0slh2RmZhXIsnqqRdJH0/bZwO8BPwZ2AKXVTG3As2l7B9CaVkQtpDjhvS8NYR2XtDzNV9w4qEypruuAF9O8x/PASknNaQJ8ZYqZmVkVNGTYZw6wNa2AmgZsi4i/k7QH2CZpHXAYuB4gIg5K2ga8BvQDt0bEyVTXzcAjwNnAzvQAeAh4TFIXxR5Ga6qrV9K9wP603z0R0TuWAzYzs9FT8QN9/SgUCtHR0VHtZpiZ5YqkAxFRGGk/XxFuZmaZOWmYmVlmThpmZpaZk4aZmWXmpGFmZpk5aZiZWWZOGmZmlpmThpmZZeakYWZmmTlpmJlZZk4aZmaWmZOGmZll5qRhZmaZOWmYmVlmThpmZpaZk4aZmWXmpGFmZpk5aZiZ1bDu90/wxQf30H38RLWbAjhpmJnVtM27Otn/Zi+bX+isdlMAaKh2A8zM7NddetdO+vpP/evrx/ce5vG9h2lsmMahTVdXrV3uaZiZ1aDdt1/FNVfMpWl68c900/RprLliLrvvuKqq7XLSMDOrQbNnNjGjsYG+/lM0Nkyjr/8UMxobmD2jqart8vCUmVmNeveDPm5YdglfXnoxT+w7TE8NTIYrIobfQboIeBT4GHAK2BIRfylpFvA0sAB4E/hiRBxLZTYC64CTwJ9ExPMpvgR4BDgbeA64LSJCUmN6jyXA/wO+FBFvpjJtwF2pOZsiYutw7S0UCtHR0ZH9X8DMzJB0ICIKI+2XZXiqH/jPEfFJYDlwq6TFwAZgV0QsAnal16SftQKXAauBb0s6K9X1ALAeWJQeq1N8HXAsIj4OfBO4P9U1C7gbWAYsBe6W1JyhzWZmNgFGTBoRcTQifpi2jwOvA/OANUDpU/9WYG3aXgM8FRF9EfEG0AUslTQHmBkRe6LYvXl0UJlSXduBFZIErALaI6I39WLaOZ1ozMxsklU0ES5pAXAlsBe4MCKOQjGxALPTbvOAt8qKHUmxeWl7cHxAmYjoB94Dzh+mLjMzq4LMSUPSbwB/A3w9It4fbtchYjFMfLRlytu2XlKHpI6enp5hmmZmZmORKWlImk4xYXw3Iv42hd9JQ06k5+4UPwJcVFZ8PvB2is8fIj6gjKQG4Dygd5i6BoiILRFRiIhCS0tLlkMyM7NRGDFppLmFh4DXI+K/l/1oB9CWttuAZ8virZIaJS2kOOG9Lw1hHZe0PNV546AypbquA15M8x7PAyslNacJ8JUpZmZmVZBlye3ngN3AKxSX3AJ8g+K8xjbgYuAwcH1E9KYydwI3UVx59fWI2JniBU4vud0JfC0tuW0CHqM4X9ILtEbET1OZm9L7AfxZRDw8Qnt7gJ9lPP5acwHwbrUbMc7q7Zjq7Xig/o6p3o4HJueYLomIEYdqRkwaNnkkdWRZJ50n9XZM9XY8UH/HVG/HA7V1TL6NiJmZZeakYWZmmTlp1JYt1W7ABKi3Y6q344H6O6Z6Ox6ooWPynIaZmWXmnoaZmWXmpDHBJP2VpG5Jr5bFfkvSHkmvSPqfkmam+EckPZziP5L0u2Vlvi/pkKSX02P2EG834SRdJOn/SHpd0kFJt6X4LEntkjrTc3NZmY2SulL7V5XFl6Rj7ZK0OV2/k+fjyeU5knR+2v8DSd8aVFfuztEIx5PXc/T7kg6kc3FA0ufL6prccxQRfkzgA/gd4DPAq2Wx/cC/S9s3Afem7VuBh9P2bOAAMC29/j5QqIHjmQN8Jm3PAP4ZWAz8ObAhxTcA96ftxcCPgEZgIfAT4Kz0s33Ab1O8XcxO4OqcH09ez9G5wOeArwLfGlRXHs/RcMeT13N0JTA3bV8O/Lxa58g9jQkWEX9P8YLFcpcCf5+224E/TNuLKd5mnojoBn4J1MTa7JKYnLseT5rxOp7JbfXwKj2miPhVRPxfYMA3/OT1HJ3peGrJKI7ppYgo3ULpINCk4l03Jv0cOWlUx6vANWn7ek7fX+tHwBpJDSregmUJA++99XDqUv+XagwTDKaJu+txVYzxeEryeI7OJK/naCR5P0d/CLwUEX1U4Rw5aVTHTRS/zOoAxa7pv6T4X1E86R3AXwD/SPFWLAA3RMSngH+bHl+Z1BYPoom96/GkG4fjgfyeozNWMUQsD+doOLk+R5Iuo/gldX9cCg2x24SeIyeNKoiIH0fEyohYAjxJcVyciOiPiD+NiCsiYg3wUaAz/ezn6fk48ARVHBLRxN/1eFKN0/Hk+RydSV7P0Rnl+RxJmg88A9wYET9J4Uk/R04aVVBasSFpGsXvP/8f6fU5ks5N278P9EfEa2m46oIUnw78AcUhrmq0fTLuejxpxut4cn6OhpTjc3SmenJ7jiR9FPhfwMaI+IfSzlU5RxM5y+5HQLEncRT4kOKngnXAbRRXS/wzcB+nL7JcAByiOCn2AsW7TkJxNcgB4J8oToL9JWnFThWO53MUu7//BLycHv+e4jct7qLYM9oFzCorcyfF3tQhylZ2UJzkfzX97Fulf4c8Hk8dnKM3KS7Y+CD9P12c83P0a8eT53NE8cPlr8r2fRmYXY1z5CvCzcwsMw9PmZlZZk4aZmaWmZOGmZll5qRhZmaZOWmYmVlmThpmZpaZk4aZmWXmpGFmZpn9f/2Yi6a8X2tzAAAAAElFTkSuQmCC\n",
+      "text/plain": [
+       "<Figure size 432x288 with 1 Axes>"
+      ]
+     },
+     "metadata": {
+      "needs_background": "light"
+     },
+     "output_type": "display_data"
+    }
+   ],
+   "source": [
+    "yearly_incidence.plot(style='*')"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 24,
+   "metadata": {},
+   "outputs": [
+    {
+     "data": {
+      "text/plain": [
+       "2020    221186\n",
+       "2021    376290\n",
+       "2002    516689\n",
+       "2018    542312\n",
+       "2017    551041\n",
+       "1996    564901\n",
+       "2019    584066\n",
+       "2015    604382\n",
+       "2000    617597\n",
+       "2001    619041\n",
+       "2012    624573\n",
+       "2005    628464\n",
+       "2006    632833\n",
+       "2011    642368\n",
+       "1993    643387\n",
+       "1995    652478\n",
+       "1994    661409\n",
+       "1998    677775\n",
+       "1997    683434\n",
+       "2014    685769\n",
+       "2013    698332\n",
+       "2007    717352\n",
+       "2008    749478\n",
+       "1999    756456\n",
+       "2003    758363\n",
+       "2004    777388\n",
+       "2016    782114\n",
+       "2010    829911\n",
+       "1992    832939\n",
+       "2009    842373\n",
+       "dtype: int64"
+      ]
+     },
+     "execution_count": 24,
+     "metadata": {},
+     "output_type": "execute_result"
+    }
+   ],
+   "source": [
+    "yearly_incidence.sort_values()"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 25,
+   "metadata": {},
+   "outputs": [
+    {
+     "data": {
+      "text/plain": [
+       "<matplotlib.axes._subplots.AxesSubplot at 0x7f2bbde21e48>"
+      ]
+     },
+     "execution_count": 25,
+     "metadata": {},
+     "output_type": "execute_result"
+    },
+    {
+     "data": {
+      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAW8AAAEICAYAAACQzXX2AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAEaZJREFUeJzt3XuQZGV5x/Hvwy7owgCiC0MFMEPQGImrImPwUsFZsCwV1BJNvCCBlGZNeSNmUxZJqWgSlcRgaSxIakWUEuOoYEoB46Wio+IF2VXLFQlqBBGIRGJAl6Cw8uSP9wyZXXd2zlx6pp/N91M1Nae7T/d5nj7dv3P67dPdkZlIkmrZa6ULkCTNn+EtSQUZ3pJUkOEtSQUZ3pJUkOEtSQUZ3pJUkOEtSQUZ3pJU0OpB3fDatWtzbGxsh/PuvPNO9ttvv0EtclnYw3Cwh+FgD0tvy5Ytt2XmwXPNN7DwHhsbY/PmzTucNzU1xcTExKAWuSzsYTjYw3Cwh6UXET/oM5/DJpJUkOEtSQUZ3pJUkOEtSQUZ3pJUkOEtSQUZ3pJUkOEtSQUN7EM6kn7V2FlX3De9cd12zphxepBuOOekZVmOlo973pJUkOEtSQUZ3pJUkOEtSQUZ3pJUkOEtSQUZ3pJUkOEtSQUZ3pJUkOEtSQUZ3pJUkOEtSQUZ3pJUkOEtSQUZ3pJUkOEtSQUZ3pJUkOEtSQUZ3pJUUO/wjohXR8Q1EfGtiPhARNx/kIVJkmbXK7wj4jDgVcB4Zj4CWAU8f5CFSZJmN59hk9XAmohYDewL3DKYkiRJc4nM7DdjxJnAm4C7gE9l5qm7mGcDsAFgdHT02MnJyR0u37ZtGyMjI4uteUXZw3Co2sPWm++4b3p0Ddx61/Isd91hBw7kdquuh5mGrYf169dvyczxuebrFd4RcRBwKfA84Hbgw8AlmXnxbNcZHx/PzZs373De1NQUExMTcy5vmNnDcKjaw9hZV9w3vXHdds7dunpZlnvDOScN5HarroeZhq2HiOgV3n2HTZ4MXJ+ZP87Me4CPAE9YTIGSpIXrG943Ao+LiH0jIoATgWsHV5YkaXd6hXdmXgVcAnwN2Npdb9MA65Ik7UbvAbfMPBs4e4C1SJJ68hOWklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklRQ7/COiAdExCUR8W8RcW1EPH6QhUmSZrd6HvO+A/hEZj43IvYB9h1QTZKkOfQK74g4ADgeOAMgM+8G7h5cWZKk3YnMnHumiEcDm4BvA48CtgBnZuadO823AdgAMDo6euzk5OQOt7Nt2zZGRkaWpvIVYg/DoWoPW2++477p0TVw610rWMwS6NPDusMOXJ5iFmjYHkvr16/fkpnjc83XN7zHga8AT8zMqyLiHcBPM/N1s11nfHw8N2/evMN5U1NTTExMzLm8YWYPw6FqD2NnXXHf9MZ12zl363xGLodPnx5uOOekZapmYYbtsRQRvcK77xuWNwE3ZeZV3elLgMcstDhJ0uL0Cu/M/BHww4h4WHfWibQhFEnSCpjPa7ZXAu/vjjT5PvCHgylJkjSX3uGdmd8A5hyHkSQNnp+wlKSCDG9JKsjwlqSCDG9JKsjwlqSCDG9JKsjwlqSCDG9JKsjwlqSCDG9JKsjwlqSCDG9JKsjwlqSCDG9JKsjwlqSCDG9JKsjwlqSCDG9JKsjwlqSCDG9JKsjwlqSCDG9JKsjwlqSCDG9JKsjwlqSCDG9JKsjwlqSCDG9JKsjwlqSCDG9JKsjwlqSCDG9JKsjwlqSCDG9JKsjwlqSCDG9JKsjwlqSC5hXeEbEqIr4eEZcPqiBJ0tzmu+d9JnDtIAqRJPXXO7wj4nDgJOCCwZUjSeojMrPfjBGXAG8B9gf+LDNP3sU8G4ANAKOjo8dOTk7ucPm2bdsYGRlZbM0ryh6Gw2J62HrzHUtczcKMroFb71rpKhbHHnZt3WEHLvi669ev35KZ43PNt7rPjUXEycB/ZuaWiJiYbb7M3ARsAhgfH8+JiR1nnZqaYufzqrGH4bCYHs4464qlLWaBNq7bzrlbez0Fh5Y97NoNp04s6e3tSt9hkycCz4yIG4BJ4ISIuHhgVUmSdqtXeGfmn2fm4Zk5Bjwf+ExmvmiglUmSZuVx3pJU0LwHejJzCpha8kokSb255y1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBfUK74g4IiI+GxHXRsQ1EXHmoAuTJM1udc/5tgMbM/NrEbE/sCUiPp2Z3x5gbZKkWfTa887M/8jMr3XTPwOuBQ4bZGGSpNlFZs7vChFjwOeBR2TmT3e6bAOwAWB0dPTYycnJHa67bds2RkZGFlHuyrOH4bCYHrbefMcSV7Mwo2vg1rtWuorFsYddW3fYgQu+7vr167dk5vhc880rvCNiBPgc8KbM/Mju5h0fH8/NmzfvcN7U1BQTExO9lzeM7GE4LKaHsbOuWNpiFmjjuu2cu7XvyOVwsoddu+GckxZ83YjoFd69jzaJiL2BS4H3zxXckqTB6nu0SQDvBq7NzLcNtiRJ0lz67nk/ETgNOCEivtH9PX2AdUmSdqPXQE9mXgnEgGuRJPXkJywlqSDDW5IKMrwlqSDDW5IKMrwlqSDDW5IKMrwlqSDDW5IKMrwlqSDDW5IKMrwlqSDDW5IKMrwlqSDDW5IKMrwlqSDDW5IKMrwlqaCh/Nnnlfx178X86vNiLGfPG9dt54wZy6vY8849SP/fuOctSQUZ3pJUkOEtSQUZ3pJUkOEtSQUZ3pJUkOEtSQUZ3pJUkOEtSQUZ3pJUkOEtSQUZ3pJUkOEtSQUZ3pJUkOEtSQUZ3pJUkOEtSQUZ3pJUkOEtSQX1Du+IeGpEXBcR34uIswZZlCRp93qFd0SsAs4DngYcDbwgIo4eZGGSpNn13fP+HeB7mfn9zLwbmASeNbiyJEm7E5k590wRzwWempkv6U6fBhyXma/Yab4NwIbu5MOA63a6qbXAbYsteoXZw3Cwh+FgD0vv1zPz4LlmWt3zxmIX5/1K6mfmJmDTrDcSsTkzx3sucyjZw3Cwh+FgDyun77DJTcARM04fDtyy9OVIkvroG95XAw+NiCMjYh/g+cDHBleWJGl3eg2bZOb2iHgF8ElgFXBhZl6zgOXNOqRSiD0MB3sYDvawQnq9YSlJGi5+wlKSCjK8Jakgw1uSChr68I6IEyLiyJWuYzGq91C9/ml7Qh/2MByGoYehfcOy++6USeB24F7g7Mz83MpWNT/Ve6he/7Q9oQ97GA7D1MPQ7HlHxOERccCMs54HXJqZx9PurBdExONXprp+qvcwn/ojYlefuh0K1dcD2MOwGOYeVjy8I+LhEfFx4ErgLyNi+guvfg7s201/iPbdA48bxtCo3sNC6s8hfMlWfT2APQyLCj2sSHhHxH4zTj4auCkzx4DPAH/Xnf8T4BcRsX9m/gT4DjAKjC1jqbOq3sNu6v8sv1r/SFf/dxmS+qdVXw9gD9jDgixbeEfEQRHx3oi4GjgnIg7utlaPBL7Y7c19DLg9Ik6i3Sn7A+u6m/gu7du/7l6umndWvYfd1L9uRv0f3UX9j+xu4jsrWf+06usB7KG7CXtYhOXc8z4e2A48nfYthX8BHNDVcOiMl+EXAS8Evgr8jPYDEGTml4ETgJ8uY807q97DbPWvokb906qvB7AHe1iszFzSP1oQvBT4HO27vdd2538IeFU3fSRwTnf5Y2njSqu6y0aAH3e3cxhwLfAK4D3A+cC+S13zntZD9fr3pD7swR4G9TeIPe+TgWcCbwQeD/xtd/6ngSd00z8EvgA8LTOvpm3x1gNk5jbgKuCxmXkzcBptTOlHwGsz838GUPOe1kP1+vekPuzBHgai748x7GD6aIOIeCztpcQXgCsy8xfAbwLfz8zPRMT1wFsj4inAFuDZEbE2M2+LiO8Cd0bEg4F3Ai+KiENo3xX+X7SXJ2TmZmDzIvvc43qoXv+e1kfXiz34nF42897znnEHHQ9cSDt05snAW7pZ7gW+ExFrMvN62suLR9LGiW6hHScJ8EvaS5C9gEtpX8t4KnAssCkz711wV3P3sKrr4Um0lz2leujqyoiYoOg66Po4oHofEfGg7v8TgPcW7eGQiHhQRIzTxnYr9rBP5ef0gvQZW6Ed1/jHwD8BfwTsDfwJ8PLu8oOAbwLH0O6Ec4Cx7rKTaXfA2m56K3AgbZD/48A+M5az16DGh4D9gJfQVshG2psS1XrYH7iC9n3qAK+uVP+Mx9LpwL/SPuxQqg/+71PJv0cbE52i/V5rtcfSfsAZtGGDO4CTCvawN/Ay4DLgH4CHAGdW6mExf3PueUfEocDlwATwPtpg/im0caPtAJn538BHgVd1D4ZDgId3N/F54EnA3Zl5OfBu4BLgPNpW/p7pZeXg9lT3o4XFCcC7gKcAz6G9KXFvhR46a4D7AUdFxFrgKNpeQon6I2Jv4BrgucBbM/M53UXHzFj+UPeRmRkRBwK/D7w9Mycy8zranlmJHiLiIbQhhROB1wI3AzfS3rAr83gCXk57Tr+d9pu6p3Snf1moh4XrsXVbQ/ul+OnTZ9AG+08Hvjrj/F8DbummX0776OhB3fUvAx48Y961y72VAh4wY/o1tBV6arEeTgfeCrwOeDHwDODqKvV3y/0IcOpO5z0PuKpKH7S9vb/qpqf3xE+p0gMtoO834/SFtA3qs6r00C3zMuAPuukXA6/sHktlntOL+esz5v1z4KvdgevQXoYcl5kX0fYADwXIzFuAayLiuMw8D/ge8EHajxdfmZk3Tt9gZt7WY7lLKjNvj4gDIuK9tGGTtbSVd1REjA5zDzPu+72Af6ftvZ6QmZcBRw57/Tu5EDg7Is6NiKmIeD3wFdpvpB7S1TbsfdwG/G5EnApsiYiLaHt7v9W9IhrqHjLzl9nexJsesw/aFy1dRlsPJZ7TwL8Ap0fEh4E3AI8CvkXr4eCurmHvYeHmsZWb3sO4CDizm34f8Dfd9AOBC+i2ZLTxqEcA91/pLdROfbyM9vJwE20M/EvA62kP4KHuAfgwba/pQNrxqq+lPVhfV2wdfJL2YYgjaE+kM4EvF1oPD+3q/XvaXtwLgbfRjgN+DW0jO9Q97NTP14FTuumLKz2nu+fwhbSdsTcAZwPXdY+vUuthvn+9jzbJzIyIw4FDaVs8ujuLiLicNqa8KrstWWbek5nfysyf913GcsjM87O923w+bZzyH2kvoT7GEPcQESO0Pb53AZ+gvVN+HPAC4KCIuIwhrn8nz87MN2fmD4E30w7jeicF1kPnRtrHoVdnG1e9vDvvctobshV6ICKmn/9X0t7sA3hTu6jM4+loYCrbXvN7aDs3H6DOY2nh5rmVewbtC1r2pm3xnkq7s14IPGalt0Tz7OUI4FPAg7rTLwKOWem6dlPv/Wmvei6gvdE0AXxqxuVDXf9u+vp12ptJD6zUB+3j0Vu66QfQXkEcU6mHrtZ9aRvOU3Y6/9Rh76HLnj8FLuhOr6XtWB5ZbT0s5G9eP8YQEV8EfgO4gXZs5Bsz85u9b2CFdUcJnEjb2BxNGzo5LzPv2e0Vh1D3IYJTgMnM/NFK1zMfEXE/2ob/NOC3aYd5nZ+Z21e0sHmKiDfRng/H0HYEzs62J15KRFwHvD4zPzj9OY6VrqmviDiK9jy+m7Yu/hn462yfiNyj9Q7v7jCvs4HrgYuze8OjkohYTft+g1/Qeij30ikiVgH3VnqC7UpEvJR2mOb7Kq6HaRHxMOAHFXuY8YG7R9PeBN9e8XHV7cg8FPhSZt610vUsl6H9GTRJ0uxW/Jd0JEnzZ3hLUkGGtyQVZHhLUkGGtyQVZHhLUkGGtyQV9L8Giii4LM1vYwAAAABJRU5ErkJggg==\n",
+      "text/plain": [
+       "<Figure size 432x288 with 1 Axes>"
+      ]
+     },
+     "metadata": {
+      "needs_background": "light"
+     },
+     "output_type": "display_data"
+    }
+   ],
+   "source": [
+    "yearly_incidence.hist(xrot=20)"
+   ]
   }
  ],
  "metadata": {