Skip to content

M
mooc-rr
Commit 200ce142 authored by Eric ARMBRUSTER's avatar Eric ARMBRUSTER
Browse files
Options

Merge branch 'master' of... 

Merge branch 'master' of https://app-learninglab.inria.fr/moocrr/gitlab/9d26ff392e3a3b9ec48275a74df5354e/mooc-rr
parents f5ddcefa 51ea5dcc
Hide whitespace changes
Inline Side-by-side
Showing with 4813 additions and 41 deletions
module2/challenger.ipynb 0 → 100644
View file @ 200ce142
{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# Risk Analysis of the Space Shuttle: Pre-Challenger Prediction of Failure"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"In this document we reperform some of the analysis provided in \n",
"*Risk Analysis of the Space Shuttle: Pre-Challenger Prediction of Failure* by *Siddhartha R. Dalal, Edward B. Fowlkes, Bruce Hoadley* published in *Journal of the American Statistical Association*, Vol. 84, No. 408 (Dec., 1989), pp. 945-957 and available at http://www.jstor.org/stable/2290069. \n",
"\n",
"On the fourth page of this article, they indicate that the maximum likelihood estimates of the logistic regression using only temperature are: $\\hat{\\alpha}=5.085$ and $\\hat{\\beta}=-0.1156$ and their asymptotic standard errors are $s_{\\hat{\\alpha}}=3.052$ and $s_{\\hat{\\beta}}=0.047$. The Goodness of fit indicated for this model was $G^2=18.086$ with 21 degrees of freedom. Our goal is to reproduce the computation behind these values and the Figure 4 of this article, possibly in a nicer looking way."
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Technical information on the computer on which the analysis is run"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"We will be using the python3 language using the pandas, statsmodels, numpy, matplotlib and seaborn libraries."
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"3.6.4 |Anaconda, Inc.| (default, Jan 16 2018, 18:10:19) \n",
"[GCC 7.2.0]\n",
"uname_result(system='Linux', node='3a716011d2b6', release='4.4.0-116-generic', version='#140-Ubuntu SMP Mon Feb 12 21:23:04 UTC 2018', machine='x86_64', processor='x86_64')\n",
"IPython 6.4.0\n",
"IPython.core.release 6.4.0\n",
"PIL 5.2.0\n",
"PIL.Image 5.2.0\n",
"PIL._version 5.2.0\n",
"_csv 1.0\n",
"_ctypes 1.1.0\n",
"_curses b'2.2'\n",
"decimal 1.70\n",
"argparse 1.1\n",
"backcall 0.1.0\n",
"cffi 1.11.5\n",
"csv 1.0\n",
"ctypes 1.1.0\n",
"cycler 0.10.0\n",
"dateutil 2.7.3\n",
"decimal 1.70\n",
"decorator 4.3.0\n",
"distutils 3.6.4\n",
"ipaddress 1.0\n",
"ipykernel 4.8.2\n",
"ipykernel._version 4.8.2\n",
"ipython_genutils 0.2.0\n",
"ipython_genutils._version 0.2.0\n",
"ipywidgets 7.2.1\n",
"ipywidgets._version 7.2.1\n",
"jedi 0.12.1\n",
"json 2.0.9\n",
"jupyter_client 5.2.3\n",
"jupyter_client._version 5.2.3\n",
"jupyter_core 4.4.0\n",
"jupyter_core.version 4.4.0\n",
"kiwisolver 1.0.1\n",
"logging 0.5.1.2\n",
"matplotlib 2.2.2\n",
"matplotlib.backends.backend_agg 2.2.2\n",
"numpy 1.13.3\n",
"numpy.core 1.13.3\n",
"numpy.core.multiarray 3.1\n",
"numpy.core.umath b'0.4.0'\n",
"numpy.lib 1.13.3\n",
"numpy.linalg._umath_linalg b'0.1.5'\n",
"numpy.matlib 1.13.3\n",
"optparse 1.5.3\n",
"pandas 0.22.0\n",
"_libjson 1.33\n",
"parso 0.3.0\n",
"patsy 0.5.0\n",
"patsy.version 0.5.0\n",
"pexpect 4.6.0\n",
"pickleshare 0.7.4\n",
"platform 1.0.8\n",
"prompt_toolkit 1.0.15\n",
"ptyprocess 0.6.0\n",
"pygments 2.2.0\n",
"pyparsing 2.2.0\n",
"pytz 2018.5\n",
"re 2.2.1\n",
"scipy 1.1.0\n",
"scipy._lib.decorator 4.0.5\n",
"scipy._lib.six 1.2.0\n",
"scipy.fftpack._fftpack b'$Revision: $'\n",
"scipy.fftpack.convolve b'$Revision: $'\n",
"scipy.integrate._dop b'$Revision: $'\n",
"scipy.integrate._ode $Id$\n",
"scipy.integrate._odepack 1.9 \n",
"scipy.integrate._quadpack 1.13 \n",
"scipy.integrate.lsoda b'$Revision: $'\n",
"scipy.integrate.vode b'$Revision: $'\n",
"scipy.interpolate._fitpack 1.7 \n",
"scipy.interpolate.dfitpack b'$Revision: $'\n",
"scipy.linalg 0.4.9\n",
"scipy.linalg._fblas b'$Revision: $'\n",
"scipy.linalg._flapack b'$Revision: $'\n",
"scipy.linalg._flinalg b'$Revision: $'\n",
"scipy.ndimage 2.0\n",
"scipy.optimize._cobyla b'$Revision: $'\n",
"scipy.optimize._lbfgsb b'$Revision: $'\n",
"scipy.optimize._minpack 1.10 \n",
"scipy.optimize._nnls b'$Revision: $'\n",
"scipy.optimize._slsqp b'$Revision: $'\n",
"scipy.optimize.minpack2 b'$Revision: $'\n",
"scipy.signal.spline 0.2\n",
"scipy.sparse.linalg.eigen.arpack._arpack b'$Revision: $'\n",
"scipy.sparse.linalg.isolve._iterative b'$Revision: $'\n",
"scipy.special.specfun b'$Revision: $'\n",
"scipy.stats.mvn b'$Revision: $'\n",
"scipy.stats.statlib b'$Revision: $'\n",
"seaborn 0.8.1\n",
"seaborn.external.husl 2.1.0\n",
"seaborn.external.six 1.10.0\n",
"six 1.11.0\n",
"statsmodels 0.9.0\n",
"statsmodels.__init__ 0.9.0\n",
"traitlets 4.3.2\n",
"traitlets._version 4.3.2\n",
"urllib.request 3.6\n",
"zlib 1.0\n",
"zmq 17.1.0\n",
"zmq.sugar 17.1.0\n",
"zmq.sugar.version 17.1.0\n"
]
}
],
"source": [
"def print_imported_modules():\n",
" import sys\n",
" for name, val in sorted(sys.modules.items()):\n",
" if(hasattr(val, '__version__')): \n",
" print(val.__name__, val.__version__)\n",
"# else:\n",
"# print(val.__name__, \"(unknown version)\")\n",
"def print_sys_info():\n",
" import sys\n",
" import platform\n",
" print(sys.version)\n",
" print(platform.uname())\n",
"\n",
"import numpy as np\n",
"import pandas as pd\n",
"import matplotlib.pyplot as plt\n",
"import statsmodels.api as sm\n",
"import seaborn as sns\n",
"\n",
"print_sys_info()\n",
"print_imported_modules()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Loading and inspecting data\n",
"Let's start by reading data."
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"<div>\n",
"<style scoped>\n",
" .dataframe tbody tr th:only-of-type {\n",
" vertical-align: middle;\n",
" }\n",
"\n",
" .dataframe tbody tr th {\n",
" vertical-align: top;\n",
" }\n",
"\n",
" .dataframe thead th {\n",
" text-align: right;\n",
" }\n",
"</style>\n",
"<table border=\"1\" class=\"dataframe\">\n",
" <thead>\n",
" <tr style=\"text-align: right;\">\n",
" <th></th>\n",
" <th>Date</th>\n",
" <th>Count</th>\n",
" <th>Temperature</th>\n",
" <th>Pressure</th>\n",
" <th>Malfunction</th>\n",
" </tr>\n",
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>0</th>\n",
" <td>4/12/81</td>\n",
" <td>6</td>\n",
" <td>66</td>\n",
" <td>50</td>\n",
" <td>0</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1</th>\n",
" <td>11/12/81</td>\n",
" <td>6</td>\n",
" <td>70</td>\n",
" <td>50</td>\n",
" <td>1</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2</th>\n",
" <td>3/22/82</td>\n",
" <td>6</td>\n",
" <td>69</td>\n",
" <td>50</td>\n",
" <td>0</td>\n",
" </tr>\n",
" <tr>\n",
" <th>3</th>\n",
" <td>11/11/82</td>\n",
" <td>6</td>\n",
" <td>68</td>\n",
" <td>50</td>\n",
" <td>0</td>\n",
" </tr>\n",
" <tr>\n",
" <th>4</th>\n",
" <td>4/04/83</td>\n",
" <td>6</td>\n",
" <td>67</td>\n",
" <td>50</td>\n",
" <td>0</td>\n",
" </tr>\n",
" <tr>\n",
" <th>5</th>\n",
" <td>6/18/82</td>\n",
" <td>6</td>\n",
" <td>72</td>\n",
" <td>50</td>\n",
" <td>0</td>\n",
" </tr>\n",
" <tr>\n",
" <th>6</th>\n",
" <td>8/30/83</td>\n",
" <td>6</td>\n",
" <td>73</td>\n",
" <td>100</td>\n",
" <td>0</td>\n",
" </tr>\n",
" <tr>\n",
" <th>7</th>\n",
" <td>11/28/83</td>\n",
" <td>6</td>\n",
" <td>70</td>\n",
" <td>100</td>\n",
" <td>0</td>\n",
" </tr>\n",
" <tr>\n",
" <th>8</th>\n",
" <td>2/03/84</td>\n",
" <td>6</td>\n",
" <td>57</td>\n",
" <td>200</td>\n",
" <td>1</td>\n",
" </tr>\n",
" <tr>\n",
" <th>9</th>\n",
" <td>4/06/84</td>\n",
" <td>6</td>\n",
" <td>63</td>\n",
" <td>200</td>\n",
" <td>1</td>\n",
" </tr>\n",
" <tr>\n",
" <th>10</th>\n",
" <td>8/30/84</td>\n",
" <td>6</td>\n",
" <td>70</td>\n",
" <td>200</td>\n",
" <td>1</td>\n",
" </tr>\n",
" <tr>\n",
" <th>11</th>\n",
" <td>10/05/84</td>\n",
" <td>6</td>\n",
" <td>78</td>\n",
" <td>200</td>\n",
" <td>0</td>\n",
" </tr>\n",
" <tr>\n",
" <th>12</th>\n",
" <td>11/08/84</td>\n",
" <td>6</td>\n",
" <td>67</td>\n",
" <td>200</td>\n",
" <td>0</td>\n",
" </tr>\n",
" <tr>\n",
" <th>13</th>\n",
" <td>1/24/85</td>\n",
" <td>6</td>\n",
" <td>53</td>\n",
" <td>200</td>\n",
" <td>2</td>\n",
" </tr>\n",
" <tr>\n",
" <th>14</th>\n",
" <td>4/12/85</td>\n",
" <td>6</td>\n",
" <td>67</td>\n",
" <td>200</td>\n",
" <td>0</td>\n",
" </tr>\n",
" <tr>\n",
" <th>15</th>\n",
" <td>4/29/85</td>\n",
" <td>6</td>\n",
" <td>75</td>\n",
" <td>200</td>\n",
" <td>0</td>\n",
" </tr>\n",
" <tr>\n",
" <th>16</th>\n",
" <td>6/17/85</td>\n",
" <td>6</td>\n",
" <td>70</td>\n",
" <td>200</td>\n",
" <td>0</td>\n",
" </tr>\n",
" <tr>\n",
" <th>17</th>\n",
" <td>7/2903/85</td>\n",
" <td>6</td>\n",
" <td>81</td>\n",
" <td>200</td>\n",
" <td>0</td>\n",
" </tr>\n",
" <tr>\n",
" <th>18</th>\n",
" <td>8/27/85</td>\n",
" <td>6</td>\n",
" <td>76</td>\n",
" <td>200</td>\n",
" <td>0</td>\n",
" </tr>\n",
" <tr>\n",
" <th>19</th>\n",
" <td>10/03/85</td>\n",
" <td>6</td>\n",
" <td>79</td>\n",
" <td>200</td>\n",
" <td>0</td>\n",
" </tr>\n",
" <tr>\n",
" <th>20</th>\n",
" <td>10/30/85</td>\n",
" <td>6</td>\n",
" <td>75</td>\n",
" <td>200</td>\n",
" <td>2</td>\n",
" </tr>\n",
" <tr>\n",
" <th>21</th>\n",
" <td>11/26/85</td>\n",
" <td>6</td>\n",
" <td>76</td>\n",
" <td>200</td>\n",
" <td>0</td>\n",
" </tr>\n",
" <tr>\n",
" <th>22</th>\n",
" <td>1/12/86</td>\n",
" <td>6</td>\n",
" <td>58</td>\n",
" <td>200</td>\n",
" <td>1</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n",
"</div>"
],
"text/plain": [
" Date Count Temperature Pressure Malfunction\n",
"0 4/12/81 6 66 50 0\n",
"1 11/12/81 6 70 50 1\n",
"2 3/22/82 6 69 50 0\n",
"3 11/11/82 6 68 50 0\n",
"4 4/04/83 6 67 50 0\n",
"5 6/18/82 6 72 50 0\n",
"6 8/30/83 6 73 100 0\n",
"7 11/28/83 6 70 100 0\n",
"8 2/03/84 6 57 200 1\n",
"9 4/06/84 6 63 200 1\n",
"10 8/30/84 6 70 200 1\n",
"11 10/05/84 6 78 200 0\n",
"12 11/08/84 6 67 200 0\n",
"13 1/24/85 6 53 200 2\n",
"14 4/12/85 6 67 200 0\n",
"15 4/29/85 6 75 200 0\n",
"16 6/17/85 6 70 200 0\n",
"17 7/2903/85 6 81 200 0\n",
"18 8/27/85 6 76 200 0\n",
"19 10/03/85 6 79 200 0\n",
"20 10/30/85 6 75 200 2\n",
"21 11/26/85 6 76 200 0\n",
"22 1/12/86 6 58 200 1"
]
},
"execution_count": 2,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"data = pd.read_csv(\"https://app-learninglab.inria.fr/gitlab/moocrr-session1/moocrr-reproducibility-study/raw/master/data/shuttle.csv\")\n",
"data"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"We know from our previous experience on this data set that filtering data is a really bad idea. We will therefore process it as such."
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAEKCAYAAAD9xUlFAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAGBNJREFUeJzt3XuQnXWd5/H3t5MACYmAxMkwBAQGhpUCJkDLReaSCFqBKpN1AYUtwXEGM1uSskRHYWZdhmGdqpVRmXFlFGRxhC2NXEbIzmaWixAdprgFiOEmTA8gdEDAGCANIemkv/vHefrxpOnLOZ1++vQ5vF9VqZznOb9++vvtp09/+rn070RmIkkSQFerC5AkTR2GgiSpZChIkkqGgiSpZChIkkqGgiSpVFkoRMTVEfFSRDwywvMREV+PiJ6IWBcRR1dViySpMVUeKfwDsHiU508BDin+LQO+WWEtkqQGVBYKmfkT4FejDFkKXJM19wB7RsQ+VdUjSRrb9BZ+7n2B5+qWe4t1LwwdGBHLqB1NMHPmzGP222+/SSmwUQMDA3R1dd7lmU7tCzq3N/tqP5PV25NPPvnLzHzXWONaGQoxzLph59zIzCuBKwG6u7tzzZo1VdbVtNWrV7Nw4cJWlzHhOrUv6Nze7Kv9TFZvEfHzRsa1Mnp7gfpf+ecDz7eoFkkSrQ2FlcA5xV1IxwOvZuZbTh1JkiZPZaePIuL7wEJgbkT0An8JzADIzG8Bq4BTgR7gDeATVdUiSWpMZaGQmWeN8XwC51X1+SVJzevMy/mSpHExFCRJJUNBklQyFCRJJUNBklQyFCRJJUNBklQyFCRJJUNBklQyFCRJJUNBklQyFCRJJUNBklQyFCRJJUNBklQyFCRJJUNBklQyFCRJJUNBklQyFCRJJUNBklQyFCRJJUNBklQyFCRJJUNBklQyFCRJJUNBklQyFCRJJUNBklQyFCRJJUNBklQyFCRJJUNBklQyFCRJJUNBklSqNBQiYnFEPBERPRFx4TDP7x8Rd0bEQxGxLiJOrbIeSdLoKguFiJgGXA6cAhwGnBURhw0Z9kXgusw8CjgT+Puq6pEkja3KI4VjgZ7MfCoztwIrgKVDxiTwjuLxHsDzFdYjSRpDZGY1G444HVicmecWy2cDx2Xm8rox+wC3AnsBuwMnZ+YDw2xrGbAMYN68ecesWLGikprHq6+vj9mzZ7e6jAnXqX1B5/ZmX+1nsnpbtGjRA5nZPda46RXWEMOsG5pAZwH/kJlfjYgTgGsj4vDMHNjhgzKvBK4E6O7uzoULF1ZR77itXr2aqVbTROjUvqBze7Ov9jPVeqvy9FEvsF/d8nzeenroT4DrADLzbmA3YG6FNUmSRlFlKNwPHBIRB0bELtQuJK8cMuZZ4CSAiHgPtVB4ucKaJEmjqCwUMnMbsBy4BXic2l1Gj0bEJRGxpBj2OeCTEfFT4PvAH2VVFzkkSWOq8poCmbkKWDVk3UV1jx8DTqyyBklS4/yLZklSyVCQJJUMBUlSyVCQJJUMBUlSyVCQJJUMBUlSyVCQJJUMBUlSyVCQJJUMBUlSyVCQJJUMBUlSyVCQJJUMBUlSyVCQJJUMBUlSyVCQJJUMBUlSyVCQJJUMBUlSyVCQJJUMBUlSyVCQJJUMBUlSyVCQJJUMBUlSyVCQJJUMBUlSyVCQJJUMBUlSyVCQJJUMBUlSyVCQJJUqDYWIWBwRT0RET0RcOMKYj0TEYxHxaER8r8p6JEmjm97IoIg4PDMfaWbDETENuBz4ANAL3B8RKzPzsboxhwB/DpyYmRsj4jea+RySpInV6JHCtyLivoj4VETs2eDHHAv0ZOZTmbkVWAEsHTLmk8DlmbkRIDNfanDbkqQKRGY2NrD2W/0fA2cA9wHfyczbRhl/OrA4M88tls8GjsvM5XVjbgKeBE4EpgEXZ+b/G2Zby4BlAPPmzTtmxYoVjXU3Sfr6+pg9e3ary5hwndoXdG5v9tV+Jqu3RYsWPZCZ3WMOzMyG/1H7wX0asB54HPgZ8J9GGHsGcFXd8tnA/xwy5p+AHwIzgAOpnWbac7QajjnmmJxq7rzzzlaXUIlO7Suzc3uzr/YzWb0Ba7KBn/MNnT6KiCMj4rIiCN4PfCgz31M8vmyED+sF9qtbng88P8yYmzOzPzOfBp4ADmmkJknSxGv0msI3gAeB383M8zLzQYDMfB744ggfcz9wSEQcGBG7AGcCK4eMuQlYBBARc4HfAZ5qrgVJ0kRp6O4j4FRgc2ZuB4iILmC3zHwjM68d7gMyc1tELAduoXba6erMfDQiLqF2GLOyeO6DEfEYsB34fGZu2MmeJEnj1Ggo3A6cDPQVy7OAW4H3jfZBmbkKWDVk3UV1jxP4bPFPktRijZ4+2i0zBwOB4vGsakqSJLVKo6HwekQcPbgQEccAm6spSZLUKo2ePvoMcH1EDN49tA/w0WpKkiS1SkOhkJn3R8R/AA4FAvhZZvZXWpkkadI1eqQA8F7ggOJjjooIMvOaSqqSJLVEoxPiXQv8NrCW2q2jAAkYCpLUQRo9UugGDituIZUkdahG7z56BPjNKguRJLVeo0cKc4HHIuI+YMvgysxcUklVkqSWaDQULq6yCEnS1NDoLak/joh3A4dk5u0RMYvafEaSpA7S6NTZnwRuAK4oVu1LbYZTSVIHafRC83nU3h3tNYDM/DfA91OWpA7TaChsydr7LAMQEdOp/Z2CJKmDNBoKP46IvwBmRsQHgOuB/1NdWZKkVmg0FC4EXgYeBv6U2nskjPSOa5KkNtXo3UcDwLeLf5KkDtXo3EdPM8w1hMw8aMIrkiS1TDNzHw3aDTgDeOfElyNJaqWGrilk5oa6f+sz82+B91dcmyRpkjV6+ujousUuakcOcyqpSJLUMo2ePvpq3eNtwDPARya8GklSSzV699GiqguRJLVeo6ePPjva85n5tYkpR5LUSs3cffReYGWx/CHgJ8BzVRQlSWqNZt5k5+jM3AQQERcD12fmuVUVJkmafI1Oc7E/sLVueStwwIRXI0lqqUaPFK4F7ouIH1L7y+YPA9dUVpUkqSUavfvoryPin4HfL1Z9IjMfqq4sSVIrNHr6CGAW8Fpm/h3QGxEHVlSTJKlFGn07zr8ELgD+vFg1A/jfVRUlSWqNRo8UPgwsAV4HyMzncZoLSeo4jYbC1sxMiumzI2L36kqSJLVKo6FwXURcAewZEZ8Ebsc33JGkjtPo3UdfKd6b+TXgUOCizLyt0sokSZNuzCOFiJgWEbdn5m2Z+fnM/LNGAyEiFkfEExHRExEXjjLu9IjIiOgeaYwkqXpjhkJmbgfeiIg9mtlwREwDLgdOAQ4DzoqIw4YZNwf4NHBvM9uXJE28Rv+i+U3g4Yi4jeIOJIDM/PQoH3Ms0JOZTwFExApgKfDYkHH/HbgU+LNGi5YkVaPRUPi/xb9m7MuOs6j2AsfVD4iIo4D9MvOfImLEUIiIZcAygHnz5rF69eomS6lWX1/flKtpInRqX9C5vdlX+5lqvY0aChGxf2Y+m5nfHce2Y5h1WbftLuAy4I/G2lBmXglcCdDd3Z0LFy4cRznVWb16NVOtponQqX1B5/ZmX+1nqvU21jWFmwYfRMSNTW67F9ivbnk+8Hzd8hzgcGB1RDwDHA+s9GKzJLXOWKFQ/9v+QU1u+37gkIg4MCJ2Ac7k12/SQ2a+mplzM/OAzDwAuAdYkplrmvw8kqQJMlYo5AiPx5SZ24DlwC3A48B1mfloRFwSEUuaK1OSNBnGutD8uxHxGrUjhpnFY4rlzMx3jPbBmbkKWDVk3UUjjF3YUMWSpMqMGgqZOW2yCpEktV4z76cgSepwhoIkqWQoSJJKhoIkqfS2CYUNfVv46XOvsKFvS6tLkdSEDX1b2Ny/3dfuJHlbhMLNa9dz4pfv4GNX3cuJX76DlWvXt7okSQ0YfO0+/fLrvnYnSceHwoa+LVxw4zre7B9g05ZtvNk/wBduXOdvHdIUV//a3Z7pa3eSdHwo9G7czIyuHduc0dVF78bNLapIUiN87bZGx4fC/L1m0j8wsMO6/oEB5u81s0UVSWqEr93W6PhQ2Hv2rlx62pHsNqOLObtOZ7cZXVx62pHsPXvXVpcmaRT1r91pEb52J0mjb7LT1pYs2JcTD55L78bNzN9rpt9UUpsYfO3ed/dd/OuS3/O1OwneFqEAtd86/IaS2s/es3dl5oxpvn4nScefPpIkNc5QkCSVDAVJUslQkCSVDAVJUslQkCSVDAVJUslQkCSVDAVJUslQkCSVDAVJUslQkCSVDAVJUslQkCSVDAVJUslQkCSVDAVJUslQkCSVDAVJUslQkCSVDAVJUqnSUIiIxRHxRET0RMSFwzz/2Yh4LCLWRcSPIuLdVdYjSRpdZaEQEdOAy4FTgMOAsyLisCHDHgK6M/NI4Abg0qrqkSSNrcojhWOBnsx8KjO3AiuApfUDMvPOzHyjWLwHmF9hPZKkMURmVrPhiNOBxZl5brF8NnBcZi4fYfw3gF9k5peGeW4ZsAxg3rx5x6xYsaKSmserr6+P2bNnt7qMCdepfUHn9mZf7Weyelu0aNEDmdk91rjpFdYQw6wbNoEi4mNAN/CHwz2fmVcCVwJ0d3fnwoULJ6jEibF69WqmWk0ToVP7gs7tzb7az1TrrcpQ6AX2q1ueDzw/dFBEnAz8V+APM3NLhfVIksZQ5TWF+4FDIuLAiNgFOBNYWT8gIo4CrgCWZOZLFdYiSWpAZaGQmduA5cAtwOPAdZn5aERcEhFLimF/A8wGro+ItRGxcoTNSZImQZWnj8jMVcCqIesuqnt8cpWfv51t6NtC78bNzN9rJnvP3nXCxraTTu2rKj0vbmLjG/30vLiJg+fNaXU5alOVhoLG5+a167ngxnXM6Oqif2CAS087kiUL9t3pse2kU/uqykU3Pcw19zzL547YxvmX/YRzTtifS5Ye0eqy1Iac5mKK2dC3hQtuXMeb/QNs2rKNN/sH+MKN69jQ99Zr8M2MbSed2ldVel7cxDX3PLvDumvufpaeFze1qCK1M0NhiunduJkZXTvulhldXfRu3LxTY9tJp/ZVlbXPvdLUemk0hsIUM3+vmfQPDOywrn9ggPl7zdypse2kU/uqyoL99mxqvTQaQ2GK2Xv2rlx62pHsNqOLObtOZ7cZXVx62pHDXmhtZmw76dS+qnLwvDmcc8L+O6w754T9vdiscfFC8xS0ZMG+nHjw3IbuvGlmbDvp1L6qcsnSIzjn+AN4+IF7uP384w0EjZuhMEXtPXvXhn8QNjO2nXRqX1U5eN4cemfNMBC0Uzx9JEkqGQqSpJKhIEkqGQqSpJKhIEkqGQqSpJKhIEkqGQqSpJKhIEkqGQqSpJKhIEkqGQqSpJKhIEkqGQqSpJKhIEkqGQqSpJKhIEkqGQqSpJKhIEkqGQqSpJKhIEkqGQqSpJKhIEkqGQqSpJKhIEkqGQqSpJKhIEkqGQqSpFKloRARiyPiiYjoiYgLh3l+14j4QfH8vRFxQJX1SM3a0LeFnz73Chv6tow6bs3TG/jarU+w5ukNE7bNZsf2vLiJjW/00/PipjHHNqOqepv5/Jv7tzf8NbhhzXMd9zWocrtDTa9qwxExDbgc+ADQC9wfESsz87G6YX8CbMzMgyPiTODLwEerqklqxs1r13PBjeuY0dVF/8AAl552JEsW7PuWcR+76h7u6qmFwdfv6OH3D96ba889fqe22ezYi256mGvueZbPHbGN8y/7CeecsD+XLD1inJ1XX2+zn//T7+nn/C/f0dDXYFCnfA2q3O5wqjxSOBboycynMnMrsAJYOmTMUuC7xeMbgJMiIiqsSWrIhr4tXHDjOt7sH2DTlm282T/AF25c95bf0tY8vaEMhEH/0rNh2COGRrfZ7NieFzft8MMQ4Jq7n93p35arqnc8n3975tvya1DldkcSmVnNhiNOBxZn5rnF8tnAcZm5vG7MI8WY3mL534sxvxyyrWXAsmLxUOCJSooev7nAL8cc1X46tS8Yo7eYMXPW9L32+Z3o6po2uC4HBrZv2/jCk9m/+Y3BddPmzP2tabvvuc/Qj9/++isvbN/0y+fHs81mx3bN2mPv6e941wEA2994lWmz9gBg22svPzPwxqtjn8/aya9Bs2PH8/kH+2rka1CvTb4GE/K92IB3Z+a7xhpU2ekjYLjf+IcmUCNjyMwrgSsnoqgqRMSazOxudR0TrVP7gs7tLSLWbHv1JftqI1Pte7HK00e9wH51y/OB50caExHTgT2AX1VYkyRpFFWGwv3AIRFxYETsApwJrBwyZiXw8eLx6cAdWdX5LEnSmCo7fZSZ2yJiOXALMA24OjMfjYhLgDWZuRL4X8C1EdFD7QjhzKrqqdiUPbW1kzq1L+jc3uyr/Uyp3iq70CxJaj/+RbMkqWQoSJJKhsI4RMQzEfFwRKyNiDXFuosjYn2xbm1EnNrqOpsVEXtGxA0R8bOIeDwiToiId0bEbRHxb8X/e7W6zmaN0Fcn7K9D6+pfGxGvRcRn2n2fjdJXJ+yz8yPi0Yh4JCK+HxG7FTfj3Fvsrx8UN+a0rkavKTQvIp4Buuv/yC4iLgb6MvMrraprZ0XEd4F/ycyrim/MWcBfAL/KzP9RzF+1V2Ze0NJCmzRCX5+hzfdXvWJamfXAccB5tPk+GzSkr0/QxvssIvYF7gIOy8zNEXEdsAo4FfjHzFwREd8CfpqZ32xVnR4pCICIeAfwB9TuCCMzt2bmK+w4Fcl3gf/YmgrHZ5S+Os1JwL9n5s9p8302RH1fnWA6MLP4u6xZwAvA+6lN8wNTYH8ZCuOTwK0R8UAxBceg5RGxLiKubrdDduAg4GXgOxHxUERcFRG7A/My8wWA4v/faGWR4zBSX9De+2uoM4HvF4/bfZ/Vq+8L2nifZeZ64CvAs9TC4FXgAeCVzNxWDOsFqpnprkGGwvicmJlHA6cA50XEHwDfBH4bWEBth3+1hfWNx3TgaOCbmXkU8DrwlunO29BIfbX7/ioVp8SWANe3upaJNExfbb3PihBbChwI/BawO7WfIUO19Jy+oTAOmfl88f9LwA+BYzPzxczcnpkDwLepzRLbTnqB3sy8t1i+gdoP0xcjYh+A4v+XWlTfeA3bVwfsr3qnAA9m5ovFcrvvs0E79NUB++xk4OnMfDkz+4F/BN4H7FmcToLhpwOaVIZCkyJi94iYM/gY+CDwyOCLsPBh4JFW1DdemfkL4LmIOLRYdRLwGDtORfJx4OYWlDduI/XV7vtriLPY8RRLW++zOjv01QH77Fng+IiYFRHBr19jd1Kb5gemwP7y7qMmRcRB1I4OoHZq4nuZ+dcRcS21w9oEngH+dPC8bruIiAXAVcAuwFPU7vboAq4D9qf2TX1GZrbVpIUj9PV12nx/AUTELOA54KDMfLVYtzftv8+G66sTXmN/Re2NxLYBDwHnUruGsAJ4Z7HuY5lZ7durjVajoSBJGuTpI0lSyVCQJJUMBUlSyVCQJJUMBUlSqbJ3XpMmW3Er5o+Kxd8EtlOb4gJqf2C4tSWFjSIi/hhYVfw9hdRy3pKqjjSVZq2NiGmZuX2E5+4Clmfm2ia2N71urhxpQnn6SG8LEfHxiLivmIf/7yOiKyKmR8QrEfE3EfFgRNwSEcdFxI8j4qnB+foj4tyI+GHx/BMR8cUGt/uliLgPODYi/ioi7i/m0f9W1HyU2h9j/aD4+F0iojci9iy2fXxE3F48/lJEXBERt1Gb3G96RHyt+NzrIuLcyf+qqhMZCup4EXE4tWkR3peZC6idNj2zeHoP4NZigsOtwMXUph84A7ikbjPHFh9zNPCfI2JBA9t9MDOPzcy7gb/LzPcCRxTPLc7MHwBrgY9m5oIGTm8dBXwoM88GlgEvZeaxwHupTcy4/3i+PlI9ryno7eBkaj8419SmnGEmtSkUADZn5m3F44eBVzNzW0Q8DBxQt41bMnMjQETcBPwetdfPSNvdyq+nQwE4KSI+D+wGzKU2ZfI/N9nHzZn5ZvH4g8B7IqI+hA6hNq2FNG6Ggt4OArg6M//bDitrM1PW/3Y+AGype1z/+hh68S3H2O7mLC7YFfP4fIPa7KzrI+JL1MJhONv49RH80DGvD+npU5n5I6QJ5OkjvR3cDnwkIuZC7S6lcZxq+WDU3ut5FrU58f+1ie3OpBYyvyxm2D2t7rlNwJy65WeAY4rH9eOGugX41OCUy1F7X+OZTfYkvYVHCup4mflwMTvl7RHRBfQD/4Xm5q2/C/getTd5uXbwbqFGtpuZG6L2PtGPAD8H7q17+jvAVRGxmdp1i4uBb0fEL4D7RqnnCmqzoK4tTl29RC2spJ3iLanSGIo7ew7PzM+0uhapap4+kiSVPFKQJJU8UpAklQwFSVLJUJAklQwFSVLJUJAklf4/GKF1l7kqzyEAAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"%matplotlib inline\n",
"pd.set_option('mode.chained_assignment',None) # this removes a useless warning from pandas\n",
"import matplotlib.pyplot as plt\n",
"\n",
"data[\"Frequency\"]=data.Malfunction/data.Count\n",
"data.plot(x=\"Temperature\",y=\"Frequency\",kind=\"scatter\",ylim=[0,1])\n",
"plt.grid(True)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Logistic regression\n",
"\n",
"Let's assume O-rings independently fail with the same probability which solely depends on temperature. A logistic regression should allow us to estimate the influence of temperature."
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"<table class=\"simpletable\">\n",
"<caption>Generalized Linear Model Regression Results</caption>\n",
"<tr>\n",
" <th>Dep. Variable:</th> <td>Frequency</td> <th> No. Observations: </th> <td> 23</td> \n",
"</tr>\n",
"<tr>\n",
" <th>Model:</th> <td>GLM</td> <th> Df Residuals: </th> <td> 21</td> \n",
"</tr>\n",
"<tr>\n",
" <th>Model Family:</th> <td>Binomial</td> <th> Df Model: </th> <td> 1</td> \n",
"</tr>\n",
"<tr>\n",
" <th>Link Function:</th> <td>logit</td> <th> Scale: </th> <td> 1.0000</td> \n",
"</tr>\n",
"<tr>\n",
" <th>Method:</th> <td>IRLS</td> <th> Log-Likelihood: </th> <td> -3.9210</td> \n",
"</tr>\n",
"<tr>\n",
" <th>Date:</th> <td>Wed, 24 Oct 2018</td> <th> Deviance: </th> <td> 3.0144</td> \n",
"</tr>\n",
"<tr>\n",
" <th>Time:</th> <td>11:05:55</td> <th> Pearson chi2: </th> <td> 5.00</td> \n",
"</tr>\n",
"<tr>\n",
" <th>No. Iterations:</th> <td>6</td> <th> Covariance Type: </th> <td>nonrobust</td>\n",
"</tr>\n",
"</table>\n",
"<table class=\"simpletable\">\n",
"<tr>\n",
" <td></td> <th>coef</th> <th>std err</th> <th>z</th> <th>P>|z|</th> <th>[0.025</th> <th>0.975]</th> \n",
"</tr>\n",
"<tr>\n",
" <th>Intercept</th> <td> 5.0850</td> <td> 7.477</td> <td> 0.680</td> <td> 0.496</td> <td> -9.570</td> <td> 19.740</td>\n",
"</tr>\n",
"<tr>\n",
" <th>Temperature</th> <td> -0.1156</td> <td> 0.115</td> <td> -1.004</td> <td> 0.316</td> <td> -0.341</td> <td> 0.110</td>\n",
"</tr>\n",
"</table>"
],
"text/plain": [
"<class 'statsmodels.iolib.summary.Summary'>\n",
"\"\"\"\n",
" Generalized Linear Model Regression Results \n",
"==============================================================================\n",
"Dep. Variable: Frequency No. Observations: 23\n",
"Model: GLM Df Residuals: 21\n",
"Model Family: Binomial Df Model: 1\n",
"Link Function: logit Scale: 1.0000\n",
"Method: IRLS Log-Likelihood: -3.9210\n",
"Date: Wed, 24 Oct 2018 Deviance: 3.0144\n",
"Time: 11:05:55 Pearson chi2: 5.00\n",
"No. Iterations: 6 Covariance Type: nonrobust\n",
"===============================================================================\n",
" coef std err z P>|z| [0.025 0.975]\n",
"-------------------------------------------------------------------------------\n",
"Intercept 5.0850 7.477 0.680 0.496 -9.570 19.740\n",
"Temperature -0.1156 0.115 -1.004 0.316 -0.341 0.110\n",
"===============================================================================\n",
"\"\"\""
]
},
"execution_count": 4,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"import statsmodels.api as sm\n",
"\n",
"data[\"Success\"]=data.Count-data.Malfunction\n",
"data[\"Intercept\"]=1\n",
"\n",
"logmodel=sm.GLM(data['Frequency'], data[['Intercept','Temperature']], \n",
" family=sm.families.Binomial(sm.families.links.logit)).fit()\n",
"\n",
"logmodel.summary()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"The maximum likelyhood estimator of the intercept and of Temperature are thus $\\hat{\\alpha}=5.0849$ and $\\hat{\\beta}=-0.1156$. This **corresponds** to the values from the article of Dalal *et al.* The standard errors are $s_{\\hat{\\alpha}} = 7.477$ and $s_{\\hat{\\beta}} = 0.115$, which is **different** from the $3.052$ and $0.04702$ reported by Dallal *et al.* The deviance is $3.01444$ with 21 degrees of freedom. I cannot find any value similar to the Goodness of fit ($G^2=18.086$) reported by Dalal *et al.* There seems to be something wrong. Oh I know, I haven't indicated that my observations are actually the result of 6 observations for each rocket launch. Let's indicate these weights (since the weights are always the same throughout all experiments, it does not change the estimates of the fit but it does influence the variance estimates)."
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"<table class=\"simpletable\">\n",
"<caption>Generalized Linear Model Regression Results</caption>\n",
"<tr>\n",
" <th>Dep. Variable:</th> <td>Frequency</td> <th> No. Observations: </th> <td> 23</td> \n",
"</tr>\n",
"<tr>\n",
" <th>Model:</th> <td>GLM</td> <th> Df Residuals: </th> <td> 21</td> \n",
"</tr>\n",
"<tr>\n",
" <th>Model Family:</th> <td>Binomial</td> <th> Df Model: </th> <td> 1</td> \n",
"</tr>\n",
"<tr>\n",
" <th>Link Function:</th> <td>logit</td> <th> Scale: </th> <td> 1.0000</td> \n",
"</tr>\n",
"<tr>\n",
" <th>Method:</th> <td>IRLS</td> <th> Log-Likelihood: </th> <td> -23.526</td> \n",
"</tr>\n",
"<tr>\n",
" <th>Date:</th> <td>Wed, 24 Oct 2018</td> <th> Deviance: </th> <td> 18.086</td> \n",
"</tr>\n",
"<tr>\n",
" <th>Time:</th> <td>11:05:55</td> <th> Pearson chi2: </th> <td> 30.0</td> \n",
"</tr>\n",
"<tr>\n",
" <th>No. Iterations:</th> <td>6</td> <th> Covariance Type: </th> <td>nonrobust</td>\n",
"</tr>\n",
"</table>\n",
"<table class=\"simpletable\">\n",
"<tr>\n",
" <td></td> <th>coef</th> <th>std err</th> <th>z</th> <th>P>|z|</th> <th>[0.025</th> <th>0.975]</th> \n",
"</tr>\n",
"<tr>\n",
" <th>Intercept</th> <td> 5.0850</td> <td> 3.052</td> <td> 1.666</td> <td> 0.096</td> <td> -0.898</td> <td> 11.068</td>\n",
"</tr>\n",
"<tr>\n",
" <th>Temperature</th> <td> -0.1156</td> <td> 0.047</td> <td> -2.458</td> <td> 0.014</td> <td> -0.208</td> <td> -0.023</td>\n",
"</tr>\n",
"</table>"
],
"text/plain": [
"<class 'statsmodels.iolib.summary.Summary'>\n",
"\"\"\"\n",
" Generalized Linear Model Regression Results \n",
"==============================================================================\n",
"Dep. Variable: Frequency No. Observations: 23\n",
"Model: GLM Df Residuals: 21\n",
"Model Family: Binomial Df Model: 1\n",
"Link Function: logit Scale: 1.0000\n",
"Method: IRLS Log-Likelihood: -23.526\n",
"Date: Wed, 24 Oct 2018 Deviance: 18.086\n",
"Time: 11:05:55 Pearson chi2: 30.0\n",
"No. Iterations: 6 Covariance Type: nonrobust\n",
"===============================================================================\n",
" coef std err z P>|z| [0.025 0.975]\n",
"-------------------------------------------------------------------------------\n",
"Intercept 5.0850 3.052 1.666 0.096 -0.898 11.068\n",
"Temperature -0.1156 0.047 -2.458 0.014 -0.208 -0.023\n",
"===============================================================================\n",
"\"\"\""
]
},
"execution_count": 5,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"logmodel=sm.GLM(data['Frequency'], data[['Intercept','Temperature']], \n",
" family=sm.families.Binomial(sm.families.links.logit),\n",
" var_weights=data['Count']).fit()\n",
"\n",
"logmodel.summary()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Good, now I have recovered the asymptotic standard errors $s_{\\hat{\\alpha}}=3.052$ and $s_{\\hat{\\beta}}=0.047$.\n",
"The Goodness of fit (Deviance) indicated for this model is $G^2=18.086$ with 21 degrees of freedom (Df Residuals).\n",
"\n",
"**I have therefore managed to fully replicate the results of the Dalal *et al.* article**."
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Predicting failure probability\n",
"The temperature when launching the shuttle was 31°F. Let's try to estimate the failure probability for such temperature using our model.:"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAEKCAYAAADpfBXhAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzt3Xl8VOXZ//HPNUs2sgABwhI2NYDIngUQa8EqoFXcUEDEpSD2qUutlVb6WLVWuzz0+blXoYBrFakVROsjCIoLIgQEWWVHSNiXhITsyfX7YwYMMZAhmWSWXO/XK6/MOXOfc647J/nOyZkz9xFVxRhjTHhxBLoAY4wx/mfhbowxYcjC3RhjwpCFuzHGhCELd2OMCUMW7sYYE4ZqDHcRmSkiB0Rk3WmeFxF5RkS2isgaEenn/zKNMcacDV+O3F8Ghp/h+cuBFO/XROCFupdljDGmLmoMd1X9DDhyhiZXA6+qx1dAUxFp468CjTHGnD2XH9bRDthdaTrLO29v1YYiMhHP0T3R0dGp7du3r9UGKyoqcDjC4+0C60vwCZd+gPUlWNWlL5s3bz6kqi1rauePcJdq5lU7poGqTgOmAaSlpemKFStqtcHFixczePDgWi0bbKwvwSdc+gHWl2BVl76IyHe+tPPHy2AWUPkQPBnY44f1GmOMqSV/hPs84BbvVTMDgFxV/cEpGWOMMQ2nxtMyIvImMBhoISJZwCOAG0BVXwQ+AK4AtgIFwO31Vawxxhjf1BjuqjqmhucVuMtvFRljQkJpaSlZWVkUFRU1yPYSEhLYuHFjg2yrvvnSl6ioKJKTk3G73bXahj/eUDXGNEJZWVnExcXRqVMnRKq7rsK/8vLyiIuLq/ftNISa+qKqHD58mKysLDp37lyrbYTHdUXGmAZXVFREYmJigwR7YyMiJCYm1um/Igt3Y0ytWbDXn7r+bC3cjTEmDNk5d2NMyHI6nfTs2fPk9Ny5c+nUqVPgCgoiFu7GmJAVHR3N6tWrT/t8WVkZLlfjjDk7LWOMCSsvv/wyN9xwA1dddRVDhw4FYMqUKaSnp9OrVy8eeeSRk22feOIJunbtyqWXXsqYMWP429/+BsDgwYM5MTzKoUOHTv43UF5ezqRJk06ua+rUqcD3wwmMHDmSbt26MXbsWDxXiUNmZiYXXnghvXv3JiMjg7y8PIYNG3bKi9KgQYNYs2aNX38OjfMlzRjjV394bz0b9hzz6zq7t43nkasuOGObwsJC+vTpA0Dnzp2ZM2cOAEuXLmXNmjU0b96cBQsWsGXLFpYvX46qMmLECD777DOaNGnCrFmzWLVqFWVlZfTr14/U1NQzbm/GjBkkJCSQmZlJcXExgwYNOvkCsmrVKtavX0/btm0ZNGgQS5YsISMjg1GjRvHWW2+Rnp7OsWPHiI6O5pZbbuHll1/mqaeeYvPmzRQXF9OrVy8//NS+Z+FujAlZpzstc9lll9G8eXMAFixYwIIFC+jbty8A+fn5bNmyhby8PK699lpiYmIAGDFiRI3bW7BgAWvWrOHtt98GIDc3ly1bthAREUFGRgbJyckA9OnTh507d5KQkECbNm1IT08HID4+HoBrr72WQYMGMWXKFGbOnMltt91Wtx9ENSzcjTF1VtMRdkNr0qTJyceqyuTJk7nzzjtPafPUU0+d9nJDl8tFRUUFwCnXmqsqzz77LMOGDTul/eLFi4mMjDw57XQ6KSsrQ1Wr3UZMTAyXXXYZ7777LrNnz6a2I+SeiZ1zN8aEtWHDhjFz5kzy8/MByM7O5sCBA1x88cXMmTOHwsJC8vLyeO+9904u06lTJ1auXAlw8ij9xLpeeOEFSktLAdi8eTPHjx8/7ba7devGnj17yMzMBDyfTC0rKwNgwoQJ3HvvvaSnp5/8L8Of7MjdGBPWhg4dysaNGxk4cCAAsbGxvP766/Tr149Ro0bRp08fOnbsyI9+9KOTyzzwwAPceOONvPbaa1xyySUn50+YMIGdO3fSr18/VJWWLVsyd+7c0247IiKCt956i3vuuYfCwkKio6NZuHAhAKmpqcTHx3P77fU01qKqBuQrNTVVa+uTTz6p9bLBxvoSfMKlH6r125cNGzbU27qrc+zYsXpd/yOPPKJTpkyp122ccOzYMc3OztaUlBQtLy8/bbvqfsbACvUhY+20jDHGNLA33niD/v3788QTT9TbrQPttIwxxgCPPvpog23rpptu+sEbvP5mR+7GmFpTrfZ2ycYP6vqztXA3xtRKVFQUhw8ftoCvB+odzz0qKqrW67DTMsaYWklOTiYrK4uDBw82yPaKiorqFHbBxJe+nLgTU21ZuBtjasXtdtf6LkG1sXjx4pOfMg11DdEXOy1jjDFhyMLdGGPCkIW7McaEIQt3Y4wJQxbuxhgThizcjTEmDFm4G2NMGLJwN8aYMGThbowxYcjC3RhjwlDIhfuBY0V8mlVqgxUZY8wZhFy4v75sFy+tK2HCKys4mFcc6HKMMSYohVy43/eTFG7qFsHnWw8x/KnPWLhhf6BLMsaYoBNy4e5wCEM7ufnPPReRFB/FhFdX8Mi76ygqLQ90acYYEzRCLtxPSEmKY85dFzL+os68svQ7rnl+CVsP5Ae6LGOMCQohG+4AkS4nv7+yOy/dns6BvGJGPPcF767ODnRZxhgTcD6Fu4gMF5FNIrJVRB6s5vkOIvKJiKwSkTUicoX/Sz29IV1b8Z97L+KCtvH8ctZqfjdnLcVldprGGNN41RjuIuIEngcuB7oDY0Ske5VmDwGzVbUvMBr4u78LrUmbhGjevGMAP//xubyxbBc3vriU7JzChi7DGGOCgi9H7hnAVlXdrqolwCzg6iptFIj3Pk4A9vivRN+5nA4evLwbU8elsv3gca585nO+3HooEKUYY0xASU0fBhKRkcBwVZ3gnR4H9FfVuyu1aQMsAJoBTYBLVXVlNeuaCEwESEpKSp01a1atis7Pzyc2NvaMbfYdr+CZVUXsO66M6hrB0I4uRKRW26tPvvQlVIRLX8KlH2B9CVZ16cuQIUNWqmpajQ1V9YxfwA3A9ErT44Bnq7S5H/i19/FAYAPgONN6U1NTtbY++eQTn9rlFZXqHa9kasffvq/3v7Vai0rLar3N+uJrX0JBuPQlXPqhan0JVnXpC7BCa8htVfXptEwW0L7SdDI/PO0yHpjtfbFYCkQBLXxYd72KjXTx4s2p3HdpCv/+Ooux/1jGoXz7VKsxJvz5Eu6ZQIqIdBaRCDxvmM6r0mYX8BMAETkfT7gf9GehteVwCPdd2oXnb+rHuj25XP3cEjbtywt0WcYYU69qDHdVLQPuBuYDG/FcFbNeRB4TkRHeZr8G7hCRb4A3gdu8/z4EjZ/2asPsOwdSWl7ByBe+5PMtQfHaY4wx9cKn69xV9QNV7aKq56rqE955D6vqPO/jDao6SFV7q2ofVV1Qn0XXVq/kpsy9axDtmkVz+0uZvJW5K9AlGWNMvQjpT6jWRtum0fzr5wO58LwW/Pbfa3lq4WYbPtgYE3YaXbgDxEW5mXFrGiNTk3lq4RYmv7OWsvKKQJdljDF+4wp0AYHidjqYMrIXbROieObjrRzKL+a5m/oR5XYGujRjjKmzRnnkfoKIcP/Qrvzx6gtY9O0BbpmxnNzC0kCXZYwxddaow/2EcQM78czovqzafZRRU5faHZ6MMSHPwt3rqt5tmXFrOjsPH2fUVBt0zBgT2izcK7m4S0teG9+fg3nF3PjiUnYeOh7okowxplYs3KtI79ScNycOoLC0nBunLrW7OxljQpKFezV6tEtg1sQBVCiMnraUb/cdC3RJxhhzVizcT6NLUhyz7xyAy+Fg9LSvWL8nN9AlGWOMzyzcz+CclrG8decAYtxOxk5fxrpsC3hjTGiwcK9Bx8QmzJo40ALeGBNSLNx90CExhlkTB9IkwsnNM5axca+dgzfGBDcLdx91SIzhzYkDiHI5uXn6MrbstzHhjTHBy8L9LHRMbMIbd/TH6RDG/GMZ2w/aZZLGmOBk4X6WzmkZyxt39EdVGTt9GbuPFAS6JGOM+QEL91o4r1Ucr43vT0FJOWOnL2NfblGgSzLGmFNYuNdS97bxvPKzDI4cL+HmGcs4crwk0CUZY8xJFu510Kd9U6bfmsbuIwXcOnM5eUU2XLAxJjhYuNfRgHMSeeHmfmzce4wJr6ygqLQ80CUZY4yFuz9c0i2J/72xN8t3HuHuN1bZLfuMMQFn4e4nV/dpx6NXXcDCjfuZ/M5au+m2MSagGu09VOvDrRd24vDxEp5ZtIXE2EgevLxboEsyxjRSFu5+9qtLUzhyvJgXP91Gq7hIfnZR50CXZIxphCzc/UxE+MOIHhzKK+GP/9lAy7hIrurdNtBlGWMaGTvnXg+cDuGp0X1I79ic+2ev5suthwJdkjGmkbFwrydRbif/uCWNzi2acOdrK+1uTsaYBmXhXo8SYty8dHsGMZFObn8pk725hYEuyRjTSFi417N2TaN56bYM8orKuP2lTPsUqzGmQVi4N4DubeN54eZ+bD2Qzy/++TWl9iEnY0w9s3BvID9Kacmfru3J51sO8dCcdfYhJ2NMvbJLIRvQjent2X20gGc/3kqHxBjuGnJeoEsyxoQpC/cGdv9lXdh1pIAp8zfRMTGG2EAXZIwJS3ZapoGJCH+9vhdpHZtx/+xv2HrURpE0xvifT+EuIsNFZJOIbBWRB0/T5kYR2SAi60XkDf+WGV6i3E6m3ZJGm4Qonl5VZLfqM8b4XY3hLiJO4HngcqA7MEZEuldpkwJMBgap6gXAffVQa1hp3iSCmbelU14B41/J5JhdImmM8SNfjtwzgK2qul1VS4BZwNVV2twBPK+qRwFU9YB/ywxP57aM5e6+UWw/eNzGgTfG+JXUdEmeiIwEhqvqBO/0OKC/qt5dqc1cYDMwCHACj6rqh9WsayIwESApKSl11qxZtSo6Pz+f2NjweCsyPz+flUcjeWl9CZd2cHFz98hAl1Rr4bJfwqUfYH0JVnXpy5AhQ1aqalpN7Xy5WkaqmVf1FcEFpACDgWTgcxHpoao5pyykOg2YBpCWlqaDBw/2YfM/tHjxYmq7bLBZvHgxj1w5GOf7G5j+xQ4G9+vGzQM6BrqsWgmX/RIu/QDrS7BqiL74clomC2hfaToZ2FNNm3dVtVRVdwCb8IS98dHkK85nSNeWPDJvvY0iaYypM1/CPRNIEZHOIhIBjAbmVWkzFxgCICItgC7Adn8WGu6cDuGZMX05t2UTfv76SnYcOh7okowxIazGcFfVMuBuYD6wEZitqutF5DERGeFtNh84LCIbgE+ASap6uL6KDldxUW6m35KO0yGMfyWT3EK7gsYYUzs+Xeeuqh+oahdVPVdVn/DOe1hV53kfq6rer6rdVbWnqtbunVJDh8QYXrg5lV2HC7jnTbuCxhhTO/YJ1SA04JxE/nhNDz7bfJA//9+3gS7HGBOCbGyZIDUmowOb9uUx44sddGsdxw1p7WteyBhjvOzIPYg99NPzGXReIv89Zx0rvzsa6HKMMSHEwj2IuZwOnr+pH22aRnHnayvtNn3GGJ9ZuAe5pjERTL8ljaLScia+upKiUhtF0hhTMwv3EJCSFMdTo/qwbk8uv/33GruLkzGmRhbuIeLS7kk8MLQr767ew9TP7PNhxpgzs3APIb8YfC4/7dWGv374LYs32cCbxpjTs3APISLClJG96NY6nnveXMX2g/mBLskYE6Qs3ENMTISLaeNScTsd3PHqCrvJhzGmWhbuIah98xiev6kfOw8XcP9bq6mosDdYjTGnsnAPUQPPTeThK7uzcOMBnly4OdDlGGOCjA0/EMJuGdiR9XtyefbjrZzfJp4rerYJdEnGmCBhR+4hTET44zU96NuhKQ/86xu+3Xcs0CUZY4KEhXuIi3Q5efHmVGIjXdzx6gpyCkoCXZIxJghYuIeBpPgoXhyXyv7cYu5+w8aAN8ZYuIeNfh2a8fg1Pfhi6yH+YmPAG9Po2RuqYeTG9Pas25PL9C92cEG7eK7tmxzokowxAWJH7mHm91d2p3/n5jz477WsycoJdDnGmACxcA8zbqeDv4/tR4vYSO58bSUH84oDXZIxJgAs3MNQYmwkU8elcrSghF/8cyUlZfYGqzGNjYV7mOrRLoG/Xt+LzJ1HefS99YEuxxjTwOwN1TB2dZ92bNh7jKmfbueCtvGM7d8x0CUZYxqIHbmHud8M68aPu7TkkXfXs3zHkUCXY4xpIBbuYc7pEJ4Z05f2zWP4r9dXkp1jN9k2pjGwcG8EEqLd/OOWNErKKpj46goKS+wm28aEOwv3RuK8VrE8PaYPG/YeY9Lb39hNto0Jcxbujcgl3ZKYNKwr76/Zy98Xbwt0OcaYemTh3sj814/P5arebfnbgk0s2rg/0OUYY+qJhXsjIyL8z/W9uKBtPL+ctZot+/MCXZIxph5YuDdC0RFOpo1LI8rtZIKNAW9MWLJwb6TaNo1m6rhU9uYU8Yt/fk2pjQFvTFixcG/EUjs240/X9eTLbYf54/sbAl2OMcaPbPiBRm5kajKb9+cx7bPtpCTFMW6ADVFgTDiwI3fDb4d3Y0jXljw6bz1fbj0U6HKMMX7gU7iLyHAR2SQiW0XkwTO0GykiKiJp/ivR1LcTQxSc06IJ//XPr9lx6HigSzLG1FGN4S4iTuB54HKgOzBGRLpX0y4OuBdY5u8iTf2Li3Iz49Z0nA5h/MuZ5BaUBrokY0wd+HLkngFsVdXtqloCzAKurqbdH4H/AYr8WJ9pQB0SY3jx5lR2Hy3grjfsChpjQpnUNMaIiIwEhqvqBO/0OKC/qt5dqU1f4CFVvV5EFgMPqOqKatY1EZgIkJSUlDpr1qxaFZ2fn09sbGytlg02wdiXz7NKmbGuhMHtXdzaPQIR8Wm5YOxLbYRLP8D6Eqzq0pchQ4asVNUaT337crVMdX/ZJ18RRMQBPAncVtOKVHUaMA0gLS1NBw8e7MPmf2jx4sXUdtlgE4x9GQxEfPgtLyzexo96d2H8RZ19Wi4Y+1Ib4dIPsL4Eq4boiy+nZbKA9pWmk4E9labjgB7AYhHZCQwA5tmbqqFt0tCuDL+gNY//ZwMLN9gYNMaEGl/CPRNIEZHOIhIBjAbmnXhSVXNVtYWqdlLVTsBXwIjqTsuY0OFwCE+O6kOPtgncO2sV67JzA12SMeYs1BjuqloG3A3MBzYCs1V1vYg8JiIj6rtAEzjREU5m3JpG02g341/JZG+u3cXJmFDh03XuqvqBqnZR1XNV9QnvvIdVdV41bQfbUXv4aBUfxczb0zleXM7tL2WSV2SXSBoTCuwTqqZG3VrH8/zYfmw5kM9db6yySySNCQEW7sYnP+7Skieu6cFnmw/y8Lvr7DZ9xgQ5GzjM+Gx0Rgd2Hy3g+U+2kdwshruGnBfokowxp2Hhbs7Kry/rStbRQqbM30SbhCiu65cc6JKMMdWwcDdnxeEQpozszcG8Yn7z9hpaxUVxUUqLQJdljKnCzrmbsxbhcvDiuFTOaxXLz19fyfo9tb8Gfu6qbAb95WM6P/gfBv3lY+auyvZjpaa+2f4LXhbuplbio9y8fHsG8VEubnspk91HCs56HXNXZTP5nbVk5xSiQHZOIZPfWWsBESJs/wU3C3dTa60Tonh1fAYlZRXcMnM5x0rO7gqaKfM3UVhafsq8wtJypszf5M8yTT2x/RfcLNxNnZzXKo6Zt6WxN7eQJ1cUkV9c5vOye3Kq/8Tr6eab4GL7L7hZuJs6S+3YnL+P7cd3eRVMfHUFxWXlNS8EtG0afVbzTXCx/RfcLNyNX1zSLYnxPSL4ctth7pu1mvKKmk/RTBrWlWi385R50W4nk4Z1ra8yjR/Z/gtuFu7Gbwa1c/P7K7vzf+v2MfmdNTV+ivWavu3483U9adc0GgHaNY3mz9f15Jq+7RqmYFMntv+Cm13nbvxq/EWdyS0s5ZlFW4iPcvPfPz3/jHdyuqZvOwuDEGb7L3hZuBu/+9WlKRwrLGX6FzuIi3Lzy0tTAl2SMY2OhbvxOxHh4Su7k19cxpMLNxMT4eSOi88JdFnGNCoW7qZeOBzCX6/vRWFpOU98sJEot4NxAzsFuixjGg0Ld1NvnA7hyRv7UFRSzu/fXU+Ey8Go9A6BLsuYRsGuljH1KsLl4Pmx/bi4S0sefGct/16ZFeiSjGkULNxNvYtyO5k2LpULz01k0tvf8O5qG3vEmPpm4W4aRJTbyfRb0sno3JxfvbXaBpcypp5ZuJsGEx3hZOZt6fTvnMj9s1czZ5WdojGmvli4mwYVE+Fi5m3pDDgnkftnf8PsFbsDXZIxYcnC3TS46AgnM25N56LzWvCbt9fw+lffBbokY8KOhbsJiOgIJ/+4JY2fdGvFQ3PXMeOLHYEuyZiwYuFuAibK7eSFm1O5vEdr/vj+Bp5euKXGwcaMMb6xcDcBFeFy8OyYvlzfL5knF27mif9stIA3xg/sE6om4FxOB1NG9iIuysX0L3aQW1jKn6/rictpxx7G1JaFuwkKDofwyFXdSYh28/SiLRwtKOW5m/oSVeVmEMYY39ihkQkaIsKvLuvCH0ZcwKJv9zNuxjJyCkoCXZYxIcnC3QSdWy/sxDOj+/LN7lxGvriUrKMFgS7JmJBj4W6C0lW92/LKzzLYf6yI6/7+JeuycwNdkjEhxcLdBK2B5yby9s8vxOkQbpy6lI+/3R/okowJGRbuJqh1bR3H3LsG0blFEya8soJXvtwZ6JKMCQkW7iboJcVHMfvOgVzSrRWPzFvPw++uo6y8ItBlGRPUfAp3ERkuIptEZKuIPFjN8/eLyAYRWSMii0Sko/9LNY1Zk0gXU8elMfHic3h16Xfc9lImuQWlgS7LmKBVY7iLiBN4Hrgc6A6MEZHuVZqtAtJUtRfwNvA//i7UGKdD+N0V5zNlZC+W7TjMiOe/YPP+vECXZUxQ8uXIPQPYqqrbVbUEmAVcXbmBqn6iqieuV/sKSPZvmcZ874a09syaOICCknKufX4JH67bG+iSjAk6UtM4HiIyEhiuqhO80+OA/qp692naPwfsU9XHq3luIjARICkpKXXWrFm1Kjo/P5/Y2NhaLRtsrC+1d7SogmdXFbM9t4IrOru5PsWN0yF1Xq/tk+BkffEYMmTISlVNq7Ghqp7xC7gBmF5pehzw7Gna3oznyD2ypvWmpqZqbX3yySe1XjbYWF/qpqi0TCe/s0Y7/vZ9HT11qR7MK6rzOm2fBCfriwewQmvIV1X16bRMFtC+0nQysKdqIxG5FPhvYISqFvuwXmPqLNLl5E/X9uRvN/Tm611HueLpz1m67XCgyzIm4HwJ90wgRUQ6i0gEMBqYV7mBiPQFpuIJ9gP+L9OYMxuZmszcuwYRG+li7PSveGbRFsorbOhg03jVGO6qWgbcDcwHNgKzVXW9iDwmIiO8zaYAscC/RGS1iMw7zeqMqTfnt4ln3j0XMaJ3W/7fR5sZO/0r9uYWBrosYwLCpyF/VfUD4IMq8x6u9PhSP9dlTK0s3LCf5TuOALBs+xF+8r+fMjq9PfPX72dPTiFtm0YzaVhXrunbzu/bnrsqmynzN9X7dnzx0Ny1vLlsN/f1KGX85A8Y0789j1/TMyC1mMCw8dxN2Ji7KpvJ76ylsLQcAAUKS8qZuWTnyTbZOYVMfmctgF+Dt+q262s7vnho7lpe/2rXyely1ZPTFvCNhw0/YMLGlPmbTobrCdWddS8sLWfK/E31vu362I4v3ly2+6zmm/Bk4W7Cxp4c38+vZ59F27ps+2xq8pfy03x25XTzTXiycDdho23TaJ/bOh3CF1sO1fu2z6Ymf3FK9R/kOt18E54s3E3YmDSsK9FV7rnqdghu56mhFuF00LxJBDfPWMYD//qGI8frfiu/6rYd7XYyaVjXOq/7bI3p3/6s5pvwZG+omrBx4o3LqlesVDdveI/WPLNoC9M+286ijfv53RXnMzI1Ganl0e3pth2Iq2VOvGl64hy7U8SulmmELNxNWLmmb7tqA7W6eb8Z3o0Rfdry0Jx1THp7DbNX7OYPI3r4fduB8Pg1PXn8mp4sXryYbWMHB7ocEwB2WsY0at1axzP7zoH89fqebDt4nCuf/ZzXNhSTU1D3UzXGBJKFu2n0HA5hVHoHPvn1YG4e0JGPd5Xx4ymLeWnJDkrtjk8mRFm4G+OVEOPmsat78NigaHq0i+cP721g2JOf8eG6fSdGPTUmZFi4G1NF+zgHr4/vz/Rb0nA4hJ+/vpKRLy5l2XYbbdKEDgt3Y6ohIlzaPYkPf/kj/nxdT3YfKWDUtK+4ZeZy1mTlBLo8Y2pk4W7MGbicDsZkdODTSUP43RXdWJOVw4jnljD+5UwLeRPULNyN8UF0hJOJF5/L578ZwgNDu7Diu6OMeG4Jt720nMydRwJdnjE/YOFuzFmIi3Jz9yUpfPHbIUwa1pW1Wbnc8OJSbnxxKQs37KfCbhBigoSFuzG1EBfl5q4h5/HFby/h4Su7k51TyIRXV3DZk5/yxrJdFJaU17wSY+qRhbsxdRAd4eRnF3Vm8aTBPD26D1FuJ7+bs5aBf1nEX/7vW3YfKQh0iaaRsuEHjPEDt9PB1X3aMaJ3WzJ3HuWlJTuY9tk2pn62jUu6tmLsgA78uEsrnA4bmdE0DAt3Y/xIRMjo3JyMzs3Zk1PIm8t38eby3Sx6eQVtE6K4Ia09I1OTad88JtClmjBn4W5MPWnbNJpfD+3KPZeksGjjft5YvotnPt7C04u2MPCcRK5PTWZ4j9bERtqfofE/+60ypp5FuBxc3rMNl/dsQ9bRAt75Opu3V2bxwL++4aG5a7mse2tG9G7LxV1aEOly1rxCY3xg4W5MA0puFsO9P0nhnkvO4+tdR5mzKpv31+zlvW/2EBflYmj31lzeozUXpbQgym1Bb2rPwt2YABARUjs2J7Vjcx656gKWbD3Ee9/s5aMN+/j311nERrr4cdeWDO2exOCurUiIdge6ZBNiLNyNCTC308Hgrq0Y3LUVJWU9+XLbIT5ct4+FGw/wnzV7cTqE9E7NuKSbp01Kq9ha3zHKNB4W7sYEkQjX90FfUaGs2p3Dx9/uZ9HGA/zpg2/50wff0iYhih+ltOCilJZceG4iLWIjA122CUIW7sbIwl0aAAANK0lEQVQEKYdDSO3YjNSOzZg0rBt7cgr5bPNBPt18kA/X7WP2iiwAurWOY8A5iQw4J5H0Ts1ItLA3WLgbEzLaNo1mdEYHRmd0oLxCWZudy5Kth1i67TCzMnfx8pc7ATivVSxp3heF8uMVqKqdxmmELNyNCUFOh9CnfVP6tG/KXUPOo7isnLVZuSzfeYTMHUf4YO1eZmXuBuAvKz+id3JTerdvSu/kBHomJ9AqLirAPTD1zcLdmDAQ6XKS1qk5aZ2aw2CoqFC2HcznjQVfURiTxOrdOTz38RZODFrZKi6SHu0S6N4mnu5t4zm/TTwdmsfY8AhhxMLdmDDkcAgpSXH8uL2bwYN7AVBQUsb6PcdYk5XL+j25rMvO5dPNByn3Jn6U20FKqzi6JMXRJSmWlKRYzmsZR7tm0Rb6IcjC3ZhGIibCRXqn5qR3an5yXlFpOVsP5LNh7zE27ctj0748Pt9ykH9/nXWyTaTLQecWTejcogmdWjShc2ITOibG0KlFE1rGRuKw4A9KFu7GNGJRbic92iXQo13CKfNzCkrYeiCfbQfz2Xognx2HjrNpfx4fbdhPWaUbkkS6HLRvHkNys2jaN4uhXbNo2jWNPvm9RWykHfUHiIW7MeYHmsZEfH8Ov5Ky8gr25BSx4/Bxdh0pYPeRAr47fJyso4Ws2pVDbmHpKe1dDiEpPorWCVG0jo8iKT6KpPhIWsVH0iouilZxkbSIjaRpjNuu6PEzC3djjM9cTgcdEmPokFj9kMV5RaVk5xSSfbSQPblF7M0pZF9uEfuOFbFx7zE+2XSAgmruUuV2ColNIkmMjSAxNpLEJhE0r/TVLCaC746U02ZfHs1i3MRHu23snRpYuBtj/CYuyk231m66tY4/bZv84jL2HyviwLFiDuQVcSi/hEP5xRzMK+bIcc/j7QfzOXK85AcvBH9e/tnJx1FuBwnRbuKj3J7v0W7io1zERbmJ836PjXIRF+miSaSLWO9Xk0gnsZEuYiJdxLidYfuegU/hLiLDgacBJzBdVf9S5flI4FUgFTgMjFLVnf4t1ZjwNXdVNlPmb2JPTiFtm0YzaVhX/rViF0u2HTnZZtC5zbkhrcMP2gE/mLfiuyO8uWw39/UoZfzkDxjTvz2PX9PTp+1Wt75r+rbzue4T2y5XxSnyg23HRrqIbRnL2qzcGrf96FUp/KhLC44cL+HTpSvokHI+RwtK+WrbYT7dfJD9x4q9p4JiKCwtZ+uBMo4VlZJXVHbyKqCaRLudxEQ4iY448d1FtNtBTISLaLeTSLeDaLeTKLeTKLeDKNf3jyNdnucjXd7HLgeRbgcRTicRLsf3X87vv7udgmr930i9xnAXESfwPHAZkAVkisg8Vd1Qqdl44Kiqnicio4G/AqPqo2Bjws3cVdlMfmcthaWeo9TsnELue2v1D9ot2XbklLDPzilk0tvfgEKpN8iycwq5/63VVFRarlyV17/aBXBKyFa33Un/+gYESsu/X9/kd9YC/CDgq1ve39t+ZN56/nxdT67p246DiU4G92rL3FXZfPztgZPLFpVWkHW08GQ7AFWlqLSCvOJS8ovKyC/2fhWVUVBSzvGSMo4Xl3G8uJzjxWUUlJZTWFJOQUkZhaUVFJaUcTCvmELv/KLScgpLPd99fM04o3HdIxhS99WckS9H7hnAVlXdDiAis4CrgcrhfjXwqPfx28BzIiLaEC9PxoS4KfM3nQyqs3UiCCurqKYdwJvLdp8SsNVtt7Sa5CosLWfK/E0/CPfqlm+IbVe3bNV2IkK092i8VdxpiqoFVaW0XCkqK6e4tIKi0nJKyisoLq2guKyckrIKissqKCmr8MwvK6e0TCku98wrLa+gtKyCuPxd/ivqNKSm/BWRkcBwVZ3gnR4H9FfVuyu1Wedtk+Wd3uZtc6jKuiYCE72TXYFNtay7BXCoxlahwfoSfBq0HxGtz0utr3WXF+TijPn+MseSfVtX1na7lZet6/K1XLYFcOhMy1atMYjV5Xeso6q2rKmRL0fu1b3bUPUVwZc2qOo0YJoP2zxzQSIrVDWtrusJBtaX4BMu/QBPX8pyD4RNX8Jpv9R3Xxw+tMkC2leaTgb2nK6NiLiABOAIxhhjAsKXcM8EUkSks4hEAKOBeVXazANu9T4eCXxs59uNMSZwajwto6plInI3MB/PpZAzVXW9iDwGrFDVecAM4DUR2YrniH10fRaNH07tBBHrS/AJl36A9SVY1XtfanxD1RhjTOjx5bSMMcaYEGPhbowxYSjow11EokRkuYh8IyLrReQP3vmdRWSZiGwRkbe8b/YGPRFxisgqEXnfOx2q/dgpImtFZLWIrPDOay4iH3n78pGINAt0nb4QkaYi8raIfCsiG0VkYCj2RUS6evfHia9jInJfiPblV96/93Ui8qY3B0L1b+WX3n6sF5H7vPPqfZ8EfbgDxcAlqtob6AMMF5EBeIY4eFJVU4CjeIZACAW/BDZWmg7VfgAMUdU+la7XfRBY5O3LIu90KHga+FBVuwG98eyfkOuLqm7y7o8+eMZ5KgDmEGJ9EZF2wL1Amqr2wHMhx4lhTULqb0VEegB34Pmkf2/gShFJoSH2iaqGzBcQA3wN9Mfz6S6Xd/5AYH6g6/Oh/mTvjrwEeB/Ph79Crh/eWncCLarM2wS08T5uA2wKdJ0+9CMe2IH34oJQ7kuV+ocCS0KxL0A7YDfQHM8Vfe8Dw0LxbwW4Ac9giyemfw/8piH2SSgcuZ84lbEaOAB8BGwDclS1zNskC88vRLB7Cs+OPTEERyKh2Q/wfAJ5gYis9A4rAZCkqnsBvN9bBaw6350DHARe8p4umy4iTQjNvlQ2GnjT+zik+qKq2cDfgF3AXiAXWElo/q2sAy4WkUQRiQGuwPOBz3rfJyER7qparp5/NZPx/HtzfnXNGraqsyMiVwIHVLXy2Bc+DdsQpAapaj/gcuAuEbk40AXVkgvoB7ygqn2B4wT5aYuaeM9FjwD+FehaasN7/vlqoDPQFmiC5/esqqD/W1HVjXhOJ30EfAh8A5SdcSE/CYlwP0FVc4DFwACgqXeoA6h+SIRgMwgYISI7gVl4Ts08Rej1AwBV3eP9fgDPed0MYL+ItAHwfj8QuAp9lgVkqeoy7/TbeMI+FPtywuXA16q63zsdan25FNihqgdVtRR4B7iQ0P1bmaGq/VT1Yjwf8txCA+yToA93EWkpIk29j6Px7PiNwCd4hjoAz9AH7wamQt+o6mRVTVbVTnj+Zf5YVccSYv0AEJEmIhJ34jGe87vrOHUYipDoi6ruA3aLSFfvrJ/gGc465PpSyRi+PyUDodeXXcAAEYkREeH7fRJyfysAItLK+70DcB2efVPv+yToP6EqIr2AV/C8Y+4AZqvqYyJyDp4j4ObAKuBmVS0OXKW+E5HBwAOqemUo9sNb8xzvpAt4Q1WfEJFEYDbQAc8f6A2qGvQDyIlIH2A6EAFsB27H+7tG6PUlBs+bkeeoaq53XsjtF+8lz6PwnMJYBUzAc449pP5WAETkczzvr5UC96vqoobYJ0Ef7sYYY85e0J+WMcYYc/Ys3I0xJgxZuBtjTBiycDfGmDBk4W6MMWHIlxtkG9OgvJeJLfJOtgbK8QwRAJChqiUBKewMRORnwAfe6+aNCTi7FNIENRF5FMhX1b8FQS1OVS0/zXNfAHer6uqzWJ+r0lgpxviVnZYxIUVEbhXP+P6rReTvIuIQEZeI5IjIFBH5WkTmi0h/EflURLaLyBXeZSeIyBzv85tE5CEf1/u4iCwHMkTkDyKS6R2f+0XxGIVnOOq3vMtHiEhWpU9WDxCRhd7Hj4vIVBH5CM9gZS4R+X/eba8RkQkN/1M14cjC3YQM79jY1wIXegeSc/H9zdgTgAXewcxKgEfxfGz9BuCxSqvJ8C7TD7hJRPr4sN6vVTVDVZcCT6tqOtDT+9xwVX0LWA2MUs946jWdNuoLXKWq44CJeAaUywDS8QzC1qE2Px9jKrNz7iaUXIonAFd4hhwhGs9H7QEKVfUj7+O1QK6qlonIWqBTpXXMV9WjACIyF7gIz9/B6dZbwvdDLQD8REQmAVFACzxD0f7fWfbjXVUt8j4eCpwvIpVfTFLwfCTdmFqzcDehRICZqvr7U2Z6RgqsfLRcgecOXiceV/49r/omk9aw3kL1vjHlHbflOaCfqmaLyON4Qr46ZXz/n3HVNser9OkXqroIY/zITsuYULIQuFFEWoDnqppanMIYKp57psbgGTN8yVmsNxrPi8Uh76iY11d6Lg+IqzS9E8+t7qjSrqr5wC9ODGUrnvugRp9ln4z5ATtyNyFDVdd6RwtcKCIOPKPs/ZyzG9f7C+AN4FzgtRNXt/iyXlU9LCKv4Bne+DtgWaWnXwKmi0ghnvP6jwL/EJF9wPIz1DMVz8iAq72nhA7gedExpk7sUkjTaHivROmhqvcFuhZj6pudljHGmDBkR+7GGBOG7MjdGGPCkIW7McaEIQt3Y4wJQxbuxhgThizcjTEmDP1/pGenSMj5bcYAAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"%matplotlib inline\n",
"data_pred = pd.DataFrame({'Temperature': np.linspace(start=30, stop=90, num=121), 'Intercept': 1})\n",
"data_pred['Frequency'] = logmodel.predict(data_pred)\n",
"data_pred.plot(x=\"Temperature\",y=\"Frequency\",kind=\"line\",ylim=[0,1])\n",
"plt.scatter(x=data[\"Temperature\"],y=data[\"Frequency\"])\n",
"plt.grid(True)"
]
},
{
"cell_type": "markdown",
"metadata": {
"hideCode": false,
"hidePrompt": false,
"scrolled": true
},
"source": [
"This figure is very similar to the Figure 4 of Dalal *et al.* **I have managed to replicate the Figure 4 of the Dalal *et al.* article.**"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Computing and plotting uncertainty"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Following the documentation of [Seaborn](https://seaborn.pydata.org/generated/seaborn.regplot.html), I use regplot."
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYgAAAEKCAYAAAAIO8L1AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzt3Xl8VPW9+P/XObNkmewhC0vYA4RFXBAQV6LIvrhVkSoWUWsv9ddqq3bzWnurtbetxfZ+Vdy3LooLQhSXoFAVURSJ7GsgAbJvk8w+5/z+mGRISIBJyGSWvJ+Ppsk5c+bk8zHMvOezvT+Krus6QgghxAnUUBdACCFEeJIAIYQQokMSIIQQQnRIAoQQQogOSYAQQgjRIQkQQgghOhS0APGLX/yCCy64gDlz5nT4uK7r/M///A/Tpk1j7ty5bN++PVhFEUII0QVBCxBXX301zzzzzEkf37BhA8XFxXzwwQf87ne/48EHHwxWUYQQQnRB0ALE+eefT3Jy8kkfLywsZMGCBSiKwtlnn01DQwMVFRXBKo4QQohOMobqF5eXl5Odne0/zs7Opry8nMzMzFM+z2pzoQAoCoDv5+OHzT+3ekzxHbe9TkFVfAfN3/zPEUII4ROyANFRho9A3qQbbW6qqhu7vTz+YNJcDkUBVVH8P7c+pzYfq4qCqoKqKv5rz1RGRiKVldYzvk+4iub6RXPdQOoX6TIyEjv9nJAFiOzsbMrKyvzHZWVlp209BJPe/H++73rrswFTFDAoii9gNH8ZmoOHQVUwGBQMqkwcE0JEhpAFiPz8fF555RVmz57N1q1bSUxMDGmA6A66Dh5dB+3kgUWB5mChYlAVjAYVg0HBZFBRVenmEkKEj6AFiLvvvpsvv/yS2tpaLrnkEn784x/j8XgAWLhwIZdeeinr169n2rRpxMXF8fDDDwerKGFFBzyajkfztntMVUAxGWmwuTAZVMwmVVocQoiQUSIt3fexqqagjEGEi7Q0CzU1Tf5jVVWIMaqYTQZiTIaIb2VEcz9vNNcNpH6RLqLGIERgNE3H7vJid/laHGajSqzZQKzZGPHBQggR3iRARBiXR8Pl0bDa3JhNBuJjjMSYDaEulhAiCkmAiFA64HR7cbq9GFSF+FgjcTFGVFnPIYToJjICGgW8mo7V5qaqzo7N4Ql1cYQQUUICRBTRdGiwuaiqs+N0t58lJYQQnSEBIgp5NJ1aq5P6JhdaZE1SE0KEEQkQUczu9FBV75DWhBCiSyRARDmtuTVhtblCXRQhRISRANFLNDk81DQ40E6RBkQIIVqTANGLuDwaVQ0O3B7pchJCnJ4EiF5G03RqGpzYnTIdVghxahIgeiEdqG9y0dDk6nBfDiGEAAkQvZrN6aHW6pSpsEKIDkmA6OVcHo2aegdeTQt1UYQQYUYChMCj6VQ3OHF7JEgIIY6TACGA5sFrq8xwEkIcJwFC+Ok61FqdeLzSkhBCSIAQJ9B0qJOBayEEEiBEBzyaTp3VKVNghejlJECIDrXsWieE6L0kQIiTsjk9OFyy4lqI3koChDilhiaXrJEQopeSACFOSdOhvlFShQvRG0mAEKfl8mg02mU8QojeRgKECEiT3S3rI4ToZSRAiIDoILOahOhlJECIgDndXpwuScUhRG8hAUJ0itUme0gI0VtIgBCd4tF0bLIbnRC9QsQFiD+89BV7SupCXYxerdHuRtOkFSFEtDOGugCddeBIPQeO1HP28D7MumAQCXGmUBep19F1aHK4SYw3h7ooQoggirgWhNL8/dt9Vfz1ta18u7dK+sRDwOb0SMZXIaJcxAWI+24+n+y0eMD3JvXax/t45YM9WG2y2rcn6TrYHDIWIUQ0i7gAMbR/Mv919VimTcjBoPraEzsP1fLX17fy7T5pTfQkm8MtrQgholjEBQgAg6oy9dz+LLt6HP0zLADYnV5eW7ePfxXulU+2PUTTwS4zmoSIWkENEBs2bGD69OlMmzaNFStWtHv86NGj3HTTTSxYsIC5c+eyfv36Tt0/Ky2eH84fy/SJx1sT3x2o4fGVW9lbKjOdekKTwyOtNiGiVNAChNfr5aGHHuKZZ56hoKCANWvWsG/fvjbXPPHEE8ycOZO3336bxx57jN/+9red/j0GVeHSs/vzo6vG+scmGmxunn93FwUbiyV/UJBpmo7dKaurhYhGQQsQRUVFDBo0iJycHMxmM7Nnz6awsLDNNYqi0NjYCIDVaiUzM7PLv69vuoUfXTWWi8/q65/p9Nl3ZTy5ajtVdfYu31ecns0hOZqEiEZBWwdRXl5Odna2/zgrK4uioqI21yxbtoxbb72VV155BbvdzvPPPx/QvdPSLCd9bNGs0UwYk81zq3dQ3+jkaFUT//fWNhZOH8nksX27Vpkedqr6havE5FhizYH9c8rISAxyaUInmusGUr/eJmgBoqN+aUVR2hwXFBRw1VVXsWTJErZs2cK9997LmjVrUNVTN2xqappO+XhGYgzLrh7Lm+sPsPNQLU63lxfW7GDb3krmXjgEkzF8x+bT0iynrV84arI6SE2MOe11GRmJVFZae6BEPS+a6wZSv0jXleAXtHfK7OxsysrK/Mfl5eXtupBWrlzJzJkzATjnnHNwOp3U1tZ2y++3xJr4/pUjmDtlsH8Ae/PuSp54e5t0OQWB0+2VrUmFiDJBCxDjxo2juLiYkpISXC4XBQUF5Ofnt7mmb9++bNy4EYD9+/fjdDpJS0vrtjIoisIFY7P54fwx/k+3ZTU2/v7Wd2w7WNNtv0f4yPRiIaJL0AKE0WjkgQceYOnSpcyaNYuZM2eSm5vL8uXL/YPV999/P6+99hrz5s3j7rvv5g9/+EO7bqju0D8jgWVXj2P04FQAXG6Nf3y4h/e/PCxJ57qR3SlTXoWIJooeYa/oY1VNVFU3dum5uq7z6XfHeH/TYVriwvD+yVx/+XAsseGR9C9SxyBaJFvMxMWcfGgrmvt5o7luIPWLdGE1BhGOFEXh4rP68YPZeVhifW9i+47U8//e2kZZjS3EpYsOsrJaiOjRqwJEi2H9kvmvq8cxoDlNR63VyZNvb2NHsYxLnCmXR5PFiUJEiV4ZIABSEmK4be4YzsntA/je2F75YA/rvimVfvQz5JB9q4WICr02QACYjCrXXjaMmZMH0jI2/tHmUl77eB9uj3wK7irpZhIiOvTqAAHHxyUWzxhFjMkAwNZ91TxbsINGu6SQ6AqvpuNySytCiEjX6wNEixE5KfxwwfH1EofLG3ni7W1U1Mqiuq6wSzeTEBFPAkQrWanx3LlgLAOzEoDmwetV2zh4rCHEJYs8TpesiRAi0kmAOEFCnIlbZ4/mrGHpgG/A9bmCnWzdVxXikkUWTfel3xBCRC4JEB0wGVW+lz+cS8/uB/j61P+9bh//2Xo0xCWLLLJPhBCRTQLESaiKwvSJA5l/0RD/DKf3Nh3m3S8OyT7MAXK5vZLKRIgIJgHiNCaNzuL7V47EZPD9p/q06BgrP94vmUsDoAMOl0x5FSJSSYAIQN6gVJbMziMuxjcN9tt9Vby0drdM5QyATdZECBGxJEAEaFB2IrfPG0OyxQzA3tJ6nnt3pywKOw2PV8ftkUAqRCSSANEJWanx3DF/DBkpsYBvrcTTq3dgtblCXLLwJvtECBGZJEB0UkpCDLfPG0P/Pr5Ef2U1Nla8s4NaqzPEJQtfDpdXBvaFiEASILrAEmvi1jl5DOnry69e3eBgxTvbqaqXVdcd0QGHdMUJEXEkQHRRrNnILTPzGDUwBYD6JhdPv7ODctlXokMyWC1E5JEAcQZMRpUbp41g7FDfPtpWu5unV+/gaFXk7ggXLDJYLUTkkQBxhowGlevzc/37SticHp5Zs4PSiq5tixrNbLKyWoiIIgGiGxhUhWsuG8bEvEzANyj7bMFODpdH7/62XSEJ/ISILBIguomqKMy/aAgXjMkGfInqnn93F4fKJEi0kAR+QkQWCRDdSFEU5kwZxIVjWwWJ93ZKkGhFFhYKETkkQHQzRVGYdcEgLj6rLwAutyZBohWHjEMIETEkQASBoijMmDSQS8YfDxIvvLdLxiQATdelm0mICCEBIkiU5nThLS2JljGJkgoJEg7ZjlSIiCABIohaWhIXjWsbJI708nUSTkkBLkREkAARZIqiMHPyQKY0D1y3bGFa1otXXGs6kipdiAggAaIHKIrC7AsGMWl0FuCbyfPsmh1U1Pbe3E3SzSRE+JMA0UMURWHuhYM5b2QGAE0OD88W7KC63hHikoWGQ1oQQoQ9CRA9SFUUrrp4KGcP96XlsNrcPFuwg7rG3pcqXNN06WYSIsxJgOhhanNajjFDfAn+6hpdPFuws1duOiStCCHCmwSIEDCoCtfnD2dEji9VeHW9g+cKdmJzuENcsp4l4xBChDcJECFiNKgsmjaCIX2TACivtfPCe7t61cY60s0kRHiTABFCJqPKzdNHkpOZAEBpZRNPvFmE26OFuGQ9R7qZhAhfAQWIW2+9lY8//rjTqZo3bNjA9OnTmTZtGitWrOjwmnfffZdZs2Yxe/Zs7rnnnk7dPxrEmA0snjGKrNQ4AHYfquVfhXvxar0jSEg3kxDhK6AAcf311/Piiy9yxRVXsGLFCmpra0/7HK/Xy0MPPcQzzzxDQUEBa9asYd++fW2uKS4uZsWKFfzzn/+koKCAX/7yl12rRYSLjzXyg9l5pCXFALDzUC1vrj+A1gv2TtA02WlOiHAVUIC48soreeGFF3j66aepqKhgzpw53HvvvWzbtu2kzykqKmLQoEHk5ORgNpuZPXs2hYWFba557bXXWLRoEcnJyQCkp6efQVUiW1K8mVtn55Gc4AsSW/ZW8e7GQ71igx27tCKECEvGrjzJZDIRExPDfffdx8UXX8z999/f7pry8nKys7P9x1lZWRQVFbW5pri4GIAbbrgBTdNYtmwZl1xyyWl/f1qapSvFDntpaRb+vxvO5s+vfE2Tw8Pn28rokxbPrClDQl20bnXi38+gKmSkR8ffNCMjMdRFCCqpX+8SUID44IMPeOWVV6iurubGG2+koKAAi8WCx+Phyiuv7DBAdPTJV1GUNsder5dDhw7x8ssvU1ZWxqJFi1izZg1JSUmnLE9NTfQmu+vXJ4GbZ4zk2TU7cXk03tlwADSNyaOzT//kCJCWZunw76e53JiMhhCUqPtkZCRSWRm92XqlfpGtK8EvoACxcuVKbrvtNi6++OK2TzYa+fWvf93hc7KzsykrK/Mfl5eXk5mZ2eaarKwszj77bEwmEzk5OQwZMoTi4mLOOuusztYjquRkJrLoyhG8tHY3Xk1n9afFxMeYOGtY9HbB2V3eiA8QQkSbgMYgnnrqqXbBoUV+fn6H58eNG0dxcTElJSW4XC4KCgraXXvFFVewadMmAGpqaiguLiYnJ6cz5Y9auQNS+F7+cBRAB17/eB/7jtSHulhBI7OZhAg/AQWIG2+8kfr6429OdXV1LFq06JTPMRqNPPDAAyxdupRZs2Yxc+ZMcnNzWb58uX+w+uKLLyYlJYVZs2axePFi7r33XlJTU8+gOtFl3NB05l40GACvpvPKB7s5UtkY2kIFiabJTnNChBtFD2CazPz581m1atVpz/WEY1VNVFVH55skdNxHX/h1KYVflwJgiTVyx/wx9EmOC0XxztjJxiAA4mKMJFvMPVyi7tMb+rClfpGrK2MQAbUgNE3DZju+wU1TUxNer3za6yn55/b37yXR5PDw/Lu7ojK5n9Pl6RXTeoWIFAEFiDlz5rBkyRJWrVrFqlWruPXWW5k3b16wyyaaKYrC3CmDGducAbbW6vTlbYqyrTs1HelmEiKMBDSL6Y477iAzM5N169ah6zo33HADCxYsCHbZRCuqqnDd1OE0OXZx8FgDx6ptvPLBHm6ZOQqjIXpSajlcXmLNXVqeI4ToZgGNQYST3jgG0ZrD5eHp1Ts4Vu3r8hs3NI3rL89FPWGNSbg6Xf0UICM1LmLq01pv6MOW+kWuoK2DqK6u5uWXX6akpASP53i3xvLlyzv9C8WZiTUbWTxzFE++vY26RhffHaghMf4Qsy8Y1G4hYiTSAafLS1yMtCKECLWAXoU//vGPGTZsGBdccAEGgyxmCrWkeDNLZuXx5Krt2Jy+lBxJFjOXjO8X6qJ1C4cECCHCQkCvwoaGBn73u98FuyyiE/qkxLF45kieWb0Tt1dj7abDJMabOCc3I9RFO2NOtxevpmFQo2dsRYhIFNArMDc3l/Ly8mCXRXRSTmYiC6/IRW3uWXrjkwPsLa0LbaG6id0ps5mECLWAWxDz5s3jnHPOISYmxn9exiBCb9SgVBZcPJQ3N/j2j3j1wz3cPncM/fpEdnZUm9NDQpwp1MUQolcLKEDMmTOHOXPmBLssoosmjMqkwebio82luNwaL763izvmjyEtKTbUResyTdNxurzEmGXMS4hQCShAXHXVVcEuhzhDU8/pT32ji692VWC1u3mhOUhYYiP3U7jN6ZEAIUQIBTQGUVxczMKFC/3ZWLdv387f/va3oBZMdI6iKMy7aAh5g3zJDqvqHby0djeuCN7Os2WwWggRGgEFiAcffJA777yTxETfQou8vDzWrl0b1IKJzjOoCtdfPpyBWQkAlFQ08u/CfXi1iFoL2YYMVgsROgEFCKvVyiWXXOJfiKWqKiZT5HZdRDOz0cBN00eSnuwbf9h5qJbVnx2M2CR4Nmd05ZsSIpIEFCAMBgNut9sfIMrLy1FljnrYssSa+MHMUViaZwF9ubOCT7YcDXGpukbT9KhLSihEpAh4w6Bly5ZRW1vL3/72N2688UaWLFkS7LKJM5CWFMstM0ZiNvr+xB9uLuGbPZUhLlXX2BwSIIQIhYBmMS1YsIABAwbw8ccfY7fbefTRR5kwYUKwyybOUP+MBG6cNoKX1u5C0+HN9QdIjDeROyAl1EXrFJdHw+3RMBml1SpETwo44c2ECRMkKESgETkpXHXJUN5Yf3wh3W1zx9A/whbS2Zweko2Ru9ucEJEooABxzTXXdJgpdOXKld1eINH9zhuZSX3T8YV0L0XgQjqH00NinAlVjfyMtUJEioACxH333ef/2el0UlBQQGZmZtAKJbpfpC+k05H0G0L0tIACxMSJE9scX3TRRTJIHWFaFtJZbW52Ha6lqt7By+/vZsnsPMzGyFitLAFCiJ7VpVG/xsZGSkpKurssIsgMqsINVwwnJ9O3kO5weWQtpJMpr0L0rE6PQWiaRmlpKT/4wQ+CWjARHGajgZtnjOTJVduprnf4F9LNv2hIROxIZ3N4ZM9qIXpIp8cgDAYDAwYMICsrK2iFEsHVspDuyVXbabS7+XJnBUkWM/nnDgh10U7L5dHweDWMBpnyKkSwdWkMQkS+tKRYFs8cxdOrt+Nya3y0uZSkeDMTRoX/5AObw0OSRaa8ChFsAQWIyZMnd9j9oOs6iqKwcePGbi+YCL7+fSwsmjaCF9/bjabrvP2fAyTEmxg1MDXURTslu8tDQrwJNQK6xISIZAEFiIULF1JXV8f111+Pruu88cYbZGVlMWvWrGCXTwRZ7oAUrrlsKK9/vB9Nh39+tJelc/LIyUwMddFOStfB4fQSHytjEUIEU0AduV999RX//d//zahRo8jLy+PXv/4169evp3///vTv3z/YZRRBdk5uBjMmDgTA7dF48b3dVNbZQ1yqU7M53aEughBRL6AAUVFRQU1Njf+4pqaGysrITPwmOnbx+L5MGZsN+NYbPP/uThqaXCEu1cl5vDout+wVIUQwBdRGX7x4MfPnz2fq1KkArF+/njvuuCOoBRM9S1EUZl0wCKvNzXcHqqlrdPHi2l3cNnd02E4rtTk9mE2RschPiEgU0Ct/0aJFnHfeeXz11Vfous6iRYsYOXJksMsmepiqKFw3dRhNDjcHjjZwrNrGy+/v4ZaZo8Iyk6rT5duS1CB7kwgRFAG/sgYMGMC5557LzTffLMEhihkNKt+/cgR90+MBOHisgdfW7UMLw9XWOrIlqRDBFFCAWL9+PbNnz+bHP/4xAN999x0//OEPg1owETqxZiO3zBxFWmIMANuLa3gnTLcttTncYVkuIaJBQAHi8ccfZ+XKlSQlJQEwbtw4Dh8+HNSCidBKjDfzg9l5bbYtLfy6NMSlak/TweGSVoQQwRBwF1NGRkabY7NZVrJGu/SkWH4wcxQxzQPB6745wsZtZSEuVXuyJakQwRFQgLBYLFRVVflXU2/atInExNMvpNqwYQPTp09n2rRprFix4qTXrV27lpEjR/Ldd98FWGzRU/r1sfD96SMwNG/Us/rzYr7dVxXiUrXl9moy5VWIIAhoFtM999zDbbfdRmlpKTfddBPFxcU88cQTp3yO1+vloYce4vnnnycrK4trr72W/Px8hg8f3ua6xsZGXn75ZcaPH9/1WoigGtYvmesvz+WfH+1B12Hlx/uJMxsYGUYpOWTKqxDdL6AWxPjx43nppZf405/+xNKlSykoKGDs2LGnfE5RURGDBg0iJycHs9nM7NmzKSwsbHfd8uXLWbp0KTExMV2rgegRY4ekseDioQBous4/PtzLoTJriEt1XMuUVyFE9zltC8Lr9fK9732PN954g0svvTTgG5eXl5Odne0/zsrKoqioqM01O3bsoKysjKlTp/Lcc88FfO+0NEvA10aicK3f9ClDQFV465P9uL0aL72/m3tuPJcBWZ3L2xSs+sXFmUhOCO0HjYyM8M1h1R2kfr3LaQOEwWAgNTUVp9PZqU/5HU09bJ0RVtM0HnnkER555JGA79mipqap08+JFGlplrCu34TcPlTV2PhP0THsTg+P/WsLd8wbTZ/kuICeH8z61QL2lNiQLZzLyEiksjJ8WlXdTeoX2boS/AIagxg8eDCLFi1i+vTpxMfH+88vWrTopM/Jzs6mrOz4jJfy8nIyM4/vNdDU1MSePXu4+eabAaisrOTOO+/kiSeeYNy4cZ2uiOgZiqIwY9JA7C4vm3dV0GR381zBTu6YNybkn951oMkue0UI0V0CChBNTU3k5uZy4MCBgG88btw4iouLKSkpISsri4KCAv785z/7H09MTGTTpk3+45tuuol7771XgkMEUBSFBRcNweHysO1ADXWNLp4t2Mnt88aQ0LxuIlTsTg+WOKOk3xCiG5wyQPzhD3/g/vvv55FHHuGzzz7jwgsvDPzGRiMPPPAAS5cuxev1cs0115Cbm8vy5csZO3Ysl19++RkXXoSOqip8b+pwnK7d7C2tp6rewfPv7mTpnNHExYQuuZ8ONNo9JEsrQogzpuinyFNw1VVX8dZbb7X7OZSOVTVRVd0Y6mIETbiPQZzI5fHywru7KG6e0ZSTmcCSWXnEmDuectoT9VOA9OTYHt+3ujf0YUv9IldXxiBO+QpqHTsk343oiNlo4OYZI+mf4ZuZVFLRyMsf7MbtCd2UU99YhGwoJMSZOmWAcLlc7N+/n3379rX5ueVLCPAl9/vBzFFkpfpmMh042sCrH+7G4w1dkLC7vLg9srpaiDNxyi6m/Pz8kz9RUTpc+BZs0sUUvqw2FytW76C63gFA3qBUbpyW22bAuCfrZzKopCfH9sjvgt7RRSH1i1zdPs113bp1XS6M6H0S483cOjuPp1fvoNbqZOehWv69bh/X5+f6czn1JLdXw+bwEB8bnjviCRHuZC6g6FYpCTHcOjvPP4to24EaVn4Sug2HGu2usNzsSIhIIB+tRLdLS4rl1jl5PP3ODqx2N1v3VaMqCtdcOqzL99xbWsfmXRXUWp2kJsYwYVQmuQNSTvs8TQer3S3TXiPItoPVfFp0jMo6OxkpcVx0Vl/GDkkPdbF6JWlBiKDokxzHrXNG+xfObdlbxZsbDqB1YTbc3tI63v+yhOoGJ5oO1Q1O3v+yhL2ldQE93+70hHRWlQjctoPVvLH+AOW1djQdymvtvLH+ANsOVoe6aL2SBAgRNJmpcdw6Jw9L8xjAN3sqeeW9nZ0OEpt3VXTqfEdsTtlUKBJ8WnSsU+dFcEmAEEGVlRrPrXNG+weKPy86xpvrD3RqXKDW6uzU+Y44nB4Zi4gAlXX2k5x39HBJBEiAED0gOy2epa2CxDd7Knlj/f6A37BTEztOAniy8x3RkVZEJMhI6TgrcEZKz01XFsdJgBA9oiVIJMYfH5NY+cl+vAEEiQmjMjt1/mRsTo9kBAhzF53Vt1PnRXBJgBA9Jjstnp8uPBdL88D1t/uq+Pe6vafdCS53QArTJ+aQnhSDqkB6UgzTJ+YENIupNU3TcbhkdXU4GzsknWsuHUpWahyqopCVGsc1lw6VWUwhItNcRY/ql5HAbXNG8+wa3xTYbQdq8Hj2svCKXEzGk39eyR2Q0umA0BGbwxPSbLPi9MYOSZeAECakBSF6XGZqHLfNG+1fm7DrcC2vfLAbVw/kTnJ7NcnRJESAJECIkOiTHMft80b7B5r3ltbzwru7cLiCP5Bsc8hgtRCBkAAhQiY1MZbb542hT3NCveIyK8+s2UljkFN1O1zeLi3YE6K3kQAhQirZYua2uaPpm+7b6/xoVRNPr95OfWPgaxw6S8e3LkIIcWoSIETIJcabWTpnNAOzEgDfoqin3tlO1UkWTXUHWRMhxOlJgBBhIS7GyJJZeeQOSAagrtHFk+9sp7QiOHt/eLw6LrcMVgtxKhIgRNgwmwzcNH0k44amAb7B5GfW7Ag4KV9n2aUVIcQpSYAQYcVoULk+P5fJo7MAcHk0Xlq7m2/3VXX773K4vJKfSYhTkAAhwo6qKsy9cDBXTBgAgFfTeW3dPtZ/e6RbU2Xo0CPTaoWIVBIgRFhSFIX8cwdw1cVDaNmt9P0vS3jns+Ju/dQvayKEODkJECKsnZ+Xxfenj/Sn4di0o5xXPtiDs5sGmD2aTpMjuOsuhIhUEiBE2Bs1MJXb5o72J/nbdbiWFe9031qJRpsbj1d2nBPiRBIgREQYkJHAnfPH+PcLOFZt44m3t3G0qumM760DDU2uM76PENFGAoSIGGlJsfxw/hiG9/etlWiwuXnqne1sO1hzxvd2eTSZ9irECSRAiIgSF2Nk8cyRnN+8WZDbo/GPD/dQ+HXpGedXstpcMu1ViFYkQIiIY1BVFlw8hFmTB6E0z3Aq/LqUf32bgCyaAAAfmklEQVS094xWR2s61Fgdp93ASIjeQgKEiEiKonDRWX1ZPGMUsWYDANsO1vDkqu1UN3R9g3uPV6emwSmD1kIgAUJEuBE5KfxowVh/yvCyGhv/9+Z37CnpenoOr6ZTY3Xi9kiQEL2bBAgR8fqkxHHngrGMGpgK+FJovPjeLtZ90/VxCU3TqbU6JEiIXk0ChIgKcTFGvj99BJefNwAF39TVjzaX8vLa3di6uBBO05EgIXo1CRAiaqiKwuXnDeCmGSP94xK7S+r4+5vfUVJh7dI9JUiI3iyoAWLDhg1Mnz6dadOmsWLFinaPP//888yaNYu5c+eyePFijhw5EsziiF5i1MBUll09jv59LIBvb4kV7+zgP0VHu9TldDxIyP4RoncJWoDwer089NBDPPPMMxQUFLBmzRr27dvX5pq8vDzeeOMNVq9ezfTp0/nf//3fYBVH9DJpSb79ric1pw33ajrvfXGYl9bu6tKe15oONQ3OLndXCRGJghYgioqKGDRoEDk5OZjNZmbPnk1hYWGbayZPnkxcnC91wtlnn01ZWVmwiiN6IZNRZf5FQ7g+fzgxJl+X056Sev62sqhLmxDp+FZv1zU6z3hRnhCRwBisG5eXl5Odne0/zsrKoqio6KTXr1y5kksuuSSge6elWc64fOFM6te9pk60MHZEJs+u2kbxsQasdjfPv7uL/Ak5XHXZMExGQ+dvalBJS47FYGj7GSsjI7GbSh2epH69S9ACREcbuygty15PsGrVKrZt28Yrr7wS0L1ras48QVu4SkuzSP2CwAAsmTWKjzaX8J+tx9CBdZtL2La/iuvzh9M3vfNBq6q6kdSEGH8q8oyMRCoruzYYHgmkfpGtK8EvaF1M2dnZbbqMysvLyczMbHfd559/zpNPPskTTzyB2WwOVnGEwGhQmTFpEEvm5JFs8f1bq6i18//e2sbH3xzB28k8TJqmU2N1dNveFEKEm6AFiHHjxlFcXExJSQkul4uCggLy8/PbXLNjxw4eeOABnnjiCdLT04NVFCHaGNYvmbuuPYuzhvn+zXk1nQ83l/Dkqm2U19g6dS9dhzqrU3amE1FJ0btzk98TrF+/nocffhiv18s111zDnXfeyfLlyxk7diyXX345t9xyC3v27CEjIwOAvn378uSTT57ynseqmqiqbgxWkUNOuph61tZ9VbzzWbE/1bdBVZh6bn8uGd8Po6Fzn59y+qXgsjtP2pUa6XpDF0y016+zghoggkECRGQLx/pt3V/FexsP0WA7PoU1NTGGtMQYnG4vqYkxTBiVSe6AlFPeJy3NQmODneQEMwa1c8Fl28FqPi06RmWdnYyUOC46qy9jh4RHq3rNxmI+2XKEJocHS6yRy87pz5wLBoe6WN1OAkR7QRukFiIS7C2t4z9bj2GJM6GqKg1NzuaFcU5qrU4ssUY8Xp33vywBOG2QcHk0quocxJoNxMea/APYp7LtYDVvrD/gPy6vtfuPQx0k1mwsZs1nxYBvkkmjze0/jsYgIdqSVBuiV9u8qwLwvfnFxxrJSI1HVY93ETU5PFTU2bE7PXy1szyge+qA3eWlusFBTYMDu9PT4ay+Fp8WHevU+Z70yZaOsxuc7LyILtKCEL1ardXZ5tigKhhU8IUIBa+mN2d29Q1EV9TZyWzeFzsQLo+Gy+PCaleIjzESF2No1/1UWWfv8LmVdV3f16K7nGzVeVMXVqOLyCMtCNGrpSbGtDtnUFVMRgMZqXEkxJn8551uL4+/XsS7XxzC4ercrCVN02m0u6mqc1DX6Gyz813GSQJORkpsp35HMLSuf2uWk5wX0UUChOjVJoxqvzYnPtaIJdaIqigkWcxkpMRhNvleKpqu82nRMf787618tbO803tY6/j2q6ixOqmqt2NzeLhwXHaH1150Vt9O16e7XXZO/06dF9HF8OCDDz4Y6kJ0RqPNjc3uCnUxgiYuzow9ipvv4Va/9KRYUhNjqLM6cbq8pCXFkH/uAEYOTPWfy0iJZfrEHEYPTqOkohGHy4vbo7HrcB3bD9aQlhRLenJsp+um6b5WSWKciczUOOqsThwuL5mpccyYNDDkA9Tg27EPBY5WN+HxaljiTFw5cWBUDlBbLDHYbNH73mKxtG8tn45Mcw0z4TgNtDtFev3cHo0NW4+yYevRNntEDO+fzHVXjCAxpgt5nVoxqAoxZgOxJgNm05ndq7v1hmmg0V6/zpIAEWYi/Q30dKKlfg1NLj7aXMLXuytp/QIaMySNaefndGog+2QMqkKs2UCs2RjQdNlg6w1voNFev86SWUxCdEGSxczVlw7jgrHZvP/lYfaU1AOw/WANO4prOHt4H6ae058+ZxAovJpOk8NDk8ODUVWIjTESYzKERbAQvYMECCHOQN90C7fMzOPgsQYKvznCgSP16Dps2VvFt/uqOHt4Hy47p/9JZyoFytM8C6rR7kZVwGT0BQqTQcVsUqM2vYcIrYjrYrI53JRXWNF0HV33TR/Udd9cdQ3QNZ2IqtAJoqUL5mSiuX6pqfFs/PYIH31dytGq43VUgDFD07js7P7069P9e2EogNlkIMakYjYZOp1DKlC9oQsm2uvXWRHXgoiPNZ10bnYLrTlg6Hrbn0GnZVairjcHEp3m78cDi64f389C03R/MIrkwCOCT1EURg1KZeTAFHYdqqXw61KOVtvQgW0Hath2oIYROclcdFY/hvVL6rZP/Tq+2VC+tONuVFUhxqhiMqoYDSpGo4oqLQzRBREXIAKhKgqqoftfELquNwec4z+3BKE2rZnmABMtrRrROYqikDc4jVGDUtldUscnW45wuNw3sWJPST17Surpmx7PRWf1ZdzQ9G7/xK9pOnaXF7vr+GI8o6pgMhkwtwocQpxOxHUxARHZDGwTULTWgcUXUDRdR9d0UtMsVFU1tgs+0SKau5hOVjdd1zl4zMr6b4+wt7S+zWOJcSYmjs5iYl4mifE9t2GWquBvXZgMKkaDgsFw6pZGb+iCifb6dVZUtiDCkaIoGFpefKeY3p6eHIfWQRqH1l1lHbZemgOMt/lLk1ZL2FAUhaH9khjaL4myGhufFh1j674qvJqO1e6m8OtSPtlyhLFD05iYl8Xg7MSgDzprekueKK3NeV8uKl+wMKgKJoOvxdE6gaHoPaQFEWa681OMV9Pwen0BQ9eh5T2n9diM/7um4+2B1kpvbEF0pKHJxaad5Xy5s6Jd4ruMlDgm5mVydm4fLLHhkfPIoCpkZyVRX2fDoCq+1odBiarZU9KCaE8CRJgJ9T9STdPxahoeb6uWSOtxlVY/n+5fTstbR+vLJEC05fFqfLe/mi92lFNS0XYBqEH1DXpPGJnB8AEpGEL8Kb6j+rUEC4OqoCi+1pKqKBgMiv+xSBHq116wSReTOGOqqqCqBkwB/stoGUdRUGj+H0CbT5atWzKJ8WZsjQ48Xg1NP359i+MzyaJr7OVkjAaVc0ZkcM6IDI5WNfHlznK+3VeFy63h1XS2H6xh+8EaEuJMjB+Wztm5fejXxxI2n9x9XZrekz6ugL/LymhoaXm0bX20fPBQUKQrK8xICyLM9IZPMYHW7/isMd/MMa318YlrYPTjU5hDpbtaR063l20Hqvl6dyXFZe3/W/VJjuWsYemMG5ZOVmr8Gf++QHV3609VlXYfBBSlZRzEN+6hNh+3tExU1fe8zm7pGoje8NrrLGlBiLDVMrDfmV6K1uMq/rGWDoLKic/xhtGgfozJwHkjMzlvZCZVdXb/quyWzY2q6h2s++YI6745QnZaPGOGpDFmSBpZqXFh07IIREep0nUdPF4dj/fkrRLwtUzU5gF1RfG1PFqCi6oo/haKqhzv+hKdJy2IMNMbPsWEa/1axl80zdd1ounNs8K8Gpqm4zlNEyWY4yuarnO43Mq3e6vYdrAGm6P9TLf0pFjyBqeSNyiVgVmJ3T5mEcnjR0rz/yn4goVC+4Wv6ekJ1NQ0+rpLOT6pQ2kJMs0/t7mPAgqK/1q1+XG1+RcqSvO6rDDoOpNB6igQzm+g3SGS6+fxajQ5PDicng5bGz31BurVNA4cbaBofzU7imuxO9sHi7gYIyNzUhiRk0JuTnK3zIaK5AARiGDWTwEU1dci9o3zKc2tneaWsnq8FdQScDrT6mm9iPd4l+vxljPA8MGd319EupiECJDRoJJsMZMQZ8Tm8GBzekIykG5QVXIHpJA7IIUFF2scPGb1Z5G12nxTZu1OD9/u83VNKUD/DAvDB6QwvH8SA7MSI2p2UTTQ8WVU0NDh1L1nfv5WT3OLp3Xw8K+HIrBMDV1tTUqAEKKTDKpKYrwZS5wJh9OXjtsbohFyg6oyvH8yw/snM/fCwRytamLXoVp2HarlaLUN8L05lVY2UVrZxCdbjmAyqgzOTvQv3uvXJyHkU2hFe/5ccW0+hfTsvzMJEEJ0kaooxMeaiI814dU0UpJj8ThcuD0aTo/W6f2qu6M8AzISGJCRwBUTcmiwudhbUseekjr2HanH7vR9dHV7NPaW1vvTfpiNKgOzEhmUncjgvonkZCSE3W52IjQkQAjRDQyqSqzZSHyrvn6PV/NlWW3ew7qn2xhJ8Wb/bChN0zla3cS+0nr2HanncLkVj9dXIpdHY98R33nwDbBmpcUzMCuRnMwEBmQm0Cc5todLL8KBBAghgqRlUZgl1oSm6TjdXmwOD26vdvondzNVPd66uOyc/ni8GiUVjRw42kBxWQOHyxv9e2xrOhyrtnGs2samHeWAb+rt4L5JZKbE0q+PhX59LKQnx0oa8SgnAUKIHqCqCnExRuJijDjdXprs7naJ8nqS0aAypG8SQ/omAb6ZUceqbBSXWTlcYaWkvJH6Jpf/eqfby+7Dtew+fPweJqNKdlo8fdPjyU6LJystnqzUeOJj5W0lWshfUogeFmMyEGMy4PZo2F0eHC5vj49XnMigqgxo7k6CvgDUNzqbB7cbKa1s5GiVrc2UWrfH1wo5MYdUYryJzNQ4MlLiyEzxfe+TEkdSvEkWrEUYCRBChIjJqGIymkmKB5fbi8PlxeHyhDxlSIvkhBiSE2IYMyQN8G2puv9QDUeqmjha1URZtY1jNTYaWrU0AKw2N1abm/1HGtqcN5tU+iTHkZ4UQ3pSLOnJsaQlxZKWGEOixSzdVWFIAoQQYcBsMmA2GUiymP0D2w536FsWrSmK4ntDT4pl3NDji66aHG7Ka2yU19gpq7FRUWunos7mnzXVwuXWONocXE5kNCikJMSQmhjT5ntygpmUhBiSLKag5F8SpyYBQogw09IFlYSvG8fp9uJyh2YmVCAssSaG9ktmaL9k/zld12m0u6mss1NZ56Cq3k5VnYOqBge1DU60E1YYerw6VfUOquodHf4OBUiIM5GUYCYp3kySxfc9Md7U/OX72RJrCou0FtFCAoQQYczUvIc0cSY0Xcft1nB5fMHC7dXCNiW6oijNb9rmNoEDfHmu6hqdVNc7qLU6qWlwUNPgpK7RSa3Via2D1CE6YLW7sdrdHOHk6TAUID7WSEKcCUucyfc91oQlzogl1kR8rJH4WN/PcTFG4mOMvv++okMSIISIEKqiEGM2EGM+vojN49VwuX3Bwu32hlVW2pMxqIpvDCKp47UVTpeX2kYn9Y1O6hpd1Dc6qW9y+b8amlz+Kbkn0oEmh291O7X2gMpjMqrExRhJjDdjMirEmY3ExRiIMxuJMRuINRuJNRuav1rOGfwtPZPp1Ht5RzIJEEJEsJa1Fq15vC0bNGn+XQEjaa/yGLOB7DTf1NmO6LpvTUmDzY21yYXV7qbR5sZqc9Fod7f5arJ72nVnncjt0XB7XO0G2zvDbFKJMRowmw3EGFX/mJLZqGI2qZiMLT/7vptafxlUjCf+3LwXeEcbLPWkoAaIDRs28Pvf/x5N07juuuu4/fbb2zzucrm499572b59OykpKTz22GMMGDAgmEUSImptO1jNp0XHqKyzk5ESx0Vn9WXskHSeLdjBVzsrcHs1TAaVCaMymDAqk8++O0ZlnYP0pFgm5mUyIieVnYdq+HJnBTUNDlITY5gwKpPSyka+3FGOzeUl3mxg4ugspp7T8et0b2kdm3dVUGt1+p+fOyDlpOc7c4+Pt5T6yuH0EB9jDKgcNQ0Oki1mRg1KxeH2sv1ADQ1NLmLMBjJS4og1+xIv2l0e3B4Nq82F3enF6Q4wo14zl9vXkuOE/cW7U8uOfAaDirFlq1fD8f3BW2//alBV/7avBoOKyaAwOjez078zaOm+vV4v06dP5/nnnycrK4trr72Wv/zlLwwfPtx/zauvvsru3bt56KGHKCgo4MMPP+Svf/3rae8dqemiAxHJ6bADEc31C2Xdth2s5o31B9qdt8Qa2Vlc2+acjq+fPiMlrs3580Zm8PXuyuPX6b6xgia7u3njHcWfOG7axBwuPzfHn2Za12F3SS1rN5W0K8OYIalsP1jb7vz0iTntgsTe0jre/7L9PbLT4vhuf3W785ed279dkOjoHg6XBwWIMbf9TNy6DK3Tfe8uqeX9TSX+LXVbdjQcOzSNFEsMDpcviLRMIHA0p1NpOdd6ckFLSpNQW/3n+Z1+TtBaEEVFRQwaNIicnBwAZs+eTWFhYZsAsW7dOpYtWwbA9OnTeeihh3z7G0dpf54QwfJp0bEOz+861P6NGcDewYZDn2w5QmK82X+sKApNdl9Kc0Vt28XxxfZyrr10eJvnv7F+v3/At/Xnzs27Kvz39Z/VoWh/FRNGHv9Uq+uwdV9Vh7OQOgoOAF/trGDW5MG0Hq3fsreSE99CWga+Y2PavuVt2VPJ2CG+KbuxZl/3D8DWvVW+T+An7JpubXIx/8IhHZblZHxpVjTcHi8uj2+syOXRfD/7v7y+caRW57ya7p+M4PF/6Xg8zT9rvp+9mt7mca3luLlL8UwELUCUl5eTnZ3tP87KyqKoqKjdNX37+lZtGo1GEhMTqa2tJS0t7ZT37srOSJFE6he5QlW32kZXh7NxNJ12b5bovjfqE69vcnhIO2HguKX/viU4tHy3OTzt6nqyMticXtKT22eHtdo9DByQ2u5crLn9tR5Nx2xsf97h8rbbCKfJ4SXuhECga4BCu/M2l5cRQ/v4j9OTfa2qJueOdte2XD9yWEa789EqaAGio56rE1sGgVwjhDi9v/zk0lAXoVvKEE33iAZBmwCcnZ1NWVmZ/7i8vJzMzMx21xw75msaezwerFYrKSkdD1wJIYToWUELEOPGjaO4uJiSkhJcLhcFBQXk5+e3uSY/P5+33noLgPfff5/JkydLC0IIIcJE0GYxAaxfv56HH34Yr9fLNddcw5133sny5csZO3Ysl19+OU6nk5///Ofs3LmT5ORkHnvsMf+gthBCiNAKaoAQQggRuSQJiRBCiA5JgBBCCNGhsM7F5HQ6WbRoES6Xy78y+6677qKkpIS7776b+vp6Ro8ezR//+EfMZvPpbxiGWsZnsrKyeOqpp6Kqbvn5+VgsFlRVxWAw8Oabb1JXV8dPf/pTjhw5Qv/+/fnrX/9KcnLy6W8WhhoaGvj1r3/Nnj17UBSFhx9+mCFDhkRF/Q4cOMBPf/pT/3FJSQl33XUXCxYsiIr6vfDCC7z++usoisKIESN45JFHqKioiJrX3osvvsjrr7+Orutcd9113HLLLV177elhTNM0vbGxUdd1XXe5XPq1116rb9myRb/rrrv0NWvW6Lqu67/5zW/0V199NZTFPCPPPfecfvfdd+u33367rut6VNVt6tSpenV1dZtzjz76qP7UU0/puq7rTz31lP7HP/4xFEXrFvfee6/+2muv6bqu606nU6+vr4+q+rXweDz6lClT9NLS0qioX1lZmT516lTdbrfruu57zb3xxhtR89rbvXu3Pnv2bN1ms+lut1tfvHixfvDgwS797cK6i0lRFCwWC+BbJ+HxeFAUhS+++ILp06cDcNVVV1FYWBjKYnZZWVkZn3zyCddeey3gWzgYLXU7mcLCQhYsWADAggUL+Oijj0Jcoq5pbGzkq6++8v/tzGYzSUlJUVO/1jZu3EhOTg79+/ePmvp5vV4cDgcejweHw0FGRkbUvPb279/P+PHjiYuLw2g0cv755/Phhx926W8X1gECfH/I+fPnM2XKFKZMmUJOTg5JSUkYjb7esezsbMrLy0Ncyq55+OGH+fnPf47avJVibW1t1NStxa233srVV1/Nv//9bwCqq6v9CyYzMzOpqakJZfG6rKSkhLS0NH7xi1+wYMECfvWrX2Gz2aKmfq0VFBQwZ84cIDr+fllZWSxZsoSpU6dy0UUXkZCQwJgxY6LmtTdixAg2b95MbW0tdrudDRs2UFZW1qW/XdgHCIPBwKpVq1i/fj1FRUUcONA+Y2UkLq77+OOPSUtLY+zYsae8LhLr1uKf//wnb731Fk8//TSvvvoqX331VaiL1G08Hg87duxg4cKFvP3228TFxbFixYpQF6vbuVwu1q1bx4wZM0JdlG5TX19PYWEhhYWF/Oc///G/iZ4oUl97w4YNY+nSpSxZsoSlS5cycuRIDIb2eawCEfYBokVSUhKTJk3i22+/paGhAY/Hl52xrKysXQqPSPDNN9+wbt068vPzufvuu/niiy/4/e9/HxV1a5GVlQVAeno606ZNo6ioiPT0dCoqKgCoqKg4bWLGcJWdnU12djbjx48HYMaMGezYsSNq6tdiw4YNjBkzhj59fAntoqF+n3/+OQMGDCAtLQ2TycSVV17Jli1bouq1d9111/HWW2/x6quvkpKSwqBBg7r0twvrAFFTU0NDQwMADoeDzz//nGHDhjFp0iTef/99AN566612KTwiwT333MOGDRtYt24df/nLX5g8eTJ//vOfo6JuADabjcbGRv/Pn332Gbm5ueTn5/P2228D8Pbbb3P55ZeHsphdlpGRQXZ2tr9Fu3HjRoYNGxY19WtRUFDA7Nmz/cfRUL9+/fqxdetW7HY7uq6zceNGhg8fHjWvPfB1BQIcPXqUDz74gDlz5nTpbxfWK6l37drF/fffj9frRdd1ZsyYwbJlyygpKeGnP/0p9fX15OXl8ac//Slip6MBbNq0ieeee84/zTUa6lZSUsJ//dd/Ab5xpDlz5nDnnXdSW1vLT37yE44dO0bfvn1Zvnx5xCZo3LlzJ7/61a9wu93k5OTwyCOPoGla1NTPbrdz2WWX8dFHH5GY6EvtHS1/v8cff5x3330Xo9FIXl4ev//97ykvL4+K1x7AjTfeSF1dHUajkV/84hdccMEFXfrbhXWAEEIIETph3cUkhBAidCRACCGE6JAECCGEEB2SACGEEKJDEiCEEEJ0KKyzuQpxKtdddx0ulwu3201xcTG5ubkAjB49mkceeSTEpQvM9u3bKSkpiaqVyiJ6yDRXEfFKS0u55ppr2LRpU6iL0o7H4/Hn9+nI66+/zueff85jjz3W7fcW4kzJvy4RlVauXMm//vUvvF4vSUlJ/Pa3v2Xw4MG8/vrrrF27FovFwp49e+jbty+//OUvefTRRykpKWH8+PE8+uijKIrCz372M+Li4jh8+DBlZWVMmjSJ3/zmN5hMJqxWKw8//DB79+7F6XQyZcoU7rvvPlRVZeHChUycOJEtW7YQHx/P448/7l8k6HQ6GT9+PL/97W9paGjg//7v/2hqamL+/PlMmjSJRYsWceONN/LZZ58BcOjQIf/xoUOHWLhwIddffz1ffPEFV199NfPnz+cvf/kLmzdvxuVykZeXx4MPPkhcXFyI/wIiKgQpJbkQPaakpESfOHGi//iLL77Q77jjDt3pdOq6ruuFhYX6okWLdF3X9ddee02fOHGiXlZWpuu6ri9ZskRfsGCBbrVadZfLpc+aNUv/4osvdF3X9XvuuUefP3++3tTUpLtcLv3mm2/W//GPf+i6ruv33Xefvnr1al3Xdd3r9ep33XWXvnLlSl3Xdf2GG27Qf/SjH+kej8f/eF1dnf/nu+++27+PxGuvvab/5Cc/8Ze9uLhYnzJlSofHxcXF+ogRI/S1a9f6H3/88cf9Of51XdcfeeQRffny5Wf2H1SIZtKCEFFn3bp17Nixg+uuuw7w7bPR1NTkf/y8887zJxIcPXo0DoeDhIQEAEaOHMnhw4eZNGkSALNmzSI+Ph7w5dD/5JNPWLhwIR9//DHbt2/n6aefBny5wgYOHOj/HXPnzvVn0NQ0jRUrVvDpp5+iaRp1dXVd3oUtPj7ev2dBS13tdjsFBQWAL/vqmDFjunRvIU4kAUJEHV3X+d73vseyZcs6fDwmJsb/s6qq7Y5bMnp2dN+WFNCapvHUU0/Rr1+/Dq9tCSoAq1atoqioiH/84x9YLBb+/ve/c+zYsQ6fZzAY0DTNf+x0Ok9635Yy/e53v+P888/v8H5CnAmZ5iqiTkvWypYNX7xeL9u2bevSvd577z3sdjtut5vVq1f7Wxb5+fmsWLECr9cL+DIPl5SUdHgPq9VKamoqFouF+vp6/6d9AIvFgtVq9R9nZmbicDj891qzZs1p6/rcc8/5A0ljYyP79+/vUl2FOJEECBF1Jk+ezLJly7jjjjuYN28ec+fO5ZNPPunSvc477zzuvPNO5syZQ05Ojn+L0d/85jdomsb8+fOZO3cut912G5WVlR3e46qrrqKuro45c+Zw9913t/m0f+GFF2K1Wpk3bx4PP/wwZrOZ+++/n8WLF3PTTTdhMplOWb4f/vCHDBs2jGuvvZa5c+eyaNEiDh482KW6CnEimeYqxEn87Gc/47zzzmPhwoWhLooQISEtCCGEEB2SFoQQQogOSQtCCCFEhyRACCGE6JAECCGEEB2SACGEEKJDEiCEEEJ06P8H8uC1Bmw5yFwAAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"sns.set(color_codes=True)\n",
"plt.xlim(30,90)\n",
"plt.ylim(0,1)\n",
"sns.regplot(x='Temperature', y='Frequency', data=data, logistic=True)\n",
"plt.show()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"**I think I have managed to correctly compute and plot the uncertainty of my prediction.** Although the shaded area seems very similar to [the one obtained by with R](https://app-learninglab.inria.fr/gitlab/moocrr-session1/moocrr-reproducibility-study/raw/5c9dbef11b4d7638b7ddf2ea71026e7bf00fcfb0/challenger.pdf), I can spot a few differences (e.g., the blue point for temperature 63 is outside)... Could this be a numerical error ? Or a difference in the statistical method ? It is not clear which one is \"right\"."
]
}
],
"metadata": {
"celltoolbar": "Hide code",
"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.6.6"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
module3/exo1/analyse-syndrome-grippal.ipynb
View file @ 200ce142
This source diff could not be displayed because it is too large. You can view the blob instead.
module3/exo1/incidence-PAY-3.csv 0 → 100644
View file @ 200ce142
# @source="rseau Sentinelles, INSERM, Sorbonne Universit, http://www.sentiweb.fr", @meta={"period":[198444,202015],"geo":["PAY","1"],"geo_ref":"insee","indicator":"3","type":"all","conf_int":true,"compact":false}, @date=2020-04-23T16:15:17+02:00
week,indicator,inc,inc_low,inc_up,inc100,inc100_low,inc100_up,geo_insee,geo_name
202015,3,0,0,0,0,0,0,FR,France
202014,3,0,0,0,0,0,0,FR,France
202013,3,0,0,0,0,0,0,FR,France
202012,3,8321,5873,10769,13,9,17,FR,France
202011,3,101704,93652,109756,154,142,166,FR,France
202010,3,104977,96650,113304,159,146,172,FR,France
202009,3,110696,102066,119326,168,155,181,FR,France
202008,3,143753,133984,153522,218,203,233,FR,France
202007,3,183610,172812,194408,279,263,295,FR,France
202006,3,206669,195481,217857,314,297,331,FR,France
202005,3,187957,177445,198469,285,269,301,FR,France
202004,3,122331,113492,131170,186,173,199,FR,France
202003,3,78413,71330,85496,119,108,130,FR,France
202002,3,53614,47654,59574,81,72,90,FR,France
202001,3,36850,31608,42092,56,48,64,FR,France
201952,3,28135,23220,33050,43,36,50,FR,France
201951,3,29786,25042,34530,45,38,52,FR,France
201950,3,34223,29156,39290,52,44,60,FR,France
201949,3,25662,21414,29910,39,33,45,FR,France
201948,3,22367,18055,26679,34,27,41,FR,France
201947,3,18669,14759,22579,28,22,34,FR,France
201946,3,16030,12567,19493,24,19,29,FR,France
201945,3,10138,7160,13116,15,10,20,FR,France
201944,3,7822,5010,10634,12,8,16,FR,France
201943,3,9487,6448,12526,14,9,19,FR,France
201942,3,7747,5243,10251,12,8,16,FR,France
201941,3,7122,4720,9524,11,7,15,FR,France
201940,3,8505,5784,11226,13,9,17,FR,France
201939,3,7091,4462,9720,11,7,15,FR,France
201938,3,4897,2891,6903,7,4,10,FR,France
201937,3,3172,1367,4977,5,2,8,FR,France
201936,3,2295,728,3862,3,1,5,FR,France
201935,3,1010,2,2018,2,0,4,FR,France
201934,3,1672,279,3065,3,1,5,FR,France
201933,3,1593,68,3118,2,0,4,FR,France
201932,3,1592,71,3113,2,0,4,FR,France
201931,3,2337,695,3979,4,2,6,FR,France
201930,3,1944,415,3473,3,1,5,FR,France
201929,3,2096,420,3772,3,0,6,FR,France
201928,3,1195,125,2265,2,0,4,FR,France
201927,3,2466,848,4084,4,2,6,FR,France
201926,3,1798,402,3194,3,1,5,FR,France
201925,3,1544,389,2699,2,0,4,FR,France
201924,3,876,17,1735,1,0,2,FR,France
201923,3,2004,748,3260,3,1,5,FR,France
201922,3,2181,602,3760,3,1,5,FR,France
201921,3,2585,1102,4068,4,2,6,FR,France
201920,3,2891,1202,4580,4,1,7,FR,France
201919,3,3221,1509,4933,5,2,8,FR,France
201918,3,1385,278,2492,2,0,4,FR,France
201917,3,3166,1348,4984,5,2,8,FR,France
201916,3,3842,1973,5711,6,3,9,FR,France
201915,3,5932,3669,8195,9,6,12,FR,France
201914,3,12675,9167,16183,19,14,24,FR,France
201913,3,16513,12530,20496,25,19,31,FR,France
201912,3,22184,17786,26582,34,27,41,FR,France
201911,3,31553,26556,36550,48,40,56,FR,France
201910,3,49742,43268,56216,76,66,86,FR,France
201909,3,88354,79564,97144,134,121,147,FR,France
201908,3,172604,160024,185184,262,243,281,FR,France
201907,3,307338,291220,323456,467,443,491,FR,France
201906,3,394286,376782,411790,599,572,626,FR,France
201905,3,355785,339295,372275,540,515,565,FR,France
201904,3,241090,227261,254919,366,345,387,FR,France
201903,3,147063,135890,158236,223,206,240,FR,France
201902,3,75548,67632,83464,115,103,127,FR,France
201901,3,50295,43525,57065,76,66,86,FR,France
201852,3,37903,31375,44431,58,48,68,FR,France
201851,3,39259,32977,45541,60,50,70,FR,France
201850,3,27781,22638,32924,42,34,50,FR,France
201849,3,19738,15481,23995,30,24,36,FR,France
201848,3,19501,15275,23727,30,24,36,FR,France
201847,3,15949,12105,19793,24,18,30,FR,France
201846,3,11278,7957,14599,17,12,22,FR,France
201845,3,11065,7791,14339,17,12,22,FR,France
201844,3,6586,3875,9297,10,6,14,FR,France
201843,3,6550,3988,9112,10,6,14,FR,France
201842,3,7787,5129,10445,12,8,16,FR,France
201841,3,8048,5098,10998,12,8,16,FR,France
201840,3,7409,4717,10101,11,7,15,FR,France
201839,3,7174,4235,10113,11,7,15,FR,France
201838,3,7349,4399,10299,11,7,15,FR,France
201837,3,4915,2386,7444,7,3,11,FR,France
201836,3,3215,1349,5081,5,2,8,FR,France
201835,3,1506,239,2773,2,0,4,FR,France
201834,3,1368,116,2620,2,0,4,FR,France
201833,3,1962,5,3919,3,0,6,FR,France
201832,3,1839,183,3495,3,0,6,FR,France
201831,3,2048,242,3854,3,0,6,FR,France
201830,3,1951,202,3700,3,0,6,FR,France
201829,3,1951,252,3650,3,0,6,FR,France
201828,3,1654,52,3256,3,1,5,FR,France
201827,3,3266,1145,5387,5,2,8,FR,France
201826,3,3758,1493,6023,6,3,9,FR,France
201825,3,4580,2220,6940,7,3,11,FR,France
201824,3,3223,1351,5095,5,2,8,FR,France
201823,3,1207,136,2278,2,0,4,FR,France
201822,3,3202,1330,5074,5,2,8,FR,France
201821,3,2537,763,4311,4,1,7,FR,France
201820,3,2694,967,4421,4,1,7,FR,France
201819,3,1025,0,2098,2,0,4,FR,France
201818,3,3541,1416,5666,5,2,8,FR,France
201817,3,2573,1003,4143,4,2,6,FR,France
201816,3,4818,2724,6912,7,4,10,FR,France
201815,3,16311,12168,20454,25,19,31,FR,France
201814,3,22666,18092,27240,35,28,42,FR,France
201813,3,32680,25536,39824,50,39,61,FR,France
201812,3,52040,44590,59490,79,68,90,FR,France
201811,3,65159,56506,73812,99,86,112,FR,France
201810,3,93512,83989,103035,142,128,156,FR,France
201809,3,108545,98645,118445,165,150,180,FR,France
201808,3,131870,120950,142790,201,184,218,FR,France
201807,3,141138,130177,152099,215,198,232,FR,France
201806,3,138810,128164,149456,211,195,227,FR,France
201805,3,157839,146646,169032,240,223,257,FR,France
201804,3,175483,163835,187131,267,249,285,FR,France
201803,3,178329,166640,190018,271,253,289,FR,France
201802,3,198079,185663,210495,302,283,321,FR,France
201801,3,235493,220941,250045,358,336,380,FR,France
201752,3,299114,279359,318869,459,429,489,FR,France
201751,3,229448,214366,244530,352,329,375,FR,France
201750,3,129535,118787,140283,199,183,215,FR,France
201749,3,63789,56120,71458,98,86,110,FR,France
201748,3,37568,31541,43595,58,49,67,FR,France
201747,3,24972,19923,30021,38,30,46,FR,France
201746,3,13398,9948,16848,21,16,26,FR,France
201745,3,14035,8091,19979,22,13,31,FR,France
201744,3,8544,5318,11770,13,8,18,FR,France
201743,3,7300,4402,10198,11,7,15,FR,France
201742,3,9966,6503,13429,15,10,20,FR,France
201741,3,12434,8718,16150,19,13,25,FR,France
201740,3,12725,8690,16760,20,14,26,FR,France
201739,3,14949,10534,19364,23,16,30,FR,France
201738,3,11463,7438,15488,18,12,24,FR,France
201737,3,9716,6364,13068,15,10,20,FR,France
201736,3,2815,917,4713,4,1,7,FR,France
201735,3,2794,850,4738,4,1,7,FR,France
201734,3,2497,879,4115,4,2,6,FR,France
201733,3,2406,766,4046,4,1,7,FR,France
201732,3,2667,879,4455,4,1,7,FR,France
201731,3,3256,1158,5354,5,2,8,FR,France
201730,3,3759,1299,6219,6,2,10,FR,France
201729,3,5014,1989,8039,8,3,13,FR,France
201728,3,5271,2576,7966,8,4,12,FR,France
201727,3,3924,1432,6416,6,2,10,FR,France
201726,3,3171,1166,5176,5,2,8,FR,France
201725,3,837,0,1721,1,0,2,FR,France
201724,3,1566,248,2884,2,0,4,FR,France
201723,3,1664,203,3125,3,1,5,FR,France
201722,3,1305,92,2518,2,0,4,FR,France
201721,3,971,0,2046,1,0,3,FR,France
201720,3,2686,793,4579,4,1,7,FR,France
201719,3,3461,1490,5432,5,2,8,FR,France
201718,3,2102,515,3689,3,1,5,FR,France
201717,3,2071,428,3714,3,0,6,FR,France
201716,3,1380,222,2538,2,0,4,FR,France
201715,3,479,0,1242,1,0,2,FR,France
201714,3,1110,0,2549,2,0,4,FR,France
201713,3,7594,3808,11380,12,6,18,FR,France
201712,3,8780,4834,12726,13,7,19,FR,France
201711,3,7814,4329,11299,12,7,17,FR,France
201710,3,11802,7964,15640,18,12,24,FR,France
201709,3,13111,9099,17123,20,14,26,FR,France
201708,3,29545,23136,35954,45,35,55,FR,France
201707,3,59590,49764,69416,91,76,106,FR,France
201706,3,93628,82560,104696,144,127,161,FR,France
201705,3,193677,179255,208099,297,275,319,FR,France
201704,3,256428,240618,272238,394,370,418,FR,France
201703,3,267276,251345,283207,410,386,434,FR,France
201702,3,260588,245070,276106,400,376,424,FR,France
201701,3,255535,239743,271327,392,368,416,FR,France
201652,3,224489,207799,241179,345,319,371,FR,France
201651,3,187704,172655,202753,288,265,311,FR,France
201650,3,126570,115081,138059,194,176,212,FR,France
201649,3,76390,67374,85406,117,103,131,FR,France
201648,3,40502,33949,47055,62,52,72,FR,France
201647,3,22270,17612,26928,34,27,41,FR,France
201646,3,20816,16125,25507,32,25,39,FR,France
201645,3,15957,11448,20466,24,17,31,FR,France
201644,3,15145,10357,19933,23,16,30,FR,France
201643,3,9378,5830,12926,14,9,19,FR,France
201642,3,12920,9044,16796,20,14,26,FR,France
201641,3,13909,9709,18109,21,15,27,FR,France
201640,3,13379,8969,17789,21,14,28,FR,France
201639,3,10598,6814,14382,16,10,22,FR,France
201638,3,5969,3073,8865,9,5,13,FR,France
201637,3,4065,1538,6592,6,2,10,FR,France
201636,3,2372,473,4271,4,1,7,FR,France
201635,3,3420,1350,5490,5,2,8,FR,France
201634,3,2233,0,4643,3,0,7,FR,France
201633,3,1567,0,3938,2,0,6,FR,France
201632,3,2100,93,4107,3,0,6,FR,France
201631,3,3691,0,9120,6,0,14,FR,France
201630,3,3209,644,5774,5,1,9,FR,France
201629,3,1469,0,3073,2,0,4,FR,France
201628,3,2725,357,5093,4,0,8,FR,France
201627,3,5905,2343,9467,9,4,14,FR,France
201626,3,4186,1636,6736,6,2,10,FR,France
201625,3,4611,1560,7662,7,2,12,FR,France
201624,3,4977,2404,7550,8,4,12,FR,France
201623,3,4893,2416,7370,8,4,12,FR,France
201622,3,5667,3015,8319,9,5,13,FR,France
201621,3,5603,2475,8731,9,4,14,FR,France
201620,3,3908,1679,6137,6,3,9,FR,France
201619,3,4738,2308,7168,7,3,11,FR,France
201618,3,8204,4394,12014,13,7,19,FR,France
201617,3,13385,8822,17948,21,14,28,FR,France
201616,3,25670,19632,31708,39,30,48,FR,France
201615,3,48441,39824,57058,74,61,87,FR,France
201614,3,92603,81261,103945,142,125,159,FR,France
201613,3,170654,156387,184921,262,240,284,FR,France
201612,3,272174,254338,290010,418,391,445,FR,France
201611,3,304543,286453,322633,467,439,495,FR,France
201610,3,261822,245799,277845,402,377,427,FR,France
201609,3,229943,214533,245353,353,329,377,FR,France
201608,3,195067,180874,209260,299,277,321,FR,France
201607,3,207359,193223,221495,318,296,340,FR,France
201606,3,196143,182672,209614,301,280,322,FR,France
201605,3,178963,166296,191630,275,256,294,FR,France
201604,3,148473,136591,160355,228,210,246,FR,France
201603,3,75277,66630,83924,116,103,129,FR,France
201602,3,44617,37821,51413,68,58,78,FR,France
201601,3,42263,35556,48970,65,55,75,FR,France
201553,3,24835,19106,30564,38,29,47,FR,France
201552,3,19495,14319,24671,30,22,38,FR,France
201551,3,29080,23014,35146,45,36,54,FR,France
201550,3,25043,19812,30274,39,31,47,FR,France
201549,3,21479,16686,26272,33,26,40,FR,France
201548,3,16428,12273,20583,25,19,31,FR,France
201547,3,14959,10856,19062,23,17,29,FR,France
201546,3,12599,8237,16961,20,13,27,FR,France
201545,3,13577,9192,17962,21,14,28,FR,France
201544,3,14596,10025,19167,23,16,30,FR,France
201543,3,16749,11745,21753,26,18,34,FR,France
201542,3,18675,13613,23737,29,21,37,FR,France
201541,3,13585,9830,17340,21,15,27,FR,France
201540,3,13100,9349,16851,20,14,26,FR,France
201539,3,8595,5671,11519,13,8,18,FR,France
201538,3,8402,4929,11875,13,8,18,FR,France
201537,3,8817,4089,13545,14,7,21,FR,France
201536,3,1922,440,3404,3,1,5,FR,France
201535,3,229,0,644,0,0,1,FR,France
201534,3,760,0,1618,1,0,2,FR,France
201533,3,1537,208,2866,2,0,4,FR,France
201532,3,1594,66,3122,2,0,4,FR,France
201531,3,2845,475,5215,4,0,8,FR,France
201530,3,2447,353,4541,4,1,7,FR,France
201529,3,2389,332,4446,4,1,7,FR,France
201528,3,2089,105,4073,3,0,6,FR,France
201527,3,2257,393,4121,3,0,6,FR,France
201526,3,2536,774,4298,4,1,7,FR,France
201525,3,2607,680,4534,4,1,7,FR,France
201524,3,2849,1000,4698,4,1,7,FR,France
201523,3,2702,857,4547,4,1,7,FR,France
201522,3,4262,1989,6535,7,3,11,FR,France
201521,3,5547,2875,8219,9,5,13,FR,France
201520,3,4594,2040,7148,7,3,11,FR,France
201519,3,6307,3256,9358,10,5,15,FR,France
201518,3,3664,1544,5784,6,3,9,FR,France
201517,3,6382,3523,9241,10,6,14,FR,France
201516,3,12098,7913,16283,19,13,25,FR,France
201515,3,19458,14192,24724,30,22,38,FR,France
201514,3,30940,24617,37263,48,38,58,FR,France
201513,3,49563,41986,57140,77,65,89,FR,France
201512,3,70024,61079,78969,109,95,123,FR,France
201511,3,99295,88838,109752,154,138,170,FR,France
201510,3,143931,131682,156180,223,204,242,FR,France
201509,3,238947,222902,254992,370,345,395,FR,France
201508,3,409972,388670,431274,635,602,668,FR,France
201507,3,517750,494301,541199,802,766,838,FR,France
201506,3,533317,510627,556007,827,792,862,FR,France
201505,3,457911,436964,478858,710,678,742,FR,France
201504,3,292161,275491,308831,453,427,479,FR,France
201503,3,153931,141290,166572,239,219,259,FR,France
201502,3,104096,93358,114834,161,144,178,FR,France
201501,3,70962,60505,81419,110,94,126,FR,France
201452,3,63168,52705,73631,98,82,114,FR,France
201451,3,45346,38002,52690,71,60,82,FR,France
201450,3,42182,35049,49315,66,55,77,FR,France
201449,3,33537,27281,39793,52,42,62,FR,France
201448,3,28156,21896,34416,44,34,54,FR,France
201447,3,18158,13186,23130,28,20,36,FR,France
201446,3,14764,10174,19354,23,16,30,FR,France
201445,3,16172,11449,20895,25,18,32,FR,France
201444,3,15500,10464,20536,24,16,32,FR,France
201443,3,16011,10653,21369,25,17,33,FR,France
201442,3,21651,15828,27474,34,25,43,FR,France
201441,3,22486,16295,28677,35,25,45,FR,France
201440,3,16526,11583,21469,26,18,34,FR,France
201439,3,8583,5054,12112,13,8,18,FR,France
201438,3,10567,6496,14638,16,10,22,FR,France
201437,3,7714,3780,11648,12,6,18,FR,France
201436,3,5018,1949,8087,8,3,13,FR,France
201435,3,2268,355,4181,4,1,7,FR,France
201434,3,2367,250,4484,4,1,7,FR,France
201433,3,2358,0,4875,4,0,8,FR,France
201432,3,3229,447,6011,5,1,9,FR,France
201431,3,4143,349,7937,6,0,12,FR,France
201430,3,4631,902,8360,7,1,13,FR,France
201429,3,5122,1258,8986,8,2,14,FR,France
201428,3,4145,1200,7090,6,1,11,FR,France
201427,3,4823,1997,7649,8,4,12,FR,France
201426,3,3858,1235,6481,6,2,10,FR,France
201425,3,4564,1276,7852,7,2,12,FR,France
201424,3,2007,249,3765,3,0,6,FR,France
201423,3,5164,1853,8475,8,3,13,FR,France
201422,3,6907,1856,11958,11,3,19,FR,France
201421,3,4523,1039,8007,7,2,12,FR,France
201420,3,5881,2432,9330,9,4,14,FR,France
201419,3,3588,1277,5899,6,2,10,FR,France
201418,3,4708,1654,7762,7,2,12,FR,France
201417,3,2989,625,5353,5,1,9,FR,France
201416,3,7057,3381,10733,11,5,17,FR,France
201415,3,10800,6762,14838,17,11,23,FR,France
201414,3,19217,13490,24944,30,21,39,FR,France
201413,3,21577,15435,27719,34,24,44,FR,France
201412,3,23524,17260,29788,37,27,47,FR,France
201411,3,44513,35939,53087,69,56,82,FR,France
201410,3,86105,73497,98713,134,114,154,FR,France
201409,3,147035,131889,162181,229,205,253,FR,France
201408,3,188419,171950,204888,293,267,319,FR,France
201407,3,208726,191283,226169,325,298,352,FR,France
201406,3,160403,145280,175526,250,226,274,FR,France
201405,3,120261,107176,133346,187,167,207,FR,France
201404,3,85280,74414,96146,133,116,150,FR,France
201403,3,64159,54089,74229,100,84,116,FR,France
201402,3,46850,38202,55498,73,60,86,FR,France
201401,3,35624,27451,43797,55,42,68,FR,France
201352,3,30142,22721,37563,47,35,59,FR,France
201351,3,29787,23215,36359,47,37,57,FR,France
201350,3,24055,18214,29896,38,29,47,FR,France
201349,3,18071,13058,23084,28,20,36,FR,France
201348,3,17606,12195,23017,28,20,36,FR,France
201347,3,22636,16343,28929,35,25,45,FR,France
201346,3,13340,8367,18313,21,13,29,FR,France
201345,3,10596,6057,15135,17,10,24,FR,France
201344,3,10193,5622,14764,16,9,23,FR,France
201343,3,15780,9497,22063,25,15,35,FR,France
201342,3,14830,9545,20115,23,15,31,FR,France
201341,3,12391,7265,17517,19,11,27,FR,France
201340,3,9794,5140,14448,15,8,22,FR,France
201339,3,5371,1704,9038,8,2,14,FR,France
201338,3,9324,4198,14450,15,7,23,FR,France
201337,3,5027,1752,8302,8,3,13,FR,France
201336,3,4639,1425,7853,7,2,12,FR,France
201335,3,1291,0,3112,2,0,5,FR,France
201334,3,1889,0,4912,3,0,8,FR,France
201333,3,4740,348,9132,7,0,14,FR,France
201332,3,4709,549,8869,7,0,14,FR,France
201331,3,2270,0,5811,4,0,10,FR,France
201330,3,2169,0,4630,3,0,7,FR,France
201329,3,1696,0,3558,3,0,6,FR,France
201328,3,983,0,2168,2,0,4,FR,France
201327,3,4278,1327,7229,7,2,12,FR,France
201326,3,3488,839,6137,5,1,9,FR,France
201325,3,2735,0,5826,4,0,9,FR,France
201324,3,7245,2997,11493,11,4,18,FR,France
201323,3,16094,8918,23270,25,14,36,FR,France
201322,3,7245,2923,11567,11,4,18,FR,France
201321,3,5974,2401,9547,9,3,15,FR,France
201320,3,3831,1245,6417,6,2,10,FR,France
201319,3,1166,0,2680,2,0,4,FR,France
201318,3,3415,757,6073,5,1,9,FR,France
201317,3,7941,2950,12932,12,4,20,FR,France
201316,3,9520,4758,14282,15,8,22,FR,France
201315,3,19357,13311,25403,30,21,39,FR,France
201314,3,36855,27737,45973,58,44,72,FR,France
201313,3,43236,34048,52424,68,54,82,FR,France
201312,3,60822,50414,71230,95,79,111,FR,France
201311,3,84308,72456,96160,132,113,151,FR,France
201310,3,165558,148471,182645,259,232,286,FR,France
201309,3,259984,239067,280901,407,374,440,FR,France
201308,3,351626,327537,375715,550,512,588,FR,France
201307,3,470076,443682,496470,736,695,777,FR,France
201306,3,471647,445112,498182,738,696,780,FR,France
201305,3,492026,464729,519323,770,727,813,FR,France
201304,3,391525,367182,415868,613,575,651,FR,France
201303,3,224661,207134,242188,352,325,379,FR,France
201302,3,184159,168343,199975,288,263,313,FR,France
201301,3,163491,147340,179642,256,231,281,FR,France
201252,3,159290,141978,176602,251,224,278,FR,France
201251,3,113901,99835,127967,179,157,201,FR,France
201250,3,79918,68086,91750,126,107,145,FR,France
201249,3,47537,38711,56363,75,61,89,FR,France
201248,3,43006,32863,53149,68,52,84,FR,France
201247,3,39223,29140,49306,62,46,78,FR,France
201246,3,26436,18286,34586,42,29,55,FR,France
201245,3,18893,12799,24987,30,20,40,FR,France
201244,3,13106,8784,17428,21,14,28,FR,France
201243,3,17791,11284,24298,28,18,38,FR,France
201242,3,13096,8369,17823,21,14,28,FR,France
201241,3,18675,12010,25340,29,19,39,FR,France
201240,3,17873,11545,24201,28,18,38,FR,France
201239,3,20257,14180,26334,32,22,42,FR,France
201238,3,14715,9178,20252,23,14,32,FR,France
201237,3,7775,3509,12041,12,5,19,FR,France
201236,3,6384,2420,10348,10,4,16,FR,France
201235,3,8277,2206,14348,13,3,23,FR,France
201234,3,3858,0,8641,6,0,14,FR,France
201233,3,6492,840,12144,10,1,19,FR,France
201232,3,6381,0,14034,10,0,22,FR,France
201231,3,2696,0,5763,4,0,9,FR,France
201230,3,1834,0,4194,3,0,7,FR,France
201229,3,6329,1424,11234,10,2,18,FR,France
201228,3,4628,1091,8165,7,1,13,FR,France
201227,3,4162,1487,6837,7,3,11,FR,France
201226,3,4819,1097,8541,8,2,14,FR,France
201225,3,2669,530,4808,4,1,7,FR,France
201224,3,5977,2335,9619,9,3,15,FR,France
201223,3,2361,576,4146,4,1,7,FR,France
201222,3,5772,2435,9109,9,4,14,FR,France
201221,3,3722,821,6623,6,1,11,FR,France
201220,3,2731,897,4565,4,1,7,FR,France
201219,3,3267,1317,5217,5,2,8,FR,France
201218,3,2328,569,4087,4,1,7,FR,France
201217,3,6982,2701,11263,11,4,18,FR,France
201216,3,11455,6874,16036,18,11,25,FR,France
201215,3,20437,13435,27439,32,21,43,FR,France
201214,3,36740,27799,45681,58,44,72,FR,France
201213,3,62917,50171,75663,99,79,119,FR,France
201212,3,84738,71566,97910,133,112,154,FR,France
201211,3,128945,113778,144112,203,179,227,FR,France
201210,3,174704,157337,192071,275,248,302,FR,France
201209,3,225648,206810,244486,355,325,385,FR,France
201208,3,287293,264901,309685,452,417,487,FR,France
201207,3,263557,242711,284403,415,382,448,FR,France
201206,3,175123,158942,191304,276,251,301,FR,France
201205,3,106110,92524,119696,167,146,188,FR,France
201204,3,67778,56784,78772,107,90,124,FR,France
201203,3,49440,40981,57899,78,65,91,FR,France
201202,3,40507,32652,48362,64,52,76,FR,France
201201,3,45104,35640,54568,71,56,86,FR,France
201152,3,36532,28367,44697,58,45,71,FR,France
201151,3,25835,19002,32668,41,30,52,FR,France
201150,3,24155,18145,30165,38,28,48,FR,France
201149,3,27209,20376,34042,43,32,54,FR,France
201148,3,32358,24451,40265,51,38,64,FR,France
201147,3,17627,11780,23474,28,19,37,FR,France
201146,3,22368,15730,29006,35,24,46,FR,France
201145,3,18104,11995,24213,29,19,39,FR,France
201144,3,14050,8471,19629,22,13,31,FR,France
201143,3,21522,14251,28793,34,22,46,FR,France
201142,3,16932,11363,22501,27,18,36,FR,France
201141,3,7486,4024,10948,12,7,17,FR,France
201140,3,14053,8659,19447,22,13,31,FR,France
201139,3,17021,11062,22980,27,18,36,FR,France
201138,3,14007,8465,19549,22,13,31,FR,France
201137,3,8365,4174,12556,13,6,20,FR,France
201136,3,4606,1714,7498,7,2,12,FR,France
201135,3,4178,1634,6722,7,3,11,FR,France
201134,3,3312,683,5941,5,1,9,FR,France
201133,3,2421,0,5277,4,0,9,FR,France
201132,3,2590,119,5061,4,0,8,FR,France
201131,3,2409,320,4498,4,1,7,FR,France
201130,3,3972,1395,6549,6,2,10,FR,France
201129,3,2265,717,3813,4,2,6,FR,France
201128,3,1977,255,3699,3,0,6,FR,France
201127,3,3524,551,6497,6,1,11,FR,France
201126,3,1155,0,2685,2,0,4,FR,France
201125,3,3791,1006,6576,6,2,10,FR,France
201124,3,3507,1149,5865,6,2,10,FR,France
201123,3,3406,1170,5642,5,1,9,FR,France
201122,3,3122,1125,5119,5,2,8,FR,France
201121,3,2694,769,4619,4,1,7,FR,France
201120,3,3951,1443,6459,6,2,10,FR,France
201119,3,5416,1995,8837,9,4,14,FR,France
201118,3,5306,1815,8797,8,2,14,FR,France
201117,3,2659,310,5008,4,0,8,FR,France
201116,3,5061,1988,8134,8,3,13,FR,France
201115,3,6557,3038,10076,10,4,16,FR,France
201114,3,7352,3784,10920,12,6,18,FR,France
201113,3,8749,5213,12285,14,8,20,FR,France
201112,3,12491,7845,17137,20,13,27,FR,France
201111,3,14619,10009,19229,23,16,30,FR,France
201110,3,33628,26073,41183,53,41,65,FR,France
201109,3,54522,46057,62987,86,73,99,FR,France
201108,3,78392,67885,88899,124,107,141,FR,France
201107,3,149219,134798,163640,236,213,259,FR,France
201106,3,216579,200236,232922,343,317,369,FR,France
201105,3,278192,260235,296149,440,412,468,FR,France
201104,3,278218,259626,296810,440,411,469,FR,France
201103,3,275882,257313,294451,437,408,466,FR,France
201102,3,308949,288320,329578,489,456,522,FR,France
201101,3,309933,288190,331676,490,456,524,FR,France
201052,3,228105,208046,248164,363,331,395,FR,France
201051,3,158957,142576,175338,253,227,279,FR,France
201050,3,104874,91639,118109,167,146,188,FR,France
201049,3,54971,46418,63524,88,74,102,FR,France
201048,3,36982,29304,44660,59,47,71,FR,France
201047,3,26598,20123,33073,42,32,52,FR,France
201046,3,15465,10915,20015,25,18,32,FR,France
201045,3,9040,5772,12308,14,9,19,FR,France
201044,3,15833,10255,21411,25,16,34,FR,France
201043,3,16002,11001,21003,25,17,33,FR,France
201042,3,18709,13376,24042,30,22,38,FR,France
201041,3,17352,12182,22522,28,20,36,FR,France
201040,3,18331,12824,23838,29,20,38,FR,France
201039,3,11978,7913,16043,19,13,25,FR,France
201038,3,8936,4906,12966,14,8,20,FR,France
201037,3,5318,2551,8085,8,4,12,FR,France
201036,3,4954,2058,7850,8,3,13,FR,France
201035,3,4279,1604,6954,7,3,11,FR,France
201034,3,4926,566,9286,8,1,15,FR,France
201033,3,2950,0,6507,5,0,11,FR,France
201032,3,671,0,1522,1,0,2,FR,France
201031,3,2732,509,4955,4,0,8,FR,France
201030,3,3989,1068,6910,6,1,11,FR,France
201029,3,4123,1531,6715,7,3,11,FR,France
201028,3,1209,0,2640,2,0,4,FR,France
201027,3,2494,433,4555,4,1,7,FR,France
201026,3,5808,2454,9162,9,4,14,FR,France
201025,3,8019,4395,11643,13,7,19,FR,France
201024,3,3584,1657,5511,6,3,9,FR,France
201023,3,4016,1905,6127,6,3,9,FR,France
201022,3,6251,3664,8838,10,6,14,FR,France
201021,3,3860,1636,6084,6,2,10,FR,France
201020,3,2470,442,4498,4,1,7,FR,France
201019,3,2273,636,3910,4,1,7,FR,France
201018,3,4095,1639,6551,7,3,11,FR,France
201017,3,2216,817,3615,4,2,6,FR,France
201016,3,3471,1085,5857,6,2,10,FR,France
201015,3,3990,1376,6604,6,2,10,FR,France
201014,3,6188,3191,9185,10,5,15,FR,France
201013,3,6705,3665,9745,11,6,16,FR,France
201012,3,6727,3582,9872,11,6,16,FR,France
201011,3,8692,5002,12382,14,8,20,FR,France
201010,3,8844,5465,12223,14,9,19,FR,France
201009,3,12081,7334,16828,19,11,27,FR,France
201008,3,15433,10422,20444,25,17,33,FR,France
201007,3,19690,11728,27652,31,18,44,FR,France
201006,3,27858,18883,36833,44,30,58,FR,France
201005,3,22197,17019,27375,35,27,43,FR,France
201004,3,33817,26151,41483,54,42,66,FR,France
201003,3,36236,27936,44536,58,45,71,FR,France
201002,3,40691,32891,48491,65,53,77,FR,France
201001,3,66574,56164,76984,106,89,123,FR,France
200953,3,96141,81818,110464,153,130,176,FR,France
200952,3,172872,153238,192506,275,244,306,FR,France
200951,3,247809,229140,266478,395,365,425,FR,France
200950,3,379424,355974,402874,604,567,641,FR,France
200949,3,473718,447756,499680,754,713,795,FR,France
200948,3,467551,441506,493596,744,703,785,FR,France
200947,3,378924,355752,402096,603,566,640,FR,France
200946,3,223847,203918,243776,356,324,388,FR,France
200945,3,147046,131155,162937,234,209,259,FR,France
200944,3,170231,151298,189164,271,241,301,FR,France
200943,3,131408,114883,147933,209,183,235,FR,France
200942,3,99387,84491,114283,158,134,182,FR,France
200941,3,100802,86381,115223,160,137,183,FR,France
200940,3,111587,97628,125546,178,156,200,FR,France
200939,3,121267,106153,136381,193,169,217,FR,France
200938,3,137924,121242,154606,220,193,247,FR,France
200937,3,101332,87522,115142,161,139,183,FR,France
200936,3,46024,36435,55613,73,58,88,FR,France
200935,3,31660,23669,39651,50,37,63,FR,France
200934,3,24985,17699,32271,40,28,52,FR,France
200933,3,24435,16308,32562,39,26,52,FR,France
200932,3,23151,14624,31678,37,23,51,FR,France
200931,3,30255,9818,50692,48,15,81,FR,France
200930,3,13868,7861,19875,22,12,32,FR,France
200929,3,10704,5314,16094,17,8,26,FR,France
200928,3,15516,273,30759,25,1,49,FR,France
200927,3,10354,2760,17948,16,4,28,FR,France
200926,3,9932,5257,14607,16,9,23,FR,France
200925,3,7619,4026,11212,12,6,18,FR,France
200924,3,7794,3409,12179,12,5,19,FR,France
200923,3,4342,1975,6709,7,3,11,FR,France
200922,3,6390,3020,9760,10,5,15,FR,France
200921,3,7118,3704,10532,11,6,16,FR,France
200920,3,9482,5804,13160,15,9,21,FR,France
200919,3,7538,4299,10777,12,7,17,FR,France
200918,3,9627,5121,14133,15,8,22,FR,France
200917,3,9537,5267,13807,15,8,22,FR,France
200916,3,14666,9218,20114,23,14,32,FR,France
200915,3,17533,12137,22929,28,19,37,FR,France
200914,3,26324,19150,33498,42,31,53,FR,France
200913,3,27737,20781,34693,44,33,55,FR,France
200912,3,38811,29935,47687,62,48,76,FR,France
200911,3,40410,30705,50115,64,49,79,FR,France
200910,3,53395,41943,64847,85,67,103,FR,France
200909,3,70558,57400,83716,112,91,133,FR,France
200908,3,109100,94840,123360,174,151,197,FR,France
200907,3,141531,125101,157961,225,199,251,FR,France
200906,3,276274,252759,299789,440,403,477,FR,France
200905,3,401868,375032,428704,640,597,683,FR,France
200904,3,545313,512892,577734,868,816,920,FR,France
200903,3,520957,487271,554643,829,775,883,FR,France
200902,3,287145,262869,311421,457,418,496,FR,France
200901,3,200349,180878,219820,319,288,350,FR,France
200852,3,198418,115769,281067,323,188,458,FR,France
200851,3,112588,96656,128520,184,158,210,FR,France
200850,3,51746,43319,60173,84,70,98,FR,France
200849,3,34950,27527,42373,57,45,69,FR,France
200848,3,25604,18563,32645,42,31,53,FR,France
200847,3,18706,10838,26574,30,17,43,FR,France
200846,3,10899,5461,16337,18,9,27,FR,France
200845,3,12472,7082,17862,20,11,29,FR,France
200844,3,9772,5661,13883,16,9,23,FR,France
200843,3,9345,5541,13149,15,9,21,FR,France
200842,3,9087,5572,12602,15,9,21,FR,France
200841,3,9925,5637,14213,16,9,23,FR,France
200840,3,7869,4535,11203,13,8,18,FR,France
200839,3,2195,506,3884,4,1,7,FR,France
200838,3,3679,1481,5877,6,2,10,FR,France
200837,3,2429,561,4297,4,1,7,FR,France
200836,3,2512,492,4532,4,1,7,FR,France
200835,3,1901,0,4869,3,0,8,FR,France
200834,3,2552,0,6509,4,0,10,FR,France
200833,3,7229,0,16106,12,0,26,FR,France
200832,3,4378,107,8649,7,0,14,FR,France
200831,3,3972,621,7323,6,1,11,FR,France
200830,3,3821,643,6999,6,1,11,FR,France
200829,3,2055,0,4225,3,0,7,FR,France
200828,3,3257,208,6306,5,0,10,FR,France
200827,3,2862,474,5250,5,1,9,FR,France
200826,3,897,0,2062,1,0,3,FR,France
200825,3,2006,360,3652,3,0,6,FR,France
200824,3,3217,823,5611,5,1,9,FR,France
200823,3,2906,1006,4806,5,2,8,FR,France
200822,3,1792,451,3133,3,1,5,FR,France
200821,3,3642,925,6359,6,2,10,FR,France
200820,3,6452,2790,10114,11,5,17,FR,France
200819,3,5173,1485,8861,8,2,14,FR,France
200818,3,5315,1770,8860,9,3,15,FR,France
200817,3,10150,5835,14465,17,10,24,FR,France
200816,3,11876,7054,16698,19,11,27,FR,France
200815,3,22871,16253,29489,37,26,48,FR,France
200814,3,45610,27907,63313,74,45,103,FR,France
200813,3,38600,29300,47900,63,48,78,FR,France
200812,3,53169,43299,63039,87,71,103,FR,France
200811,3,62014,50430,73598,101,82,120,FR,France
200810,3,90342,75729,104955,147,123,171,FR,France
200809,3,119696,104205,135187,195,170,220,FR,France
200808,3,198688,180749,216627,324,295,353,FR,France
200807,3,286498,264535,308461,467,431,503,FR,France
200806,3,377097,350714,403480,615,572,658,FR,France
200805,3,344050,320436,367664,561,523,599,FR,France
200804,3,290196,268193,312199,473,437,509,FR,France
200803,3,260268,233385,287151,424,380,468,FR,France
200802,3,160872,145311,176433,262,237,287,FR,France
200801,3,102562,88547,116577,167,144,190,FR,France
200752,3,86742,72295,101189,141,117,165,FR,France
200751,3,61954,51764,72144,101,84,118,FR,France
200750,3,45366,36923,53809,74,60,88,FR,France
200749,3,38351,26177,50525,63,43,83,FR,France
200748,3,48174,35418,60930,79,58,100,FR,France
200747,3,30716,23351,38081,50,38,62,FR,France
200746,3,18204,12267,24141,30,20,40,FR,France
200745,3,15521,10195,20847,25,16,34,FR,France
200744,3,13904,9072,18736,23,15,31,FR,France
200743,3,16212,10778,21646,26,17,35,FR,France
200742,3,20197,13955,26439,33,23,43,FR,France
200741,3,12987,8464,17510,21,14,28,FR,France
200740,3,13517,2983,24051,22,5,39,FR,France
200739,3,15985,7519,24451,26,12,40,FR,France
200738,3,4769,837,8701,8,2,14,FR,France
200737,3,1449,300,2598,2,0,4,FR,France
200736,3,1100,0,2424,2,0,4,FR,France
200735,3,1177,0,2872,2,0,5,FR,France
200734,3,3353,654,6052,5,1,9,FR,France
200733,3,4189,360,8018,7,1,13,FR,France
200732,3,132,0,688,0,0,1,FR,France
200731,3,1965,0,4461,3,0,7,FR,France
200730,3,2158,0,4583,4,0,8,FR,France
200729,3,3406,240,6572,6,1,11,FR,France
200728,3,1459,0,3467,2,0,5,FR,France
200727,3,1453,0,3382,2,0,5,FR,France
200726,3,157,0,513,0,0,1,FR,France
200725,3,3699,0,7778,6,0,13,FR,France
200724,3,2854,481,5227,5,1,9,FR,France
200723,3,5343,2164,8522,9,4,14,FR,France
200722,3,7689,4224,11154,13,7,19,FR,France
200721,3,4319,1659,6979,7,3,11,FR,France
200720,3,6400,2685,10115,10,4,16,FR,France
200719,3,3012,703,5321,5,1,9,FR,France
200718,3,2232,348,4116,4,1,7,FR,France
200717,3,4705,1736,7674,8,3,13,FR,France
200716,3,3817,1438,6196,6,2,10,FR,France
200715,3,4757,717,8797,8,1,15,FR,France
200714,3,11377,6957,15797,19,12,26,FR,France
200713,3,21271,13093,29449,35,22,48,FR,France
200712,3,25433,17717,33149,41,28,54,FR,France
200711,3,30764,23356,38172,50,38,62,FR,France
200710,3,49469,40368,58570,81,66,96,FR,France
200709,3,138901,117018,160784,226,190,262,FR,France
200708,3,316747,284988,348506,516,464,568,FR,France
200707,3,417698,386505,448891,681,630,732,FR,France
200706,3,500016,466668,533364,815,761,869,FR,France
200705,3,375457,346822,404092,612,565,659,FR,France
200704,3,209878,189332,230424,342,309,375,FR,France
200703,3,126161,110596,141726,206,181,231,FR,France
200702,3,82849,69493,96205,135,113,157,FR,France
200701,3,70842,58613,83071,115,95,135,FR,France
200652,3,78190,64191,92189,128,105,151,FR,France
200651,3,62734,50210,75258,103,82,124,FR,France
200650,3,35668,27280,44056,58,44,72,FR,France
200649,3,18773,13164,24382,31,22,40,FR,France
200648,3,16439,11225,21653,27,18,36,FR,France
200647,3,16104,9343,22865,26,15,37,FR,France
200646,3,14032,5590,22474,23,9,37,FR,France
200645,3,12285,7745,16825,20,13,27,FR,France
200644,3,8031,4172,11890,13,7,19,FR,France
200643,3,21157,12483,29831,35,21,49,FR,France
200642,3,16673,8249,25097,27,13,41,FR,France
200641,3,13920,9081,18759,23,15,31,FR,France
200640,3,11950,7132,16768,20,12,28,FR,France
200639,3,8945,4333,13557,15,7,23,FR,France
200638,3,2532,138,4926,4,0,8,FR,France
200637,3,1786,233,3339,3,0,6,FR,France
200636,3,1107,0,2346,2,0,4,FR,France
200635,3,2170,0,4760,4,0,8,FR,France
200634,3,1572,0,3756,3,0,7,FR,France
200633,3,0,0,0,0,0,0,FR,France
200632,3,0,0,0,0,0,0,FR,France
200631,3,1773,0,3892,3,0,6,FR,France
200630,3,665,0,1614,1,0,3,FR,France
200629,3,1339,0,3260,2,0,5,FR,France
200628,3,843,0,2245,1,0,3,FR,France
200627,3,1351,0,3031,2,0,5,FR,France
200626,3,4564,1506,7622,7,2,12,FR,France
200625,3,6259,2072,10446,10,3,17,FR,France
200624,3,4279,964,7594,7,2,12,FR,France
200623,3,4532,1489,7575,7,2,12,FR,France
200622,3,5807,2667,8947,10,5,15,FR,France
200621,3,4661,1761,7561,8,3,13,FR,France
200620,3,8859,4754,12964,15,8,22,FR,France
200619,3,5083,2011,8155,8,3,13,FR,France
200618,3,4454,1632,7276,7,2,12,FR,France
200617,3,17762,11802,23722,29,19,39,FR,France
200616,3,22096,15420,28772,36,25,47,FR,France
200615,3,26969,19631,34307,44,32,56,FR,France
200614,3,40321,31022,49620,66,51,81,FR,France
200613,3,77053,63531,90575,126,104,148,FR,France
200612,3,108004,93005,123003,177,152,202,FR,France
200611,3,122437,106388,138486,201,175,227,FR,France
200610,3,124749,109104,140394,204,178,230,FR,France
200609,3,164150,146099,182201,269,239,299,FR,France
200608,3,199088,178830,219346,326,293,359,FR,France
200607,3,243908,221071,266745,400,363,437,FR,France
200606,3,256758,235190,278326,421,386,456,FR,France
200605,3,184928,166901,202955,303,273,333,FR,France
200604,3,104620,90909,118331,171,149,193,FR,France
200603,3,66255,55727,76783,109,92,126,FR,France
200602,3,39211,31418,47004,64,51,77,FR,France
200601,3,34603,26269,42937,57,43,71,FR,France
200552,3,32132,24389,39875,53,40,66,FR,France
200551,3,40344,28397,52291,67,47,87,FR,France
200550,3,36880,28635,45125,61,47,75,FR,France
200549,3,27589,20576,34602,46,34,58,FR,France
200548,3,24860,17839,31881,41,29,53,FR,France
200547,3,26884,19810,33958,44,32,56,FR,France
200546,3,23967,17155,30779,40,29,51,FR,France
200545,3,14609,8709,20509,24,14,34,FR,France
200544,3,20481,13625,27337,34,23,45,FR,France
200543,3,18045,11916,24174,30,20,40,FR,France
200542,3,23653,15209,32097,39,25,53,FR,France
200541,3,26466,18270,34662,44,30,58,FR,France
200540,3,26647,18088,35206,44,30,58,FR,France
200539,3,23994,15381,32607,40,26,54,FR,France
200538,3,21881,13979,29783,36,23,49,FR,France
200537,3,10347,5384,15310,17,9,25,FR,France
200536,3,7082,2806,11358,12,5,19,FR,France
200535,3,2448,0,4987,4,0,8,FR,France
200534,3,5281,539,10023,9,1,17,FR,France
200533,3,2815,0,6252,5,0,11,FR,France
200532,3,1572,0,3690,3,0,6,FR,France
200531,3,3767,0,10955,6,0,18,FR,France
200530,3,2203,0,5072,4,0,9,FR,France
200529,3,3965,739,7191,7,2,12,FR,France
200528,3,1988,96,3880,3,0,6,FR,France
200527,3,2305,0,4713,4,0,8,FR,France
200526,3,4720,1385,8055,8,2,14,FR,France
200525,3,2354,0,4747,4,0,8,FR,France
200524,3,5892,2029,9755,10,4,16,FR,France
200523,3,5246,1722,8770,9,3,15,FR,France
200522,3,3935,1031,6839,6,1,11,FR,France
200521,3,4284,474,8094,7,1,13,FR,France
200520,3,7841,3611,12071,13,6,20,FR,France
200519,3,7638,3346,11930,13,6,20,FR,France
200518,3,2755,388,5122,5,1,9,FR,France
200517,3,7166,3530,10802,12,6,18,FR,France
200516,3,9741,5344,14138,16,9,23,FR,France
200515,3,14881,9018,20744,25,15,35,FR,France
200514,3,21988,14834,29142,36,24,48,FR,France
200513,3,42971,33053,52889,71,55,87,FR,France
200512,3,90040,75374,104706,149,125,173,FR,France
200511,3,183032,161407,204657,302,266,338,FR,France
200510,3,222526,199990,245062,367,330,404,FR,France
200509,3,251784,229464,274104,416,379,453,FR,France
200508,3,277423,254825,300021,458,421,495,FR,France
200507,3,459444,430026,488862,759,710,808,FR,France
200506,3,568610,536207,601013,939,885,993,FR,France
200505,3,568285,534912,601658,938,883,993,FR,France
200504,3,316856,292011,341701,523,482,564,FR,France
200503,3,154208,137061,171355,255,227,283,FR,France
200502,3,88853,75576,102130,147,125,169,FR,France
200501,3,55801,45339,66263,92,75,109,FR,France
200453,3,40695,31926,49464,68,53,83,FR,France
200452,3,32687,24316,41058,54,40,68,FR,France
200451,3,46209,36017,56401,77,60,94,FR,France
200450,3,42358,31865,52851,70,53,87,FR,France
200449,3,27092,17910,36274,45,30,60,FR,France
200448,3,30522,21236,39808,51,36,66,FR,France
200447,3,21148,13386,28910,35,22,48,FR,France
200446,3,13934,8335,19533,23,14,32,FR,France
200445,3,15658,9703,21613,26,16,36,FR,France
200444,3,23166,14852,31480,38,24,52,FR,France
200443,3,16780,10545,23015,28,18,38,FR,France
200442,3,20480,13909,27051,34,23,45,FR,France
200441,3,15852,9789,21915,26,16,36,FR,France
200440,3,25639,17291,33987,43,29,57,FR,France
200439,3,22680,13222,32138,38,22,54,FR,France
200438,3,17483,9417,25549,29,16,42,FR,France
200437,3,9103,4768,13438,15,8,22,FR,France
200436,3,7753,3596,11910,13,6,20,FR,France
200435,3,4408,775,8041,7,1,13,FR,France
200434,3,4501,610,8392,7,1,13,FR,France
200433,3,1697,0,4046,3,0,7,FR,France
200432,3,3496,0,8181,6,0,14,FR,France
200431,3,2949,20,5878,5,0,10,FR,France
200430,3,3557,642,6472,6,1,11,FR,France
200429,3,5164,1162,9166,9,2,16,FR,France
200428,3,6358,2224,10492,11,4,18,FR,France
200427,3,7486,3040,11932,12,5,19,FR,France
200426,3,7113,0,14709,12,0,25,FR,France
200425,3,3275,508,6042,5,0,10,FR,France
200424,3,6737,2394,11080,11,4,18,FR,France
200423,3,9792,807,18777,16,1,31,FR,France
200422,3,8077,3482,12672,13,5,21,FR,France
200421,3,8073,3588,12558,13,6,20,FR,France
200420,3,7885,3454,12316,13,6,20,FR,France
200419,3,7235,3101,11369,12,5,19,FR,France
200418,3,5108,1830,8386,8,3,13,FR,France
200417,3,3647,523,6771,6,1,11,FR,France
200416,3,5265,1359,9171,9,3,15,FR,France
200415,3,8272,3655,12889,14,6,22,FR,France
200414,3,15657,8888,22426,26,15,37,FR,France
200413,3,8316,4067,12565,14,7,21,FR,France
200412,3,12482,6899,18065,21,12,30,FR,France
200411,3,13964,7879,20049,23,13,33,FR,France
200410,3,22531,14164,30898,37,23,51,FR,France
200409,3,12618,7283,17953,21,12,30,FR,France
200408,3,17738,11055,24421,29,18,40,FR,France
200407,3,25484,17228,33740,42,28,56,FR,France
200406,3,27452,19505,35399,46,33,59,FR,France
200405,3,35110,25869,44351,58,43,73,FR,France
200404,3,44571,35266,53876,74,59,89,FR,France
200403,3,60482,50389,70575,100,83,117,FR,France
200402,3,95234,81407,109061,158,135,181,FR,France
200401,3,109116,84565,133667,181,140,222,FR,France
200352,3,72239,5520,138958,121,10,232,FR,France
200351,3,445464,389214,501714,744,650,838,FR,France
200350,3,522427,491781,553073,873,822,924,FR,France
200349,3,555492,524706,586278,928,877,979,FR,France
200348,3,514668,485364,543972,860,811,909,FR,France
200347,3,332927,305972,359882,556,511,601,FR,France
200346,3,146101,126840,165362,244,212,276,FR,France
200345,3,95565,80900,110230,160,135,185,FR,France
200344,3,54703,43193,66213,91,72,110,FR,France
200343,3,45834,36611,55057,77,62,92,FR,France
200342,3,39010,28885,49135,65,48,82,FR,France
200341,3,29579,21478,37680,49,35,63,FR,France
200340,3,30862,19484,42240,52,33,71,FR,France
200339,3,26592,18927,34257,44,31,57,FR,France
200338,3,14811,8777,20845,25,15,35,FR,France
200337,3,14735,7495,21975,25,13,37,FR,France
200336,3,10951,3732,18170,18,6,30,FR,France
200335,3,1504,0,3349,3,0,6,FR,France
200334,3,2194,0,4889,4,0,9,FR,France
200333,3,1504,0,3921,3,0,7,FR,France
200332,3,1619,0,4072,3,0,7,FR,France
200331,3,5164,0,10601,9,0,18,FR,France
200330,3,3183,28,6338,5,0,10,FR,France
200329,3,6801,450,13152,11,0,22,FR,France
200328,3,3975,609,7341,7,1,13,FR,France
200327,3,2293,0,4605,4,0,8,FR,France
200326,3,5333,2138,8528,9,4,14,FR,France
200325,3,2602,0,5228,4,0,8,FR,France
200324,3,8429,0,19483,14,0,32,FR,France
200323,3,4821,1479,8163,8,2,14,FR,France
200322,3,2573,158,4988,4,0,8,FR,France
200321,3,5563,2465,8661,9,4,14,FR,France
200320,3,7554,3748,11360,13,7,19,FR,France
200319,3,8325,3777,12873,14,6,22,FR,France
200318,3,12268,6392,18144,20,10,30,FR,France
200317,3,46253,36234,56272,77,60,94,FR,France
200316,3,50759,39835,61683,85,67,103,FR,France
200315,3,69269,55778,82760,116,93,139,FR,France
200314,3,70084,55374,84794,117,92,142,FR,France
200313,3,70335,56370,84300,118,95,141,FR,France
200312,3,88582,73568,103596,148,123,173,FR,France
200311,3,99729,83353,116105,167,140,194,FR,France
200310,3,121516,104361,138671,203,174,232,FR,France
200309,3,177255,151800,202710,296,253,339,FR,France
200308,3,222638,194269,251007,372,325,419,FR,France
200307,3,237104,207851,266357,396,347,445,FR,France
200306,3,206955,177973,235937,346,298,394,FR,France
200305,3,151864,125877,177851,254,211,297,FR,France
200304,3,95356,76211,114501,159,127,191,FR,France
200303,3,50983,38264,63702,85,64,106,FR,France
200302,3,31447,21706,41188,53,37,69,FR,France
200301,3,31228,19202,43254,52,32,72,FR,France
200252,3,47016,33692,60340,79,57,101,FR,France
200251,3,47629,36912,58346,80,62,98,FR,France
200250,3,35221,26217,44225,59,44,74,FR,France
200249,3,24656,16900,32412,41,28,54,FR,France
200248,3,19799,12868,26730,33,21,45,FR,France
200247,3,16449,9214,23684,28,16,40,FR,France
200246,3,12753,5779,19727,21,9,33,FR,France
200245,3,19304,8475,30133,32,14,50,FR,France
200244,3,18082,10671,25493,30,18,42,FR,France
200243,3,14326,7548,21104,24,13,35,FR,France
200242,3,13774,7800,19748,23,13,33,FR,France
200241,3,11036,5650,16422,19,10,28,FR,France
200240,3,8137,3375,12899,14,6,22,FR,France
200239,3,15989,9892,22086,27,17,37,FR,France
200238,3,7067,3299,10835,12,6,18,FR,France
200237,3,12409,4981,19837,21,9,33,FR,France
200236,3,1679,0,4072,3,0,7,FR,France
200235,3,2126,0,6301,4,0,11,FR,France
200234,3,2843,0,7483,5,0,13,FR,France
200233,3,2670,0,5462,4,0,9,FR,France
200232,3,1087,0,3003,2,0,5,FR,France
200231,3,5455,399,10511,9,1,17,FR,France
200230,3,3601,0,7836,6,0,13,FR,France
200229,3,5576,379,10773,9,0,18,FR,France
200228,3,6602,0,14289,11,0,24,FR,France
200227,3,1882,0,6151,3,0,10,FR,France
200226,3,4523,0,12320,8,0,21,FR,France
200225,3,4765,1076,8454,8,2,14,FR,France
200224,3,12063,7174,16952,20,12,28,FR,France
200223,3,15468,9806,21130,26,16,36,FR,France
200222,3,4854,1696,8012,8,3,13,FR,France
200221,3,5901,2212,9590,10,4,16,FR,France
200220,3,3803,1215,6391,6,2,10,FR,France
200219,3,6139,2598,9680,10,4,16,FR,France
200218,3,6083,2646,9520,10,4,16,FR,France
200217,3,7223,2678,11768,12,4,20,FR,France
200216,3,17762,9887,25637,30,17,43,FR,France
200215,3,21654,13662,29646,36,23,49,FR,France
200214,3,23877,15709,32045,40,26,54,FR,France
200213,3,27497,16289,38705,46,27,65,FR,France
200212,3,33510,20938,46082,56,35,77,FR,France
200211,3,40426,29171,51681,68,49,87,FR,France
200210,3,42915,24537,61293,72,41,103,FR,France
200209,3,63735,48368,79102,107,81,133,FR,France
200208,3,117204,95548,138860,197,161,233,FR,France
200207,3,170251,137065,203437,286,230,342,FR,France
200206,3,325084,296111,354057,546,497,595,FR,France
200205,3,456186,422559,489813,767,710,824,FR,France
200204,3,504715,470859,538571,848,791,905,FR,France
200203,3,383720,353019,414421,645,593,697,FR,France
200202,3,244028,217752,270304,410,366,454,FR,France
200201,3,115422,95696,135148,194,161,227,FR,France
200152,3,57999,43225,72773,98,73,123,FR,France
200151,3,45961,36017,55905,78,61,95,FR,France
200150,3,31309,22878,39740,53,39,67,FR,France
200149,3,34792,23268,46316,59,40,78,FR,France
200148,3,43283,20164,66402,73,34,112,FR,France
200147,3,27455,20021,34889,46,33,59,FR,France
200146,3,21127,14294,27960,36,24,48,FR,France
200145,3,20766,13979,27553,35,24,46,FR,France
200144,3,11410,6593,16227,19,11,27,FR,France
200143,3,23311,15652,30970,39,26,52,FR,France
200142,3,18963,11652,26274,32,20,44,FR,France
200141,3,20241,14032,26450,34,24,44,FR,France
200140,3,18113,9971,26255,31,17,45,FR,France
200139,3,25407,7785,43029,43,13,73,FR,France
200138,3,9736,4985,14487,16,8,24,FR,France
200137,3,8025,3527,12523,14,6,22,FR,France
200136,3,4663,654,8672,8,1,15,FR,France
200135,3,4399,0,9810,7,0,16,FR,France
200134,3,2990,0,6419,5,0,11,FR,France
200133,3,0,0,0,0,0,0,FR,France
200132,3,0,0,0,0,0,0,FR,France
200131,3,18999,4286,33712,32,7,57,FR,France
200130,3,2144,0,5670,4,0,10,FR,France
200129,3,2408,0,5215,4,0,9,FR,France
200128,3,3555,553,6557,6,1,11,FR,France
200127,3,4739,1122,8356,8,2,14,FR,France
200126,3,3335,466,6204,6,1,11,FR,France
200125,3,2643,527,4759,4,0,8,FR,France
200124,3,4784,2158,7410,8,4,12,FR,France
200123,3,6299,3161,9437,11,6,16,FR,France
200122,3,4544,2006,7082,8,4,12,FR,France
200121,3,8925,4094,13756,15,7,23,FR,France
200120,3,15269,8883,21655,26,15,37,FR,France
200119,3,17369,10575,24163,29,18,40,FR,France
200118,3,21043,14262,27824,36,25,47,FR,France
200117,3,25787,19109,32465,44,33,55,FR,France
200116,3,24025,16833,31217,41,29,53,FR,France
200115,3,24188,17340,31036,41,29,53,FR,France
200114,3,40429,31774,49084,68,53,83,FR,France
200113,3,48644,39517,57771,82,67,97,FR,France
200112,3,45411,35780,55042,77,61,93,FR,France
200111,3,49809,40534,59084,84,68,100,FR,France
200110,3,46649,37240,56058,79,63,95,FR,France
200109,3,58623,47346,69900,99,80,118,FR,France
200108,3,77423,63237,91609,131,107,155,FR,France
200107,3,126970,107898,146042,215,183,247,FR,France
200106,3,225969,198142,253796,382,335,429,FR,France
200105,3,278674,249486,307862,471,422,520,FR,France
200104,3,209193,189129,229257,354,320,388,FR,France
200103,3,139823,124243,155403,236,210,262,FR,France
200102,3,105356,92392,118320,178,156,200,FR,France
200101,3,98930,84599,113261,167,143,191,FR,France
200052,3,143525,124203,162847,244,211,277,FR,France
200051,3,125224,108719,141729,213,185,241,FR,France
200050,3,99177,84857,113497,169,145,193,FR,France
200049,3,74170,54196,94144,126,92,160,FR,France
200048,3,43642,34064,53220,74,58,90,FR,France
200047,3,64848,11322,118374,110,19,201,FR,France
200046,3,20598,14053,27143,35,24,46,FR,France
200045,3,25375,17229,33521,43,29,57,FR,France
200044,3,16330,9762,22898,28,17,39,FR,France
200043,3,36951,27693,46209,63,47,79,FR,France
200042,3,41200,30772,51628,70,52,88,FR,France
200041,3,40022,30939,49105,68,53,83,FR,France
200040,3,27237,19986,34488,46,34,58,FR,France
200039,3,19338,12649,26027,33,22,44,FR,France
200038,3,7417,3193,11641,13,6,20,FR,France
200037,3,2243,0,4722,4,0,8,FR,France
200036,3,2990,0,6416,5,0,11,FR,France
200035,3,5590,0,12003,10,0,21,FR,France
200034,3,5188,0,12418,9,0,21,FR,France
200033,3,1690,0,3859,3,0,7,FR,France
200032,3,1405,0,4239,2,0,7,FR,France
200031,3,2159,0,4744,4,0,8,FR,France
200030,3,218,0,783,0,0,1,FR,France
200029,3,1532,0,3751,3,0,7,FR,France
200028,3,1071,0,2881,2,0,5,FR,France
200027,3,2427,150,4704,4,0,8,FR,France
200026,3,1997,0,4996,3,0,8,FR,France
200025,3,4239,221,8257,7,0,14,FR,France
200024,3,5411,1122,9700,9,2,16,FR,France
200023,3,8757,4101,13413,15,7,23,FR,France
200022,3,4966,1694,8238,8,2,14,FR,France
200021,3,4271,1389,7153,7,2,12,FR,France
200020,3,4520,1478,7562,8,3,13,FR,France
200019,3,13867,8548,19186,24,15,33,FR,France
200018,3,5386,2261,8511,9,4,14,FR,France
200017,3,7018,1581,12455,12,3,21,FR,France
200016,3,5328,1433,9223,9,2,16,FR,France
200015,3,12020,6436,17604,20,11,29,FR,France
200014,3,10942,5691,16193,19,10,28,FR,France
200013,3,8572,3443,13701,15,6,24,FR,France
200012,3,7081,3097,11065,12,5,19,FR,France
200011,3,12172,7447,16897,21,13,29,FR,France
200010,3,9091,5450,12732,15,9,21,FR,France
200009,3,14631,9200,20062,25,16,34,FR,France
200008,3,23647,17172,30122,40,29,51,FR,France
200007,3,59173,48440,69906,101,83,119,FR,France
200006,3,146517,129961,163073,249,221,277,FR,France
200005,3,261687,239983,283391,445,408,482,FR,France
200004,3,374822,349785,399859,637,594,680,FR,France
200003,3,484233,456706,511760,824,777,871,FR,France
200002,3,521069,492185,549953,886,837,935,FR,France
200001,3,541945,511326,572564,922,870,974,FR,France
199952,3,387655,360652,414658,663,617,709,FR,France
199951,3,272075,250906,293244,465,429,501,FR,France
199950,3,181113,164571,197655,310,282,338,FR,France
199949,3,112328,99058,125598,192,169,215,FR,France
199948,3,58710,48651,68769,100,83,117,FR,France
199947,3,35912,23462,48362,61,40,82,FR,France
199946,3,8601,1503,15699,15,3,27,FR,France
199945,3,17781,10817,24745,30,18,42,FR,France
199944,3,18238,12870,23606,31,22,40,FR,France
199943,3,18135,12601,23669,31,22,40,FR,France
199942,3,19183,13258,25108,33,23,43,FR,France
199941,3,25171,18689,31653,43,32,54,FR,France
199940,3,24185,16184,32186,41,27,55,FR,France
199939,3,28553,20776,36330,49,36,62,FR,France
199938,3,12430,7759,17101,21,13,29,FR,France
199937,3,7543,4283,10803,13,7,19,FR,France
199936,3,8191,2827,13555,14,5,23,FR,France
199935,3,3334,1188,5480,6,2,10,FR,France
199934,3,7600,0,27730,13,0,47,FR,France
199933,3,9757,354,19160,17,1,33,FR,France
199932,3,7149,0,15077,12,0,26,FR,France
199931,3,2640,0,5463,5,0,10,FR,France
199930,3,1478,0,3153,3,0,6,FR,France
199929,3,3983,498,7468,7,1,13,FR,France
199928,3,3731,712,6750,6,1,11,FR,France
199927,3,5004,1597,8411,9,3,15,FR,France
199926,3,3459,699,6219,6,1,11,FR,France
199925,3,5302,2057,8547,9,3,15,FR,France
199924,3,4327,1646,7008,7,2,12,FR,France
199923,3,5370,2374,8366,9,4,14,FR,France
199922,3,4404,1643,7165,8,3,13,FR,France
199921,3,4180,1485,6875,7,2,12,FR,France
199920,3,6073,2775,9371,10,4,16,FR,France
199919,3,6241,3220,9262,11,6,16,FR,France
199918,3,3057,678,5436,5,1,9,FR,France
199917,3,2904,701,5107,5,1,9,FR,France
199916,3,7539,3360,11718,13,6,20,FR,France
199915,3,11535,6814,16256,20,12,28,FR,France
199914,3,18730,13131,24329,32,22,42,FR,France
199913,3,38453,29560,47346,66,51,81,FR,France
199912,3,58090,48072,68108,99,82,116,FR,France
199911,3,91329,78686,103972,156,134,178,FR,France
199910,3,147976,131058,164894,253,224,282,FR,France
199909,3,207450,188710,226190,355,323,387,FR,France
199908,3,313293,292280,334306,536,500,572,FR,France
199907,3,524136,496702,551570,896,849,943,FR,France
199906,3,493908,466605,521211,844,797,891,FR,France
199905,3,476884,449696,504072,815,769,861,FR,France
199904,3,368514,345131,391897,630,590,670,FR,France
199903,3,242388,224208,260568,414,383,445,FR,France
199902,3,160657,145209,176105,275,249,301,FR,France
199901,3,130133,115125,145141,222,196,248,FR,France
199853,3,107611,92301,122921,185,159,211,FR,France
199852,3,72373,61032,83714,124,105,143,FR,France
199851,3,78179,67156,89202,134,115,153,FR,France
199850,3,54775,45478,64072,94,78,110,FR,France
199849,3,33385,26332,40438,57,45,69,FR,France
199848,3,22636,16668,28604,39,29,49,FR,France
199847,3,14505,9746,19264,25,17,33,FR,France
199846,3,13804,9049,18559,24,16,32,FR,France
199845,3,13021,8113,17929,22,14,30,FR,France
199844,3,12631,7410,17852,22,13,31,FR,France
199843,3,12180,7125,17235,21,12,30,FR,France
199842,3,18885,12726,25044,32,21,43,FR,France
199841,3,24196,17506,30886,42,31,53,FR,France
199840,3,20235,14114,26356,35,25,45,FR,France
199839,3,17075,11538,22612,29,20,38,FR,France
199838,3,10837,6336,15338,19,11,27,FR,France
199837,3,6968,3227,10709,12,6,18,FR,France
199836,3,3033,580,5486,5,1,9,FR,France
199835,3,2491,0,5022,4,0,8,FR,France
199834,3,3853,589,7117,7,1,13,FR,France
199833,3,2934,0,6598,5,0,11,FR,France
199832,3,1308,0,2894,2,0,5,FR,France
199831,3,12147,0,29237,21,0,50,FR,France
199830,3,5309,1628,8990,9,3,15,FR,France
199829,3,6486,2465,10507,11,4,18,FR,France
199828,3,9577,4464,14690,16,7,25,FR,France
199827,3,12358,7158,17558,21,12,30,FR,France
199826,3,10415,6547,14283,18,11,25,FR,France
199825,3,18246,13050,23442,31,22,40,FR,France
199824,3,16263,11126,21400,28,19,37,FR,France
199823,3,8553,4911,12195,15,9,21,FR,France
199822,3,5181,2669,7693,9,5,13,FR,France
199821,3,4751,2416,7086,8,4,12,FR,France
199820,3,12129,8099,16159,21,14,28,FR,France
199819,3,24118,17906,30330,41,30,52,FR,France
199818,3,31612,24512,38712,54,42,66,FR,France
199817,3,61148,51894,70402,105,89,121,FR,France
199816,3,123195,109878,136512,211,188,234,FR,France
199815,3,171603,156715,186491,294,268,320,FR,France
199814,3,319092,298398,339786,547,512,582,FR,France
199813,3,298212,278513,317911,512,478,546,FR,France
199812,3,263585,246356,280814,452,422,482,FR,France
199811,3,220578,204894,236262,378,351,405,FR,France
199810,3,202354,186304,218404,347,319,375,FR,France
199809,3,227553,209970,245136,390,360,420,FR,France
199808,3,271362,252005,290719,465,432,498,FR,France
199807,3,190187,173694,206680,326,298,354,FR,France
199806,3,88217,76829,99605,151,131,171,FR,France
199805,3,45180,36970,53390,77,63,91,FR,France
199804,3,30005,23272,36738,51,39,63,FR,France
199803,3,23900,18221,29579,41,31,51,FR,France
199802,3,31143,23994,38292,53,41,65,FR,France
199801,3,24809,17786,31832,43,31,55,FR,France
199752,3,29039,21984,36094,50,38,62,FR,France
199751,3,34924,27737,42111,60,48,72,FR,France
199750,3,33527,26388,40666,58,46,70,FR,France
199749,3,15978,11147,20809,27,19,35,FR,France
199748,3,11507,7744,15270,20,14,26,FR,France
199747,3,11059,7151,14967,19,12,26,FR,France
199746,3,8073,4613,11533,14,8,20,FR,France
199745,3,18202,13113,23291,31,22,40,FR,France
199744,3,19770,14405,25135,34,25,43,FR,France
199743,3,16707,11945,21469,29,21,37,FR,France
199742,3,15340,10351,20329,26,17,35,FR,France
199741,3,16933,11941,21925,29,20,38,FR,France
199740,3,13292,8709,17875,23,15,31,FR,France
199739,3,10189,6279,14099,18,11,25,FR,France
199738,3,3231,1061,5401,6,2,10,FR,France
199737,3,4542,1949,7135,8,4,12,FR,France
199736,3,5131,2205,8057,9,4,14,FR,France
199735,3,1842,0,3740,3,0,6,FR,France
199734,3,3226,0,6717,6,0,12,FR,France
199733,3,1799,0,4368,3,0,7,FR,France
199732,3,2086,0,4411,4,0,8,FR,France
199731,3,1386,0,3332,2,0,5,FR,France
199730,3,2116,0,4919,4,0,9,FR,France
199729,3,1973,0,4347,3,0,7,FR,France
199728,3,1560,216,2904,3,1,5,FR,France
199727,3,2540,754,4326,4,1,7,FR,France
199726,3,5582,2699,8465,10,5,15,FR,France
199725,3,4012,1090,6934,7,2,12,FR,France
199724,3,4696,1824,7568,8,3,13,FR,France
199723,3,4709,2136,7282,8,4,12,FR,France
199722,3,2283,665,3901,4,1,7,FR,France
199721,3,2108,664,3552,4,2,6,FR,France
199720,3,3257,1519,4995,6,3,9,FR,France
199719,3,4646,2436,6856,8,4,12,FR,France
199718,3,5246,2764,7728,9,5,13,FR,France
199717,3,14236,9906,18566,24,17,31,FR,France
199716,3,12169,8111,16227,21,14,28,FR,France
199715,3,9516,6192,12840,16,10,22,FR,France
199714,3,13602,9659,17545,23,16,30,FR,France
199713,3,18032,13345,22719,31,23,39,FR,France
199712,3,17199,12707,21691,30,22,38,FR,France
199711,3,19988,15419,24557,34,26,42,FR,France
199710,3,21903,17338,26468,38,30,46,FR,France
199709,3,29860,24261,35459,51,41,61,FR,France
199708,3,55445,47857,63033,95,82,108,FR,France
199707,3,70876,62518,79234,122,108,136,FR,France
199706,3,76906,68152,85660,132,117,147,FR,France
199705,3,95900,86095,105705,165,148,182,FR,France
199704,3,96896,87079,106713,167,150,184,FR,France
199703,3,143844,131917,155771,248,227,269,FR,France
199702,3,211578,197395,225761,364,340,388,FR,France
199701,3,280440,263783,297097,483,454,512,FR,France
199652,3,388201,369211,407191,670,637,703,FR,France
199651,3,640627,617388,663866,1106,1066,1146,FR,France
199650,3,535205,514568,555842,924,888,960,FR,France
199649,3,411498,393335,429661,710,679,741,FR,France
199648,3,195568,182563,208573,338,316,360,FR,France
199647,3,71980,63619,80341,124,110,138,FR,France
199646,3,32138,26358,37918,55,45,65,FR,France
199645,3,25551,20136,30966,44,35,53,FR,France
199644,3,10753,7230,14276,19,13,25,FR,France
199643,3,12971,9311,16631,22,16,28,FR,France
199642,3,11929,8472,15386,21,15,27,FR,France
199641,3,9918,6674,13162,17,11,23,FR,France
199640,3,10765,7067,14463,19,13,25,FR,France
199639,3,8377,5183,11571,14,8,20,FR,France
199638,3,9470,5965,12975,16,10,22,FR,France
199637,3,2925,829,5021,5,1,9,FR,France
199636,3,1122,23,2221,2,0,4,FR,France
199635,3,1674,93,3255,3,0,6,FR,France
199634,3,1273,122,2424,2,0,4,FR,France
199633,3,540,0,1246,1,0,2,FR,France
199632,3,1453,176,2730,3,1,5,FR,France
199631,3,3010,752,5268,5,1,9,FR,France
199630,3,2992,706,5278,5,1,9,FR,France
199629,3,1559,144,2974,3,1,5,FR,France
199628,3,3109,541,5677,5,1,9,FR,France
199627,3,4456,1595,7317,8,3,13,FR,France
199626,3,2871,1138,4604,5,2,8,FR,France
199625,3,2308,856,3760,4,1,7,FR,France
199624,3,2490,1067,3913,4,2,6,FR,France
199623,3,2819,1262,4376,5,2,8,FR,France
199622,3,4005,1984,6026,7,4,10,FR,France
199621,3,7042,4291,9793,12,7,17,FR,France
199620,3,10619,7365,13873,18,12,24,FR,France
199619,3,11898,8065,15731,21,14,28,FR,France
199618,3,18338,13352,23324,32,23,41,FR,France
199617,3,24954,19100,30808,43,33,53,FR,France
199616,3,29313,23868,34758,51,42,60,FR,France
199615,3,51703,43644,59762,89,75,103,FR,France
199614,3,44414,36774,52054,77,64,90,FR,France
199613,3,30622,25062,36182,53,43,63,FR,France
199612,3,31818,25852,37784,55,45,65,FR,France
199611,3,24427,19499,29355,42,33,51,FR,France
199610,3,20910,16198,25622,36,28,44,FR,France
199609,3,19594,15228,23960,34,26,42,FR,France
199608,3,19046,14850,23242,33,26,40,FR,France
199607,3,23041,18338,27744,40,32,48,FR,France
199606,3,21810,17153,26467,38,30,46,FR,France
199605,3,25275,20410,30140,44,36,52,FR,France
199604,3,43359,36948,49770,75,64,86,FR,France
199603,3,54691,48016,61366,94,82,106,FR,France
199602,3,105606,96106,115106,182,166,198,FR,France
199601,3,223830,209274,238386,386,361,411,FR,France
199552,3,397122,378252,415992,688,655,721,FR,France
199551,3,749988,725434,774542,1299,1256,1342,FR,France
199550,3,604926,583100,626752,1047,1009,1085,FR,France
199549,3,375500,358200,392800,650,620,680,FR,France
199548,3,206446,193457,219435,357,335,379,FR,France
199547,3,120472,110526,130418,209,192,226,FR,France
199546,3,62380,55042,69718,108,95,121,FR,France
199545,3,29791,24523,35059,52,43,61,FR,France
199544,3,13168,9820,16516,23,17,29,FR,France
199543,3,16188,12389,19987,28,21,35,FR,France
199542,3,15178,11387,18969,26,19,33,FR,France
199541,3,14980,10851,19109,26,19,33,FR,France
199540,3,16245,11610,20880,28,20,36,FR,France
199539,3,18499,13844,23154,32,24,40,FR,France
199538,3,13123,9410,16836,23,17,29,FR,France
199537,3,4895,2666,7124,8,4,12,FR,France
199536,3,3659,1574,5744,6,2,10,FR,France
199535,3,2496,644,4348,4,1,7,FR,France
199534,3,1887,214,3560,3,0,6,FR,France
199533,3,905,0,2007,2,0,4,FR,France
199532,3,700,0,1699,1,0,3,FR,France
199531,3,1946,92,3800,3,0,6,FR,France
199530,3,1898,46,3750,3,0,6,FR,France
199529,3,2148,687,3609,4,1,7,FR,France
199528,3,3312,1430,5194,6,3,9,FR,France
199527,3,4853,2683,7023,8,4,12,FR,France
199526,3,4986,2911,7061,9,5,13,FR,France
199525,3,4571,2520,6622,8,4,12,FR,France
199524,3,8058,5054,11062,14,9,19,FR,France
199523,3,13820,9801,17839,24,17,31,FR,France
199522,3,14124,9999,18249,24,17,31,FR,France
199521,3,11931,8655,15207,21,15,27,FR,France
199520,3,18846,14311,23381,33,25,41,FR,France
199519,3,26261,21488,31034,45,37,53,FR,France
199518,3,39231,33258,45204,68,58,78,FR,France
199517,3,88640,79343,97937,153,137,169,FR,France
199516,3,140488,129439,151537,243,224,262,FR,France
199515,3,181007,169110,192904,313,292,334,FR,France
199514,3,248704,235144,262264,431,408,454,FR,France
199513,3,209839,197191,222487,363,341,385,FR,France
199512,3,122434,112497,132371,212,195,229,FR,France
199511,3,81730,73346,90114,142,127,157,FR,France
199510,3,49387,42703,56071,86,74,98,FR,France
199509,3,39696,33259,46133,69,58,80,FR,France
199508,3,38008,32225,43791,66,56,76,FR,France
199507,3,30550,25479,35621,53,44,62,FR,France
199506,3,37351,32007,42695,65,56,74,FR,France
199505,3,27876,23218,32534,48,40,56,FR,France
199504,3,27643,22818,32468,48,40,56,FR,France
199503,3,24844,20391,29297,43,35,51,FR,France
199502,3,23480,19249,27711,41,34,48,FR,France
199501,3,26556,21916,31196,46,38,54,FR,France
199452,3,22036,17495,26577,38,30,46,FR,France
199451,3,26912,21750,32074,47,38,56,FR,France
199450,3,28831,23470,34192,50,41,59,FR,France
199449,3,21161,16559,25763,37,29,45,FR,France
199448,3,18487,14456,22518,32,25,39,FR,France
199447,3,15206,11578,18834,26,20,32,FR,France
199446,3,14252,10525,17979,25,19,31,FR,France
199445,3,12152,8448,15856,21,15,27,FR,France
199444,3,11307,7688,14926,20,14,26,FR,France
199443,3,11254,8002,14506,20,14,26,FR,France
199442,3,14406,10653,18159,25,18,32,FR,France
199441,3,14073,10341,17805,24,18,30,FR,France
199440,3,15962,11850,20074,28,21,35,FR,France
199439,3,15878,11643,20113,28,21,35,FR,France
199438,3,15366,11100,19632,27,20,34,FR,France
199437,3,5949,3561,8337,10,6,14,FR,France
199436,3,2806,1133,4479,5,2,8,FR,France
199435,3,4556,2358,6754,8,4,12,FR,France
199434,3,3459,1444,5474,6,2,10,FR,France
199433,3,4272,977,7567,7,1,13,FR,France
199432,3,4786,1613,7959,8,2,14,FR,France
199431,3,5027,1580,8474,9,3,15,FR,France
199430,3,2100,520,3680,4,1,7,FR,France
199429,3,1188,0,2384,2,0,4,FR,France
199428,3,1581,316,2846,3,1,5,FR,France
199427,3,713,0,1459,1,0,2,FR,France
199426,3,990,80,1900,2,0,4,FR,France
199425,3,5431,3222,7640,9,5,13,FR,France
199424,3,6288,3723,8853,11,7,15,FR,France
199423,3,6440,1052,11828,11,2,20,FR,France
199422,3,845,0,2449,1,0,4,FR,France
199421,3,2197,395,3999,4,1,7,FR,France
199420,3,2217,852,3582,4,2,6,FR,France
199419,3,1483,406,2560,3,1,5,FR,France
199418,3,3705,1940,5470,6,3,9,FR,France
199417,3,2288,1064,3512,4,2,6,FR,France
199416,3,5156,3176,7136,9,6,12,FR,France
199415,3,6952,4678,9226,12,8,16,FR,France
199414,3,4671,2784,6558,8,5,11,FR,France
199413,3,3343,1641,5045,6,3,9,FR,France
199412,3,4305,2412,6198,7,4,10,FR,France
199411,3,5399,3327,7471,9,5,13,FR,France
199410,3,4548,2529,6567,8,4,12,FR,France
199409,3,5635,3251,8019,10,6,14,FR,France
199408,3,9539,6652,12426,17,12,22,FR,France
199407,3,14639,11188,18090,25,19,31,FR,France
199406,3,17308,13395,21221,30,23,37,FR,France
199405,3,18489,14290,22688,32,25,39,FR,France
199404,3,23712,19048,28376,41,33,49,FR,France
199403,3,29882,24817,34947,52,43,61,FR,France
199402,3,47402,41415,53389,82,72,92,FR,France
199401,3,118119,108182,128056,205,188,222,FR,France
199352,3,233137,218608,247666,406,381,431,FR,France
199351,3,377472,359792,395152,658,627,689,FR,France
199350,3,627639,607097,648181,1094,1058,1130,FR,France
199349,3,898045,873419,922671,1565,1522,1608,FR,France
199348,3,527582,507752,547412,920,885,955,FR,France
199347,3,241371,226769,255973,421,396,446,FR,France
199346,3,76034,68003,84065,133,119,147,FR,France
199345,3,32289,27022,37556,56,47,65,FR,France
199344,3,17832,14015,21649,31,24,38,FR,France
199343,3,15098,11333,18863,26,19,33,FR,France
199342,3,19187,14893,23481,33,26,40,FR,France
199341,3,18285,14384,22186,32,25,39,FR,France
199340,3,13805,10348,17262,24,18,30,FR,France
199339,3,14783,11173,18393,26,20,32,FR,France
199338,3,10492,7329,13655,18,12,24,FR,France
199337,3,7023,4503,9543,12,8,16,FR,France
199336,3,4361,2445,6277,8,5,11,FR,France
199335,3,5589,3066,8112,10,6,14,FR,France
199334,3,2984,1209,4759,5,2,8,FR,France
199333,3,3217,758,5676,6,2,10,FR,France
199332,3,3755,1051,6459,7,2,12,FR,France
199331,3,4163,1768,6558,7,3,11,FR,France
199330,3,4055,1907,6203,7,3,11,FR,France
199329,3,3714,1668,5760,6,2,10,FR,France
199328,3,1774,484,3064,3,1,5,FR,France
199327,3,2574,918,4230,4,1,7,FR,France
199326,3,3048,1450,4646,5,2,8,FR,France
199325,3,8348,5260,11436,15,10,20,FR,France
199324,3,6624,3747,9501,12,7,17,FR,France
199323,3,7424,4302,10546,13,8,18,FR,France
199322,3,3432,1849,5015,6,3,9,FR,France
199321,3,5942,3670,8214,10,6,14,FR,France
199320,3,5137,3157,7117,9,6,12,FR,France
199319,3,4593,2725,6461,8,5,11,FR,France
199318,3,5027,3075,6979,9,6,12,FR,France
199317,3,8111,5138,11084,14,9,19,FR,France
199316,3,24932,19468,30396,43,33,53,FR,France
199315,3,22181,17147,27215,39,30,48,FR,France
199314,3,38266,30865,45667,67,54,80,FR,France
199313,3,62987,54506,71468,110,95,125,FR,France
199312,3,63493,55956,71030,111,98,124,FR,France
199311,3,80548,72628,88468,140,126,154,FR,France
199310,3,127910,117688,138132,223,205,241,FR,France
199309,3,184150,171378,196922,321,299,343,FR,France
199308,3,222932,209051,236813,389,365,413,FR,France
199307,3,248494,234936,262052,433,409,457,FR,France
199306,3,286939,271830,302048,500,474,526,FR,France
199305,3,255890,241516,270264,446,421,471,FR,France
199304,3,193107,181261,204953,337,316,358,FR,France
199303,3,147854,137256,158452,258,240,276,FR,France
199302,3,90269,82164,98374,157,143,171,FR,France
199301,3,75288,67244,83332,131,117,145,FR,France
199253,3,56622,49107,64137,99,86,112,FR,France
199252,3,49588,42568,56608,87,75,99,FR,France
199251,3,41156,35024,47288,72,61,83,FR,France
199250,3,40100,34699,45501,70,61,79,FR,France
199249,3,39146,33847,44445,69,60,78,FR,France
199248,3,32392,27499,37285,57,48,66,FR,France
199247,3,26276,21797,30755,46,38,54,FR,France
199246,3,22758,18462,27054,40,32,48,FR,France
199245,3,22023,17804,26242,39,32,46,FR,France
199244,3,27619,22945,32293,48,40,56,FR,France
199243,3,24969,20704,29234,44,37,51,FR,France
199242,3,29980,25261,34699,52,44,60,FR,France
199241,3,24091,19953,28229,42,35,49,FR,France
199240,3,17593,13885,21301,31,25,37,FR,France
199239,3,13535,10301,16769,24,18,30,FR,France
199238,3,9837,7106,12568,17,12,22,FR,France
199237,3,7480,5130,9830,13,9,17,FR,France
199236,3,5155,2958,7352,9,5,13,FR,France
199235,3,6383,3687,9079,11,6,16,FR,France
199234,3,2953,922,4984,5,1,9,FR,France
199233,3,4046,1549,6543,7,3,11,FR,France
199232,3,4963,2215,7711,9,4,14,FR,France
199231,3,4233,1733,6733,7,3,11,FR,France
199230,3,3091,1145,5037,5,2,8,FR,France
199229,3,3586,1583,5589,6,2,10,FR,France
199228,3,4179,2133,6225,7,3,11,FR,France
199227,3,5568,3134,8002,10,6,14,FR,France
199226,3,5886,3229,8543,10,5,15,FR,France
199225,3,5304,2805,7803,9,5,13,FR,France
199224,3,5928,3482,8374,10,6,14,FR,France
199223,3,3914,1855,5973,7,3,11,FR,France
199222,3,4899,2641,7157,9,5,13,FR,France
199221,3,5059,2889,7229,9,5,13,FR,France
199220,3,4308,2226,6390,8,4,12,FR,France
199219,3,2982,1181,4783,5,2,8,FR,France
199218,3,3887,1777,5997,7,3,11,FR,France
199217,3,3614,1216,6012,6,2,10,FR,France
199216,3,5460,2934,7986,10,6,14,FR,France
199215,3,8642,5029,12255,15,9,21,FR,France
199214,3,11346,7536,15156,20,13,27,FR,France
199213,3,11693,8096,15290,20,14,26,FR,France
199212,3,11752,7834,15670,21,14,28,FR,France
199211,3,11689,7575,15803,20,13,27,FR,France
199210,3,19224,12971,25477,34,23,45,FR,France
199209,3,26629,20275,32983,47,36,58,FR,France
199208,3,37667,30887,44447,66,54,78,FR,France
199207,3,57312,48777,65847,100,85,115,FR,France
199206,3,112450,99989,124911,197,175,219,FR,France
199205,3,165474,151342,179606,290,265,315,FR,France
199204,3,147959,133595,162323,259,234,284,FR,France
199203,3,137893,122944,152842,241,215,267,FR,France
199202,3,185852,166558,205146,325,291,359,FR,France
199201,3,247688,222574,272802,434,390,478,FR,France
199152,3,349815,325499,374131,615,572,658,FR,France
199151,3,378413,355928,400898,666,626,706,FR,France
199150,3,198721,182878,214564,350,322,378,FR,France
199149,3,93903,82992,104814,165,146,184,FR,France
199148,3,45173,37384,52962,79,65,93,FR,France
199147,3,35691,28224,43158,63,50,76,FR,France
199146,3,23685,17531,29839,42,31,53,FR,France
199145,3,19565,13551,25579,34,23,45,FR,France
199144,3,21571,15477,27665,38,27,49,FR,France
199143,3,28842,21958,35726,51,39,63,FR,France
199142,3,27483,20385,34581,48,36,60,FR,France
199141,3,28185,20655,35715,50,37,63,FR,France
199140,3,26740,16049,37431,47,28,66,FR,France
199139,3,14373,8734,20012,25,15,35,FR,France
199138,3,7460,3700,11220,13,6,20,FR,France
199137,3,2429,434,4424,4,0,8,FR,France
199136,3,2380,0,4882,4,0,8,FR,France
199135,3,2372,0,5025,4,0,9,FR,France
199134,3,2099,0,5104,4,0,9,FR,France
199133,3,1654,0,4102,3,0,7,FR,France
199132,3,625,0,2190,1,0,4,FR,France
199131,3,2464,0,5830,4,0,10,FR,France
199130,3,2382,0,6857,4,0,12,FR,France
199129,3,3837,40,7634,7,0,14,FR,France
199128,3,3297,152,6442,6,0,12,FR,France
199127,3,3164,934,5394,6,2,10,FR,France
199126,3,6376,3324,9428,11,6,16,FR,France
199125,3,8218,4776,11660,14,8,20,FR,France
199124,3,8223,4660,11786,14,8,20,FR,France
199123,3,6801,3596,10006,12,6,18,FR,France
199122,3,6050,2682,9418,11,5,17,FR,France
199121,3,6391,3385,9397,11,6,16,FR,France
199120,3,8633,4683,12583,15,8,22,FR,France
199119,3,7785,3986,11584,14,7,21,FR,France
199118,3,10856,6457,15255,19,11,27,FR,France
199117,3,15603,10950,20256,27,19,35,FR,France
199116,3,15665,11184,20146,28,20,36,FR,France
199115,3,15062,10715,19409,26,18,34,FR,France
199114,3,27862,21583,34141,49,38,60,FR,France
199113,3,32179,25829,38529,57,46,68,FR,France
199112,3,27873,22312,33434,49,39,59,FR,France
199111,3,46186,38876,53496,81,68,94,FR,France
199110,3,88750,77784,99716,156,137,175,FR,France
199109,3,139093,124326,153860,245,219,271,FR,France
199108,3,216369,198630,234108,381,350,412,FR,France
199107,3,208836,191789,225883,367,337,397,FR,France
199106,3,134737,121280,148194,237,213,261,FR,France
199105,3,73564,62716,84412,129,110,148,FR,France
199104,3,48253,39031,57475,85,69,101,FR,France
199103,3,29382,22258,36506,52,39,65,FR,France
199102,3,29759,22044,37474,52,38,66,FR,France
199101,3,42927,33706,52148,76,60,92,FR,France
199052,3,41262,31607,50917,73,56,90,FR,France
199051,3,46924,37544,56304,83,66,100,FR,France
199050,3,44125,36137,52113,78,64,92,FR,France
199049,3,38379,30603,46155,68,54,82,FR,France
199048,3,29571,22315,36827,52,39,65,FR,France
199047,3,23358,15890,30826,41,28,54,FR,France
199046,3,20246,7930,32562,36,14,58,FR,France
199045,3,17571,11591,23551,31,20,42,FR,France
199044,3,16182,10305,22059,29,19,39,FR,France
199043,3,17579,11588,23570,31,20,42,FR,France
199042,3,13011,7875,18147,23,14,32,FR,France
199041,3,19723,13330,26116,35,24,46,FR,France
199040,3,20433,14614,26252,36,26,46,FR,France
199039,3,15618,10315,20921,28,19,37,FR,France
199038,3,11146,6197,16095,20,11,29,FR,France
199037,3,4051,1022,7080,7,2,12,FR,France
199036,3,2542,124,4960,4,0,8,FR,France
199035,3,934,0,2681,2,0,5,FR,France
199034,3,201,0,595,0,0,1,FR,France
199033,3,0,0,0,0,0,0,FR,France
199032,3,1399,0,3849,2,0,6,FR,France
199031,3,881,0,2206,2,0,4,FR,France
199030,3,425,0,1262,1,0,2,FR,France
199029,3,1077,0,2643,2,0,5,FR,France
199028,3,959,0,2252,2,0,4,FR,France
199027,3,2948,793,5103,5,1,9,FR,France
199026,3,7236,2651,11821,13,5,21,FR,France
199025,3,7879,3136,12622,14,6,22,FR,France
199024,3,5727,2460,8994,10,4,16,FR,France
199023,3,3655,1043,6267,6,1,11,FR,France
199022,3,4196,1292,7100,7,2,12,FR,France
199021,3,1004,0,2367,2,0,4,FR,France
199020,3,4506,1771,7241,8,3,13,FR,France
199019,3,5897,2440,9354,10,4,16,FR,France
199018,3,1738,141,3335,3,0,6,FR,France
199017,3,5278,2197,8359,9,4,14,FR,France
199016,3,8129,3754,12504,14,6,22,FR,France
199015,3,11707,5831,17583,21,11,31,FR,France
199014,3,12862,8177,17547,23,15,31,FR,France
199013,3,21879,15088,28670,39,27,51,FR,France
199012,3,17035,11135,22935,30,20,40,FR,France
199011,3,20018,14311,25725,35,25,45,FR,France
199010,3,16294,11045,21543,29,20,38,FR,France
199009,3,21793,15460,28126,39,28,50,FR,France
199008,3,26110,19342,32878,46,34,58,FR,France
199007,3,61742,49720,73764,109,88,130,FR,France
199006,3,399984,353345,446623,707,625,789,FR,France
199005,3,161765,146354,177176,286,259,313,FR,France
199004,3,287771,267698,307844,509,474,544,FR,France
199003,3,315466,296550,334382,558,525,591,FR,France
199002,3,345060,326087,364033,610,576,644,FR,France
199001,3,483600,460913,506287,855,815,895,FR,France
198952,3,581149,556023,606275,1035,990,1080,FR,France
198951,3,821540,793846,849234,1463,1414,1512,FR,France
198950,3,749283,723633,774933,1334,1288,1380,FR,France
198949,3,358043,340053,376033,638,606,670,FR,France
198948,3,119087,108253,129921,212,193,231,FR,France
198947,3,40963,33897,48029,73,60,86,FR,France
198946,3,23547,17952,29142,42,32,52,FR,France
198945,3,22711,16999,28423,40,30,50,FR,France
198944,3,23405,17475,29335,42,31,53,FR,France
198943,3,29549,23287,35811,53,42,64,FR,France
198942,3,39410,32172,46648,70,57,83,FR,France
198941,3,43763,35245,52281,78,63,93,FR,France
198940,3,31453,24288,38618,56,43,69,FR,France
198939,3,24622,18425,30819,44,33,55,FR,France
198938,3,18296,13483,23109,33,24,42,FR,France
198937,3,12022,8082,15962,21,14,28,FR,France
198936,3,4769,2299,7239,8,4,12,FR,France
198935,3,4661,1696,7626,8,3,13,FR,France
198934,3,6039,2185,9893,11,4,18,FR,France
198933,3,7940,3858,12022,14,7,21,FR,France
198932,3,7114,2361,11867,13,5,21,FR,France
198931,3,2721,0,5478,5,0,10,FR,France
198930,3,3168,561,5775,6,1,11,FR,France
198929,3,3909,1504,6314,7,3,11,FR,France
198928,3,1284,105,2463,2,0,4,FR,France
198927,3,3641,1350,5932,6,2,10,FR,France
198926,3,5753,2885,8621,10,5,15,FR,France
198925,3,6299,3454,9144,11,6,16,FR,France
198924,3,6337,3537,9137,11,6,16,FR,France
198923,3,11779,7984,15574,21,14,28,FR,France
198922,3,9104,4764,13444,16,8,24,FR,France
198921,3,9554,4316,14792,17,8,26,FR,France
198920,3,6093,1023,11163,11,2,20,FR,France
198919,3,0,,,0,,,FR,France
198918,3,13452,1365,25539,24,2,46,FR,France
198917,3,10949,7147,14751,19,12,26,FR,France
198916,3,10908,7494,14322,19,13,25,FR,France
198915,3,9667,6209,13125,17,11,23,FR,France
198914,3,15230,10469,19991,27,19,35,FR,France
198913,3,10564,6767,14361,19,12,26,FR,France
198912,3,13419,9542,17296,24,17,31,FR,France
198911,3,16633,12140,21126,30,22,38,FR,France
198910,3,12387,8750,16024,22,16,28,FR,France
198909,3,12895,9183,16607,23,16,30,FR,France
198908,3,19746,15258,24234,35,27,43,FR,France
198907,3,27277,21760,32794,49,39,59,FR,France
198906,3,47844,40320,55368,85,72,98,FR,France
198905,3,52840,45665,60015,94,81,107,FR,France
198904,3,54942,47750,62134,98,85,111,FR,France
198903,3,66010,58065,73955,118,104,132,FR,France
198902,3,112284,102003,122565,200,182,218,FR,France
198901,3,202266,188305,216227,360,335,385,FR,France
198852,3,467971,447033,488909,837,800,874,FR,France
198851,3,872748,846468,899028,1562,1515,1609,FR,France
198850,3,1001824,974799,1028849,1793,1745,1841,FR,France
198849,3,966523,940068,992978,1729,1682,1776,FR,France
198848,3,637811,614785,660837,1141,1100,1182,FR,France
198847,3,246062,231234,260890,440,413,467,FR,France
198846,3,92484,83160,101808,165,148,182,FR,France
198845,3,53940,46148,61732,97,83,111,FR,France
198844,3,42106,35336,48876,75,63,87,FR,France
198843,3,38252,31909,44595,68,57,79,FR,France
198842,3,45495,38467,52523,81,68,94,FR,France
198841,3,46336,39344,53328,83,70,96,FR,France
198840,3,39882,33507,46257,71,60,82,FR,France
198839,3,43730,36755,50705,78,66,90,FR,France
198838,3,35243,29367,41119,63,52,74,FR,France
198837,3,15537,11911,19163,28,22,34,FR,France
198836,3,7539,4840,10238,13,8,18,FR,France
198835,3,7797,4762,10832,14,9,19,FR,France
198834,3,4950,2005,7895,9,4,14,FR,France
198833,3,7937,3861,12013,14,7,21,FR,France
198832,3,7033,2595,11471,13,5,21,FR,France
198831,3,8758,3591,13925,16,7,25,FR,France
198830,3,4205,1224,7186,8,3,13,FR,France
198829,3,6103,1013,11193,11,2,20,FR,France
198828,3,8562,3483,13641,15,6,24,FR,France
198827,3,4690,2196,7184,8,4,12,FR,France
198826,3,7065,3905,10225,13,7,19,FR,France
198825,3,7657,4664,10650,14,9,19,FR,France
198824,3,6676,3978,9374,12,7,17,FR,France
198823,3,8695,5755,11635,16,11,21,FR,France
198822,3,11566,8282,14850,21,15,27,FR,France
198821,3,11107,7742,14472,20,14,26,FR,France
198820,3,14118,10223,18013,25,18,32,FR,France
198819,3,18260,13847,22673,33,25,41,FR,France
198818,3,20415,16281,24549,37,30,44,FR,France
198817,3,21750,17436,26064,39,31,47,FR,France
198816,3,33730,27496,39964,60,49,71,FR,France
198815,3,70616,61754,79478,126,110,142,FR,France
198814,3,123868,113211,134525,222,203,241,FR,France
198813,3,193710,181519,205901,347,325,369,FR,France
198812,3,283445,268834,298056,507,481,533,FR,France
198811,3,316175,300858,331492,566,539,593,FR,France
198810,3,235142,222004,248280,421,397,445,FR,France
198809,3,135564,125263,145865,243,225,261,FR,France
198808,3,87753,78449,97057,157,140,174,FR,France
198807,3,77099,68423,85775,138,122,154,FR,France
198806,3,69349,61502,77196,124,110,138,FR,France
198805,3,72705,64952,80458,130,116,144,FR,France
198804,3,72818,64944,80692,130,116,144,FR,France
198803,3,68869,60922,76816,123,109,137,FR,France
198802,3,59743,51907,67579,107,93,121,FR,France
198801,3,60929,52573,69285,109,94,124,FR,France
198753,3,47257,38170,56344,85,69,101,FR,France
198752,3,11935,2353,21517,21,4,38,FR,France
198751,3,47744,34266,61222,86,62,110,FR,France
198750,3,57767,50181,65353,104,90,118,FR,France
198749,3,49380,42461,56299,89,77,101,FR,France
198748,3,41871,35115,48627,75,63,87,FR,France
198747,3,29688,23465,35911,53,42,64,FR,France
198746,3,37216,30199,44233,67,54,80,FR,France
198745,3,35456,29024,41888,64,52,76,FR,France
198744,3,41024,33829,48219,74,61,87,FR,France
198743,3,46356,38370,54342,83,69,97,FR,France
198742,3,38021,30236,45806,68,54,82,FR,France
198741,3,42021,32016,52026,76,58,94,FR,France
198740,3,24641,15930,33352,44,28,60,FR,France
198739,3,32328,23548,41108,58,42,74,FR,France
198738,3,16782,11870,21694,30,21,39,FR,France
198737,3,11208,6542,15874,20,12,28,FR,France
198736,3,7906,4014,11798,14,7,21,FR,France
198735,3,6149,2970,9328,11,5,17,FR,France
198734,3,6631,2734,10528,12,5,19,FR,France
198733,3,6404,1834,10974,12,4,20,FR,France
198732,3,7215,2353,12077,13,4,22,FR,France
198731,3,8233,3655,12811,15,7,23,FR,France
198730,3,4642,1301,7983,8,2,14,FR,France
198729,3,4389,913,7865,8,2,14,FR,France
198728,3,2702,138,5266,5,0,10,FR,France
198727,3,10171,5904,14438,18,10,26,FR,France
198726,3,12341,7824,16858,22,14,30,FR,France
198725,3,23413,15490,31336,42,28,56,FR,France
198724,3,11289,7282,15296,20,13,27,FR,France
198723,3,15857,10943,20771,29,20,38,FR,France
198722,3,14887,9743,20031,27,18,36,FR,France
198721,3,17619,12649,22589,32,23,41,FR,France
198720,3,10344,6454,14234,19,12,26,FR,France
198719,3,12539,7943,17135,23,15,31,FR,France
198718,3,16095,10330,21860,29,19,39,FR,France
198717,3,21256,14850,27662,38,26,50,FR,France
198716,3,32416,25100,39732,58,45,71,FR,France
198715,3,44292,36456,52128,80,66,94,FR,France
198714,3,50642,41866,59418,91,75,107,FR,France
198713,3,62246,52749,71743,112,95,129,FR,France
198712,3,59625,50653,68597,107,91,123,FR,France
198711,3,53561,45035,62087,96,81,111,FR,France
198710,3,93067,80519,105615,167,144,190,FR,France
198709,3,147006,133519,160493,264,240,288,FR,France
198708,3,217779,202179,233379,391,363,419,FR,France
198707,3,288180,270718,305642,518,487,549,FR,France
198706,3,296702,279051,314353,533,501,565,FR,France
198705,3,236046,219483,252609,424,394,454,FR,France
198704,3,158333,143625,173041,285,259,311,FR,France
198703,3,103841,91699,115983,187,165,209,FR,France
198702,3,74213,64014,84412,133,115,151,FR,France
198701,3,69745,59207,80283,125,106,144,FR,France
198652,3,76263,64958,87568,138,118,158,FR,France
198651,3,72920,63047,82793,132,114,150,FR,France
198650,3,73413,64019,82807,133,116,150,FR,France
198649,3,56296,48229,64363,102,87,117,FR,France
198648,3,46031,38578,53484,83,70,96,FR,France
198647,3,36740,29710,43770,66,53,79,FR,France
198646,3,32134,25327,38941,58,46,70,FR,France
198645,3,38964,31382,46546,70,56,84,FR,France
198644,3,28646,22048,35244,52,40,64,FR,France
198643,3,27532,20842,34222,50,38,62,FR,France
198642,3,31815,25124,38506,57,45,69,FR,France
198641,3,32864,25939,39789,59,47,71,FR,France
198640,3,33516,26776,40256,60,48,72,FR,France
198639,3,28665,22477,34853,52,41,63,FR,France
198638,3,26615,20479,32751,48,37,59,FR,France
198637,3,15998,10802,21194,29,20,38,FR,France
198636,3,9182,4982,13382,17,9,25,FR,France
198635,3,10587,4470,16704,19,8,30,FR,France
198634,3,5762,858,10666,10,1,19,FR,France
198633,3,3115,554,5676,6,1,11,FR,France
198632,3,1027,0,2311,2,0,4,FR,France
198631,3,2247,0,4574,4,0,8,FR,France
198630,3,2346,141,4551,4,0,8,FR,France
198629,3,3654,931,6377,7,2,12,FR,France
198628,3,2327,0,5318,4,0,9,FR,France
198627,3,7529,3353,11705,14,6,22,FR,France
198626,3,8778,5106,12450,16,9,23,FR,France
198625,3,11075,6804,15346,20,12,28,FR,France
198624,3,11160,7252,15068,20,13,27,FR,France
198623,3,14372,9896,18848,26,18,34,FR,France
198622,3,9717,5594,13840,18,11,25,FR,France
198621,3,14460,9575,19345,26,17,35,FR,France
198620,3,17133,12106,22160,31,22,40,FR,France
198619,3,22177,16008,28346,40,29,51,FR,France
198618,3,31736,24115,39357,57,43,71,FR,France
198617,3,34874,27440,42308,63,50,76,FR,France
198616,3,41196,32894,49498,74,59,89,FR,France
198615,3,42659,33537,51781,77,61,93,FR,France
198614,3,49792,39728,59856,90,72,108,FR,France
198613,3,87249,63562,110936,157,114,200,FR,France
198612,3,124777,97237,152317,225,175,275,FR,France
198611,3,168146,151035,185257,303,272,334,FR,France
198610,3,257305,231225,283385,464,417,511,FR,France
198609,3,321140,296644,345636,580,536,624,FR,France
198608,3,448786,419219,478353,810,757,863,FR,France
198607,3,490723,460150,521296,886,831,941,FR,France
198606,3,439300,411764,466836,793,743,843,FR,France
198605,3,364727,339868,389586,658,613,703,FR,France
198604,3,348509,322464,374554,629,582,676,FR,France
198603,3,271539,248654,294424,490,449,531,FR,France
198602,3,178138,159344,196932,322,288,356,FR,France
198601,3,112614,96539,128689,203,174,232,FR,France
198552,3,86452,72250,100654,157,131,183,FR,France
198551,3,114666,99538,129794,208,181,235,FR,France
198550,3,86031,73604,98458,156,133,179,FR,France
198549,3,65931,55326,76536,120,101,139,FR,France
198548,3,68291,57581,79001,124,105,143,FR,France
198547,3,66177,55221,77133,120,100,140,FR,France
198546,3,70534,58788,82280,128,107,149,FR,France
198545,3,57514,47167,67861,104,85,123,FR,France
198544,3,67356,55744,78968,122,101,143,FR,France
198543,3,87402,74208,100596,158,134,182,FR,France
198542,3,76218,64604,87832,138,117,159,FR,France
198541,3,45712,36755,54669,83,67,99,FR,France
198540,3,49739,40335,59143,90,73,107,FR,France
198539,3,39819,31234,48404,72,56,88,FR,France
198538,3,27605,20614,34596,50,37,63,FR,France
198537,3,22219,15626,28812,40,28,52,FR,France
198536,3,27986,18992,36980,51,35,67,FR,France
198535,3,30215,18546,41884,55,34,76,FR,France
198534,3,40389,24786,55992,73,45,101,FR,France
198533,3,22611,8819,36403,41,16,66,FR,France
198532,3,7488,1851,13125,14,4,24,FR,France
198531,3,16958,8479,25437,31,16,46,FR,France
198530,3,11598,5507,17689,21,10,32,FR,France
198529,3,13054,6474,19634,24,12,36,FR,France
198528,3,14588,7659,21517,26,13,39,FR,France
198527,3,19670,11761,27579,36,22,50,FR,France
198526,3,18609,12637,24581,34,23,45,FR,France
198525,3,19362,12454,26270,35,22,48,FR,France
198524,3,19855,13577,26133,36,25,47,FR,France
198523,3,19373,10010,28736,35,18,52,FR,France
198522,3,24099,17190,31008,44,31,57,FR,France
198521,3,26096,19621,32571,47,35,59,FR,France
198520,3,27896,20885,34907,51,38,64,FR,France
198519,3,43154,32821,53487,78,59,97,FR,France
198518,3,40555,29935,51175,74,55,93,FR,France
198517,3,34053,24366,43740,62,44,80,FR,France
198516,3,50362,36451,64273,91,66,116,FR,France
198515,3,63881,45538,82224,116,83,149,FR,France
198514,3,134545,114400,154690,244,207,281,FR,France
198513,3,197206,176080,218332,357,319,395,FR,France
198512,3,245240,223304,267176,445,405,485,FR,France
198511,3,276205,252399,300011,501,458,544,FR,France
198510,3,353231,326279,380183,640,591,689,FR,France
198509,3,369895,341109,398681,670,618,722,FR,France
198508,3,389886,359529,420243,707,652,762,FR,France
198507,3,471852,432599,511105,855,784,926,FR,France
198506,3,565825,518011,613639,1026,939,1113,FR,France
198505,3,637302,592795,681809,1155,1074,1236,FR,France
198504,3,424937,390794,459080,770,708,832,FR,France
198503,3,213901,174689,253113,388,317,459,FR,France
198502,3,97586,80949,114223,177,147,207,FR,France
198501,3,85489,65918,105060,155,120,190,FR,France
198452,3,84830,60602,109058,154,110,198,FR,France
198451,3,101726,80242,123210,185,146,224,FR,France
198450,3,123680,101401,145959,225,184,266,FR,France
198449,3,101073,81684,120462,184,149,219,FR,France
198448,3,78620,60634,96606,143,110,176,FR,France
198447,3,72029,54274,89784,131,99,163,FR,France
198446,3,87330,67686,106974,159,123,195,FR,France
198445,3,135223,101414,169032,246,184,308,FR,France
198444,3,68422,20056,116788,125,37,213,FR,France
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment