diff --git a/module3/exo3/Paradoxe_de_Simpson.ipynb b/module3/exo3/Paradoxe_de_Simpson.ipynb new file mode 100644 index 0000000000000000000000000000000000000000..36d4a179876ea3ead0cb0403ce7d5d549106132b --- /dev/null +++ b/module3/exo3/Paradoxe_de_Simpson.ipynb @@ -0,0 +1,1562 @@ +{ + "cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Concentration de CO2 dans l'atmosphère depuis 1958" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "Les données pour cette étude ont été téléchargées au format CSV le 8 Janvier 2020 via [le lien du MOOC](https://gitlab.inria.fr/learninglab/mooc-rr/mooc-rr-ressources/blob/master/module3/Practical_session/Subject6_smoking.csv).\n", + "Sur chaque ligne il est indiqué si la personne fume ou non, si elle est vivante ou décédée au moment de la seconde étude, et son âge lors du premier sondage." + ] + }, + { + "cell_type": "code", + "execution_count": 87, + "metadata": {}, + "outputs": [], + "source": [ + "import pandas as pd\n", + "%matplotlib inline\n", + "import matplotlib.pyplot as plt\n", + "import isoweek\n", + "import pandas as pd\n", + "import numpy as np\n" + ] + }, + { + "cell_type": "code", + "execution_count": 88, + "metadata": {}, + "outputs": [ + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + "
SmokerStatusAge
0YesAlive21.0
1YesAlive19.3
2NoDead57.5
3NoAlive47.1
4YesAlive81.4
5NoAlive36.8
6NoAlive23.8
7YesDead57.5
8YesAlive24.8
9YesAlive49.5
10YesAlive30.0
11NoDead66.0
12YesAlive49.2
13NoAlive58.4
14NoDead60.6
15NoAlive25.1
16NoAlive43.5
17NoAlive27.1
18NoAlive58.3
19YesAlive65.7
20NoDead73.2
21YesAlive38.3
22NoAlive33.4
23YesDead62.3
24NoAlive18.0
25NoAlive56.2
26YesAlive59.2
27NoAlive25.8
28NoDead36.9
29NoAlive20.2
............
1284YesDead36.0
1285YesAlive48.3
1286NoAlive63.1
1287NoAlive60.8
1288YesDead39.3
1289NoAlive36.7
1290NoAlive63.8
1291NoDead71.3
1292NoAlive57.7
1293NoAlive63.2
1294NoAlive46.6
1295YesDead82.4
1296YesAlive38.3
1297YesAlive32.7
1298NoAlive39.7
1299YesDead60.0
1300NoDead71.0
1301NoAlive20.5
1302NoAlive44.4
1303YesAlive31.2
1304YesAlive47.8
1305YesAlive60.9
1306NoDead61.4
1307YesAlive43.0
1308NoAlive42.1
1309YesAlive35.9
1310NoAlive22.3
1311YesDead62.1
1312NoDead88.6
1313NoAlive39.1
\n", + "

1314 rows × 3 columns

\n", + "
" + ], + "text/plain": [ + " Smoker Status Age\n", + "0 Yes Alive 21.0\n", + "1 Yes Alive 19.3\n", + "2 No Dead 57.5\n", + "3 No Alive 47.1\n", + "4 Yes Alive 81.4\n", + "5 No Alive 36.8\n", + "6 No Alive 23.8\n", + "7 Yes Dead 57.5\n", + "8 Yes Alive 24.8\n", + "9 Yes Alive 49.5\n", + "10 Yes Alive 30.0\n", + "11 No Dead 66.0\n", + "12 Yes Alive 49.2\n", + "13 No Alive 58.4\n", + "14 No Dead 60.6\n", + "15 No Alive 25.1\n", + "16 No Alive 43.5\n", + "17 No Alive 27.1\n", + "18 No Alive 58.3\n", + "19 Yes Alive 65.7\n", + "20 No Dead 73.2\n", + "21 Yes Alive 38.3\n", + "22 No Alive 33.4\n", + "23 Yes Dead 62.3\n", + "24 No Alive 18.0\n", + "25 No Alive 56.2\n", + "26 Yes Alive 59.2\n", + "27 No Alive 25.8\n", + "28 No Dead 36.9\n", + "29 No Alive 20.2\n", + "... ... ... ...\n", + "1284 Yes Dead 36.0\n", + "1285 Yes Alive 48.3\n", + "1286 No Alive 63.1\n", + "1287 No Alive 60.8\n", + "1288 Yes Dead 39.3\n", + "1289 No Alive 36.7\n", + "1290 No Alive 63.8\n", + "1291 No Dead 71.3\n", + "1292 No Alive 57.7\n", + "1293 No Alive 63.2\n", + "1294 No Alive 46.6\n", + "1295 Yes Dead 82.4\n", + "1296 Yes Alive 38.3\n", + "1297 Yes Alive 32.7\n", + "1298 No Alive 39.7\n", + "1299 Yes Dead 60.0\n", + "1300 No Dead 71.0\n", + "1301 No Alive 20.5\n", + "1302 No Alive 44.4\n", + "1303 Yes Alive 31.2\n", + "1304 Yes Alive 47.8\n", + "1305 Yes Alive 60.9\n", + "1306 No Dead 61.4\n", + "1307 Yes Alive 43.0\n", + "1308 No Alive 42.1\n", + "1309 Yes Alive 35.9\n", + "1310 No Alive 22.3\n", + "1311 Yes Dead 62.1\n", + "1312 No Dead 88.6\n", + "1313 No Alive 39.1\n", + "\n", + "[1314 rows x 3 columns]" + ] + }, + "execution_count": 88, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "ini_data = pd.read_csv('Subject6_smoking.csv')\n", + "ini_data" + ] + }, + { + "cell_type": "code", + "execution_count": 89, + "metadata": {}, + "outputs": [ + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + "
SmokerStatusAge
\n", + "
" + ], + "text/plain": [ + "Empty DataFrame\n", + "Columns: [Smoker, Status, Age]\n", + "Index: []" + ] + }, + "execution_count": 89, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "ini_data[ini_data.isnull().any(axis=1)]" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "Y a-t-il des points manquants dans ce jeux de données ? Non il n'y en a pas" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "Dans un tableau sont représentés le nombre total de femmes vivantes et décédées sur la période en fonction de leur habitude de tabagisme. Pour chaque groupe (fumeuses / non fumeuses) le taux de mortalité (le rapport entre le nombre de femmes décédées dans un groupe et le nombre total de femmes dans ce groupe) est aussi calculé." + ] + }, + { + "cell_type": "code", + "execution_count": 90, + "metadata": {}, + "outputs": [ + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + "
StatusAliveDeadtaux_mortalité
Smoker
No50223031.420765
Yes44313923.883162
\n", + "
" + ], + "text/plain": [ + "Status Alive Dead taux_mortalité\n", + "Smoker \n", + "No 502 230 31.420765\n", + "Yes 443 139 23.883162" + ] + }, + "execution_count": 90, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "raw_data=ini_data.copy()\n", + "raw_data=pd.crosstab(index=raw_data['Smoker'], columns=raw_data['Status'])\n", + "\n", + "taux1=raw_data.Dead[0]/(raw_data.Dead[0]+raw_data.Alive[0])\n", + "taux2=raw_data.Dead[1]/(raw_data.Dead[1]+raw_data.Alive[1])\n", + "raw_data['taux_mortalité']= [taux1*100, taux2*100]\n", + "raw_data" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "Ces résultats sont surpenants car ils montrent qu'il y a plus de femmes qui ne fumaient pas lors de la première étude qui sont mortes au moement de la deuxième (230 non fumeuses mortes) que de femmes qui fumaient (139 fumeuses mortes). Le taux de mortalité est plus élevé chez les non fumeuses (31 %) que chez les fumeuses (24 %).\n", + "Après 20 ans, le taux de mortalité chez les fumeuses est de 24%, alors que celui des non-fumeuses est de 31%. Alors, est-ce que non-fumer tue ?\n", + "\n", + "Afin de préciser cette étude les effectifs et les taux de mortalité vont être recalculés en ajoutant une nouvelle catégorie liée à la classe d'âge. Ici sont considéres les classes d'ages suivantes : 18-34 ans, 34-54 ans, 55-64 ans et plus de 65 ans." + ] + }, + { + "cell_type": "code", + "execution_count": 91, + "metadata": {}, + "outputs": [ + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + "
SmokerStatusAge
0YesAlive21.0
1YesAlive19.3
6NoAlive23.8
8YesAlive24.8
10YesAlive30.0
15NoAlive25.1
17NoAlive27.1
22NoAlive33.4
24NoAlive18.0
27NoAlive25.8
29NoAlive20.2
33NoAlive19.4
37YesAlive29.5
38YesDead33.0
44NoAlive25.3
47NoAlive18.5
49YesAlive22.1
54NoAlive28.4
58NoAlive22.9
65YesAlive33.0
67YesAlive27.9
71YesAlive26.2
76NoAlive27.6
77YesAlive31.4
79NoAlive18.9
81YesAlive25.4
84NoAlive27.3
86NoAlive32.8
91NoAlive18.3
92YesAlive20.2
............
1205NoAlive23.2
1207YesAlive31.4
1208YesAlive30.0
1213NoAlive21.4
1216YesAlive27.9
1217YesAlive29.5
1219YesAlive27.0
1223YesAlive28.3
1226YesAlive31.0
1232NoAlive28.3
1240YesAlive29.7
1247NoAlive26.0
1250NoAlive19.8
1251YesAlive27.8
1253YesAlive27.8
1255NoDead28.5
1256NoAlive26.7
1260YesAlive20.4
1263YesAlive20.9
1265NoAlive26.7
1267NoAlive33.7
1271YesAlive24.9
1272NoAlive33.0
1274NoAlive25.7
1275NoAlive19.5
1277NoAlive23.4
1297YesAlive32.7
1301NoAlive20.5
1303YesAlive31.2
1310NoAlive22.3
\n", + "

400 rows × 3 columns

\n", + "
" + ], + "text/plain": [ + " Smoker Status Age\n", + "0 Yes Alive 21.0\n", + "1 Yes Alive 19.3\n", + "6 No Alive 23.8\n", + "8 Yes Alive 24.8\n", + "10 Yes Alive 30.0\n", + "15 No Alive 25.1\n", + "17 No Alive 27.1\n", + "22 No Alive 33.4\n", + "24 No Alive 18.0\n", + "27 No Alive 25.8\n", + "29 No Alive 20.2\n", + "33 No Alive 19.4\n", + "37 Yes Alive 29.5\n", + "38 Yes Dead 33.0\n", + "44 No Alive 25.3\n", + "47 No Alive 18.5\n", + "49 Yes Alive 22.1\n", + "54 No Alive 28.4\n", + "58 No Alive 22.9\n", + "65 Yes Alive 33.0\n", + "67 Yes Alive 27.9\n", + "71 Yes Alive 26.2\n", + "76 No Alive 27.6\n", + "77 Yes Alive 31.4\n", + "79 No Alive 18.9\n", + "81 Yes Alive 25.4\n", + "84 No Alive 27.3\n", + "86 No Alive 32.8\n", + "91 No Alive 18.3\n", + "92 Yes Alive 20.2\n", + "... ... ... ...\n", + "1205 No Alive 23.2\n", + "1207 Yes Alive 31.4\n", + "1208 Yes Alive 30.0\n", + "1213 No Alive 21.4\n", + "1216 Yes Alive 27.9\n", + "1217 Yes Alive 29.5\n", + "1219 Yes Alive 27.0\n", + "1223 Yes Alive 28.3\n", + "1226 Yes Alive 31.0\n", + "1232 No Alive 28.3\n", + "1240 Yes Alive 29.7\n", + "1247 No Alive 26.0\n", + "1250 No Alive 19.8\n", + "1251 Yes Alive 27.8\n", + "1253 Yes Alive 27.8\n", + "1255 No Dead 28.5\n", + "1256 No Alive 26.7\n", + "1260 Yes Alive 20.4\n", + "1263 Yes Alive 20.9\n", + "1265 No Alive 26.7\n", + "1267 No Alive 33.7\n", + "1271 Yes Alive 24.9\n", + "1272 No Alive 33.0\n", + "1274 No Alive 25.7\n", + "1275 No Alive 19.5\n", + "1277 No Alive 23.4\n", + "1297 Yes Alive 32.7\n", + "1301 No Alive 20.5\n", + "1303 Yes Alive 31.2\n", + "1310 No Alive 22.3\n", + "\n", + "[400 rows x 3 columns]" + ] + }, + "execution_count": 91, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "# Pour ne garder que les femmes ayant entre 18 et 34 ans : enlèves colonnes > 34 ans\n", + "dfini_18_34=ini_data.drop(ini_data[ini_data.Age > 34] .index)\n", + "# Pour ne garder que les femmes ayant entre 35 et 54 ans : enlèves colonnes < ou = à 34 ans et colonnes > 54 ans\n", + "dfini_35_54=ini_data.drop(ini_data[(ini_data.Age <= 34) | (ini_data.Age >54)] .index)\n", + "# Pour ne garder que les femmes ayant entre 55 et 64 ans : enlèves colonnes < ou = à 54 ans et colonnes > 64 ans\n", + "dfini_55_64=ini_data.drop(ini_data[(ini_data.Age <= 54) | (ini_data.Age >64)] .index)\n", + "# Pour ne garder que les femmes ayant plus de 65 ans : enlèves colonnes < ou = à 64 ans\n", + "dfini_65=ini_data.drop(ini_data[ini_data.Age <= 64] .index)\n", + "\n", + "dfini_18_34" + ] + }, + { + "cell_type": "code", + "execution_count": 92, + "metadata": {}, + "outputs": [ + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + "
StatusAliveDeadtaux_mortalité
Smoker
No21362.739726
Yes17652.762431
\n", + "
" + ], + "text/plain": [ + "Status Alive Dead taux_mortalité\n", + "Smoker \n", + "No 213 6 2.739726\n", + "Yes 176 5 2.762431" + ] + }, + "execution_count": 92, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "df_18_34=dfini_18_34.copy()\n", + "df_18_34=pd.crosstab(index=df_18_34['Smoker'], columns=df_18_34['Status'])\n", + "taux11=df_18_34.Dead[0]/(df_18_34.Dead[0]+df_18_34.Alive[0])\n", + "taux22=df_18_34.Dead[1]/(df_18_34.Dead[1]+df_18_34.Alive[1])\n", + "df_18_34['taux_mortalité']= [taux11*100, taux22*100]\n", + "\n", + "df_18_34" + ] + }, + { + "cell_type": "code", + "execution_count": 93, + "metadata": {}, + "outputs": [ + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + "
StatusAliveDeadtaux_mortalité
Smoker
No180199.547739
Yes1964117.299578
\n", + "
" + ], + "text/plain": [ + "Status Alive Dead taux_mortalité\n", + "Smoker \n", + "No 180 19 9.547739\n", + "Yes 196 41 17.299578" + ] + }, + "execution_count": 93, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "df_35_54=dfini_35_54.copy()\n", + "df_35_54=pd.crosstab(index=df_35_54['Smoker'], columns=df_35_54['Status'])\n", + "taux111=df_35_54.Dead[0]/(df_35_54.Dead[0]+df_35_54.Alive[0])\n", + "taux222=df_35_54.Dead[1]/(df_35_54.Dead[1]+df_35_54.Alive[1])\n", + "df_35_54['taux_mortalité']= [taux111*100, taux222*100]\n", + "\n", + "df_35_54" + ] + }, + { + "cell_type": "code", + "execution_count": 94, + "metadata": {}, + "outputs": [ + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + "
StatusAliveDeadtaux_mortalité
Smoker
No814033.057851
Yes645144.347826
\n", + "
" + ], + "text/plain": [ + "Status Alive Dead taux_mortalité\n", + "Smoker \n", + "No 81 40 33.057851\n", + "Yes 64 51 44.347826" + ] + }, + "execution_count": 94, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "df_55_64=dfini_55_64.copy()\n", + "df_55_64=pd.crosstab(index=df_55_64['Smoker'], columns=df_55_64['Status'])\n", + "taux1111=df_55_64.Dead[0]/(df_55_64.Dead[0]+df_55_64.Alive[0])\n", + "taux2222=df_55_64.Dead[1]/(df_55_64.Dead[1]+df_55_64.Alive[1])\n", + "df_55_64['taux_mortalité']= [taux1111*100, taux2222*100]\n", + "\n", + "df_55_64" + ] + }, + { + "cell_type": "code", + "execution_count": 95, + "metadata": {}, + "outputs": [ + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + "
StatusAliveDeadtaux_mortalité
Smoker
No2816585.492228
Yes74285.714286
\n", + "
" + ], + "text/plain": [ + "Status Alive Dead taux_mortalité\n", + "Smoker \n", + "No 28 165 85.492228\n", + "Yes 7 42 85.714286" + ] + }, + "execution_count": 95, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "df_65=dfini_65.copy()\n", + "df_65=pd.crosstab(index=df_65['Smoker'], columns=df_65['Status'])\n", + "taux11111=df_65.Dead[0]/(df_65.Dead[0]+df_65.Alive[0])\n", + "taux22222=df_65.Dead[1]/(df_65.Dead[1]+df_65.Alive[1])\n", + "df_65['taux_mortalité']= [taux11111*100, taux22222*100]\n", + "\n", + "df_65" + ] + }, + { + "cell_type": "code", + "execution_count": 96, + "metadata": {}, + "outputs": [ + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYgAAAEICAYAAABF82P+AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAHRdJREFUeJzt3Xm4XFWZ7/HvLwkQSAhDEzDMiNGAXKHhgCAOYGga0QYckNCgYVDEqyDQiGijNF68BBUHBLxGEFAQiYjNJDIEUEEIBAIyhEnGSAgHIkMwBJK894+1TrNz2FVnJydVuzjn93meeqr2tNa7T+2z373WHkoRgZmZWW9D6g7AzMw6kxOEmZmVcoIwM7NSThBmZlbKCcLMzEo5QZiZWSknCGtI0jhJC+uOY3mQtJukhwvDf5W0Q4Xlhkq6UtKnWxvhG+r9mKTrJK3YznotkbSVpHl1x1E3J4gOIGle4bVY0vzC8H51x9eJJE2SdOayLh8Rm0bEzRXK+g7w3xHx82Wta1lExMXAmcBpy1qGpFNzInxJ0n2S9u01fVtJd0r6h6RbJW3R37hbRdLbe/2fhKSXC8Pb9rH8EZIub1e8A8WwugMwiIiRPZ8lPQZ8JiKurS+iziapbdttRBzVrrogrVtELMx1/xL4ZT+KexH4EPAw8B7gCkkPRsTtklYGLgFOBM4CDgd+K2mznvo7SUQ8CIwEkDQcmA+8IyJm1RrYQBcRfnXQC3gM2KXXuB2BacALwFPA94Fhedo4YGGv+W8B9s+fzwbOL0z7IXBFg7qH5enPkXYqhxXLBtYEfg48DTwJHA8MaVDWJOB84EJgHnAnsEle5tm8njsX5t8Q+B0wF3gQmNirrF/msl4C9gdeBV7LZd+a5/sccH+e52HgoEIZuwEPF4afBt4L7NWgrKVZ12bfz3AggC/mde4GvgUoTz8UuA44Hfg7cFxhXR7I464ENsjjh5JaFd25vrtIO8oq29bVwBfy5z2ARwrThuR13anBsh/Ndb0IPA58rTBtHLAQOBCYlWP7cq+/z4y87NPASf38H+n5m67fa/xawK/z9vUIcGQevwOwIMc4D3gsj/8kcHfeXh4Dji6UtRUwr+79Qd2v2gPwq9cXUp4gtgO2zTuHTfPO79A8ra8EsSrwKDABGA88A7ylQd1H5H+YdYHRwI0smSCuBH4ErAKMyf/0ExuUNQn4B7AzKfFcmOM4Og8fBswszD+NtGNdCegiJYodC2UtAHbPO7KV87gze9W5BykJCdiFdJT5zjytNEEUyu9d1tKsa7Pvp2dndhWweo7vkcL3c2jecX02L79y/q4eAN4BrACcBPwxz78ncDMwKv8t3gmsXWG7Gknace6Uh78K/LbXPNeSE0jJ8uNzXUOArfP3s1thGwxSkhue/xavAm/N02cAexe2x3f383+kUYK4FDgvf2fjSIn944Vt+/Je8++a51P+Dl8o/H2cICJ8DuLNICJujYjbImJRRPyV1Df9gYrLvgR8mnTUeS5px/V0g9k/CZwSEU9FRDfw7Z4JkjYC3g8cFRH/iIjZwKmknVkjUyPi+khdFheRdmqn5OFfAeMkrSxpLLAl6ah0QURMz7F+qlDWHyLidxGxOCLmN1jXSyPi0UiuBf5AaiUslaVd14rfz0kR8XxEPEr6LornAx6JiJ/m5eeTWg8nRsQDEfEaqfWyg6S3kFo6o8g75Yi4NyKe6WN9lGO6MSJuyKNHknaIRS+QduBl6zg117U4Iu4AppSs4/ER8UpE3EZqyb0rj38NeLukf4qIlyJiWrN4l4WkVYF/A47J39n9pL/zpxotExFXR8T9eXu5FfhtyToNak4QbwKSNs9X0syR9CLwDVJzuqobgdnAK6R/gkbWJR119Xi88Hkj0pFbt6TnJT1P6o5ap0l5cwqf5wPdkQ/P8jDAiFxvd68d/+PAeoXhYlylJO2RT7bOzfF9kKX7O/VYqnWt+P30/ruu22BaT/0nSrpf0v3AX0hdTWNILZuzgJ8AcySdIWkkzZ2ay9y/MG4eKdEUjSJ1t5St446S/iCpW9ILwAEsuY6LIuLZwvA/yOcMgImkZPGgpGmS/rVBHdcVTjp/vI916m1dUmv3qcK43ttQ7/p2knSjpGfzOu3Lsm0vA5YTxJvDT4E7gE0jYhTwTVKzGOBlYKiklQrzv6XX8keRjuJeJDW1G5kNbFAY3rDw+UnSTmWNiFg9v0ZFxNZLvTZv9BQwOp84Ldb9t8Jw78cOLzEsaQSp//n/kLpcVif17Yu+9S57ade12ffTo/fftbgjK6v/6IgYV3itHREz8tHu9yLin0k73S2BLzVaMUmTSK2oD0VE8bLNe/OyPfMNAbbI48tMIXUTbhARqwHnlKxjqYiYGRH7AGuTktXFZZfvRsQHI2Jkfv2mStkFTwHDJBUTb3Eb6r29iNSqPRtYL6/TBVXXabBwgnhzWBV4ISLmSXonqb+6x1Okk4L75Wv2/zeFo6Z86eJxpKPH/YFvSNq8QT1TgCMljZG0FnBMz4TcNXIL8G1Jq0oaImmspKXuwinxMOko+URJK0namnTUeX6TZeYAm+R/dEh99yuQzrEslrQHsFPF+pcoaxnWtdn30+MrklaTtDHphPWFTeL5f8DXJW0GIGmNniNqSdtL6spXcr1M6utfVFaIpBNI5yx2jYjne02+BlhZ0qH54OLIXN6NJeWI1Bp4LiJekfQeYO8m8fde/tO5e2kRqRsrgMVVl68id6VeAUyStIqktwNfIJ2TgPQdb1i4Am4IqfX6LPCqpA+QTsRbgRPEm8ORwGeUbtw5ncLOJf/TfYbXrw7aALgdIB+lnQecEBH3RcR9pKPbX0haoaSe04A/kY4ip5ESRtG+pBOt95NOUl5I8y6mSnK30yeBzUknjy8kXQXzpyaL/Yp0MnKupD/n7o2jgctIV2HtRboqqoolysrjlmZdG34/BVeQrgKaTmrpnFcyDwARcQGpS+ui3GV1J/AvefLqpKP350knux8nHZUvIe/0vwG8FXi00HVzVK5jPil5HJrLmgDsFSWXuObv51Dgu5JeIh04/LpR/CU+AjyQlz0J+GRZPcvBQaTv8UnSFVunRsRFedrlpAOpbkl/zf83hwJnkJLWYaQWhRX0XGpnZi1QuGZ/g/A1+/Ym4xaEmZmValmCkPQzSc9Iuqcwbk1J10h6KL+vUZj2VUkPS3qg0VUOZmbWPq1sQZxDujmp6FjStfFjgal5mHzSdALpRpzdgDMkDW1hbGZtke8LkLuX7M2oZQkiIv5IOrlXtCfpBijy+16F8b/KN0k9SrqqZbtWxWZmZn1r98P61sl3pRIRsyWtncevR7qssMcsGtzgIukQ4BCAESNGbDNu3LgWhmtmNvDcfvvtz0bE6L7m65SnuZbdnFJ6eVVETAYmA3R1dcX06dNbGZeZ2YAj6fG+52r/VUxzJI0ByO89z5CZxZJ3mq7PkneamplZm7U7QVxKukOW/H5JYfyEfBftJsBY4NY2x2ZmZgUt62KSdAHpUQdrSZpFutN3EjBF0sHAE+Tb9SPiXklTgPtIjz7+Qr7T0czMatKyBBER+zaYNL7B/N8i/ZCKmZl1AN9JbWZmpZwgzMyslBOEmZmVcoIwM7NSThBmZlbKCcLMzEo5QZiZWSknCDMzK+UEYWZmpZwgzMyslBOEmZmVcoIwM7NSThBmZlbKCcLMzEo5QZiZWSknCDMzK+UEYWZmpZwgzMyslBOEmZmVcoIwM7NSThBmZlbKCcLMzEo5QZiZWSknCDMzK+UEYWZmpZwgzMyslBOEmZmVcoIwM7NSThBmZlbKCcLMzEo5QZiZWSknCDMzK+UEYWZmpZwgzMyslBOEmZmVqiVBSDpS0r2S7pF0gaThktaUdI2kh/L7GnXEZmZmSdsThKT1gMOBrojYAhgKTACOBaZGxFhgah42M7Oa1NXFNAxYWdIwYBXgKWBP4Nw8/Vxgr5piMzMzakgQEfE34LvAE8Bs4IWIuBpYJyJm53lmA2uXLS/pEEnTJU3v7u5uV9hmZoNOHV1Ma5BaC5sA6wIjJO1fdfmImBwRXRHRNXr06FaFaWY26NXRxbQL8GhEdEfEa8DFwHuAOZLGAOT3Z2qIzczMsjoSxBPA9pJWkSRgPDATuBSYmOeZCFxSQ2xmZpYNa3eFETFN0kXAHcBCYAYwGRgJTJF0MCmJ7N3u2MzM7HVtTxAAEXE8cHyv0QtIrQkzM+sAvpPazMxKOUGYmVmppUoQkjaStFmrgjEzs85R+RyEpK8AXcBiSfMj4oCWRWVmZrVr2IKQ9HlJxelbR8TeEbEPsHXrQzMzszo162KaD/xe0ofy8FRJ10m6nvQwPTMzG8AadjFFxDmSpgBfkXQI8HXgAmDFiHiuXQGamVk9+joHsQHpyaoLgBOBV3jj/QtmZjYANUwQks4CRgArA/dFxIGSuoCzJd0YESe1K0gzM2u/ZucguiJiQkTsCewGEBHTI+LDwINtic7MzGrTrIvpGknXASsCFxYnRMRvWhqVmZnVrtlJ6qMlrQksiogX2hiTmZl1gKYnqSNibrsCMTOzzuJnMZmZWSknCDMzK1XpWUyStgM2Ls4fEb9sUUxmZtYB+kwQks4BNgfuBBbl0QE4QZiZDWBVWhDbA5tHxOJWB2NmZp2jyjmIe4G1Wh2ImZl1liotiNWAmZJuIT2TCYCI+FjLojIzs9pVSRB+5pKZ2SDUZ4KIiKmS1gfGRsT1koYDQ1sfmpmZ1anPcxCSDgIuBc7MozYELmllUGZmVr8qJ6kPJ13J9CJARDwIrNPKoMzMrH5VEsQrEfFqz4Akdy+ZmQ0CVRLETZKOAYZL2pn06O/LWxuWmZnVrUqCOAZ4Cbgf+BIwFfjPVgZlZmb1q3IV0yLgx/llZmaDhJ/mamZmpZwgzMysVOUEIWmlVgZiZmadpTRBSBpW+LydpLuBh/LwlpJ+1Kb4zMysJo1aEJ+TtEP+fCrwEeA5gIi4C9i5DbGZmVmNGiWInwD79cwTEY/3mr4IMzMb0EoTREQsjIgv5sEn80+OhqShko4AHuxPpZJWl3SRpPslzZS0g6Q1JV0j6aH8vkZ/6jAzs/6pcpL688BRpIf0zSE9l+nQftb7Q+D3ETEO2BKYCRwLTI2IsaSb8Y7tZx1mZtYPVX4P4q0RMaE4QtL25HMSS0vSKOD9wAEA+TlPr0raE9gpz3YucAPwlWWpw8zM+q9KC+KMknGn96POtwLdwNmSZkg6U9IIYJ2ImA2Q39cuW1jSIZKmS5re3d3djzDMzKyZhi2IfN5hB2C0pMMLk0YBK/Szzq2BwyJimqQfshTdSRExGZgM0NXVFf2Iw8zMmmjWghgBrEXaoY8uvF4F9u5HnbOAWRExLQ9fREoYcySNAcjvz/SjDjMz66eGLYiIuB64XtLZEfHI8qowIp6W9KSkd0TEA8B44L78mghMyu/+1Tozsxo162I6JSL+AzhF0hu6ciLiY/2o9zDgfEkrAo8AB5JaM1MkHQw8Qf9aKWZm1k/NrmK6ML+ftrwrjYg7ga6SSeOXd11mZrZsmnUx3Zrfp7YvHDMz6xTNuphmAA2vEoqIrVsSkZmZdYRmXUyfaFsUZmbWcZp1Mf21nYGYmVln6fNOaknbSrpF0guSXpG0QNKL7QjOzMzqU+VZTGcA+wO/ArYjPUNpgxbGZGZmHaDKs5iG5BvahkXEaxHxU2CXFsdlZmY1q9KCeDnf0HaXpP8LzAZGtjYsMzOrW5UWxAF5vi+SfkluLPDxFsZkZmYdoEqC2D0iXomI5yPi6xFxOPCvrQ7MzMzqVSVBHFQy7uDlHYiZmXWWZndS7wNMAN4q6eLCpFHA860OzMzM6tXsJPWtpJ8VXZ8lf0HuJWBGK4MyM7P6NbuT+lFJTwBz/cA+M7PBp+k5iIhYBLwqaVSb4jEzsw5R5T6IeaR7IK4GXu4ZGRFHtSwqMzOrXZUEcW1+mZnZINJngoiIsyQNA96WRz0cEQtbG5aZmdWtzwQh6X3AL4C/AQLeIulTEXFTq4MzM7P6VOli+j7pbur7ACRtRkoYZb8pbWZmA0SVO6lX7EkOABExE1ixdSGZmVknqNKCuEPST0itBoD98I1yZmYDXpUEcShwOHAM6RzEH4EftTIoMzOrX5WrmF4Bvp1fZmY2SFT5TerdJN0m6RlJcyX9XdLcdgRnZmb1qdLFdBrwSeBuYHFrwzEzs05RJUHMAu6MCCcHM7NBpEqCOAa4TNINwIKekRFxaquCMjOz+lVJECcArwGr4y4mM7NBo0qCWDsitml5JGZm1lGq3Ek9VdIHWx6JmZl1lCoJ4rPAtZLm+TJXM7PBo0oX01otj8LMzDpOlTupF7UjEDMz6yxVupjMzGwQqi1BSBoqaYaky/PwmpKukfRQfl+jrtjMzKxigpC0vaRP58//JGnD5VD3l4CZheFjgakRMRaYmofNzKwmVR7WdxxwPHBcHjUc+GV/KpW0PvBh4MzC6D2Bc/Pnc4G9+lOHmZn1T5UWxCeA3YGXASLib8Coftb7A9IjPIp3Zq8TEbNzHbOBtcsWlHSIpOmSpnd3d/czDDMza6RKglgQEQEEgKRV+lOhpI8Az0TE7cuyfERMjoiuiOgaPXp0f0IxM7MmqtwHcbGk04HVJB0IHAz8rB917gjsIWl3UnfVKEnnAXMkjYmI2ZLGAM/0ow4zM+unPlsQEXEycDlwKbAl8K2I+MGyVhgRX42I9SNiY2ACcF1E7J/Ln5hnmwhcsqx1mJlZ/1VpQRARVwJXtjiWScAUSQcDTwB7t7g+MzNromGCkPR38nmHMhGxZn8rj4gbgBvy5+eA8f0t08zMlo9mLYi1AJEuce0GfpGH9wP6daLazMw6X8ME0fMMJkm7RsS7C5N+JOkW4ORWB2dmZvWpcplrSNpHkgAk7dPimMzMrANUSRD/DnwaeE7Ss8CnSN1MZmY2gFV53PcjpMdimJnZIOLHfZuZWSknCDMzK+UEYWZmpao87vtsSasWhteXdHVrwzIzs7pVaUFMB26VtGt+WN/1wI9bG5aZmdWtylVMp0u6i5QYngW27vndBjMzG7iqdDHtS3q890HAecBlkrZodWBmZlavKk9z3Q94f0Q8DSBpB+B80qO/39R0guoOwTpUHN/wOZVmg0aVLqaP9Bq+WdK7G81vZmYDQ58JQtLkBpMOWc6xmJlZB6nSxTS18Hk48FHgydaEY2ZmnaJKF9OFxWFJvwCuaVlEZmbWEZblTupNgI2WdyBmZtZZqpyDKP706BBgLnBsK4MyM7P6VTkHsVbh8+KI8PV/ZmaDQJVzEIskrQZsCgzPPyxHRPy5xbGZmVmNqnQxHQwcBawH3A1sC9wC7NTSyMzMrFZVTlIfAXQBj0XE+4BtAD+LycxsgKuSIF6JiPkAklaMiHuBca0Ny8zM6tawi0nSsIhYCMyWtDpwGXCVpLnAnHYFaGZm9Wh2DuJW0qO998jDX5c0HlgNuKLlkZmZWa2aJYg3POo0IqaWzWhmZgNPswQxWtJRjSZGxPdaEI+ZmXWIZgliKDCSkpaEmZkNfM0SxOyI+GbbIjGzN5KPz6yBNjzUotllrt4yzcwGsWYJYnzbojAzs47TMEFExNx2BmJmZp1lWX4PwszMBoG2JwhJG0i6XtJMSfdK+lIev6akayQ9lN/XaHdsZmb2ujpaEAuB/4iIzYDtgS9I2pz0I0RTI2Is6Xew/aNEZmY1anuCiIjZEXFH/vwSMJP0KPE9gXPzbOcCe7U7NjMze12t5yAkbQz8MzANWCciZkNKIsDaDZY5RNJ0SdO7u7vbFaqZ2aBTW4KQNBL4DXBERLxYdbmImBwRXRHRNXr06NYFaGY2yNWSICStQEoO50fExXn0HElj8vQxwDN1xGZmZkkdVzEJOAuY2euBf5cCE/PnicAl7Y7NzMxe1+dvUrfAjsCngLsl3ZnHfQ2YBEzJv4H9BLB3DbGZmVnW9gQRETfS+DlPfryHmVmH8J3UZmZWygnCzMxKOUGYmVkpJwgzMyvlBGFmZqWcIMzMrJQThJmZlXKCMDOzUk4QZmZWygnCzMxKOUGYmVkpJwgzMyvlBGFmZqWcIMzMrJQThJmZlXKCMDOzUk4QZmZWygnCzMxKOUGYmVkpJwgzMyvlBGFmZqWcIMzMrJQThJmZlXKCMDOzUk4QZmZWygnCzMxKOUGYmVkpJwgzMyvlBGFmZqWcIMzMrJQThJmZlXKCMDOzUk4QZmZWygnCzMxKOUGYmVmpjksQknaT9ICkhyUdW3c8ZmaDVUclCElDgdOBDwGbA/tK2rzeqMzMBqeOShDAdsDDEfFIRLwK/ArYs+aYzMwGpWF1B9DLesCTheFZwLuLM0g6BDgkD86T9ECbYhvo1gKerTuITqH/Ut0h2Bt5Gy1Sv7bRjarM1GkJomyNY4mBiMnA5PaEM3hImh4RXXXHYdaIt9H267QuplnABoXh9YGnaorFzGxQ67QEcRswVtImklYEJgCX1hyTmdmg1FFdTBGxUNIXgauAocDPIuLemsMaLNxtZ53O22ibKSL6nsvMzAadTutiMjOzDuEEYWZmpZwgaiQpJJ1SGD5a0n8tp7JHS5omaYak9y2PMs2akbRI0p2F18Z1x2T901EnqQehBcDHJJ0UEcv7BqDxwP0RMXE5l2vWyPyI2KruIGz5cQuiXgtJV2Yc2XuCpI0kTZX0l/y+YR5/jqRTJf1Z0iOSPlGy7FbAt4Hd85HcypLmFaZ/QtI5hfJ+LOn6XN4HJP1M0syeefJ8u0q6WdIdkn4taWQe/5iktfLnLkk35M8fKBxJzpC0ah7/ZUm35fU6IY8bIekKSXdJukfSPsvlr2u1k3SApNMKw5dL2il/nifpZEm3S7pW0naSbsjb4R55nqGSvlPYZj6Xx+8k6fJCuadJOiB/niTpvjz/d/O40ZJ+k8u5TdKOeXzpdmqJE0T9Tgf2k7Rar/GnAT+PiHcB5wOnFqaNAd4LfASY1LvAiLgT+AZwYURsFRHz+4hhDeCDpER1GfB94J3A/5K0VU4AxwG7RMTWwHTgqD7KPBr4Qj6ifB8wX9KuwFjSM7e2AraR9H5gN+CpiNgyIrYAft9H2daZVi7sbH9bYf4RwA0RsQ3wEnAi8C/AR4Fv5nkOBl6IiG2BbYHPStqkUYGS1szLvzP/75yYJ/0Q+H4u5+PAmXn8G7bT6qs78LmLqWYR8aKknwOHs+TGuQPwsfz5F6QWQY//jojFwH2S1lkOYVwWESHpbmBORNwNIOleYGPSHe2bAzcpPf9lReDmPsq8CfiepPOBiyNiVk4QuwIz8jwjSQnjT8B3JZ0MXB4Rf1oO62Ttt7RdTK/y+sHA3cCCiHgtb4cb5/G7Au8qtJRXI20zrzYo80XgFeBMSVcAPa2MXYDN9frzi0bl1sIbttOliH/Ac4LoDD8A7gDObjJP8YaVBYXPApD0LeDDAA3+SYvLD+81rae8xb3KXkzaRhYB10TEviXlLuT1luj/lBsRk/I/6O7ALZJ2ybGeFBE/6V2IpG3yvCdJujoivtl7HntTKm4fsOS291q8fiPW/2x7EbFYUs++ScBhEXFVsVBJ7y0rN99sux3pHNwE4Iuk1vEQYIeS1vQbttOIuH/ZVnXgcRdTB4iIucAUUnO6x59JGzjAfsCNfZTxn7k7qdER3BxJm0kaQmqCL41bgB0lvQ1A0iqS3p6nPQZskz9/vGcBSZtGxN0RcTKpS2oc6Q75gwrnL9aTtLakdYF/RMR5wHeBrZcyPutcjwFbSRoiaQNS9+LSuAr4vKQVACS9XdII4HFSi2Cl3D07Pk8fCawWEb8DjiB1ZQJcTUoW5Pm2yu9l26llbkF0jlMobMCkLqefSfoy0A0c2M/yjyU1t58E7iF171QSEd35BOAFklbKo48DHgROAM6S9DVgWmGxIyTtTGp93AdcGRELJG0G3Jyb+vOA/YG3Ad+RtBh4Dfj8Mq+ldZqbgEdJXUj3kFrKS+NMUnfTHUobTTewV0Q8KWkK8BfgIV7vtlwVuETScFLro+cCkMOB0yX9hbTf+yNwKCXb6bKs5EDlR22YmVkpdzGZmVkpJwgzMyvlBGFmZqWcIMzMrJQThJmZlXKCMDOzUk4QZmZW6v8DLGSay6+Roj8AAAAASUVORK5CYII=\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + }, + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 96, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYgAAAEICAYAAABF82P+AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzt3XmYFOW5/vHvDcOmgCsaDCqQo0EERASDoh4MBo3iftyiEbegOeIaFzwmbkfzI1Gjoh4Tj2sMcfm5JzEGwd2oBHTcAAUVFSWyKQIBFHnOH1UztmN10zDT3cPM/bmuvrprf+qdnn6q3rfqLUUEZmZmdbWodABmZtY4OUGYmVkmJwgzM8vkBGFmZpmcIMzMLJMThJmZZXKCsDUmqYekFZWOoyFI2kvSjJzhtyXtVMRyLSX9VdLRpY2wcSrld0BSW0khqUuBeR6WNK4U2zcniLWCpMU5r5WSluYMH1np+BojSaMl3bSmy0fEdyLi+SLWdTnwYET8fk23VV/13de1haS7JB2eM7wHsACYKumAykXWdFVVOgBbtYhoX/NZ0kzghIgYX7mIGjdJZfteR8SZ5doWJPsWEat1xC6pBUBErCxNVBXTCTgTWA4cXOFYmiSfQTQBkgZJelHSQkkfSbqq5kcyqwpA0guSjko/3yppbM60ayT9Jc92qtLp89PqmB+sIq7NJT0kaZ6kdySdlDNttKSxku5Oz4SqJXWTdGE6/0xJu+fMv4WkRyQtkPSWpOF11vXHdF2LgMNJfjiGp+uemM53oqRpkhZJmiHpuAKx/1PSLumRada6NpT0+3S+D9K4M/+f0vjulHRfuu1/SNo2Z/oFkt5Np70uaZ+caSdJelzS9ZI+AUbVWXe++F6QdImkF4F/AZsV2v+aKjZJ/yVprqQPc89OJa0raUy6rwslPZWbiCUdK2lWuuzZOeNbSvpF+vefl/7N1y9Q7udL+ljSLOCoAvP1AE4EZgAzgSGSOuRM31HSK+m+/lHS/ZJ+njP9QEmvSvpU0jOSeubbVrMWEX6tRS+Sf4Y96ozbERgAtAS+Q/JPc1I6rQewos78LwBHpZ87AO+S/KgOAeYA38qz7dOB14DNSI7enq277px5W6bzngu0BrYG3gf+PZ0+muSHa3eSM9m70zjOSodPAabmrO9F4CqgDdCfpGphUM66lgN7kxz0tEvH3VQnpv2AboCAPYClwLbptL2AGTnz/hPYJWf9ddf1V+BaYB2gM/AyMDxPWYwGPk+33wr4OfAm0DKdfli6jhbAj4FFwMbptJOAFcBP0jJtl2f9deN7AXgH+G66zaoi9v8L4Px0/gPTONqn028GxgHfSuPYNX3vAQRwPdCW5Hv4OdA9XW4U8Ez6nWkL3AbcmqecDgA+TNfZHrgvXXeXjHl7AN8n+W59K93f0em0tsBHadlVkXy3vwB+nk4fCMwGdkj3YQTwFlBV6f/vxvaqeAB+reYfLCNBZMwzCrgz/VwwQaTDuwLzgFnAQQXW+3fgmJzh/equO2favwPT64y7GLgh/Twa+FPOtEOA+YDS4U7pj0M7YCtgWe6PI0my+G3OusbV2dY3fjQzYnwUODH9XHSCALYElgCtcsYdC/w1z3ZGA0/mDFel+zogz/zTgD3TzycBb61iP/IliP9azf1fCLTImf4Z0JckYXwBfDdjHTUJYuOcca8CB6Sf3yVN5OlwN5IDA2Ws64/ARTnDfciTIDKWPRx4Pv08FHinzvRJfJUgbgXOrzP9PeB7xf4fNpeX2yCagPT0+EqgH8kPahXw3Gqs4lmSI6p2wAMF5tsM+CBn+L0C824JdJX0ac64lkBu28nHOZ+XAnMj/W9NhwHWTbc7NyKW5sz/HskZT43cuDJJ2o/k6P3fSI7W1yE5ul1dW5Icpc6VVDOuBcmZWz618UXECkkfkewXko4HTgO2SGdpD2yctexq+tpyRez/3Ph6O8W/0lg6k3yn3smznS8jYl7d5ZQUzubAI5JyewVtAWxEclCSazNgQs5w3u+XpM2Aa4CdSc6CW5B8h2vWM6vOIrllsSVwaG5VGMmZyLfzba+5chtE0/C/wEvAdyKiI3AJSTUCJEe6LSW1yZn/W3WWP5PkCPEzkmqkfGaT/MPX2CLfjCT/kNMiYv2cV4eIOHDVu/MNHwGdJLWrs+0Pc4brdkv8tWFJ6wL/H/hvYJOIWB94nK/KqZC66/4AWAxskLNvHSOiX4F11JabpJYkP2IfSdqapKpqBLBhGteMOnGtqsvlfNNrx9dz/2eTVHN1L2LerzaeJPsPge/X+R60rZNQcrdT7PfrcpLvdq/0O38CX+3LbKDupbG56/0AuKBOTOtExP2rsXvNghNE09ABWBgRi9PGz5/kTPsImAscmTYY/ic5R0qSepEcVR6Vvi4o0GB3D3CGpM6SNgbOKRDTs+n6T1dyPXuVpD6SCv2I5jODpNriUklt0nUMB8YWWOZjoJu+OsRvR1JVMgdYmR5NDy5y+19bV0S8S1KF82tJHSS1kLSVpF0KrGNnScMktSIpt/kkSb09sJLkb9RCSUP+vxUZV2Z8eazx/kfEF8DvgWskbZp+j3ZJE92q/BYYLWlzAEmbSNo3z7z3ACdI2lpSe+CCAuvtQJKkP5O0BclBTo2ngXaSRqTfu0OB7XKm3wicIqm/Eu0l7SdpnSL2p1lxgmgaziD5x1pM0lh4d82EiPiS5OjqQpJT+s2ByQCSWgN/AC6OiCkRMYXk7OOO9IesrutIqiTeIGk0vidfQOmPyt4kVQDvkfwA3kDyg7ha0iPRQ4GeJG0DdwNnR0Sh6qG7SKpQFkj6e3rEehbwJ5If5wOAR4oM4WvrSscdAaxP0l6wII1p0wLruA84DviE5JLMgyPiy4h4ieRHdBLJkW+39PPqyIrva+q5/wCnAm+TNMbPJzkTKebs49ck1YqPK7nC7O8kVaFZMT5A8uP9DEm5/q3Aei8AdiFpN3mApHxr1rMUOIjkQodPSPb1byQXMhARz6X78zvgU5IG6h+x6jO1ZqemQdDMSkTSaJJG3BMqHUtzJekVkquc7qx0LGsTn0GYWZMjafe0OquVpBEkl38/Vum41jYlSxCSbpE0R9LrOeM2lPSYpOnp+wY5085TcqPOm5L2LFVcZtYsbAu8TlLF9J8kl29nNYxbASWrYpK0G0kj0u8jolc67tfAgogYLWkUyVUg56aNoneS3PC1GUmd5dZp/bmZmVVAyc4gIuJpksa7XPsDt6efbydpPKoZf1dELE+vEJlBkizMzKxCyn2j3KYRMRsgImZL2iQd/22SywZrzCLPTStpfeIIgHXXXXeHHj16lDBcM7OmZ/LkyfMiotOq5mssd1JnXS6XWfcVETeSXApH//79Y9Kk1b0i0MyseZNUqBeEWuW+iuljSZ0B0vc56fhZfP1Oxy4kN3iZmVmFlDtBPExyByzp+0M54w9P75LtRtI528Qyx2ZmZjlKVsUk6U6SW/k3Tvt2v5Ck18l70s7J3ifpwZOIeEPSPcAUkj5fTvYVTGZmlVWyBBERR+SZNCRrZERcBlxW3+1+8cUXzJo1i2XLltV3VVYGbdu2pUuXLrRqldWzh5lVUmNppG4ws2bNokOHDnTt2pXCfZdZpUUE8+fPZ9asWXTr1q3S4ZhZHU2uq41ly5ax0UYbOTmsBSSx0UYb+WzPrJFqcgkCcHJYi/hvZdZ4NckEYWZm9dfk2iDq0sUNe4QaF6667ypJnHnmmVx55ZUAXHHFFSxevJiLLrqo3tufO3cuw4YN4/PPP2fMmDHsuuuu9V6nmVkWn0GUQJs2bbj//vuZN6/hO4+cMGECPXr04OWXX3ZyMLOScoIogaqqKkaMGMFVV131jWnvvfceQ4YMoU+fPgwZMoT3338fgGOOOYZTTz2VnXfeme7du3Pvvfd+Y9nq6mrOOeccHnnkEfr27cvSpUtp3/6rB7Tde++9HHPMMbXr++lPf8ruu+9O9+7deeqppzjuuOPYZpttaucBGDduHDvttBP9+vXjkEMOYfHixQB07dq1NsFNmjSJwYMHA/DUU0/Rt29f+vbty/bbb8+iRYsAuPzyyxkwYAB9+vThwgsvBGDJkiXss88+bLfddvTq1Yu776590J2ZrQWcIErk5JNPZuzYsSxcuPBr40eOHMnRRx/Nq6++ypFHHsmpp55aO2327Nk8++yz/PnPf2bUqFHfWGffvn255JJLOOyww6iurqZdu3YFY/jkk094/PHHueqqq9h3330544wzeOONN3jttdeorq5m3rx5XHrppYwfP56XXnqJ/v3785vf/KbgOq+44gquv/56qqureeaZZ2jXrh3jxo1j+vTpTJw4kerqaiZPnszTTz/No48+ymabbcYrr7zC66+/zl577bUaJWhmleYEUSIdO3bk6KOPZsyYMV8b//zzz/OjH/0IgB//+Mc8++yztdMOOOAAWrRoQc+ePfn444/rHcO+++6LJHr37s2mm25K7969adGiBdtuuy0zZ87khRdeYMqUKQwaNIi+ffty++238957hfvwGjRoEGeeeSZjxozh008/paqqinHjxjFu3Di23357+vXrx7Rp05g+fTq9e/dm/PjxnHvuuTzzzDOst9569d4nMyufJt9IXUmnn346/fr149hjj807T+5lnm3atKn9XPMgp/PPP5+//OUvQFLFVGj5uvcT1KyvRYsWX1t3ixYtWLFiBS1btuQHP/gBd975zcf0VlVVsXLlym+sd9SoUeyzzz488sgjDBw4kPHjxxMRnHfeeZx44onfWM/kyZN55JFHOO+88xg6dCgXXHBB3rIws8bFZxAltOGGG3LooYdy8803147beeedueuuuwAYO3Ysu+yyS8F1XHbZZVRXV2cmB4BNN92UqVOnsnLlSh544IHVim/gwIE899xzzJgxA4B//etfvPXWW0DSBjF58mQA7rvvvtpl3n77bXr37s25555L//79mTZtGnvuuSe33HJLbfvFhx9+yJw5c/joo49YZ511OOqoozjrrLN46aWXVis+M6usJn8GUcxlqaX0s5/9jOuuu652eMyYMRx33HFcfvnldOrUiVtvvbVe6x89ejTDhg1j8803p1evXrU/0sXo1KkTt912G0cccQTLly8H4NJLL2Xrrbfmwgsv5Pjjj+eXv/wl3/ve92qXufrqq3niiSdo2bIlPXv25Ic//CFt2rRh6tSp7LTTTgC0b9+eP/zhD8yYMYOzzz6bFi1a0KpVK2644YZ67auZlVfJnkldDlkPDJo6dSrbbLNNhSKyNeG/mVl5SZocEf1XNZ+rmMzMLJMThJmZZXKCMDOzTE4QZmaWyQnCzMwyOUGYmVmmJn8fBA39QJoiLgtu2bIlvXv3rh1+8MEH6dq1a8PGYWZWYk0/QVRAu3bt8t75bGa2tnAVU5ncdtttjBw5snZ42LBhPPnkk0By5/G5557LDjvswB577MHEiRMZPHgw3bt35+GHHwbgyy+/5Oyzz67tUvt3v/sdAE8++STDhg2rXe/IkSO57bbbgKTfpJ49e9KnTx/OOussIHng0MEHH8yAAQMYMGAAzz33HJC/G28za758BlECS5cupW/fvgB069ZtlX0kLVmyhMGDB/OrX/2KAw88kJ///Oc89thjTJkyheHDh7Pffvtx8803s9566/GPf/yD5cuXM2jQIIYOHZp3nQsWLOCBBx5g2rRpSOLTTz8F4LTTTuOMM85gl1124f3332fPPfdk6tSptd14Dxo0iMWLF9O2bduGKxAzWys5QZTA6lYxtW7duvZZCb1796ZNmza0atWK3r17M3PmTCB5sM+rr75a+yChhQsXMn36dFq3bp25zo4dO9K2bVtOOOEE9tlnn9qzjPHjxzNlypTa+T777DMWLVpU2433kUceyUEHHUSXLl3WZNfNrAlxFVOZ5HafDV/vQrtVq1a13Xbnds1d0y03JN1/X3vttbU9u7777rsMHTo073qrqqqYOHEiBx98MA8++GBtAlq5ciXPP/987Xo+/PBDOnTowKhRo7jppptYunQpAwcOZNq0aaUtEDNr9JwgyqRr165UV1ezcuVKPvjgAyZOnLhay++5557ccMMNfPHFFwC89dZbLFmyhC233JIpU6awfPlyFi5cyIQJEwBYvHgxCxcuZO+99+bqq6+uPaMZOnTo13qXrRmf1Y23mRVJKv+rDJp+FVMj6a120KBBdOvWjd69e9OrVy/69eu3WsufcMIJzJw5k379+hERdOrUiQcffJDNN9+cQw89lD59+rDVVlux/fbbA7Bo0SL2339/li1bRkTUPh97zJgxnHzyyfTp04cVK1aw22678dvf/jazG2+ztZEuLs+PZ67G8SvT8Nzdt1Wc/2bWkCqSIC4q+ybrdfDr7r7NzKxenCDMzCxTk0wQa3O1WXPjv5VZ49XkEkTbtm2ZP3++f3jWAhHB/PnzfVOeWSPV5K5i6tKlC7NmzWLu3LmVDsWK0LZtW9+UZ9ZINbkE0apVK7p161bpMMzM1npNrorJzMwaRkUShKQzJL0h6XVJd0pqK2lDSY9Jmp6+b1CJ2MzMLFH2BCHp28CpQP+I6AW0BA4HRgETImIrYEI6bGZmFVKpKqYqoJ2kKmAd4CNgf+D2dPrtwAEVis3MzKhAgoiID4ErgPeB2cDCiBgHbBoRs9N5ZgObZC0vaYSkSZIm+UolM7PSqUQV0wYkZwvdgM2AdSUdVezyEXFjRPSPiP6dOnUqVZhmZs1eJaqY9gDejYi5EfEFcD+wM/CxpM4A6fucCsRmZmapSiSI94GBktZR8pScIcBU4GFgeDrPcOChCsRmZmapst8oFxEvSroXeAlYAbwM3Ai0B+6RdDxJEjmk3LGZmdlXKnIndURcCFxYZ/RykrMJMzNrBHwntZmZZXKCMDOzTE4QZmaWyQnCzMwyOUGYmVkmJwgzM8tUdIKQ1KaUgZiZWeOSN0EocaikhyR9DMyUNF/Sq5L+n6TuZYzTzMzKrNAZxBPAtsDFwGYR0TkiNiLpS6ka+I2kI8sQo5mZVUChO6n3jIjldUdGxBzgbuBuSa1LFpmZmVVU3gRRNzmkbRBHAO2AuyLik4j4vMTxmZlZhazOVUzXAOumnx8sQSxmZtaIFGqkvkNS15xRGwG3kjwOdKPShmVmZpVWqA3iYuBXkt4FLgOuAh4B2gL/XYbYzMysggq1QcwADpM0GLgPeAD4fkSsLFNsZmZWQYWqmNaTdCLQHTiY5HkN4yTtVa7gzMyscgo1Uj8ICFgfuCMibgH2BwZJciO1mVkTV6gNohNwG0mbw48AImIJ8AtJXUofmpmZVVKhBHEJyd3UXwLn506IiFmlDMrMzCqvUCP1PcA9ZYzFzMwakUKN1DdI6pFnWjtJR0s6onShmZlZJRWqYroJ+KWk7wKvAnNJ2iO2AjYmaZ+4vtQBmplZZRSqYpoMHCSpI7Aj0BlYClwTEW+UKT4zM6uQQmcQAETEZ8D4MsRiZmaNiB85amZmmZwgzMwsk59JbWZmmTIThKSqnM87SnoNmJ4Obyfp2jLFZ2ZmFZLvDOJESTuln8cAw4D5ABHxCrB7GWIzM7MKypcgfgccWTNPRLxXZ/qXpQvJzMwag8zLXCNiBTAyHfxA0o5ASGoJnAK8Vab4zMysQopppP4pcCawBfAxMBA4qZRBmZlZ5a3yRjmge0QcnjtC0kDSNgkzM2uaijmD+J+Mce6Dycysict7BpG2O+wEdJJ0as6kjkCr+mxU0voknQH2AgI4DngTuBvoCswEDo2IT+qzHTMzW3OFziDWJem1tYrk6XI1r8+BQ+q53WuARyOiB7AdMBUYBUyIiK2ACemwmZlVSKHeXJ8AnpB0a0S801AbTHuH3Q04Jt3O58DnkvYHBqez3Q48CZzbUNs1M7PVU6iK6cqI+BlwpaSoOz0iDlrDbXYnebbErZK2AyYDpwGbRsTsdN2zJW2SJ64RwAiALbbYYg1DMDOzVSl0FdPd6ft1JdhmP+CUiHhR0jWsRnVSRNwI3AjQv3//byQuMzNrGIWqmCam7xMaeJuzgFkR8WI6fC9JgvhYUuf07KEzMKeBt2tmZquhUBXTyyRXGGWKiH5rssGI+KekDyR9NyLeBIYAU9LXcGB0+v7QmqzfzMwaRqEqpv8o4XZPAcZKag28AxxLckXVPZKOB96n/ldKmZlZPRSqYnq7VBuNiGqgf8akIaXappmZrZ5V3kktaYCkFyQtlLRM0nJJn5UjODMzq5xi+mL6H+Ao4C5gR5L7FzYvYUxmZtYIFNMXU4u0MbkqIr6IiP8F9ihxXGZmVmHFnEEsSRuTX5H0S2A20L60YZmZWaUVcwZxTDrfSJInyW0FHFzCmMzMrBEoJkHsHRHLIuLTiPhFRJwK7FnqwMzMrLKKSRDHZYw7vqEDMTOzxqXQndSHAYcD3SXdnzOpI/BpqQMzM7PKKtRIPZHksaJd+PoT5BYBL5cyKDMzq7xCd1K/K+l9YEEJOuwzM7NGrmAbRER8SfIwn45lisfMzBqJYu6DWExyD8Q4YEnNyIg4s2RRmZlZxRWTIManLzMza0ZWmSAi4mZJVcC/paNmRMSK0oZlZmaVtsoEIWlX4A7gQ0DAtyT9OCKeK3VwZmZWOcVUMV1Fcjf1FABJ25AkjKznOZiZWRNRzJ3UrWuSA0BETAValy4kMzNrDIo5g3hJ0u9IzhoAjsQ3ypmZNXnFJIiTgFOBc0jaIJ4Gri1lUGZmVnnFXMW0DPh1+jIzs2aimGdS7yXpH5LmSFog6RNJC8oRnJmZVU4xVUzXAYcCrwErSxuOmZk1FsUkiFlAdUQ4OZiZNSPFJIhzgD9JehJYXjMyIsaUKigzM6u8YhLExcAXwPq4isnMrNkoJkFsEhE7lDwSMzNrVIq5k3qCpO+XPBIzM2tUikkQPwHGS1rsy1zNzJqPYqqYNi55FGZWGVL5txlR/m3aGinmTuovyxGImZk1LsVUMZmZWTPkBGFmZpmKShCSBko6Ov28kaQtShuWmZlVWjGPHP05MAj4DvB7oC3wR2CX0oZmZmaVVMwZxH8AewNLACLiQ6BjKYMyM7PKKyZBLI+IAAJA0joNsWFJLSW9LOnP6fCGkh6TND1936AhtmNmZmummARxv6TrgfUkHQuMA25pgG2fBkzNGR4FTIiIrYAJ6bCZmVXIKhNERPwK+DPwMLAdcFlEXF2fjUrqAuwD3JQzen/g9vTz7cAB9dmGmZnVTzF3UhMRfwX+2oDbvZqkG/EOOeM2jYjZ6fZmS9oka0FJI4ARAFts4YupzMxKJe8ZRE2fS/lea7pBScOAORExeU2Wj4gbI6J/RPTv1KnTmoZhZmarUOgMYmNAwIXAXOCOdPhIoD4N1YOA/STtTXLJbEdJfwA+ltQ5PXvoDMypxzbMzKye8p5BRMSXEbECGBoRYyLik4hYEBHXAgeu6QYj4ryI6BIRXYHDgccj4iiSNo7h6WzDgYfWdBtmZlZ/xVzFFJIOk5JuHyUdVqJYRgM/kDQd+EE6bGZmFVJMI/WPgGuBGyStBF4gqWaqt4h4Engy/TwfGNIQ6zUzs/orprvvd0guSTUzs2bEvbmamVkmJwgzM8vkBGFmZplWmSAk3SqpQ85wF0njShuWmZlVWjFnEJOAiZKGpp31PQHcUNqwzMys0oq5iul6Sa+QJIZ5QL+aPpPMzKzpKqaK6QiS7r2PA/4A/ElSr1IHZmZmlVXMjXJHArtFxD8BJO0EjCXp+tvMzJqoYqqYhtUZfl7S90oXklnzpItV9m1G2bdoa5NVJghJN+aZNKKBYzEzs0akmCqmCTmf25L05PpBacIxM7PGopgqprtzhyXdATxWsojMzKxRWJM7qbsBWzZ0IGZm1rgU0wbxCV+1ZbUAFgCjShmUmZlVXjFtEBvnfF4ZEb7wwcysGSimDeJLSesB3wHapg+WIyL+XuLYzMysgoqpYjoeOBP4NvAaMIDkqXKDSxqZmZlVVDGN1KcD/YGZEbErsAPgvpjMzJq4YhLEsohYCiCpdUS8AfQobVhmZlZpeauYJFVFxApgtqT1gT8Bf5O0APi4XAGamVllFGqDmEjStfd+6fAvJA0B1gP+UvLIzMysogoliG/0HBYRE7JmNDOzpqdQgugk6cx8EyPiNyWIx8zMGolCCaIl0J6MMwkzM2v6CiWI2RFxSdkiMTOzRqXQZa4+czAza8YKJYghZYvCzMwanbwJIiIWlDMQMzNrXNbkeRBmZtYMOEGYmVkmJwgzM8vkBGFmZpmcIMzMLFPZE4SkzSU9IWmqpDcknZaO31DSY5Kmp+8blDs2MzP7SiXOIFYAP4uIbYCBwMmSegKjgAkRsRUwIR02M7MKKXuCiIjZEfFS+nkRMJXkcab7A7ens90OHFDu2MzM7CsVbYOQ1BXYHngR2DQiZkOSRIBN8iwzQtIkSZPmzp1brlDNzJqdiiUISe2B+4DTI+KzYpeLiBsjon9E9O/UqVPpAjQza+YqkiAktSJJDmMj4v509MeSOqfTOwNzKhGbmZklKnEVk4Cbgal1Hjr0MDA8/TwceKjcsZmZ2VcKPQ+iVAYBPwZek1SdjvsvYDRwj6TjgfeBQyoQm5mZpcqeICLiWfI/a8JdjJuZNRKVOIMwqz+V+XlWEeXdnlkj4K42zMwskxOEmZllcoIwM7NMThBmZpbJCcLMzDI5QZiZWSYnCDMzy+QEYWZmmZwgzMwskxOEmZllcoIwM7NMThBmZpbJCcLMzDI5QZiZWSYnCDMzy+QEYWZmmfzAIKs3XVzmh/cAfnyPWen5DMLMzDI5QZiZWSYnCDMzy+QEYWZmmZwgzMwskxOEmZllcoIwM7NMThBmZpbJCcLMzDI5QZiZWSYnCDMzy+QEYWZmmZwgzMwskxOEmZllatbdfZe7m+q4qKybSzfqjrHNbM34DMLMzDI5QZiZWaZGlyAk7SXpTUkzJI2qdDxmZs1Vo0oQkloC1wM/BHoCR0jqWdmozMyap0aVIIAdgRkR8U5EfA7cBexf4ZjMzJqlxnYV07eBD3KGZwHfy51B0ghgRDq4WNKbZYqt3gQbA/PKu9HyXqlVLmUvS5djA2606ZXlWliOWxYzU2NLEFl7/LXrNCPiRuDG8oTTsCSCyTEEAAAERklEQVRNioj+lY6jKXBZNgyXY8NoquXY2KqYZgGb5wx3AT6qUCxmZs1aY0sQ/wC2ktRNUmvgcODhCsdkZtYsNaoqpohYIWkk8DegJXBLRLxR4bAa0lpZNdZIuSwbhsuxYTTJclS4KwYzM8vQ2KqYzMyskXCCMDOzTE4Qa0jSLZLmSHo9Z1xfSS9IqpY0SdKOeZb9b0mvpvONk7RZnelbSFos6axS70elSWoraaKkVyS9IenidPxFkj5My6ha0t55li84XzMry5mSXqv5/qXjiirHdN5T0m5u3pD06zrTmk051iVpfUn3SpomaaqknVanXNdmboNYQ5J2AxYDv4+IXum4ccBVEfHX9AtzTkQMzli2Y0R8ln4+FegZESflTL8PWAm8GBFXlH5vKkeSgHUjYrGkVsCzwGnAXsDiVe2/pIsKzdfMynIm0D8i5uWMu4jiynF34Hxgn4hYLmmTiJiTM73ZlGNdkm4HnomIm9KrK9cBTmcV5ZqW/cyIuK0sgZZAo7qKaW0SEU9L6lp3NNAx/bweee7hqEkOqXXJuRlQ0gHAO8CShoq1MYvkCGVxOtgqfTXIUUtzK8t6+ikwOiKWA9RJDs22HCV1BHYDjgFIuwD6XE3wbvAsrmJqWKcDl0v6ALgCOC/fjJIuS+c7ErggHbcucC5wcRlibTQktZRUDcwBHouIF9NJI9OquFskbVBgFd+Yr5mWZQDjJE1Ou6SpUUw5bg3sKulFSU9JGgDNthxzdQfmArdKelnSTWmZQPHfz7WWE0TD+ilwRkRsDpwB3Jxvxog4P51vLDAyHX0xSRXV4nzLNUUR8WVE9CW5c35HSb2AG4DvAH2B2cCVeRbPN19zLMtBEdGPpDfkk9Nq0GLLsQrYABgInA3ck1b/NcdyzFUF9ANuiIjtSc6iRpGnXCX1rmmXAE4CLslpp9ioIntQHxHh1xq+gK7A6znDC/mqXUfAZ+nnW4Fq4JGMdWxZsw7gGWBm+voUWACMrPR+lrlMLwTOylfOqyjLri7L2rK4aHXKEXgUGJwz79tAJ5cj3yJpR6gZ3hX4S75yzfgbHFPpfajPy20QDesj4N+BJ4HvA9MBIuLY3JkkbRUR09PB/YBp6Xy75sxzEUkj2HUlj7qCJHUCvoiITyW1A/YAfiWpc0TMTmc7EHgdMssy33zNqizTao8WEbEo/TyU5Oi1qHIEHiT5zj4paWugNTCvuZVjXRHxT0kfSPpuRLwJDAGm5CvXpsYJYg1JuhMYDGwsaRbJke9PgGskVQHL+Kpb8rpGS/ouyVUh75GcijZXnYHblTwsqgVwT0T8WdIdkvqS1KvPBE7Ms/yvi5yvqdsUeCBtPK0C/hgRj65GOd4C3KLksu3PgeGRHgYbpwBj0yuY3gGOBcY0h++dL3M1M7NMbqQ2M7NMThBmZpbJCcLMzDI5QZiZWSYnCDMzy+QEYWZmmZwgzMws0/8BdFkxE9+N17gAAAAASUVORK5CYII=\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "plt.title(\"Taux de mortalité après 20 ans - Total \")\n", + "x_tot=[\"Non-fumeuses\", \"Fumeuses\"]\n", + "y_tot=[raw_data.taux_mortalité[0], raw_data.taux_mortalité[1]]\n", + "plt.ylim([0,100])\n", + "plt.bar(x_tot, y_tot, color=[\"green\", \"red\"] )\n", + "plt.ylabel(\"Taux de mortalité en %\")\n", + "plt.show()\n", + "\n", + "\n", + "labels = ['18-34', '35-54', '55-64', '65+']\n", + "y_nf=[df_18_34.taux_mortalité[0], df_35_54.taux_mortalité[0], df_55_64.taux_mortalité[0], df_65.taux_mortalité[0]]\n", + "y_f=[df_18_34.taux_mortalité[1], df_35_54.taux_mortalité[1], df_55_64.taux_mortalité[1], df_65.taux_mortalité[1]]\n", + "\n", + "x = np.arange(len(labels)) # the label locations\n", + "width = 0.35 # the width of the bars\n", + "\n", + "fig, ax = plt.subplots()\n", + "rects1 = ax.bar(x - width/2, y_nf, width, color=\"green\", label='Non-fumeuses')\n", + "rects2 = ax.bar(x + width/2, y_f, width, color=\"red\",label='Fumeuses')\n", + "\n", + "ax.set_ylabel('Taux de mortalité (%)')\n", + "ax.set_title(\"Taux d emortalité par tranche d'âge\")\n", + "ax.set_xticks(x)\n", + "ax.set_ylim([0,100])\n", + "ax.set_xticklabels(labels)\n", + "ax.legend()" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "D'après les graphiques précédents, si on raisonne par classe d’âge, dans chaque tranche la mortalité chez les fumeuses a été supérieure à celle des non-fumeuses. C'est rassurant, mais comment les chiffres peuvent-ils s’inverser quand on groupe tout le monde ?\n", + "\n", + "En fait dans la population initiale, il y avait plus de femmes âgées chez les non-fumeuses que chez les fumeuses. Et même si dans chaque tranche d’âge les non-fumeuses meurent moins, cet effet est compensé par le fait que la tranche d’âge « élevée » est sur-représentée chez les non-fumeuses… qui donc en moyenne meurent plus !" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": {}, + "outputs": [], + "source": [] + } + ], + "metadata": { + "kernelspec": { + "display_name": "Python 3", + "language": "python", + "name": "python3" + }, + "language_info": { + "codemirror_mode": { + "name": "ipython", + "version": 3 + }, + "file_extension": ".py", + "mimetype": "text/x-python", + "name": "python", + "nbconvert_exporter": "python", + "pygments_lexer": "ipython3", + "version": "3.6.4" + } + }, + "nbformat": 4, + "nbformat_minor": 2 +} diff --git a/module3/exo3/Subject6_smoking.csv b/module3/exo3/Subject6_smoking.csv new file mode 100644 index 0000000000000000000000000000000000000000..6c47065bc7aa53c0c92eef463cf16384eb028aaa --- /dev/null +++ b/module3/exo3/Subject6_smoking.csv @@ -0,0 +1,1315 @@ +"Smoker","Status","Age" +"Yes","Alive",21 +"Yes","Alive",19.3 +"No","Dead",57.5 +"No","Alive",47.1 +"Yes","Alive",81.4 +"No","Alive",36.8 +"No","Alive",23.8 +"Yes","Dead",57.5 +"Yes","Alive",24.8 +"Yes","Alive",49.5 +"Yes","Alive",30 +"No","Dead",66 +"Yes","Alive",49.2 +"No","Alive",58.4 +"No","Dead",60.6 +"No","Alive",25.1 +"No","Alive",43.5 +"No","Alive",27.1 +"No","Alive",58.3 +"Yes","Alive",65.7 +"No","Dead",73.2 +"Yes","Alive",38.3 +"No","Alive",33.4 +"Yes","Dead",62.3 +"No","Alive",18 +"No","Alive",56.2 +"Yes","Alive",59.2 +"No","Alive",25.8 +"No","Dead",36.9 +"No","Alive",20.2 +"Yes","Alive",34.6 +"Yes","Alive",51.9 +"Yes","Alive",49.9 +"No","Alive",19.4 +"No","Alive",56.9 +"Yes","Alive",46.7 +"Yes","Alive",44.4 +"Yes","Alive",29.5 +"Yes","Dead",33 +"Yes","Alive",35.6 +"Yes","Alive",39.1 +"No","Dead",69.7 +"Yes","Alive",35.7 +"No","Dead",75.8 +"No","Alive",25.3 +"No","Dead",83 +"Yes","Dead",44.3 +"No","Alive",18.5 +"Yes","Alive",37.5 +"Yes","Alive",22.1 +"No","Alive",82.8 +"No","Alive",45 +"No","Dead",73.3 +"Yes","Alive",39 +"No","Alive",28.4 +"No","Dead",73.7 +"Yes","Alive",40.1 +"No","Alive",51.2 +"No","Alive",22.9 +"No","Alive",41.9 +"Yes","Alive",58.1 +"Yes","Alive",37.3 +"No","Alive",41.7 +"Yes","Dead",36.3 +"Yes","Dead",80.7 +"Yes","Alive",33 +"Yes","Alive",38.6 +"Yes","Alive",27.9 +"No","Alive",47.6 +"No","Dead",77.6 +"No","Dead",58.1 +"Yes","Alive",26.2 +"No","Alive",45.4 +"No","Alive",62.4 +"No","Alive",62.5 +"No","Alive",39.5 +"No","Alive",27.6 +"Yes","Alive",31.4 +"No","Dead",85 +"No","Alive",18.9 +"No","Alive",35.3 +"Yes","Alive",25.4 +"No","Dead",72.8 +"Yes","Alive",58.3 +"No","Alive",27.3 +"No","Dead",55.9 +"No","Alive",32.8 +"Yes","Dead",53.6 +"No","Alive",55.9 +"Yes","Alive",48 +"Yes","Alive",56.1 +"No","Alive",18.3 +"Yes","Alive",20.2 +"No","Alive",62.8 +"Yes","Alive",18.6 +"No","Alive",46.3 +"No","Alive",36 +"Yes","Dead",55.5 +"Yes","Alive",18.6 +"No","Dead",65.7 +"No","Dead",76.5 +"Yes","Dead",61 +"No","Alive",26.8 +"Yes","Alive",47.6 +"No","Dead",70.5 +"No","Dead",81.8 +"Yes","Alive",32.5 +"No","Alive",23 +"No","Alive",83.7 +"Yes","Dead",62.8 +"Yes","Alive",45.9 +"No","Alive",59.9 +"Yes","Dead",66.5 +"No","Alive",47.5 +"No","Dead",89.3 +"No","Alive",57.2 +"Yes","Alive",21.3 +"Yes","Alive",34 +"No","Alive",59.5 +"Yes","Alive",50.1 +"No","Alive",56.1 +"Yes","Alive",30.6 +"Yes","Dead",63.8 +"Yes","Alive",27.4 +"Yes","Alive",32.5 +"No","Alive",22.5 +"Yes","Alive",24.2 +"No","Alive",56.8 +"Yes","Alive",28.9 +"Yes","Dead",87.8 +"Yes","Alive",19.4 +"No","Dead",88.4 +"Yes","Dead",35.7 +"No","Alive",33.7 +"No","Dead",62.3 +"Yes","Alive",63.6 +"Yes","Dead",71.7 +"No","Alive",48.5 +"No","Alive",82 +"Yes","Dead",40.8 +"Yes","Alive",31.3 +"No","Alive",25.1 +"No","Alive",24.4 +"No","Alive",32.2 +"No","Alive",53.1 +"No","Alive",47.5 +"No","Dead",26.3 +"No","Dead",66 +"Yes","Alive",41 +"No","Dead",58.3 +"No","Dead",86.8 +"Yes","Alive",49.7 +"No","Alive",21 +"No","Dead",75.8 +"No","Alive",50.5 +"No","Dead",63.5 +"Yes","Alive",33.1 +"No","Alive",30.7 +"Yes","Dead",59.4 +"No","Alive",67.2 +"Yes","Alive",20.7 +"No","Alive",41.6 +"No","Alive",37.7 +"No","Dead",89.7 +"No","Alive",45.2 +"Yes","Dead",59.9 +"No","Alive",55.4 +"Yes","Alive",44.8 +"No","Alive",18.3 +"No","Dead",76.7 +"Yes","Dead",62.3 +"Yes","Dead",48.4 +"No","Alive",82.7 +"No","Alive",27 +"No","Alive",30.9 +"No","Dead",82.9 +"No","Alive",40.5 +"Yes","Alive",56.8 +"No","Alive",26.5 +"Yes","Alive",41.9 +"Yes","Alive",29.9 +"No","Dead",75 +"No","Dead",66.3 +"No","Dead",87 +"No","Dead",79.5 +"Yes","Alive",51.7 +"No","Alive",30 +"No","Alive",78.2 +"No","Alive",80 +"Yes","Alive",62.9 +"Yes","Dead",78.3 +"Yes","Alive",49.8 +"Yes","Alive",36.5 +"No","Dead",60.1 +"Yes","Alive",62 +"No","Alive",19 +"No","Dead",87.6 +"No","Alive",24.3 +"No","Alive",22.2 +"Yes","Dead",68.4 +"No","Alive",32.5 +"Yes","Dead",89.2 +"Yes","Alive",44.5 +"Yes","Alive",43.3 +"Yes","Dead",57.4 +"Yes","Alive",45.6 +"Yes","Alive",18.1 +"Yes","Dead",63.4 +"No","Alive",67 +"Yes","Alive",55.6 +"Yes","Alive",51.9 +"No","Alive",23.3 +"Yes","Dead",57.6 +"No","Alive",38.4 +"No","Dead",35.2 +"Yes","Alive",26.2 +"No","Alive",60.3 +"Yes","Alive",48.7 +"Yes","Alive",23.7 +"Yes","Alive",46.9 +"Yes","Alive",18 +"No","Dead",65.8 +"Yes","Alive",33 +"Yes","Dead",83.1 +"No","Dead",58.4 +"Yes","Alive",23.2 +"No","Alive",66.7 +"No","Alive",58.8 +"Yes","Alive",56.7 +"Yes","Alive",21.5 +"No","Dead",78.3 +"No","Alive",51.5 +"Yes","Alive",63.5 +"Yes","Alive",57.8 +"Yes","Alive",59.5 +"Yes","Dead",47.8 +"Yes","Alive",53.9 +"Yes","Alive",45.5 +"Yes","Alive",24.2 +"No","Alive",63.9 +"Yes","Alive",37.5 +"No","Alive",20.6 +"No","Alive",22.9 +"No","Alive",46.1 +"No","Alive",49.6 +"No","Alive",31.4 +"No","Alive",25.9 +"Yes","Alive",46.8 +"Yes","Dead",81 +"No","Dead",84.3 +"No","Alive",30.8 +"Yes","Alive",52.4 +"No","Alive",20.1 +"Yes","Dead",58.9 +"Yes","Alive",72.1 +"No","Alive",19.6 +"No","Alive",52.6 +"Yes","Alive",35 +"Yes","Dead",35.4 +"No","Dead",55.1 +"Yes","Alive",23.7 +"No","Alive",49.1 +"Yes","Alive",39.7 +"Yes","Alive",33.7 +"No","Dead",66.4 +"No","Alive",24.2 +"No","Dead",67.2 +"No","Alive",19.4 +"No","Alive",52.4 +"No","Dead",58.6 +"Yes","Alive",36.2 +"Yes","Alive",38.8 +"Yes","Alive",47.9 +"No","Alive",36.5 +"Yes","Alive",24.3 +"No","Alive",38.8 +"No","Alive",38.4 +"No","Alive",55.3 +"No","Dead",87.7 +"Yes","Dead",56.7 +"No","Alive",74.1 +"Yes","Alive",62.3 +"No","Alive",18.5 +"Yes","Dead",59.3 +"No","Alive",39.8 +"Yes","Dead",55 +"No","Alive",42.8 +"No","Alive",34.2 +"Yes","Alive",33.7 +"No","Alive",30.6 +"No","Dead",81.7 +"Yes","Alive",62 +"No","Alive",61.3 +"Yes","Alive",58.5 +"No","Alive",41.6 +"Yes","Alive",60.6 +"No","Alive",52.9 +"Yes","Alive",34 +"No","Dead",52.4 +"No","Alive",38.5 +"No","Alive",23.7 +"Yes","Alive",38.7 +"No","Alive",49.3 +"No","Alive",59.5 +"Yes","Alive",26.2 +"Yes","Dead",65.8 +"Yes","Alive",44.3 +"No","Alive",31.9 +"No","Dead",47.9 +"Yes","Alive",57.7 +"Yes","Dead",36.5 +"Yes","Alive",36.3 +"No","Dead",56.1 +"No","Alive",21.1 +"Yes","Alive",22.7 +"No","Alive",19.7 +"Yes","Dead",60.1 +"Yes","Dead",77.6 +"No","Dead",67.6 +"No","Alive",49.3 +"Yes","Alive",37 +"No","Dead",79.9 +"No","Dead",56.3 +"Yes","Alive",20.2 +"No","Alive",31.1 +"Yes","Alive",40.9 +"Yes","Dead",35.2 +"No","Alive",24.5 +"Yes","Alive",35 +"Yes","Alive",36.3 +"Yes","Dead",34.3 +"Yes","Alive",20.5 +"Yes","Alive",29 +"Yes","Dead",74.1 +"Yes","Alive",43.6 +"Yes","Alive",33 +"Yes","Dead",42.3 +"No","Dead",63.2 +"No","Alive",53.2 +"Yes","Alive",53.7 +"No","Alive",62.7 +"Yes","Alive",39 +"Yes","Alive",39.3 +"No","Dead",47 +"No","Alive",35.8 +"No","Alive",49.4 +"No","Alive",20.7 +"No","Dead",76.7 +"Yes","Alive",31.3 +"No","Alive",20.1 +"No","Alive",56.3 +"No","Alive",51.3 +"No","Dead",85.2 +"Yes","Alive",25.2 +"Yes","Alive",20.2 +"No","Alive",58.1 +"No","Alive",49.9 +"No","Dead",79.4 +"Yes","Alive",31.6 +"No","Alive",31.6 +"No","Alive",55.4 +"No","Alive",41.6 +"No","Dead",74.6 +"No","Dead",81.3 +"No","Dead",71.4 +"Yes","Alive",56.4 +"Yes","Alive",39.7 +"Yes","Alive",59.1 +"No","Alive",20.7 +"No","Alive",89.7 +"Yes","Alive",61.8 +"Yes","Alive",26.8 +"Yes","Dead",44.3 +"Yes","Alive",36.1 +"Yes","Alive",22.1 +"No","Alive",33.5 +"Yes","Alive",44.5 +"Yes","Alive",24.1 +"No","Dead",72.5 +"Yes","Dead",57.7 +"Yes","Alive",58.7 +"Yes","Alive",40.7 +"No","Dead",75.6 +"No","Dead",74.1 +"No","Alive",37 +"Yes","Alive",31.6 +"Yes","Alive",34.7 +"Yes","Dead",38.5 +"No","Alive",22 +"Yes","Alive",39.3 +"No","Dead",61.2 +"No","Alive",37.2 +"No","Alive",25.7 +"No","Dead",88.8 +"No","Dead",65.6 +"Yes","Dead",58 +"No","Alive",20.4 +"Yes","Alive",44.4 +"No","Alive",46.2 +"No","Alive",29.7 +"Yes","Alive",51.9 +"Yes","Dead",43.7 +"Yes","Dead",61.1 +"No","Dead",78 +"No","Alive",26.8 +"No","Alive",63 +"No","Dead",82.3 +"No","Alive",32.2 +"No","Dead",57.2 +"No","Alive",36.7 +"Yes","Dead",67.5 +"No","Alive",59.2 +"No","Dead",55.6 +"No","Dead",86.2 +"Yes","Dead",79.1 +"No","Dead",75.1 +"No","Alive",52 +"No","Alive",55.3 +"Yes","Alive",40.3 +"Yes","Dead",56.9 +"Yes","Alive",53.9 +"No","Alive",26.5 +"Yes","Alive",33 +"No","Alive",19.7 +"No","Alive",23.3 +"Yes","Dead",75.9 +"Yes","Dead",35.5 +"No","Alive",31.9 +"No","Alive",34.2 +"No","Dead",76.2 +"Yes","Alive",21.7 +"Yes","Alive",50.6 +"No","Alive",25.7 +"No","Alive",24.2 +"No","Alive",42.2 +"Yes","Alive",49.2 +"Yes","Alive",33.6 +"Yes","Alive",49.5 +"Yes","Alive",61.6 +"No","Alive",23 +"No","Alive",24.3 +"No","Alive",23.6 +"No","Alive",57.1 +"Yes","Alive",32.5 +"No","Dead",83.1 +"Yes","Alive",21.8 +"Yes","Alive",43.2 +"Yes","Alive",26.6 +"Yes","Alive",45.7 +"Yes","Alive",18.1 +"Yes","Dead",45.6 +"Yes","Alive",29.7 +"Yes","Dead",73.9 +"No","Alive",56.4 +"No","Alive",55.6 +"Yes","Alive",55.1 +"No","Alive",80.8 +"Yes","Alive",29.7 +"No","Alive",25.7 +"No","Alive",52.8 +"No","Dead",81.3 +"Yes","Dead",80.5 +"Yes","Alive",34.3 +"No","Dead",59 +"No","Alive",42.5 +"No","Alive",76.9 +"Yes","Alive",33.3 +"No","Alive",20.6 +"Yes","Dead",86.8 +"No","Alive",33.1 +"No","Dead",80.2 +"Yes","Alive",30.5 +"No","Alive",31.9 +"No","Alive",19.8 +"No","Dead",84.5 +"No","Alive",56 +"No","Alive",50.3 +"No","Alive",56.8 +"Yes","Dead",60.7 +"Yes","Alive",27.6 +"Yes","Alive",32.9 +"No","Alive",56.2 +"Yes","Dead",63.4 +"No","Alive",86.9 +"No","Dead",79.9 +"No","Alive",41.5 +"Yes","Alive",45.3 +"Yes","Alive",63 +"No","Dead",77.2 +"No","Dead",69.4 +"No","Alive",49 +"No","Alive",44.7 +"Yes","Alive",27.7 +"Yes","Dead",62.3 +"No","Dead",70.7 +"No","Alive",38 +"Yes","Alive",44.3 +"No","Alive",32.3 +"Yes","Alive",56.1 +"Yes","Alive",58 +"No","Dead",82.9 +"Yes","Alive",44.4 +"No","Alive",24.9 +"Yes","Alive",63.1 +"No","Alive",35.9 +"Yes","Alive",31.1 +"No","Alive",24 +"No","Dead",88.5 +"Yes","Alive",39.5 +"No","Alive",35.6 +"No","Dead",82.4 +"No","Dead",63.8 +"No","Alive",87.4 +"No","Alive",37.2 +"No","Dead",69.5 +"No","Dead",25.3 +"Yes","Alive",59.6 +"Yes","Dead",35.7 +"Yes","Dead",56.6 +"Yes","Alive",34.5 +"Yes","Alive",58.6 +"Yes","Dead",78.2 +"Yes","Alive",48.3 +"Yes","Alive",25.4 +"Yes","Alive",74.1 +"Yes","Dead",88.7 +"No","Alive",68.4 +"No","Alive",33.4 +"No","Alive",36.5 +"No","Alive",25.5 +"Yes","Alive",21.2 +"Yes","Dead",61.8 +"Yes","Alive",38 +"No","Alive",35.1 +"No","Alive",38 +"Yes","Dead",36.2 +"Yes","Dead",87.9 +"No","Dead",76.1 +"No","Alive",59.4 +"No","Alive",18.9 +"Yes","Alive",53.3 +"Yes","Dead",82.6 +"Yes","Alive",45.3 +"No","Dead",86.3 +"Yes","Dead",63.2 +"No","Dead",88.1 +"Yes","Alive",36.1 +"No","Dead",71 +"Yes","Alive",62.1 +"Yes","Dead",55.3 +"No","Alive",52.2 +"No","Alive",25.6 +"No","Alive",33 +"No","Dead",75.3 +"Yes","Alive",21.3 +"Yes","Dead",76.9 +"No","Alive",30 +"No","Dead",77.5 +"Yes","Dead",75.2 +"No","Dead",83.9 +"Yes","Alive",53 +"No","Alive",62.4 +"Yes","Alive",43.7 +"Yes","Alive",50.9 +"No","Dead",29.8 +"Yes","Alive",32.8 +"Yes","Alive",50.7 +"Yes","Dead",66.1 +"No","Alive",33.5 +"Yes","Alive",27.2 +"No","Dead",56.2 +"Yes","Alive",38.1 +"Yes","Dead",66.8 +"Yes","Dead",55.2 +"No","Alive",51.6 +"Yes","Alive",50.9 +"No","Alive",41.4 +"No","Dead",65.4 +"No","Dead",67.7 +"No","Alive",37.8 +"Yes","Alive",42.5 +"No","Alive",23.9 +"No","Alive",60.1 +"Yes","Alive",26.6 +"Yes","Alive",23.3 +"No","Dead",75.6 +"No","Dead",72.1 +"Yes","Alive",34.8 +"No","Dead",55.3 +"Yes","Alive",28.2 +"No","Dead",79.3 +"Yes","Alive",38.5 +"Yes","Alive",41 +"No","Alive",60.7 +"No","Alive",51.8 +"Yes","Alive",25.7 +"No","Dead",62.7 +"No","Alive",23.7 +"No","Alive",23.4 +"No","Alive",56.5 +"No","Alive",28.4 +"No","Alive",42.8 +"No","Dead",83.5 +"No","Alive",36.8 +"Yes","Alive",43.8 +"Yes","Alive",59 +"No","Alive",25.5 +"No","Dead",47.2 +"Yes","Alive",23.5 +"No","Alive",19.4 +"No","Dead",68.5 +"Yes","Alive",43.4 +"No","Alive",19.5 +"Yes","Alive",62.2 +"Yes","Alive",31.1 +"No","Alive",19.2 +"No","Dead",61.9 +"No","Alive",27.6 +"Yes","Alive",30.2 +"Yes","Alive",59 +"Yes","Alive",49.2 +"No","Alive",40.3 +"Yes","Alive",45.4 +"No","Alive",29.4 +"No","Alive",36.8 +"No","Alive",29.6 +"Yes","Dead",58.6 +"No","Dead",29.3 +"No","Alive",40 +"Yes","Alive",21.7 +"Yes","Alive",40.3 +"No","Dead",81.6 +"Yes","Alive",22.9 +"Yes","Alive",42.7 +"Yes","Alive",40.4 +"No","Dead",85.7 +"Yes","Alive",32.2 +"Yes","Alive",19.3 +"Yes","Alive",24.3 +"No","Alive",55.8 +"No","Alive",28.8 +"No","Alive",55.1 +"No","Alive",31.1 +"Yes","Alive",31.8 +"No","Alive",46.7 +"Yes","Alive",36.3 +"No","Alive",23.3 +"Yes","Alive",50.7 +"Yes","Alive",40.8 +"Yes","Alive",36.9 +"Yes","Dead",81.8 +"No","Alive",55.2 +"Yes","Dead",62.4 +"No","Dead",78.4 +"Yes","Alive",18 +"No","Dead",85.7 +"Yes","Alive",43 +"Yes","Dead",88.3 +"Yes","Alive",26.2 +"Yes","Alive",52.7 +"No","Dead",81.9 +"No","Alive",44.4 +"No","Alive",71.8 +"No","Alive",35.4 +"No","Dead",71.4 +"No","Dead",55.9 +"Yes","Dead",46.6 +"No","Dead",65.6 +"No","Alive",56.2 +"No","Dead",57.9 +"Yes","Alive",43.5 +"No","Alive",22.6 +"No","Alive",27.2 +"No","Alive",27.2 +"No","Dead",20.2 +"Yes","Dead",60.2 +"Yes","Alive",55.2 +"No","Alive",39.6 +"No","Alive",24.5 +"No","Alive",36.7 +"Yes","Alive",24.2 +"No","Dead",73.3 +"No","Alive",26.6 +"No","Alive",41.7 +"No","Dead",42.6 +"No","Alive",18.6 +"Yes","Alive",31.3 +"No","Alive",51.6 +"No","Alive",19 +"No","Dead",72.6 +"No","Alive",35.7 +"No","Alive",44.1 +"No","Alive",58.3 +"Yes","Dead",65.6 +"No","Alive",62.3 +"Yes","Alive",57.4 +"No","Alive",26 +"No","Dead",85.7 +"No","Dead",47.3 +"Yes","Dead",62.1 +"Yes","Dead",66.1 +"Yes","Alive",18.5 +"Yes","Alive",24.6 +"Yes","Alive",48.3 +"Yes","Alive",28.8 +"No","Alive",52.2 +"No","Dead",85.5 +"No","Dead",58.4 +"Yes","Alive",38.1 +"Yes","Alive",27.7 +"No","Alive",42.1 +"Yes","Alive",47.9 +"No","Dead",67.4 +"No","Alive",29 +"No","Alive",29.4 +"No","Alive",21.4 +"No","Alive",41.5 +"No","Alive",74 +"No","Alive",42 +"No","Dead",68.1 +"Yes","Alive",21.5 +"No","Dead",58.5 +"No","Alive",32.8 +"Yes","Alive",37.7 +"No","Alive",55.5 +"No","Dead",78.7 +"No","Alive",31 +"Yes","Dead",51.6 +"No","Dead",66.6 +"No","Dead",40 +"Yes","Alive",52.1 +"Yes","Alive",30.4 +"No","Alive",38.1 +"Yes","Alive",23.1 +"Yes","Dead",57.9 +"Yes","Alive",25.2 +"No","Dead",76.2 +"No","Alive",63.4 +"No","Alive",21 +"Yes","Alive",45.5 +"No","Alive",46.5 +"No","Alive",48.1 +"No","Alive",32.4 +"Yes","Alive",40.1 +"No","Alive",23.4 +"Yes","Alive",62.1 +"No","Alive",45.1 +"Yes","Dead",53.6 +"No","Dead",60.6 +"No","Alive",83 +"No","Alive",55.5 +"No","Alive",41.8 +"No","Dead",40.1 +"Yes","Alive",24.4 +"Yes","Dead",62.7 +"Yes","Alive",23.7 +"No","Alive",84.9 +"Yes","Dead",50.2 +"No","Alive",40 +"Yes","Alive",27.3 +"Yes","Dead",67.2 +"Yes","Alive",48.4 +"Yes","Alive",32.7 +"No","Alive",56 +"Yes","Dead",63.4 +"No","Alive",22.5 +"No","Alive",59.8 +"Yes","Alive",22.3 +"No","Alive",38 +"No","Alive",62.3 +"Yes","Alive",43.5 +"No","Alive",47.7 +"No","Alive",34.6 +"Yes","Alive",37 +"Yes","Alive",18.7 +"Yes","Alive",35.5 +"No","Dead",62.1 +"Yes","Alive",42.5 +"Yes","Dead",61.1 +"No","Alive",45.7 +"Yes","Alive",35 +"Yes","Alive",30.3 +"Yes","Alive",27.3 +"No","Alive",43.1 +"No","Alive",20.5 +"Yes","Dead",59.7 +"Yes","Alive",35.9 +"Yes","Dead",36.9 +"No","Alive",18.8 +"No","Dead",66.4 +"Yes","Alive",27.5 +"No","Dead",67.7 +"Yes","Alive",43.6 +"No","Alive",62.2 +"No","Dead",86 +"No","Dead",85.8 +"No","Alive",29.8 +"Yes","Alive",28.7 +"Yes","Alive",61.4 +"No","Alive",73.2 +"No","Alive",57.6 +"No","Alive",29.5 +"Yes","Dead",56.5 +"Yes","Alive",19.9 +"No","Alive",55.1 +"Yes","Dead",58.9 +"No","Alive",32.6 +"No","Dead",82.5 +"No","Alive",38.4 +"Yes","Alive",47.3 +"No","Dead",76.8 +"Yes","Alive",28.8 +"No","Alive",26.4 +"No","Dead",67.5 +"Yes","Alive",40.5 +"Yes","Alive",50.8 +"No","Alive",25.8 +"Yes","Alive",52.1 +"Yes","Alive",50.2 +"No","Alive",38.5 +"No","Dead",56.1 +"Yes","Alive",29.8 +"No","Alive",58.1 +"No","Dead",87.6 +"No","Dead",69.6 +"Yes","Alive",33.5 +"No","Dead",86 +"Yes","Dead",22.6 +"No","Alive",53.3 +"No","Alive",21.9 +"No","Alive",26 +"Yes","Alive",49.4 +"Yes","Alive",18 +"No","Alive",44.1 +"No","Dead",79.1 +"Yes","Alive",38.5 +"No","Alive",26.3 +"Yes","Alive",49.4 +"No","Alive",33.9 +"No","Dead",63.2 +"No","Alive",77.1 +"No","Dead",63.8 +"Yes","Alive",37.1 +"Yes","Alive",29.6 +"Yes","Alive",22.2 +"No","Dead",83.4 +"No","Dead",61.8 +"No","Alive",32.7 +"No","Alive",32.8 +"No","Alive",24.7 +"No","Alive",60.6 +"Yes","Alive",38.2 +"No","Dead",43.7 +"Yes","Alive",29.9 +"No","Alive",50.8 +"No","Alive",19.5 +"No","Alive",55.4 +"Yes","Alive",55.1 +"No","Dead",41.9 +"Yes","Alive",30.7 +"No","Alive",27.9 +"No","Alive",26 +"No","Dead",77.8 +"No","Alive",48.3 +"Yes","Dead",42.4 +"No","Alive",20.5 +"No","Dead",77.5 +"No","Alive",58.5 +"No","Alive",44.7 +"No","Alive",28.5 +"Yes","Alive",21 +"No","Alive",50.5 +"Yes","Alive",71.5 +"No","Alive",37.8 +"No","Alive",23.1 +"Yes","Dead",44.9 +"Yes","Alive",55.6 +"Yes","Alive",20.2 +"No","Alive",48.1 +"Yes","Alive",42.8 +"Yes","Alive",44.5 +"No","Dead",85.8 +"No","Alive",44.7 +"No","Alive",87.6 +"No","Alive",27 +"No","Alive",52 +"Yes","Alive",58.7 +"No","Alive",34.3 +"Yes","Alive",19.9 +"No","Alive",19.7 +"Yes","Alive",33.8 +"No","Alive",53 +"Yes","Alive",20.7 +"Yes","Alive",59.4 +"Yes","Alive",44.4 +"Yes","Alive",20.4 +"No","Dead",69 +"Yes","Alive",60.5 +"No","Alive",44 +"Yes","Alive",33.1 +"No","Alive",42.2 +"No","Alive",50.5 +"No","Alive",30.5 +"Yes","Alive",26.6 +"Yes","Alive",21 +"Yes","Alive",36.6 +"Yes","Alive",28.9 +"No","Alive",47.8 +"No","Dead",73.3 +"No","Alive",49.6 +"No","Alive",44.8 +"Yes","Alive",38.6 +"No","Dead",79.9 +"Yes","Dead",84.4 +"No","Dead",39.1 +"Yes","Alive",47.4 +"No","Alive",57.8 +"No","Alive",41.5 +"No","Alive",20.3 +"Yes","Alive",38.1 +"Yes","Alive",44.6 +"Yes","Alive",39.3 +"Yes","Alive",18.1 +"No","Alive",51.5 +"No","Alive",23.1 +"No","Alive",22.7 +"Yes","Alive",36.8 +"No","Alive",57.4 +"Yes","Alive",57.1 +"No","Alive",19.2 +"No","Dead",84.8 +"No","Alive",26.9 +"No","Dead",88.4 +"No","Dead",77.4 +"No","Dead",41.3 +"No","Alive",53.4 +"Yes","Alive",58.9 +"Yes","Dead",38.8 +"No","Dead",82.2 +"No","Alive",46.9 +"Yes","Alive",24.6 +"Yes","Alive",30.4 +"No","Alive",42.4 +"No","Dead",64 +"No","Alive",33.3 +"Yes","Alive",60.2 +"Yes","Alive",25 +"Yes","Dead",37.1 +"Yes","Alive",47.7 +"No","Dead",66.5 +"Yes","Dead",43.3 +"No","Alive",19.1 +"No","Alive",52.4 +"No","Alive",33.9 +"No","Alive",40 +"No","Alive",29.9 +"Yes","Alive",58.4 +"Yes","Alive",48.7 +"Yes","Alive",52.3 +"No","Dead",59.9 +"No","Alive",63.5 +"Yes","Alive",48.3 +"Yes","Alive",51.1 +"Yes","Dead",34.5 +"Yes","Alive",37.5 +"Yes","Alive",73.8 +"Yes","Alive",24.6 +"No","Dead",65.3 +"No","Alive",34.2 +"No","Alive",71.8 +"No","Dead",47.5 +"No","Alive",31.3 +"Yes","Dead",28.3 +"Yes","Dead",61.9 +"Yes","Dead",74.8 +"Yes","Alive",51 +"Yes","Dead",42.5 +"Yes","Alive",38.1 +"No","Alive",47.4 +"No","Alive",32.1 +"No","Dead",86.2 +"Yes","Alive",55.7 +"Yes","Alive",43.6 +"Yes","Dead",58.2 +"No","Alive",25.1 +"No","Dead",75 +"No","Alive",23.2 +"No","Alive",20.6 +"No","Dead",66.4 +"Yes","Alive",44 +"Yes","Alive",19.4 +"Yes","Dead",61 +"Yes","Alive",29.5 +"No","Alive",22.5 +"No","Alive",60 +"Yes","Alive",84.7 +"No","Dead",85.1 +"Yes","Dead",82 +"Yes","Alive",33.4 +"No","Alive",21.3 +"No","Dead",65.2 +"Yes","Dead",83.6 +"Yes","Alive",52.4 +"Yes","Alive",38.9 +"No","Alive",32.9 +"Yes","Alive",53.6 +"No","Alive",35.7 +"Yes","Alive",19.8 +"No","Dead",65.7 +"No","Alive",40.6 +"Yes","Alive",25.7 +"Yes","Dead",44.3 +"No","Alive",68.4 +"No","Alive",33.7 +"No","Alive",26.5 +"Yes","Alive",43.6 +"Yes","Dead",32.6 +"No","Alive",21 +"No","Alive",25.2 +"No","Dead",81.4 +"No","Dead",50.2 +"No","Alive",85 +"No","Alive",45.9 +"Yes","Dead",56.5 +"No","Alive",33.4 +"Yes","Alive",61.1 +"No","Alive",35.1 +"No","Alive",22.7 +"Yes","Alive",40.4 +"No","Alive",48.6 +"No","Dead",82.9 +"No","Dead",78.1 +"No","Alive",36.9 +"Yes","Alive",61.8 +"Yes","Alive",29.5 +"No","Alive",38.9 +"No","Dead",71.3 +"No","Alive",36.5 +"No","Dead",81.8 +"No","Alive",21.7 +"Yes","Dead",78.3 +"Yes","Alive",30.5 +"No","Alive",61.5 +"Yes","Alive",33.1 +"No","Alive",32.2 +"No","Alive",48.5 +"No","Alive",20.3 +"No","Alive",62.6 +"No","Alive",28.5 +"Yes","Alive",52.4 +"Yes","Dead",55.7 +"No","Alive",53.8 +"No","Alive",20.7 +"Yes","Alive",33.4 +"No","Alive",43.8 +"Yes","Alive",53.1 +"Yes","Alive",51.5 +"No","Alive",31.8 +"No","Dead",73.2 +"Yes","Alive",41.1 +"No","Dead",82 +"Yes","Alive",27 +"Yes","Alive",44.3 +"Yes","Dead",42.9 +"Yes","Dead",56.1 +"No","Dead",60.2 +"Yes","Alive",55.8 +"Yes","Alive",29.1 +"No","Alive",49.4 +"Yes","Dead",44.9 +"No","Dead",80.9 +"No","Alive",25.8 +"No","Alive",31.5 +"No","Dead",82.6 +"No","Alive",27.3 +"No","Alive",18.8 +"No","Alive",33.2 +"No","Alive",29.7 +"Yes","Dead",52.6 +"No","Dead",81.1 +"Yes","Dead",88.6 +"No","Alive",35 +"No","Dead",75.2 +"Yes","Alive",37.3 +"Yes","Alive",52.1 +"No","Dead",84.7 +"No","Dead",85 +"No","Alive",27 +"No","Dead",85 +"No","Alive",20.2 +"No","Alive",46.3 +"Yes","Alive",60 +"No","Dead",63.5 +"Yes","Dead",84.3 +"No","Alive",66.4 +"Yes","Alive",30.2 +"Yes","Alive",23.1 +"No","Alive",61.5 +"No","Alive",40.7 +"Yes","Alive",27.1 +"Yes","Alive",36.7 +"No","Alive",58.2 +"Yes","Alive",29.7 +"No","Alive",48.9 +"No","Alive",52.9 +"No","Alive",41.7 +"No","Alive",23 +"No","Alive",18.3 +"No","Dead",89.9 +"No","Alive",60.6 +"No","Alive",30.1 +"Yes","Alive",41.9 +"Yes","Alive",47 +"No","Alive",23.8 +"Yes","Dead",31.3 +"Yes","Dead",63.3 +"No","Alive",52.4 +"No","Alive",65 +"No","Dead",74.8 +"No","Alive",32.9 +"Yes","Dead",49.6 +"No","Alive",59.9 +"No","Alive",30.8 +"No","Alive",30.1 +"No","Alive",52 +"Yes","Alive",57.2 +"No","Dead",89.5 +"Yes","Alive",32.5 +"No","Alive",19.1 +"Yes","Alive",44 +"Yes","Dead",39.2 +"No","Alive",22.9 +"Yes","Alive",18 +"No","Alive",20.1 +"Yes","Alive",28 +"No","Alive",53 +"Yes","Alive",46.7 +"No","Alive",44.6 +"No","Alive",18.7 +"No","Dead",71.1 +"Yes","Alive",42.3 +"No","Alive",64 +"Yes","Dead",71 +"Yes","Alive",26.6 +"Yes","Alive",50.8 +"No","Alive",25.5 +"Yes","Alive",24 +"No","Alive",48.1 +"Yes","Alive",50.6 +"Yes","Alive",21.5 +"No","Alive",61.2 +"No","Dead",75.9 +"No","Dead",88 +"No","Dead",66.8 +"No","Alive",50.8 +"No","Alive",34.9 +"No","Dead",83.8 +"No","Alive",25 +"Yes","Dead",41.7 +"No","Alive",42.3 +"No","Alive",62.4 +"Yes","Alive",38.1 +"Yes","Alive",23.3 +"Yes","Alive",25.6 +"No","Dead",51.1 +"Yes","Alive",21.2 +"No","Dead",56.9 +"No","Alive",35 +"Yes","Dead",45 +"Yes","Alive",25.2 +"Yes","Alive",43.7 +"No","Dead",86.7 +"No","Alive",20.2 +"No","Dead",71.6 +"No","Dead",78.3 +"No","Alive",23.1 +"No","Dead",84.8 +"Yes","Alive",58.1 +"Yes","Alive",53.9 +"No","Alive",53.3 +"No","Alive",30.9 +"Yes","Alive",60.6 +"Yes","Dead",85.2 +"No","Alive",57.5 +"No","Alive",46.5 +"No","Dead",73.8 +"No","Alive",62.6 +"No","Alive",43.5 +"No","Alive",52.5 +"Yes","Alive",34.1 +"No","Alive",38.7 +"No","Alive",22.6 +"No","Alive",20 +"No","Alive",59.9 +"No","Dead",83.3 +"Yes","Alive",52.2 +"No","Dead",76.2 +"Yes","Alive",28 +"Yes","Alive",56.6 +"No","Dead",67.8 +"No","Alive",21.2 +"No","Alive",27.9 +"Yes","Alive",29.8 +"Yes","Alive",28.1 +"Yes","Alive",53.2 +"No","Alive",23.2 +"No","Alive",39.5 +"Yes","Alive",31.4 +"Yes","Alive",30 +"Yes","Alive",37.8 +"Yes","Alive",46.9 +"Yes","Alive",43.8 +"Yes","Alive",63.1 +"No","Alive",21.4 +"No","Dead",62.5 +"No","Alive",45.5 +"Yes","Alive",27.9 +"Yes","Alive",29.5 +"Yes","Alive",61 +"Yes","Alive",27 +"Yes","Alive",61.5 +"Yes","Dead",56.2 +"Yes","Dead",87.9 +"Yes","Alive",28.3 +"No","Dead",75.1 +"No","Dead",87.9 +"Yes","Alive",31 +"Yes","Alive",55.3 +"No","Alive",40.8 +"Yes","Alive",46.2 +"No","Alive",52.3 +"Yes","Alive",51.9 +"No","Alive",28.3 +"Yes","Alive",44.4 +"Yes","Dead",63.3 +"Yes","Alive",41 +"Yes","Alive",50.2 +"No","Alive",55.4 +"No","Dead",43.3 +"No","Alive",60.1 +"Yes","Alive",29.7 +"No","Dead",79 +"No","Dead",65.1 +"Yes","Alive",40.1 +"No","Alive",46 +"No","Alive",40.2 +"No","Dead",89.2 +"No","Alive",26 +"No","Alive",43.4 +"No","Alive",48.8 +"No","Alive",19.8 +"Yes","Alive",27.8 +"Yes","Alive",52.4 +"Yes","Alive",27.8 +"Yes","Alive",41 +"No","Dead",28.5 +"No","Alive",26.7 +"No","Alive",36 +"No","Dead",74.4 +"Yes","Alive",40.8 +"Yes","Alive",20.4 +"No","Dead",42.1 +"No","Alive",41.2 +"Yes","Alive",20.9 +"Yes","Alive",45.5 +"No","Alive",26.7 +"No","Alive",41.8 +"No","Alive",33.7 +"No","Alive",56.5 +"Yes","Alive",38.8 +"Yes","Alive",55.5 +"Yes","Alive",24.9 +"No","Alive",33 +"Yes","Alive",55.7 +"No","Alive",25.7 +"No","Alive",19.5 +"Yes","Alive",58.5 +"No","Alive",23.4 +"Yes","Alive",43.7 +"No","Alive",34.4 +"No","Dead",83.9 +"No","Alive",34.9 +"Yes","Alive",51.2 +"No","Dead",86.3 +"Yes","Dead",36 +"Yes","Alive",48.3 +"No","Alive",63.1 +"No","Alive",60.8 +"Yes","Dead",39.3 +"No","Alive",36.7 +"No","Alive",63.8 +"No","Dead",71.3 +"No","Alive",57.7 +"No","Alive",63.2 +"No","Alive",46.6 +"Yes","Dead",82.4 +"Yes","Alive",38.3 +"Yes","Alive",32.7 +"No","Alive",39.7 +"Yes","Dead",60 +"No","Dead",71 +"No","Alive",20.5 +"No","Alive",44.4 +"Yes","Alive",31.2 +"Yes","Alive",47.8 +"Yes","Alive",60.9 +"No","Dead",61.4 +"Yes","Alive",43 +"No","Alive",42.1 +"Yes","Alive",35.9 +"No","Alive",22.3 +"Yes","Dead",62.1 +"No","Dead",88.6 +"No","Alive",39.1 diff --git a/module3/exo3/exercice.ipynb b/module3/exo3/exercice.ipynb deleted file mode 100644 index 0bbbe371b01e359e381e43239412d77bf53fb1fb..0000000000000000000000000000000000000000 --- a/module3/exo3/exercice.ipynb +++ /dev/null @@ -1,25 +0,0 @@ -{ - "cells": [], - "metadata": { - "kernelspec": { - "display_name": "Python 3", - "language": "python", - "name": "python3" - }, - "language_info": { - "codemirror_mode": { - "name": "ipython", - "version": 3 - }, - "file_extension": ".py", - "mimetype": "text/x-python", - "name": "python", - "nbconvert_exporter": "python", - "pygments_lexer": "ipython3", - "version": "3.6.3" - } - }, - "nbformat": 4, - "nbformat_minor": 2 -} - diff --git a/module3/exo3/monthly_in_situ_co2_mlo.csv b/module3/exo3/monthly_in_situ_co2_mlo.csv new file mode 100644 index 0000000000000000000000000000000000000000..4366337f818f4eb949ee7a3bd354c8184198d1ba --- /dev/null +++ b/module3/exo3/monthly_in_situ_co2_mlo.csv @@ -0,0 +1,825 @@ +"-------------------------------------------------------------------------------------------" +" Atmospheric CO2 concentrations (ppm) derived from in situ air measurements " +" at Mauna Loa, Observatory, Hawaii: Latitude 19.5°N Longitude 155.6°W Elevation 3397m " +" " +" Source: R. F. Keeling, S. J. Walker, S. C. Piper and A. F. Bollenbacher " +" Scripps CO2 Program ( http://scrippsco2.ucsd.edu ) " +" Scripps Institution of Oceanography (SIO) " +" University of California " +" La Jolla, California USA 92093-0244 " +" " +" Status of data and correspondence: " +" " +" These data are subject to revision based on recalibration of standard gases. Questions " +" about the data should be directed to Dr. Ralph Keeling (rkeeling@ucsd.edu), Stephen Walker" +" (sjwalker@ucsd.edu) and Stephen Piper (scpiper@ucsd.edu), Scripps CO2 Program. " +" " +" Baseline data in this file through 06-Jan-2021 from archive dated 06-Jan-2021 09:38:28 " +" " +"-------------------------------------------------------------------------------------------" +" " +" Please cite as: " +" " +" C. D. Keeling, S. C. Piper, R. B. Bacastow, M. Wahlen, T. P. Whorf, M. Heimann, and " +" H. A. Meijer, Exchanges of atmospheric CO2 and 13CO2 with the terrestrial biosphere and " +" oceans from 1978 to 2000. I. Global aspects, SIO Reference Series, No. 01-06, Scripps " +" Institution of Oceanography, San Diego, 88 pages, 2001. " +" " +" If it is necessary to cite a peer-reviewed article, please cite as: " +" " +" C. D. Keeling, S. C. Piper, R. B. Bacastow, M. Wahlen, T. P. Whorf, M. Heimann, and " +" H. A. Meijer, Atmospheric CO2 and 13CO2 exchange with the terrestrial biosphere and " +" oceans from 1978 to 2000: observations and carbon cycle implications, pages 83-113, " +" in "A History of Atmospheric CO2 and its effects on Plants, Animals, and Ecosystems", " +" editors, Ehleringer, J.R., T. E. Cerling, M. D. Dearing, Springer Verlag, " +" New York, 2005. " +" " +"-------------------------------------------------------------------------------------------" +" " +" The data file below contains 10 columns. Columns 1-4 give the dates in several redundant " +" formats. Column 5 below gives monthly Mauna Loa CO2 concentrations in micro-mol CO2 per " +" mole (ppm), reported on the 2008A SIO manometric mole fraction scale. This is the " +" standard version of the data most often sought. The monthly values have been adjusted " +" to 24:00 hours on the 15th of each month. Column 6 gives the same data after a seasonal " +" adjustment to remove the quasi-regular seasonal cycle. The adjustment involves " +" subtracting from the data a 4-harmonic fit with a linear gain factor. Column 7 is a " +" smoothed version of the data generated from a stiff cubic spline function plus 4-harmonic " +" functions with linear gain. Column 8 is the same smoothed version with the seasonal " +" cycle removed. Column 9 is identical to Column 5 except that the missing values from " +" Column 5 have been filled with values from Column 7. Column 10 is identical to Column 6 " +" except missing values have been filled with values from Column 8. Missing values are " +" denoted by -99.99 " +" " +" CO2 concentrations are measured on the '08A' calibration scale " +" " + Yr, Mn, Date, Date, CO2,seasonally, fit, seasonally, CO2, seasonally + , , , , , adjusted, ,adjusted fit, filled,adjusted filled + , , Excel, , [ppm], [ppm] , [ppm], [ppm], [ppm], [ppm] +1958, 01, 21200, 1958.0411, -99.99, -99.99, -99.99, -99.99, -99.99, -99.99 +1958, 02, 21231, 1958.1260, -99.99, -99.99, -99.99, -99.99, -99.99, -99.99 +1958, 03, 21259, 1958.2027, 315.70, 314.43, 316.19, 314.90, 315.70, 314.43 +1958, 04, 21290, 1958.2877, 317.46, 315.16, 317.30, 314.98, 317.46, 315.16 +1958, 05, 21320, 1958.3699, 317.51, 314.71, 317.86, 315.06, 317.51, 314.71 +1958, 06, 21351, 1958.4548, -99.99, -99.99, 317.24, 315.14, 317.24, 315.14 +1958, 07, 21381, 1958.5370, 315.86, 315.19, 315.86, 315.21, 315.86, 315.19 +1958, 08, 21412, 1958.6219, 314.93, 316.19, 313.99, 315.29, 314.93, 316.19 +1958, 09, 21443, 1958.7068, 313.21, 316.09, 312.45, 315.35, 313.21, 316.09 +1958, 10, 21473, 1958.7890, -99.99, -99.99, 312.43, 315.40, 312.43, 315.40 +1958, 11, 21504, 1958.8740, 313.33, 315.20, 313.61, 315.46, 313.33, 315.20 +1958, 12, 21534, 1958.9562, 314.67, 315.43, 314.76, 315.51, 314.67, 315.43 +1959, 01, 21565, 1959.0411, 315.58, 315.54, 315.62, 315.57, 315.58, 315.54 +1959, 02, 21596, 1959.1260, 316.49, 315.86, 316.27, 315.63, 316.49, 315.86 +1959, 03, 21624, 1959.2027, 316.65, 315.38, 316.98, 315.69, 316.65, 315.38 +1959, 04, 21655, 1959.2877, 317.72, 315.41, 318.09, 315.77, 317.72, 315.41 +1959, 05, 21685, 1959.3699, 318.29, 315.48, 318.66, 315.85, 318.29, 315.48 +1959, 06, 21716, 1959.4548, 318.15, 316.03, 318.04, 315.94, 318.15, 316.03 +1959, 07, 21746, 1959.5370, 316.54, 315.87, 316.67, 316.03, 316.54, 315.87 +1959, 08, 21777, 1959.6219, 314.80, 316.07, 314.82, 316.12, 314.80, 316.07 +1959, 09, 21808, 1959.7068, 313.84, 316.73, 313.31, 316.22, 313.84, 316.73 +1959, 10, 21838, 1959.7890, 313.33, 316.33, 313.32, 316.31, 313.33, 316.33 +1959, 11, 21869, 1959.8740, 314.81, 316.69, 314.54, 316.39, 314.81, 316.69 +1959, 12, 21899, 1959.9562, 315.58, 316.35, 315.72, 316.47, 315.58, 316.35 +1960, 01, 21930, 1960.0410, 316.43, 316.39, 316.61, 316.56, 316.43, 316.39 +1960, 02, 21961, 1960.1257, 316.98, 316.35, 317.27, 316.64, 316.98, 316.35 +1960, 03, 21990, 1960.2049, 317.58, 316.28, 318.03, 316.71, 317.58, 316.28 +1960, 04, 22021, 1960.2896, 319.03, 316.69, 319.15, 316.79, 319.03, 316.69 +1960, 05, 22051, 1960.3716, 320.04, 317.21, 319.68, 316.86, 320.04, 317.21 +1960, 06, 22082, 1960.4563, 319.58, 317.47, 319.02, 316.93, 319.58, 317.47 +1960, 07, 22112, 1960.5383, 318.18, 317.53, 317.59, 316.97, 318.18, 317.53 +1960, 08, 22143, 1960.6230, 315.90, 317.20, 315.68, 317.01, 315.90, 317.20 +1960, 09, 22174, 1960.7077, 314.17, 317.08, 314.12, 317.04, 314.17, 317.08 +1960, 10, 22204, 1960.7896, 313.83, 316.83, 314.09, 317.07, 313.83, 316.83 +1960, 11, 22235, 1960.8743, 315.00, 316.88, 315.25, 317.11, 315.00, 316.88 +1960, 12, 22265, 1960.9563, 316.19, 316.96, 316.40, 317.15, 316.19, 316.96 +1961, 01, 22296, 1961.0411, 316.89, 316.85, 317.25, 317.20, 316.89, 316.85 +1961, 02, 22327, 1961.1260, 317.70, 317.07, 317.91, 317.26, 317.70, 317.07 +1961, 03, 22355, 1961.2027, 318.54, 317.26, 318.63, 317.33, 318.54, 317.26 +1961, 04, 22386, 1961.2877, 319.48, 317.15, 319.75, 317.41, 319.48, 317.15 +1961, 05, 22416, 1961.3699, 320.58, 317.75, 320.33, 317.50, 320.58, 317.75 +1961, 06, 22447, 1961.4548, 319.77, 317.63, 319.71, 317.59, 319.77, 317.63 +1961, 07, 22477, 1961.5370, 318.56, 317.89, 318.32, 317.68, 318.56, 317.89 +1961, 08, 22508, 1961.6219, 316.79, 318.07, 316.45, 317.76, 316.79, 318.07 +1961, 09, 22539, 1961.7068, 314.99, 317.90, 314.92, 317.85, 314.99, 317.90 +1961, 10, 22569, 1961.7890, 315.31, 318.32, 314.92, 317.92, 315.31, 318.32 +1961, 11, 22600, 1961.8740, 316.10, 317.99, 316.12, 317.99, 316.10, 317.99 +1961, 12, 22630, 1961.9562, 317.01, 317.78, 317.30, 318.06, 317.01, 317.78 +1962, 01, 22661, 1962.0411, 317.94, 317.90, 318.18, 318.13, 317.94, 317.90 +1962, 02, 22692, 1962.1260, 318.55, 317.91, 318.84, 318.19, 318.55, 317.91 +1962, 03, 22720, 1962.2027, 319.68, 318.39, 319.56, 318.26, 319.68, 318.39 +1962, 04, 22751, 1962.2877, 320.57, 318.24, 320.67, 318.32, 320.57, 318.24 +1962, 05, 22781, 1962.3699, 321.02, 318.17, 321.23, 318.39, 321.02, 318.17 +1962, 06, 22812, 1962.4548, 320.62, 318.47, 320.58, 318.45, 320.62, 318.47 +1962, 07, 22842, 1962.5370, 319.61, 318.93, 319.15, 318.50, 319.61, 318.93 +1962, 08, 22873, 1962.6219, 317.40, 318.69, 317.24, 318.56, 317.40, 318.69 +1962, 09, 22904, 1962.7068, 316.25, 319.17, 315.66, 318.60, 316.25, 319.17 +1962, 10, 22934, 1962.7890, 315.42, 318.44, 315.63, 318.64, 315.42, 318.44 +1962, 11, 22965, 1962.8740, 316.69, 318.58, 316.80, 318.68, 316.69, 318.58 +1962, 12, 22995, 1962.9562, 317.70, 318.47, 317.96, 318.72, 317.70, 318.47 +1963, 01, 23026, 1963.0411, 318.74, 318.70, 318.82, 318.76, 318.74, 318.70 +1963, 02, 23057, 1963.1260, 319.07, 318.43, 319.45, 318.81, 319.07, 318.43 +1963, 03, 23085, 1963.2027, 319.86, 318.57, 320.16, 318.86, 319.86, 318.57 +1963, 04, 23116, 1963.2877, 321.38, 319.04, 321.26, 318.91, 321.38, 319.04 +1963, 05, 23146, 1963.3699, 322.25, 319.39, 321.81, 318.96, 322.25, 319.39 +1963, 06, 23177, 1963.4548, 321.48, 319.33, 321.15, 319.01, 321.48, 319.33 +1963, 07, 23207, 1963.5370, 319.74, 319.06, 319.71, 319.06, 319.74, 319.06 +1963, 08, 23238, 1963.6219, 317.77, 319.05, 317.78, 319.10, 317.77, 319.05 +1963, 09, 23269, 1963.7068, 316.21, 319.14, 316.20, 319.14, 316.21, 319.14 +1963, 10, 23299, 1963.7890, 315.99, 319.02, 316.16, 319.18, 315.99, 319.02 +1963, 11, 23330, 1963.8740, 317.07, 318.97, 317.35, 319.23, 317.07, 318.97 +1963, 12, 23360, 1963.9562, 318.35, 319.13, 318.51, 319.27, 318.35, 319.13 +1964, 01, 23391, 1964.0410, 319.57, 319.53, 319.37, 319.32, 319.57, 319.53 +1964, 02, 23422, 1964.1257, -99.99, -99.99, 320.01, 319.37, 320.01, 319.37 +1964, 03, 23451, 1964.2049, -99.99, -99.99, 320.74, 319.41, 320.74, 319.41 +1964, 04, 23482, 1964.2896, -99.99, -99.99, 321.84, 319.45, 321.84, 319.45 +1964, 05, 23512, 1964.3716, 322.26, 319.40, 322.35, 319.49, 322.26, 319.40 +1964, 06, 23543, 1964.4563, 321.89, 319.75, 321.65, 319.53, 321.89, 319.75 +1964, 07, 23573, 1964.5383, 320.44, 319.78, 320.18, 319.56, 320.44, 319.78 +1964, 08, 23604, 1964.6230, 318.69, 320.01, 318.23, 319.58, 318.69, 320.01 +1964, 09, 23635, 1964.7077, 316.71, 319.66, 316.64, 319.60, 316.71, 319.66 +1964, 10, 23665, 1964.7896, 316.87, 319.91, 316.59, 319.62, 316.87, 319.91 +1964, 11, 23696, 1964.8743, 317.68, 319.58, 317.75, 319.64, 317.68, 319.58 +1964, 12, 23726, 1964.9563, 318.71, 319.49, 318.89, 319.66, 318.71, 319.49 +1965, 01, 23757, 1965.0411, 319.44, 319.40, 319.74, 319.68, 319.44, 319.40 +1965, 02, 23788, 1965.1260, 320.44, 319.80, 320.37, 319.72, 320.44, 319.80 +1965, 03, 23816, 1965.2027, 320.89, 319.59, 321.08, 319.77, 320.89, 319.59 +1965, 04, 23847, 1965.2877, 322.14, 319.78, 322.20, 319.83, 322.14, 319.78 +1965, 05, 23877, 1965.3699, 322.17, 319.30, 322.78, 319.91, 322.17, 319.30 +1965, 06, 23908, 1965.4548, 321.87, 319.70, 322.15, 320.00, 321.87, 319.70 +1965, 07, 23938, 1965.5370, 321.21, 320.52, 320.76, 320.10, 321.21, 320.52 +1965, 08, 23969, 1965.6219, 318.87, 320.16, 318.88, 320.21, 318.87, 320.16 +1965, 09, 24000, 1965.7068, 317.82, 320.77, 317.36, 320.32, 317.82, 320.77 +1965, 10, 24030, 1965.7890, 317.30, 320.36, 317.39, 320.43, 317.30, 320.36 +1965, 11, 24061, 1965.8740, 318.87, 320.78, 318.66, 320.55, 318.87, 320.78 +1965, 12, 24091, 1965.9562, 319.42, 320.20, 319.90, 320.66, 319.42, 320.20 +1966, 01, 24122, 1966.0411, 320.62, 320.58, 320.84, 320.78, 320.62, 320.58 +1966, 02, 24153, 1966.1260, 321.60, 320.95, 321.56, 320.90, 321.60, 320.95 +1966, 03, 24181, 1966.2027, 322.39, 321.08, 322.33, 321.01, 322.39, 321.08 +1966, 04, 24212, 1966.2877, 323.70, 321.34, 323.50, 321.12, 323.70, 321.34 +1966, 05, 24242, 1966.3699, 324.08, 321.20, 324.10, 321.22, 324.08, 321.20 +1966, 06, 24273, 1966.4548, 323.75, 321.57, 323.48, 321.32, 323.75, 321.57 +1966, 07, 24303, 1966.5370, 322.38, 321.69, 322.06, 321.40, 322.38, 321.69 +1966, 08, 24334, 1966.6219, 320.36, 321.66, 320.15, 321.48, 320.36, 321.66 +1966, 09, 24365, 1966.7068, 318.64, 321.60, 318.58, 321.55, 318.64, 321.60 +1966, 10, 24395, 1966.7890, 318.10, 321.17, 318.56, 321.62, 318.10, 321.17 +1966, 11, 24426, 1966.8740, 319.78, 321.70, 319.78, 321.68, 319.78, 321.70 +1966, 12, 24456, 1966.9562, 321.03, 321.81, 320.97, 321.74, 321.03, 321.81 +1967, 01, 24487, 1967.0411, 322.33, 322.29, 321.86, 321.80, 322.33, 322.29 +1967, 02, 24518, 1967.1260, 322.50, 321.85, 322.51, 321.86, 322.50, 321.85 +1967, 03, 24546, 1967.2027, 323.04, 321.73, 323.23, 321.90, 323.04, 321.73 +1967, 04, 24577, 1967.2877, 324.42, 322.04, 324.35, 321.96, 324.42, 322.04 +1967, 05, 24607, 1967.3699, 325.00, 322.11, 324.91, 322.01, 325.00, 322.11 +1967, 06, 24638, 1967.4548, 324.09, 321.91, 324.24, 322.07, 324.09, 321.91 +1967, 07, 24668, 1967.5370, 322.54, 321.85, 322.80, 322.14, 322.54, 321.85 +1967, 08, 24699, 1967.6219, 320.92, 322.22, 320.87, 322.20, 320.92, 322.22 +1967, 09, 24730, 1967.7068, 319.25, 322.23, 319.29, 322.27, 319.25, 322.23 +1967, 10, 24760, 1967.7890, 319.39, 322.47, 319.28, 322.35, 319.39, 322.47 +1967, 11, 24791, 1967.8740, 320.73, 322.65, 320.51, 322.42, 320.73, 322.65 +1967, 12, 24821, 1967.9562, 321.96, 322.75, 321.72, 322.49, 321.96, 322.75 +1968, 01, 24852, 1968.0410, 322.57, 322.53, 322.62, 322.57, 322.57, 322.53 +1968, 02, 24883, 1968.1257, 323.15, 322.50, 323.31, 322.65, 323.15, 322.50 +1968, 03, 24912, 1968.2049, 323.89, 322.55, 324.08, 322.73, 323.89, 322.55 +1968, 04, 24943, 1968.2896, 325.02, 322.62, 325.24, 322.83, 325.02, 322.62 +1968, 05, 24973, 1968.3716, 325.57, 322.67, 325.82, 322.92, 325.57, 322.67 +1968, 06, 25004, 1968.4563, 325.36, 323.19, 325.18, 323.03, 325.36, 323.19 +1968, 07, 25034, 1968.5383, 324.14, 323.47, 323.77, 323.14, 324.14, 323.47 +1968, 08, 25065, 1968.6230, 322.11, 323.44, 321.88, 323.25, 322.11, 323.44 +1968, 09, 25096, 1968.7077, 320.33, 323.32, 320.36, 323.37, 320.33, 323.32 +1968, 10, 25126, 1968.7896, 320.25, 323.33, 320.42, 323.49, 320.25, 323.33 +1968, 11, 25157, 1968.8743, 321.32, 323.25, 321.71, 323.62, 321.32, 323.25 +1968, 12, 25187, 1968.9563, 322.89, 323.68, 322.97, 323.75, 322.89, 323.68 +1969, 01, 25218, 1969.0411, 324.00, 323.96, 323.94, 323.89, 324.00, 323.96 +1969, 02, 25249, 1969.1260, 324.41, 323.76, 324.69, 324.03, 324.41, 323.76 +1969, 03, 25277, 1969.2027, 325.63, 324.32, 325.49, 324.16, 325.63, 324.32 +1969, 04, 25308, 1969.2877, 326.66, 324.27, 326.71, 324.30, 326.66, 324.27 +1969, 05, 25338, 1969.3699, 327.38, 324.47, 327.35, 324.44, 327.38, 324.47 +1969, 06, 25369, 1969.4548, 326.71, 324.51, 326.75, 324.57, 326.71, 324.51 +1969, 07, 25399, 1969.5370, 325.88, 325.19, 325.36, 324.69, 325.88, 325.19 +1969, 08, 25430, 1969.6219, 323.66, 324.98, 323.46, 324.80, 323.66, 324.98 +1969, 09, 25461, 1969.7068, 322.38, 325.37, 321.90, 324.90, 322.38, 325.37 +1969, 10, 25491, 1969.7890, 321.78, 324.88, 321.91, 324.99, 321.78, 324.88 +1969, 11, 25522, 1969.8740, 322.85, 324.80, 323.16, 325.08, 322.85, 324.80 +1969, 12, 25552, 1969.9562, 324.12, 324.91, 324.38, 325.16, 324.12, 324.91 +1970, 01, 25583, 1970.0411, 325.06, 325.02, 325.30, 325.24, 325.06, 325.02 +1970, 02, 25614, 1970.1260, 325.98, 325.33, 325.98, 325.32, 325.98, 325.33 +1970, 03, 25642, 1970.2027, 326.93, 325.61, 326.73, 325.40, 326.93, 325.61 +1970, 04, 25673, 1970.2877, 328.13, 325.74, 327.89, 325.48, 328.13, 325.74 +1970, 05, 25703, 1970.3699, 328.08, 325.16, 328.47, 325.55, 328.08, 325.16 +1970, 06, 25734, 1970.4548, 327.67, 325.46, 327.82, 325.63, 327.67, 325.46 +1970, 07, 25764, 1970.5370, 326.34, 325.65, 326.37, 325.70, 326.34, 325.65 +1970, 08, 25795, 1970.6219, 324.69, 326.00, 324.42, 325.77, 324.69, 326.00 +1970, 09, 25826, 1970.7068, 323.10, 326.10, 322.82, 325.83, 323.10, 326.10 +1970, 10, 25856, 1970.7890, 323.06, 326.17, 322.79, 325.88, 323.06, 326.17 +1970, 11, 25887, 1970.8740, 324.01, 325.95, 324.00, 325.93, 324.01, 325.95 +1970, 12, 25917, 1970.9562, 325.13, 325.93, 325.19, 325.97, 325.13, 325.93 +1971, 01, 25948, 1971.0411, 326.17, 326.13, 326.06, 326.01, 326.17, 326.13 +1971, 02, 25979, 1971.1260, 326.68, 326.03, 326.71, 326.05, 326.68, 326.03 +1971, 03, 26007, 1971.2027, 327.17, 325.85, 327.43, 326.09, 327.17, 325.85 +1971, 04, 26038, 1971.2877, 327.79, 325.38, 328.56, 326.14, 327.79, 325.38 +1971, 05, 26068, 1971.3699, 328.93, 326.00, 329.12, 326.19, 328.93, 326.00 +1971, 06, 26099, 1971.4548, 328.57, 326.36, 328.46, 326.26, 328.57, 326.36 +1971, 07, 26129, 1971.5370, 327.36, 326.66, 327.00, 326.33, 327.36, 326.66 +1971, 08, 26160, 1971.6219, 325.43, 326.75, 325.05, 326.40, 325.43, 326.75 +1971, 09, 26191, 1971.7068, 323.36, 326.37, 323.45, 326.47, 323.36, 326.37 +1971, 10, 26221, 1971.7890, 323.56, 326.68, 323.44, 326.54, 323.56, 326.68 +1971, 11, 26252, 1971.8740, 324.80, 326.75, 324.69, 326.62, 324.80, 326.75 +1971, 12, 26282, 1971.9562, 326.01, 326.81, 325.92, 326.70, 326.01, 326.81 +1972, 01, 26313, 1972.0410, 326.77, 326.73, 326.84, 326.78, 326.77, 326.73 +1972, 02, 26344, 1972.1257, 327.63, 326.98, 327.54, 326.88, 327.63, 326.98 +1972, 03, 26373, 1972.2049, 327.75, 326.39, 328.35, 326.98, 327.75, 326.39 +1972, 04, 26404, 1972.2896, 329.72, 327.29, 329.55, 327.10, 329.72, 327.29 +1972, 05, 26434, 1972.3716, 330.07, 327.13, 330.16, 327.23, 330.07, 327.13 +1972, 06, 26465, 1972.4563, 329.09, 326.89, 329.56, 327.38, 329.09, 326.89 +1972, 07, 26495, 1972.5383, 328.04, 327.37, 328.18, 327.54, 328.04, 327.37 +1972, 08, 26526, 1972.6230, 326.32, 327.67, 326.33, 327.72, 326.32, 327.67 +1972, 09, 26557, 1972.7077, 324.84, 327.87, 324.86, 327.91, 324.84, 327.87 +1972, 10, 26587, 1972.7896, 325.20, 328.33, 324.99, 328.10, 325.20, 328.33 +1972, 11, 26618, 1972.8743, 326.50, 328.46, 326.37, 328.31, 326.50, 328.46 +1972, 12, 26648, 1972.9563, 327.55, 328.35, 327.72, 328.51, 327.55, 328.35 +1973, 01, 26679, 1973.0411, 328.55, 328.50, 328.77, 328.72, 328.55, 328.50 +1973, 02, 26710, 1973.1260, 329.56, 328.90, 329.59, 328.92, 329.56, 328.90 +1973, 03, 26738, 1973.2027, 330.30, 328.96, 330.46, 329.11, 330.30, 328.96 +1973, 04, 26769, 1973.2877, 331.50, 329.08, 331.74, 329.30, 331.50, 329.08 +1973, 05, 26799, 1973.3699, 332.48, 329.53, 332.43, 329.48, 332.48, 329.53 +1973, 06, 26830, 1973.4548, 332.07, 329.84, 331.85, 329.64, 332.07, 329.84 +1973, 07, 26860, 1973.5370, 330.87, 330.16, 330.44, 329.77, 330.87, 330.16 +1973, 08, 26891, 1973.6219, 329.31, 330.64, 328.51, 329.88, 329.31, 330.64 +1973, 09, 26922, 1973.7068, 327.51, 330.55, 326.91, 329.96, 327.51, 330.55 +1973, 10, 26952, 1973.7890, 327.18, 330.32, 326.88, 330.01, 327.18, 330.32 +1973, 11, 26983, 1973.8740, 328.16, 330.13, 328.10, 330.04, 328.16, 330.13 +1973, 12, 27013, 1973.9562, 328.64, 329.44, 329.28, 330.07, 328.64, 329.44 +1974, 01, 27044, 1974.0411, 329.35, 329.31, 330.15, 330.09, 329.35, 329.31 +1974, 02, 27075, 1974.1260, 330.71, 330.05, 330.79, 330.12, 330.71, 330.05 +1974, 03, 27103, 1974.2027, 331.48, 330.14, 331.50, 330.14, 331.48, 330.14 +1974, 04, 27134, 1974.2877, 332.65, 330.22, 332.62, 330.18, 332.65, 330.22 +1974, 05, 27164, 1974.3699, 333.09, 330.13, 333.18, 330.22, 333.09, 330.13 +1974, 06, 27195, 1974.4548, 332.25, 330.01, 332.48, 330.26, 332.25, 330.01 +1974, 07, 27225, 1974.5370, 331.18, 330.47, 330.99, 330.31, 331.18, 330.47 +1974, 08, 27256, 1974.6219, 329.39, 330.73, 328.99, 330.36, 329.39, 330.73 +1974, 09, 27287, 1974.7068, 327.43, 330.47, 327.35, 330.41, 327.43, 330.47 +1974, 10, 27317, 1974.7890, 327.37, 330.52, 327.32, 330.46, 327.37, 330.52 +1974, 11, 27348, 1974.8740, 328.46, 330.43, 328.57, 330.52, 328.46, 330.43 +1974, 12, 27378, 1974.9562, 329.57, 330.38, 329.79, 330.58, 329.57, 330.38 +1975, 01, 27409, 1975.0411, 330.40, 330.36, 330.71, 330.65, 330.40, 330.36 +1975, 02, 27440, 1975.1260, 331.40, 330.74, 331.41, 330.74, 331.40, 330.74 +1975, 03, 27468, 1975.2027, 332.04, 330.69, 332.17, 330.82, 332.04, 330.69 +1975, 04, 27499, 1975.2877, 333.31, 330.87, 333.36, 330.91, 333.31, 330.87 +1975, 05, 27529, 1975.3699, 333.97, 331.00, 333.97, 331.00, 333.97, 331.00 +1975, 06, 27560, 1975.4548, 333.60, 331.36, 333.33, 331.10, 333.60, 331.36 +1975, 07, 27590, 1975.5370, 331.90, 331.19, 331.88, 331.20, 331.90, 331.19 +1975, 08, 27621, 1975.6219, 330.06, 331.39, 329.92, 331.29, 330.06, 331.39 +1975, 09, 27652, 1975.7068, 328.56, 331.61, 328.32, 331.39, 328.56, 331.61 +1975, 10, 27682, 1975.7890, 328.34, 331.50, 328.32, 331.47, 328.34, 331.50 +1975, 11, 27713, 1975.8740, 329.49, 331.47, 329.60, 331.56, 329.49, 331.47 +1975, 12, 27743, 1975.9562, 330.76, 331.57, 330.84, 331.64, 330.76, 331.57 +1976, 01, 27774, 1976.0410, 331.75, 331.70, 331.77, 331.72, 331.75, 331.70 +1976, 02, 27805, 1976.1257, 332.57, 331.90, 332.47, 331.80, 332.57, 331.90 +1976, 03, 27834, 1976.2049, 333.50, 332.12, 333.26, 331.87, 333.50, 332.12 +1976, 04, 27865, 1976.2896, 334.58, 332.12, 334.43, 331.95, 334.58, 332.12 +1976, 05, 27895, 1976.3716, 334.88, 331.90, 335.00, 332.03, 334.88, 331.90 +1976, 06, 27926, 1976.4563, 334.33, 332.10, 334.32, 332.11, 334.33, 332.10 +1976, 07, 27956, 1976.5383, 333.05, 332.37, 332.85, 332.19, 333.05, 332.37 +1976, 08, 27987, 1976.6230, 330.94, 332.31, 330.88, 332.29, 330.94, 332.31 +1976, 09, 28018, 1976.7077, 329.30, 332.38, 329.30, 332.39, 329.30, 332.38 +1976, 10, 28048, 1976.7896, 328.94, 332.11, 329.35, 332.50, 328.94, 332.11 +1976, 11, 28079, 1976.8743, 330.31, 332.29, 330.67, 332.63, 330.31, 332.29 +1976, 12, 28109, 1976.9563, 331.68, 332.49, 331.98, 332.77, 331.68, 332.49 +1977, 01, 28140, 1977.0411, 332.93, 332.88, 332.99, 332.93, 332.93, 332.88 +1977, 02, 28171, 1977.1260, 333.42, 332.75, 333.78, 333.11, 333.42, 332.75 +1977, 03, 28199, 1977.2027, 334.70, 333.35, 334.64, 333.27, 334.70, 333.35 +1977, 04, 28230, 1977.2877, 336.07, 333.62, 335.93, 333.46, 336.07, 333.62 +1977, 05, 28260, 1977.3699, 336.75, 333.76, 336.62, 333.64, 336.75, 333.76 +1977, 06, 28291, 1977.4548, 336.27, 334.01, 336.06, 333.82, 336.27, 334.01 +1977, 07, 28321, 1977.5370, 334.92, 334.21, 334.67, 333.99, 334.92, 334.21 +1977, 08, 28352, 1977.6219, 332.75, 334.10, 332.77, 334.15, 332.75, 334.10 +1977, 09, 28383, 1977.7068, 331.59, 334.66, 331.22, 334.31, 331.59, 334.66 +1977, 10, 28413, 1977.7890, 331.16, 334.34, 331.29, 334.46, 331.16, 334.34 +1977, 11, 28444, 1977.8740, 332.40, 334.40, 332.63, 334.60, 332.40, 334.40 +1977, 12, 28474, 1977.9562, 333.85, 334.66, 333.94, 334.74, 333.85, 334.66 +1978, 01, 28505, 1978.0411, 334.97, 334.92, 334.93, 334.88, 334.97, 334.92 +1978, 02, 28536, 1978.1260, 335.39, 334.72, 335.69, 335.01, 335.39, 334.72 +1978, 03, 28564, 1978.2027, 336.64, 335.28, 336.50, 335.13, 336.64, 335.28 +1978, 04, 28595, 1978.2877, 337.76, 335.30, 337.73, 335.25, 337.76, 335.30 +1978, 05, 28625, 1978.3699, 338.01, 335.01, 338.37, 335.37, 338.01, 335.01 +1978, 06, 28656, 1978.4548, 337.90, 335.63, 337.74, 335.49, 337.90, 335.63 +1978, 07, 28686, 1978.5370, 336.54, 335.82, 336.28, 335.60, 336.54, 335.82 +1978, 08, 28717, 1978.6219, 334.68, 336.04, 334.31, 335.70, 334.68, 336.04 +1978, 09, 28748, 1978.7068, 332.76, 335.84, 332.70, 335.80, 332.76, 335.84 +1978, 10, 28778, 1978.7890, 332.55, 335.74, 332.71, 335.89, 332.55, 335.74 +1978, 11, 28809, 1978.8740, 333.92, 335.92, 334.01, 335.99, 333.92, 335.92 +1978, 12, 28839, 1978.9562, 334.95, 335.77, 335.28, 336.08, 334.95, 335.77 +1979, 01, 28870, 1979.0411, 336.23, 336.18, 336.25, 336.19, 336.23, 336.18 +1979, 02, 28901, 1979.1260, 336.76, 336.09, 336.99, 336.31, 336.76, 336.09 +1979, 03, 28929, 1979.2027, 337.96, 336.60, 337.79, 336.41, 337.96, 336.60 +1979, 04, 28960, 1979.2877, 338.88, 336.41, 339.03, 336.54, 338.88, 336.41 +1979, 05, 28990, 1979.3699, 339.47, 336.46, 339.68, 336.67, 339.47, 336.46 +1979, 06, 29021, 1979.4548, 339.29, 337.01, 339.06, 336.80, 339.29, 337.01 +1979, 07, 29051, 1979.5370, 337.73, 337.01, 337.62, 336.94, 337.73, 337.01 +1979, 08, 29082, 1979.6219, 336.09, 337.44, 335.68, 337.08, 336.09, 337.44 +1979, 09, 29113, 1979.7068, 333.92, 337.01, 334.11, 337.22, 333.92, 337.01 +1979, 10, 29143, 1979.7890, 333.86, 337.06, 334.18, 337.37, 333.86, 337.06 +1979, 11, 29174, 1979.8740, 335.29, 337.30, 335.54, 337.52, 335.29, 337.30 +1979, 12, 29204, 1979.9562, 336.73, 337.55, 336.88, 337.68, 336.73, 337.55 +1980, 01, 29235, 1980.0410, 338.01, 337.97, 337.91, 337.85, 338.01, 337.97 +1980, 02, 29266, 1980.1257, 338.36, 337.69, 338.70, 338.02, 338.36, 337.69 +1980, 03, 29295, 1980.2049, 340.07, 338.68, 339.58, 338.17, 340.07, 338.68 +1980, 04, 29326, 1980.2896, 340.76, 338.26, 340.85, 338.33, 340.76, 338.26 +1980, 05, 29356, 1980.3716, 341.47, 338.45, 341.50, 338.48, 341.47, 338.45 +1980, 06, 29387, 1980.4563, 341.17, 338.92, 340.87, 338.63, 341.17, 338.92 +1980, 07, 29417, 1980.5383, 339.56, 338.87, 339.43, 338.77, 339.56, 338.87 +1980, 08, 29448, 1980.6230, 337.60, 338.99, 337.47, 338.90, 337.60, 338.99 +1980, 09, 29479, 1980.7077, 335.88, 338.99, 335.89, 339.02, 335.88, 338.99 +1980, 10, 29509, 1980.7896, 336.02, 339.23, 335.93, 339.13, 336.02, 339.23 +1980, 11, 29540, 1980.8743, 337.10, 339.11, 337.25, 339.24, 337.10, 339.11 +1980, 12, 29570, 1980.9563, 338.21, 339.03, 338.54, 339.34, 338.21, 339.03 +1981, 01, 29601, 1981.0411, 339.24, 339.19, 339.50, 339.45, 339.24, 339.19 +1981, 02, 29632, 1981.1260, 340.48, 339.80, 340.23, 339.55, 340.48, 339.80 +1981, 03, 29660, 1981.2027, 341.38, 340.01, 341.02, 339.63, 341.38, 340.01 +1981, 04, 29691, 1981.2877, 342.50, 340.02, 342.23, 339.73, 342.50, 340.02 +1981, 05, 29721, 1981.3699, 342.91, 339.88, 342.84, 339.81, 342.91, 339.88 +1981, 06, 29752, 1981.4548, 342.25, 339.96, 342.16, 339.89, 342.25, 339.96 +1981, 07, 29782, 1981.5370, 340.49, 339.76, 340.67, 339.98, 340.49, 339.76 +1981, 08, 29813, 1981.6219, 338.43, 339.80, 338.67, 340.07, 338.43, 339.80 +1981, 09, 29844, 1981.7068, 336.69, 339.81, 337.03, 340.16, 336.69, 339.81 +1981, 10, 29874, 1981.7890, 336.86, 340.08, 337.05, 340.26, 336.86, 340.08 +1981, 11, 29905, 1981.8740, 338.36, 340.38, 338.38, 340.37, 338.36, 340.38 +1981, 12, 29935, 1981.9562, 339.61, 340.44, 339.67, 340.48, 339.61, 340.44 +1982, 01, 29966, 1982.0411, 340.75, 340.71, 340.65, 340.59, 340.75, 340.71 +1982, 02, 29997, 1982.1260, 341.61, 340.94, 341.39, 340.70, 341.61, 340.94 +1982, 03, 30025, 1982.2027, 342.70, 341.32, 342.18, 340.79, 342.70, 341.32 +1982, 04, 30056, 1982.2877, 343.57, 341.07, 343.39, 340.88, 343.57, 341.07 +1982, 05, 30086, 1982.3699, 344.14, 341.10, 344.00, 340.96, 344.14, 341.10 +1982, 06, 30117, 1982.4548, 343.35, 341.05, 343.32, 341.04, 343.35, 341.05 +1982, 07, 30147, 1982.5370, 342.06, 341.33, 341.81, 341.12, 342.06, 341.33 +1982, 08, 30178, 1982.6219, 339.81, 341.18, 339.79, 341.20, 339.81, 341.18 +1982, 09, 30209, 1982.7068, 337.98, 341.10, 338.14, 341.28, 337.98, 341.10 +1982, 10, 30239, 1982.7890, 337.86, 341.10, 338.16, 341.38, 337.86, 341.10 +1982, 11, 30270, 1982.8740, 339.26, 341.29, 339.49, 341.49, 339.26, 341.29 +1982, 12, 30300, 1982.9562, 340.49, 341.32, 340.81, 341.62, 340.49, 341.32 +1983, 01, 30331, 1983.0411, 341.38, 341.33, 341.84, 341.78, 341.38, 341.33 +1983, 02, 30362, 1983.1260, 342.52, 341.84, 342.64, 341.95, 342.52, 341.84 +1983, 03, 30390, 1983.2027, 343.10, 341.72, 343.51, 342.12, 343.10, 341.72 +1983, 04, 30421, 1983.2877, 344.94, 342.44, 344.83, 342.31, 344.94, 342.44 +1983, 05, 30451, 1983.3699, 345.76, 342.71, 345.54, 342.50, 345.76, 342.71 +1983, 06, 30482, 1983.4548, 345.32, 343.01, 344.97, 342.69, 345.32, 343.01 +1983, 07, 30512, 1983.5370, 343.98, 343.25, 343.56, 342.86, 343.98, 343.25 +1983, 08, 30543, 1983.6219, 342.38, 343.76, 341.62, 343.03, 342.38, 343.76 +1983, 09, 30574, 1983.7068, 339.87, 343.00, 340.04, 343.19, 339.87, 343.00 +1983, 10, 30604, 1983.7890, 339.99, 343.24, 340.11, 343.34, 339.99, 343.24 +1983, 11, 30635, 1983.8740, 341.15, 343.19, 341.47, 343.48, 341.15, 343.19 +1983, 12, 30665, 1983.9562, 342.99, 343.82, 342.81, 343.62, 342.99, 343.82 +1984, 01, 30696, 1984.0410, 343.70, 343.65, 343.82, 343.76, 343.70, 343.65 +1984, 02, 30727, 1984.1257, 344.50, 343.82, 344.58, 343.89, 344.50, 343.82 +1984, 03, 30756, 1984.2049, 345.28, 343.87, 345.44, 344.01, 345.28, 343.87 +1984, 04, 30787, 1984.2896, 347.05, 344.52, 346.69, 344.14, 347.05, 344.52 +1984, 05, 30817, 1984.3716, 347.43, 344.38, 347.32, 344.27, 347.43, 344.38 +1984, 06, 30848, 1984.4563, 346.80, 344.51, 346.65, 344.39, 346.80, 344.51 +1984, 07, 30878, 1984.5383, 345.39, 344.69, 345.17, 344.50, 345.39, 344.69 +1984, 08, 30909, 1984.6230, 343.28, 344.69, 343.18, 344.62, 343.28, 344.69 +1984, 09, 30940, 1984.7077, 341.07, 344.23, 341.57, 344.74, 341.07, 344.23 +1984, 10, 30970, 1984.7896, 341.35, 344.60, 341.63, 344.86, 341.35, 344.60 +1984, 11, 31001, 1984.8743, 342.98, 345.01, 342.98, 344.99, 342.98, 345.01 +1984, 12, 31031, 1984.9563, 344.22, 345.05, 344.30, 345.12, 344.22, 345.05 +1985, 01, 31062, 1985.0411, 344.97, 344.92, 345.31, 345.25, 344.97, 344.92 +1985, 02, 31093, 1985.1260, 345.99, 345.31, 346.08, 345.39, 345.99, 345.31 +1985, 03, 31121, 1985.2027, 347.42, 346.04, 346.91, 345.50, 347.42, 346.04 +1985, 04, 31152, 1985.2877, 348.35, 345.83, 348.16, 345.62, 348.35, 345.83 +1985, 05, 31182, 1985.3699, 348.93, 345.86, 348.80, 345.73, 348.93, 345.86 +1985, 06, 31213, 1985.4548, 348.25, 345.94, 348.13, 345.83, 348.25, 345.94 +1985, 07, 31243, 1985.5370, 346.56, 345.83, 346.63, 345.93, 346.56, 345.83 +1985, 08, 31274, 1985.6219, 344.67, 346.06, 344.60, 346.02, 344.67, 346.06 +1985, 09, 31305, 1985.7068, 343.09, 346.24, 342.94, 346.11, 343.09, 346.24 +1985, 10, 31335, 1985.7890, 342.80, 346.07, 342.94, 346.19, 342.80, 346.07 +1985, 11, 31366, 1985.8740, 344.24, 346.29, 344.26, 346.28, 344.24, 346.29 +1985, 12, 31396, 1985.9562, 345.56, 346.39, 345.55, 346.37, 345.56, 346.39 +1986, 01, 31427, 1986.0411, 346.30, 346.25, 346.53, 346.47, 346.30, 346.25 +1986, 02, 31458, 1986.1260, 346.95, 346.26, 347.28, 346.58, 346.95, 346.26 +1986, 03, 31486, 1986.2027, 347.85, 346.46, 348.10, 346.69, 347.85, 346.46 +1986, 04, 31517, 1986.2877, 349.55, 347.03, 349.35, 346.81, 349.55, 347.03 +1986, 05, 31547, 1986.3699, 350.22, 347.14, 350.01, 346.93, 350.22, 347.14 +1986, 06, 31578, 1986.4548, 349.55, 347.23, 349.37, 347.06, 349.55, 347.23 +1986, 07, 31608, 1986.5370, 347.94, 347.20, 347.89, 347.19, 347.94, 347.20 +1986, 08, 31639, 1986.6219, 345.90, 347.29, 345.89, 347.31, 345.90, 347.29 +1986, 09, 31670, 1986.7068, 344.85, 348.01, 344.26, 347.44, 344.85, 348.01 +1986, 10, 31700, 1986.7890, 344.17, 347.45, 344.30, 347.56, 344.17, 347.45 +1986, 11, 31731, 1986.8740, 345.66, 347.71, 345.67, 347.70, 345.66, 347.71 +1986, 12, 31761, 1986.9562, 346.90, 347.74, 347.01, 347.83, 346.90, 347.74 +1987, 01, 31792, 1987.0411, 348.02, 347.98, 348.03, 347.98, 348.02, 347.98 +1987, 02, 31823, 1987.1260, 348.48, 347.79, 348.83, 348.13, 348.48, 347.79 +1987, 03, 31851, 1987.2027, 349.42, 348.02, 349.70, 348.29, 349.42, 348.02 +1987, 04, 31882, 1987.2877, 350.98, 348.45, 351.02, 348.47, 350.98, 348.45 +1987, 05, 31912, 1987.3699, 351.85, 348.76, 351.75, 348.66, 351.85, 348.76 +1987, 06, 31943, 1987.4548, 351.26, 348.92, 351.17, 348.86, 351.26, 348.92 +1987, 07, 31973, 1987.5370, 349.51, 348.77, 349.77, 349.06, 349.51, 348.77 +1987, 08, 32004, 1987.6219, 348.10, 349.49, 347.85, 349.28, 348.10, 349.49 +1987, 09, 32035, 1987.7068, 346.45, 349.62, 346.31, 349.50, 346.45, 349.62 +1987, 10, 32065, 1987.7890, 346.36, 349.65, 346.44, 349.72, 346.36, 349.65 +1987, 11, 32096, 1987.8740, 347.81, 349.87, 347.91, 349.94, 347.81, 349.87 +1987, 12, 32126, 1987.9562, 348.96, 349.81, 349.34, 350.16, 348.96, 349.81 +1988, 01, 32157, 1988.0410, 350.43, 350.39, 350.45, 350.39, 350.43, 350.39 +1988, 02, 32188, 1988.1257, 351.73, 351.04, 351.31, 350.61, 351.73, 351.04 +1988, 03, 32217, 1988.2049, 352.22, 350.79, 352.26, 350.81, 352.22, 350.79 +1988, 04, 32248, 1988.2896, 353.59, 351.02, 353.60, 351.02, 353.59, 351.02 +1988, 05, 32278, 1988.3716, 354.22, 351.12, 354.31, 351.21, 354.22, 351.12 +1988, 06, 32309, 1988.4563, 353.80, 351.48, 353.70, 351.40, 353.80, 351.48 +1988, 07, 32339, 1988.5383, 352.38, 351.67, 352.25, 351.57, 352.38, 351.67 +1988, 08, 32370, 1988.6230, 350.43, 351.85, 350.28, 351.74, 350.43, 351.85 +1988, 09, 32401, 1988.7077, 348.73, 351.92, 348.68, 351.89, 348.73, 351.92 +1988, 10, 32431, 1988.7896, 348.88, 352.18, 348.75, 352.03, 348.88, 352.18 +1988, 11, 32462, 1988.8743, 350.07, 352.13, 350.12, 352.16, 350.07, 352.13 +1988, 12, 32492, 1988.9563, 351.34, 352.18, 351.45, 352.28, 351.34, 352.18 +1989, 01, 32523, 1989.0411, 352.76, 352.71, 352.45, 352.39, 352.76, 352.71 +1989, 02, 32554, 1989.1260, 353.07, 352.38, 353.20, 352.50, 353.07, 352.38 +1989, 03, 32582, 1989.2027, 353.68, 352.27, 354.02, 352.59, 353.68, 352.27 +1989, 04, 32613, 1989.2877, 355.42, 352.87, 355.26, 352.70, 355.42, 352.87 +1989, 05, 32643, 1989.3699, 355.67, 352.56, 355.90, 352.79, 355.67, 352.56 +1989, 06, 32674, 1989.4548, 355.12, 352.77, 355.22, 352.89, 355.12, 352.77 +1989, 07, 32704, 1989.5370, 353.90, 353.15, 353.70, 352.99, 353.90, 353.15 +1989, 08, 32735, 1989.6219, 351.67, 353.07, 351.65, 353.08, 351.67, 353.07 +1989, 09, 32766, 1989.7068, 349.81, 353.00, 349.97, 353.18, 349.81, 353.00 +1989, 10, 32796, 1989.7890, 349.99, 353.29, 349.98, 353.28, 349.99, 353.29 +1989, 11, 32827, 1989.8740, 351.30, 353.37, 351.33, 353.38, 351.30, 353.37 +1989, 12, 32857, 1989.9562, 352.52, 353.37, 352.65, 353.48, 352.52, 353.37 +1990, 01, 32888, 1990.0411, 353.66, 353.62, 353.64, 353.58, 353.66, 353.62 +1990, 02, 32919, 1990.1260, 354.70, 354.00, 354.39, 353.68, 354.70, 354.00 +1990, 03, 32947, 1990.2027, 355.38, 353.97, 355.20, 353.77, 355.38, 353.97 +1990, 04, 32978, 1990.2877, 356.20, 353.64, 356.45, 353.88, 356.20, 353.64 +1990, 05, 33008, 1990.3699, 357.16, 354.04, 357.10, 353.98, 357.16, 354.04 +1990, 06, 33039, 1990.4548, 356.23, 353.87, 356.43, 354.10, 356.23, 353.87 +1990, 07, 33069, 1990.5370, 354.81, 354.07, 354.93, 354.22, 354.81, 354.07 +1990, 08, 33100, 1990.6219, 352.91, 354.32, 352.90, 354.35, 352.91, 354.32 +1990, 09, 33131, 1990.7068, 350.96, 354.17, 351.27, 354.49, 350.96, 354.17 +1990, 10, 33161, 1990.7890, 351.18, 354.50, 351.32, 354.63, 351.18, 354.50 +1990, 11, 33192, 1990.8740, 352.83, 354.91, 352.72, 354.78, 352.83, 354.91 +1990, 12, 33222, 1990.9562, 354.21, 355.06, 354.09, 354.93, 354.21, 355.06 +1991, 01, 33253, 1991.0411, 354.72, 354.68, 355.13, 355.07, 354.72, 354.68 +1991, 02, 33284, 1991.1260, 355.75, 355.05, 355.92, 355.22, 355.75, 355.05 +1991, 03, 33312, 1991.2027, 357.16, 355.74, 356.76, 355.33, 357.16, 355.74 +1991, 04, 33343, 1991.2877, 358.60, 356.03, 358.03, 355.45, 358.60, 356.03 +1991, 05, 33373, 1991.3699, 359.34, 356.21, 358.67, 355.54, 359.34, 356.21 +1991, 06, 33404, 1991.4548, 358.24, 355.88, 357.96, 355.62, 358.24, 355.88 +1991, 07, 33434, 1991.5370, 356.17, 355.42, 356.40, 355.69, 356.17, 355.42 +1991, 08, 33465, 1991.6219, 354.01, 355.43, 354.30, 355.74, 354.01, 355.43 +1991, 09, 33496, 1991.7068, 352.15, 355.36, 352.57, 355.80, 352.15, 355.36 +1991, 10, 33526, 1991.7890, 352.21, 355.54, 352.54, 355.86, 352.21, 355.54 +1991, 11, 33557, 1991.8740, 353.75, 355.83, 353.86, 355.93, 353.75, 355.83 +1991, 12, 33587, 1991.9562, 354.99, 355.84, 355.16, 355.99, 354.99, 355.84 +1992, 01, 33618, 1992.0410, 355.99, 355.94, 356.12, 356.06, 355.99, 355.94 +1992, 02, 33649, 1992.1257, 356.72, 356.02, 356.84, 356.13, 356.72, 356.02 +1992, 03, 33678, 1992.2049, 357.81, 356.36, 357.66, 356.20, 357.81, 356.36 +1992, 04, 33709, 1992.2896, 359.15, 356.55, 358.88, 356.26, 359.15, 356.55 +1992, 05, 33739, 1992.3716, 359.66, 356.53, 359.45, 356.31, 359.66, 356.53 +1992, 06, 33770, 1992.4563, 359.25, 356.91, 358.69, 356.36, 359.25, 356.91 +1992, 07, 33800, 1992.5383, 357.02, 356.30, 357.08, 356.40, 357.02, 356.30 +1992, 08, 33831, 1992.6230, 355.00, 356.44, 354.95, 356.43, 355.00, 356.44 +1992, 09, 33862, 1992.7077, 353.01, 356.24, 353.21, 356.46, 353.01, 356.24 +1992, 10, 33892, 1992.7896, 353.31, 356.64, 353.17, 356.49, 353.31, 356.64 +1992, 11, 33923, 1992.8743, 354.16, 356.25, 354.47, 356.53, 354.16, 356.25 +1992, 12, 33953, 1992.9563, 355.40, 356.25, 355.73, 356.57, 355.40, 356.25 +1993, 01, 33984, 1993.0411, 356.70, 356.66, 356.68, 356.62, 356.70, 356.66 +1993, 02, 34015, 1993.1260, 357.17, 356.46, 357.40, 356.68, 357.17, 356.46 +1993, 03, 34043, 1993.2027, 358.38, 356.95, 358.19, 356.75, 358.38, 356.95 +1993, 04, 34074, 1993.2877, 359.46, 356.88, 359.42, 356.82, 359.46, 356.88 +1993, 05, 34104, 1993.3699, 360.28, 357.13, 360.05, 356.90, 360.28, 357.13 +1993, 06, 34135, 1993.4548, 359.60, 357.22, 359.35, 356.99, 359.60, 357.22 +1993, 07, 34165, 1993.5370, 357.57, 356.82, 357.81, 357.09, 357.57, 356.82 +1993, 08, 34196, 1993.6219, 355.52, 356.94, 355.75, 357.20, 355.52, 356.94 +1993, 09, 34227, 1993.7068, 353.69, 356.93, 354.08, 357.33, 353.69, 356.93 +1993, 10, 34257, 1993.7890, 353.99, 357.34, 354.13, 357.47, 353.99, 357.34 +1993, 11, 34288, 1993.8740, 355.34, 357.43, 355.54, 357.62, 355.34, 357.43 +1993, 12, 34318, 1993.9562, 356.80, 357.66, 356.93, 357.77, 356.80, 357.66 +1994, 01, 34349, 1994.0411, 358.37, 358.32, 358.00, 357.94, 358.37, 358.32 +1994, 02, 34380, 1994.1260, 358.91, 358.21, 358.83, 358.11, 358.91, 358.21 +1994, 03, 34408, 1994.2027, 359.97, 358.54, 359.71, 358.26, 359.97, 358.54 +1994, 04, 34439, 1994.2877, 361.26, 358.67, 361.04, 358.43, 361.26, 358.67 +1994, 05, 34469, 1994.3699, 361.69, 358.53, 361.75, 358.59, 361.69, 358.53 +1994, 06, 34500, 1994.4548, 360.94, 358.56, 361.12, 358.76, 360.94, 358.56 +1994, 07, 34530, 1994.5370, 359.55, 358.79, 359.64, 358.92, 359.55, 358.79 +1994, 08, 34561, 1994.6219, 357.48, 358.90, 357.63, 359.09, 357.48, 358.90 +1994, 09, 34592, 1994.7068, 355.84, 359.09, 356.01, 359.27, 355.84, 359.09 +1994, 10, 34622, 1994.7890, 356.00, 359.36, 356.10, 359.45, 356.00, 359.36 +1994, 11, 34653, 1994.8740, 357.58, 359.69, 357.55, 359.63, 357.58, 359.69 +1994, 12, 34683, 1994.9562, 359.04, 359.90, 358.97, 359.81, 359.04, 359.90 +1995, 01, 34714, 1995.0411, 359.97, 359.92, 360.05, 359.99, 359.97, 359.92 +1995, 02, 34745, 1995.1260, 361.00, 360.30, 360.89, 360.17, 361.00, 360.30 +1995, 03, 34773, 1995.2027, 361.63, 360.20, 361.78, 360.33, 361.63, 360.20 +1995, 04, 34804, 1995.2877, 363.45, 360.85, 363.12, 360.50, 363.45, 360.85 +1995, 05, 34834, 1995.3699, 363.80, 360.63, 363.83, 360.67, 363.80, 360.63 +1995, 06, 34865, 1995.4548, 363.26, 360.87, 363.20, 360.83, 363.26, 360.87 +1995, 07, 34895, 1995.5370, 361.89, 361.14, 361.71, 360.99, 361.89, 361.14 +1995, 08, 34926, 1995.6219, 359.45, 360.88, 359.69, 361.15, 359.45, 360.88 +1995, 09, 34957, 1995.7068, 358.05, 361.31, 358.04, 361.31, 358.05, 361.31 +1995, 10, 34987, 1995.7890, 357.76, 361.13, 358.11, 361.47, 357.76, 361.13 +1995, 11, 35018, 1995.8740, 359.56, 361.68, 359.54, 361.63, 359.56, 361.68 +1995, 12, 35048, 1995.9562, 360.70, 361.56, 360.94, 361.78, 360.70, 361.56 +1996, 01, 35079, 1996.0410, 362.05, 362.00, 362.00, 361.94, 362.05, 362.00 +1996, 02, 35110, 1996.1257, 363.24, 362.54, 362.80, 362.08, 363.24, 362.54 +1996, 03, 35139, 1996.2049, 364.02, 362.56, 363.69, 362.21, 364.02, 362.56 +1996, 04, 35170, 1996.2896, 364.71, 362.08, 364.98, 362.34, 364.71, 362.08 +1996, 05, 35200, 1996.3716, 365.42, 362.24, 365.62, 362.45, 365.42, 362.24 +1996, 06, 35231, 1996.4563, 364.97, 362.60, 364.91, 362.56, 364.97, 362.60 +1996, 07, 35261, 1996.5383, 363.65, 362.92, 363.35, 362.66, 363.65, 362.92 +1996, 08, 35292, 1996.6230, 361.48, 362.94, 361.25, 362.75, 361.48, 362.94 +1996, 09, 35323, 1996.7077, 359.45, 362.72, 359.54, 362.83, 359.45, 362.72 +1996, 10, 35353, 1996.7896, 359.61, 362.98, 359.55, 362.91, 359.61, 362.98 +1996, 11, 35384, 1996.8743, 360.76, 362.87, 360.90, 362.99, 360.76, 362.87 +1996, 12, 35414, 1996.9563, 362.33, 363.19, 362.22, 363.07, 362.33, 363.19 +1997, 01, 35445, 1997.0411, 363.19, 363.14, 363.21, 363.15, 363.19, 363.14 +1997, 02, 35476, 1997.1260, 363.99, 363.28, 363.96, 363.24, 363.99, 363.28 +1997, 03, 35504, 1997.2027, 364.56, 363.12, 364.78, 363.32, 364.56, 363.12 +1997, 04, 35535, 1997.2877, 366.35, 363.74, 366.06, 363.43, 366.35, 363.74 +1997, 05, 35565, 1997.3699, 366.80, 363.61, 366.73, 363.54, 366.80, 363.61 +1997, 06, 35596, 1997.4548, 365.63, 363.22, 366.06, 363.68, 365.63, 363.22 +1997, 07, 35626, 1997.5370, 364.47, 363.71, 364.55, 363.83, 364.47, 363.71 +1997, 08, 35657, 1997.6219, 362.50, 363.94, 362.53, 364.00, 362.50, 363.94 +1997, 09, 35688, 1997.7068, 360.19, 363.46, 360.91, 364.20, 360.19, 363.46 +1997, 10, 35718, 1997.7890, 360.78, 364.17, 361.04, 364.42, 360.78, 364.17 +1997, 11, 35749, 1997.8740, 362.43, 364.56, 362.56, 364.66, 362.43, 364.56 +1997, 12, 35779, 1997.9562, 364.28, 365.14, 364.06, 364.91, 364.28, 365.14 +1998, 01, 35810, 1998.0411, 365.33, 365.28, 365.24, 365.18, 365.33, 365.28 +1998, 02, 35841, 1998.1260, 366.15, 365.44, 366.18, 365.46, 366.15, 365.44 +1998, 03, 35869, 1998.2027, 367.31, 365.87, 367.17, 365.71, 367.31, 365.87 +1998, 04, 35900, 1998.2877, 368.61, 365.99, 368.62, 365.98, 368.61, 365.99 +1998, 05, 35930, 1998.3699, 369.30, 366.11, 369.44, 366.24, 369.30, 366.11 +1998, 06, 35961, 1998.4548, 368.88, 366.46, 368.90, 366.51, 368.88, 366.46 +1998, 07, 35991, 1998.5370, 367.64, 366.87, 367.48, 366.75, 367.64, 366.87 +1998, 08, 36022, 1998.6219, 365.78, 367.22, 365.51, 366.99, 365.78, 367.22 +1998, 09, 36053, 1998.7068, 363.90, 367.19, 363.90, 367.20, 363.90, 367.19 +1998, 10, 36083, 1998.7890, 364.23, 367.64, 364.00, 367.39, 364.23, 367.64 +1998, 11, 36114, 1998.8740, 365.46, 367.59, 365.46, 367.57, 365.46, 367.59 +1998, 12, 36144, 1998.9562, 366.97, 367.84, 366.87, 367.72, 366.97, 367.84 +1999, 01, 36175, 1999.0411, 368.15, 368.10, 367.91, 367.85, 368.15, 368.10 +1999, 02, 36206, 1999.1260, 368.87, 368.16, 368.69, 367.97, 368.87, 368.16 +1999, 03, 36234, 1999.2027, 369.59, 368.14, 369.52, 368.05, 369.59, 368.14 +1999, 04, 36265, 1999.2877, 371.14, 368.51, 370.78, 368.13, 371.14, 368.51 +1999, 05, 36295, 1999.3699, 371.00, 367.80, 371.41, 368.20, 371.00, 367.80 +1999, 06, 36326, 1999.4548, 370.35, 367.93, 370.67, 368.27, 370.35, 367.93 +1999, 07, 36356, 1999.5370, 369.27, 368.50, 369.07, 368.34, 369.27, 368.50 +1999, 08, 36387, 1999.6219, 366.93, 368.37, 366.92, 368.41, 366.93, 368.37 +1999, 09, 36418, 1999.7068, 364.64, 367.93, 365.17, 368.48, 364.64, 367.93 +1999, 10, 36448, 1999.7890, 365.13, 368.54, 365.16, 368.56, 365.13, 368.54 +1999, 11, 36479, 1999.8740, 366.68, 368.81, 366.53, 368.65, 366.68, 368.81 +1999, 12, 36509, 1999.9562, 368.00, 368.88, 367.87, 368.73, 368.00, 368.88 +2000, 01, 36540, 2000.0410, 369.14, 369.09, 368.88, 368.82, 369.14, 369.09 +2000, 02, 36571, 2000.1257, 369.46, 368.75, 369.64, 368.92, 369.46, 368.75 +2000, 03, 36600, 2000.2049, 370.51, 369.03, 370.51, 369.01, 370.51, 369.03 +2000, 04, 36631, 2000.2896, 371.66, 369.00, 371.80, 369.12, 371.66, 369.00 +2000, 05, 36661, 2000.3716, 371.83, 368.61, 372.45, 369.24, 371.83, 368.61 +2000, 06, 36692, 2000.4563, 371.69, 369.29, 371.75, 369.37, 371.69, 369.29 +2000, 07, 36722, 2000.5383, 370.12, 369.38, 370.21, 369.50, 370.12, 369.38 +2000, 08, 36753, 2000.6230, 368.12, 369.60, 368.13, 369.65, 368.12, 369.60 +2000, 09, 36784, 2000.7077, 366.62, 369.94, 366.46, 369.79, 366.62, 369.94 +2000, 10, 36814, 2000.7896, 366.73, 370.15, 366.53, 369.93, 366.73, 370.15 +2000, 11, 36845, 2000.8743, 368.29, 370.43, 367.96, 370.07, 368.29, 370.43 +2000, 12, 36875, 2000.9563, 369.52, 370.40, 369.35, 370.20, 369.52, 370.40 +2001, 01, 36906, 2001.0411, 370.28, 370.23, 370.39, 370.33, 370.28, 370.23 +2001, 02, 36937, 2001.1260, 371.50, 370.78, 371.18, 370.45, 371.50, 370.78 +2001, 03, 36965, 2001.2027, 372.12, 370.66, 372.03, 370.56, 372.12, 370.66 +2001, 04, 36996, 2001.2877, 372.86, 370.21, 373.34, 370.68, 372.86, 370.21 +2001, 05, 37026, 2001.3699, 374.02, 370.79, 374.02, 370.80, 374.02, 370.79 +2001, 06, 37057, 2001.4548, 373.31, 370.87, 373.34, 370.93, 373.31, 370.87 +2001, 07, 37087, 2001.5370, 371.62, 370.85, 371.80, 371.06, 371.62, 370.85 +2001, 08, 37118, 2001.6219, 369.55, 371.01, 369.71, 371.20, 369.55, 371.01 +2001, 09, 37149, 2001.7068, 367.96, 371.27, 368.02, 371.36, 367.96, 371.27 +2001, 10, 37179, 2001.7890, 368.10, 371.53, 368.09, 371.51, 368.10, 371.53 +2001, 11, 37210, 2001.8740, 369.68, 371.83, 369.55, 371.67, 369.68, 371.83 +2001, 12, 37240, 2001.9562, 371.24, 372.12, 370.97, 371.84, 371.24, 372.12 +2002, 01, 37271, 2002.0411, 372.44, 372.39, 372.07, 372.01, 372.44, 372.39 +2002, 02, 37302, 2002.1260, 373.08, 372.36, 372.91, 372.18, 373.08, 372.36 +2002, 03, 37330, 2002.2027, 373.52, 372.05, 373.82, 372.34, 373.52, 372.05 +2002, 04, 37361, 2002.2877, 374.85, 372.20, 375.20, 372.52, 374.85, 372.20 +2002, 05, 37391, 2002.3699, 375.55, 372.31, 375.95, 372.72, 375.55, 372.31 +2002, 06, 37422, 2002.4548, 375.40, 372.96, 375.35, 372.93, 375.40, 372.96 +2002, 07, 37452, 2002.5370, 374.02, 373.25, 373.88, 373.14, 374.02, 373.25 +2002, 08, 37483, 2002.6219, 371.48, 372.94, 371.87, 373.37, 371.48, 372.94 +2002, 09, 37514, 2002.7068, 370.70, 374.03, 370.26, 373.60, 370.70, 374.03 +2002, 10, 37544, 2002.7890, 370.25, 373.69, 370.40, 373.83, 370.25, 373.69 +2002, 11, 37575, 2002.8740, 372.08, 374.24, 371.93, 374.06, 372.08, 374.24 +2002, 12, 37605, 2002.9562, 373.78, 374.66, 373.42, 374.28, 373.78, 374.66 +2003, 01, 37636, 2003.0411, 374.68, 374.63, 374.57, 374.51, 374.68, 374.63 +2003, 02, 37667, 2003.1260, 375.62, 374.90, 375.46, 374.73, 375.62, 374.90 +2003, 03, 37695, 2003.2027, 376.11, 374.64, 376.41, 374.92, 376.11, 374.64 +2003, 04, 37726, 2003.2877, 377.65, 374.99, 377.81, 375.13, 377.65, 374.99 +2003, 05, 37756, 2003.3699, 378.35, 375.11, 378.58, 375.33, 378.35, 375.11 +2003, 06, 37787, 2003.4548, 378.13, 375.68, 377.97, 375.54, 378.13, 375.68 +2003, 07, 37817, 2003.5370, 376.60, 375.83, 376.47, 375.73, 376.60, 375.83 +2003, 08, 37848, 2003.6219, 374.48, 375.95, 374.42, 375.92, 374.48, 375.95 +2003, 09, 37879, 2003.7068, 372.98, 376.31, 372.75, 376.10, 372.98, 376.31 +2003, 10, 37909, 2003.7890, 373.00, 376.45, 372.83, 376.27, 373.00, 376.45 +2003, 11, 37940, 2003.8740, 374.35, 376.51, 374.29, 376.43, 374.35, 376.51 +2003, 12, 37970, 2003.9562, 375.69, 376.57, 375.70, 376.57, 375.69, 376.57 +2004, 01, 38001, 2004.0410, 376.79, 376.74, 376.77, 376.71, 376.79, 376.74 +2004, 02, 38032, 2004.1257, 377.37, 376.64, 377.58, 376.84, 377.37, 376.64 +2004, 03, 38061, 2004.2049, 378.39, 376.89, 378.48, 376.96, 378.39, 376.89 +2004, 04, 38092, 2004.2896, 380.50, 377.80, 379.80, 377.08, 380.50, 377.80 +2004, 05, 38122, 2004.3716, 380.62, 377.36, 380.45, 377.19, 380.62, 377.36 +2004, 06, 38153, 2004.4563, 379.55, 377.12, 379.72, 377.31, 379.55, 377.12 +2004, 07, 38183, 2004.5383, 377.76, 377.02, 378.15, 377.43, 377.76, 377.02 +2004, 08, 38214, 2004.6230, 375.83, 377.33, 376.03, 377.57, 375.83, 377.33 +2004, 09, 38245, 2004.7077, 374.05, 377.40, 374.34, 377.72, 374.05, 377.40 +2004, 10, 38275, 2004.7896, 374.22, 377.68, 374.43, 377.88, 374.22, 377.68 +2004, 11, 38306, 2004.8743, 375.84, 378.01, 375.92, 378.06, 375.84, 378.01 +2004, 12, 38336, 2004.9563, 377.44, 378.33, 377.38, 378.25, 377.44, 378.33 +2005, 01, 38367, 2005.0411, 378.34, 378.29, 378.51, 378.45, 378.34, 378.29 +2005, 02, 38398, 2005.1260, 379.61, 378.88, 379.40, 378.66, 379.61, 378.88 +2005, 03, 38426, 2005.2027, 380.08, 378.61, 380.35, 378.86, 380.08, 378.61 +2005, 04, 38457, 2005.2877, 382.05, 379.37, 381.78, 379.08, 382.05, 379.37 +2005, 05, 38487, 2005.3699, 382.24, 378.97, 382.56, 379.29, 382.24, 378.97 +2005, 06, 38518, 2005.4548, 382.08, 379.61, 381.96, 379.51, 382.08, 379.61 +2005, 07, 38548, 2005.5370, 380.66, 379.88, 380.47, 379.72, 380.66, 379.88 +2005, 08, 38579, 2005.6219, 378.67, 380.14, 378.43, 379.94, 378.67, 380.14 +2005, 09, 38610, 2005.7068, 376.42, 379.77, 376.77, 380.15, 376.42, 379.77 +2005, 10, 38640, 2005.7890, 376.80, 380.27, 376.88, 380.35, 376.80, 380.27 +2005, 11, 38671, 2005.8740, 378.31, 380.48, 378.39, 380.55, 378.31, 380.48 +2005, 12, 38701, 2005.9562, 379.96, 380.85, 379.86, 380.73, 379.96, 380.85 +2006, 01, 38732, 2006.0411, 381.37, 381.32, 380.98, 380.91, 381.37, 381.32 +2006, 02, 38763, 2006.1260, 382.02, 381.29, 381.83, 381.09, 382.02, 381.29 +2006, 03, 38791, 2006.2027, 382.56, 381.08, 382.73, 381.23, 382.56, 381.08 +2006, 04, 38822, 2006.2877, 384.36, 381.68, 384.10, 381.39, 384.36, 381.68 +2006, 05, 38852, 2006.3699, 384.92, 381.65, 384.81, 381.54, 384.92, 381.65 +2006, 06, 38883, 2006.4548, 384.03, 381.55, 384.14, 381.68, 384.03, 381.55 +2006, 07, 38913, 2006.5370, 382.28, 381.50, 382.57, 381.83, 382.28, 381.50 +2006, 08, 38944, 2006.6219, 380.48, 381.95, 380.46, 381.97, 380.48, 381.95 +2006, 09, 38975, 2006.7068, 378.81, 382.17, 378.74, 382.12, 378.81, 382.17 +2006, 10, 39005, 2006.7890, 379.06, 382.54, 378.80, 382.27, 379.06, 382.54 +2006, 11, 39036, 2006.8740, 380.14, 382.32, 380.26, 382.42, 380.14, 382.32 +2006, 12, 39066, 2006.9562, 381.66, 382.55, 381.69, 382.56, 381.66, 382.55 +2007, 01, 39097, 2007.0411, 382.58, 382.53, 382.78, 382.72, 382.58, 382.53 +2007, 02, 39128, 2007.1260, 383.71, 382.98, 383.62, 382.87, 383.71, 382.98 +2007, 03, 39156, 2007.2027, 384.34, 382.85, 384.52, 383.01, 384.34, 382.85 +2007, 04, 39187, 2007.2877, 386.23, 383.53, 385.89, 383.17, 386.23, 383.53 +2007, 05, 39217, 2007.3699, 386.41, 383.13, 386.61, 383.33, 386.41, 383.13 +2007, 06, 39248, 2007.4548, 385.87, 383.39, 385.95, 383.49, 385.87, 383.39 +2007, 07, 39278, 2007.5370, 384.44, 383.66, 384.40, 383.65, 384.44, 383.66 +2007, 08, 39309, 2007.6219, 381.84, 383.33, 382.29, 383.81, 381.84, 383.33 +2007, 09, 39340, 2007.7068, 380.86, 384.24, 380.58, 383.98, 380.86, 384.24 +2007, 10, 39370, 2007.7890, 380.86, 384.36, 380.65, 384.14, 380.86, 384.36 +2007, 11, 39401, 2007.8740, 382.36, 384.55, 382.13, 384.29, 382.36, 384.55 +2007, 12, 39431, 2007.9562, 383.61, 384.51, 383.56, 384.44, 383.61, 384.51 +2008, 01, 39462, 2008.0410, 385.07, 385.02, 384.65, 384.59, 385.07, 385.02 +2008, 02, 39493, 2008.1257, 385.84, 385.11, 385.47, 384.73, 385.84, 385.11 +2008, 03, 39522, 2008.2049, 385.83, 384.31, 386.40, 384.87, 385.83, 384.31 +2008, 04, 39553, 2008.2896, 386.77, 384.04, 387.77, 385.02, 386.77, 384.04 +2008, 05, 39583, 2008.3716, 388.51, 385.22, 388.47, 385.18, 388.51, 385.22 +2008, 06, 39614, 2008.4563, 388.05, 385.59, 387.79, 385.35, 388.05, 385.59 +2008, 07, 39644, 2008.5383, 386.25, 385.50, 386.23, 385.51, 386.25, 385.50 +2008, 08, 39675, 2008.6230, 384.08, 385.60, 384.13, 385.68, 384.08, 385.60 +2008, 09, 39706, 2008.7077, 383.09, 386.48, 382.44, 385.85, 383.09, 386.48 +2008, 10, 39736, 2008.7896, 382.78, 386.28, 382.52, 386.01, 382.78, 386.28 +2008, 11, 39767, 2008.8743, 384.01, 386.20, 384.00, 386.16, 384.01, 386.20 +2008, 12, 39797, 2008.9563, 385.11, 386.00, 385.44, 386.31, 385.11, 386.00 +2009, 01, 39828, 2009.0411, 386.65, 386.60, 386.53, 386.47, 386.65, 386.60 +2009, 02, 39859, 2009.1260, 387.12, 386.39, 387.37, 386.63, 387.12, 386.39 +2009, 03, 39887, 2009.2027, 388.52, 387.02, 388.28, 386.77, 388.52, 387.02 +2009, 04, 39918, 2009.2877, 389.57, 386.86, 389.67, 386.94, 389.57, 386.86 +2009, 05, 39948, 2009.3699, 390.17, 386.86, 390.41, 387.10, 390.17, 386.86 +2009, 06, 39979, 2009.4548, 389.62, 387.12, 389.76, 387.29, 389.62, 387.12 +2009, 07, 40009, 2009.5370, 388.07, 387.28, 388.23, 387.47, 388.07, 387.28 +2009, 08, 40040, 2009.6219, 386.08, 387.57, 386.14, 387.67, 386.08, 387.57 +2009, 09, 40071, 2009.7068, 384.65, 388.05, 384.47, 387.88, 384.65, 388.05 +2009, 10, 40101, 2009.7890, 384.33, 387.85, 384.59, 388.09, 384.33, 387.85 +2009, 11, 40132, 2009.8740, 386.05, 388.25, 386.14, 388.32, 386.05, 388.25 +2009, 12, 40162, 2009.9562, 387.49, 388.39, 387.65, 388.54, 387.49, 388.39 +2010, 01, 40193, 2010.0411, 388.55, 388.50, 388.83, 388.76, 388.55, 388.50 +2010, 02, 40224, 2010.1260, 390.08, 389.34, 389.74, 388.99, 390.08, 389.34 +2010, 03, 40252, 2010.2027, 391.01, 389.51, 390.70, 389.18, 391.01, 389.51 +2010, 04, 40283, 2010.2877, 392.38, 389.66, 392.13, 389.39, 392.38, 389.66 +2010, 05, 40313, 2010.3699, 393.22, 389.90, 392.90, 389.58, 393.22, 389.90 +2010, 06, 40344, 2010.4548, 392.24, 389.73, 392.25, 389.77, 392.24, 389.73 +2010, 07, 40374, 2010.5370, 390.33, 389.54, 390.70, 389.94, 390.33, 389.54 +2010, 08, 40405, 2010.6219, 388.52, 390.01, 388.58, 390.11, 388.52, 390.01 +2010, 09, 40436, 2010.7068, 386.84, 390.25, 386.85, 390.28, 386.84, 390.25 +2010, 10, 40466, 2010.7890, 387.16, 390.69, 386.92, 390.43, 387.16, 390.69 +2010, 11, 40497, 2010.8740, 388.67, 390.88, 388.39, 390.58, 388.67, 390.88 +2010, 12, 40527, 2010.9562, 389.81, 390.71, 389.83, 390.71, 389.81, 390.71 +2011, 01, 40558, 2011.0411, 391.30, 391.25, 390.91, 390.85, 391.30, 391.25 +2011, 02, 40589, 2011.1260, 391.92, 391.18, 391.73, 390.98, 391.92, 391.18 +2011, 03, 40617, 2011.2027, 392.45, 390.95, 392.62, 391.09, 392.45, 390.95 +2011, 04, 40648, 2011.2877, 393.37, 390.64, 393.98, 391.23, 393.37, 390.64 +2011, 05, 40678, 2011.3699, 394.28, 390.96, 394.70, 391.37, 394.28, 390.96 +2011, 06, 40709, 2011.4548, 393.69, 391.18, 394.02, 391.53, 393.69, 391.18 +2011, 07, 40739, 2011.5370, 392.59, 391.80, 392.45, 391.69, 392.59, 391.80 +2011, 08, 40770, 2011.6219, 390.21, 391.71, 390.33, 391.87, 390.21, 391.71 +2011, 09, 40801, 2011.7068, 389.00, 392.42, 388.62, 392.05, 389.00, 392.42 +2011, 10, 40831, 2011.7890, 388.93, 392.47, 388.70, 392.23, 388.93, 392.47 +2011, 11, 40862, 2011.8740, 390.24, 392.46, 390.22, 392.41, 390.24, 392.46 +2011, 12, 40892, 2011.9562, 391.80, 392.71, 391.70, 392.59, 391.80, 392.71 +2012, 01, 40923, 2012.0410, 393.07, 393.02, 392.83, 392.77, 393.07, 393.02 +2012, 02, 40954, 2012.1257, 393.35, 392.61, 393.71, 392.96, 393.35, 392.61 +2012, 03, 40983, 2012.2049, 394.36, 392.82, 394.69, 393.14, 394.36, 392.82 +2012, 04, 41014, 2012.2896, 396.43, 393.67, 396.12, 393.34, 396.43, 393.67 +2012, 05, 41044, 2012.3716, 396.87, 393.54, 396.88, 393.54, 396.87, 393.54 +2012, 06, 41075, 2012.4563, 395.88, 393.38, 396.23, 393.76, 395.88, 393.38 +2012, 07, 41105, 2012.5383, 394.52, 393.76, 394.71, 393.98, 394.52, 393.76 +2012, 08, 41136, 2012.6230, 392.54, 394.07, 392.64, 394.22, 392.54, 394.07 +2012, 09, 41167, 2012.7077, 391.13, 394.56, 391.00, 394.46, 391.13, 394.56 +2012, 10, 41197, 2012.7896, 391.01, 394.55, 391.16, 394.69, 391.01, 394.55 +2012, 11, 41228, 2012.8743, 392.95, 395.17, 392.74, 394.94, 392.95, 395.17 +2012, 12, 41258, 2012.9563, 394.33, 395.24, 394.28, 395.17, 394.33, 395.24 +2013, 01, 41289, 2013.0411, 395.61, 395.55, 395.47, 395.40, 395.61, 395.55 +2013, 02, 41320, 2013.1260, 396.85, 396.10, 396.39, 395.63, 396.85, 396.10 +2013, 03, 41348, 2013.2027, 397.26, 395.75, 397.36, 395.83, 397.26, 395.75 +2013, 04, 41379, 2013.2877, 398.35, 395.61, 398.80, 396.04, 398.35, 395.61 +2013, 05, 41409, 2013.3699, 399.98, 396.63, 399.58, 396.24, 399.98, 396.63 +2013, 06, 41440, 2013.4548, 398.87, 396.34, 398.95, 396.44, 398.87, 396.34 +2013, 07, 41470, 2013.5370, 397.37, 396.57, 397.40, 396.64, 397.37, 396.57 +2013, 08, 41501, 2013.6219, 395.41, 396.92, 395.28, 396.83, 395.41, 396.92 +2013, 09, 41532, 2013.7068, 393.39, 396.83, 393.56, 397.02, 393.39, 396.83 +2013, 10, 41562, 2013.7890, 393.70, 397.25, 393.65, 397.20, 393.70, 397.25 +2013, 11, 41593, 2013.8740, 395.19, 397.42, 395.17, 397.38, 395.19, 397.42 +2013, 12, 41623, 2013.9562, 396.82, 397.73, 396.65, 397.55, 396.82, 397.73 +2014, 01, 41654, 2014.0411, 397.93, 397.87, 397.78, 397.72, 397.93, 397.87 +2014, 02, 41685, 2014.1260, 398.10, 397.35, 398.65, 397.89, 398.10, 397.35 +2014, 03, 41713, 2014.2027, 399.47, 397.95, 399.57, 398.04, 399.47, 397.95 +2014, 04, 41744, 2014.2877, 401.33, 398.57, 400.97, 398.20, 401.33, 398.57 +2014, 05, 41774, 2014.3699, 401.88, 398.53, 401.71, 398.36, 401.88, 398.53 +2014, 06, 41805, 2014.4548, 401.31, 398.77, 401.03, 398.52, 401.31, 398.77 +2014, 07, 41835, 2014.5370, 399.07, 398.27, 399.44, 398.68, 399.07, 398.27 +2014, 08, 41866, 2014.6219, 397.21, 398.72, 397.29, 398.84, 397.21, 398.72 +2014, 09, 41897, 2014.7068, 395.40, 398.85, 395.54, 399.01, 395.40, 398.85 +2014, 10, 41927, 2014.7890, 395.65, 399.22, 395.62, 399.17, 395.65, 399.22 +2014, 11, 41958, 2014.8740, 397.23, 399.46, 397.14, 399.35, 397.23, 399.46 +2014, 12, 41988, 2014.9562, 398.79, 399.70, 398.63, 399.52, 398.79, 399.70 +2015, 01, 42019, 2015.0411, 399.85, 399.80, 399.77, 399.71, 399.85, 399.80 +2015, 02, 42050, 2015.1260, 400.31, 399.56, 400.66, 399.90, 400.31, 399.56 +2015, 03, 42078, 2015.2027, 401.51, 399.99, 401.62, 400.08, 401.51, 399.99 +2015, 04, 42109, 2015.2877, 403.45, 400.69, 403.06, 400.29, 403.45, 400.69 +2015, 05, 42139, 2015.3699, 404.10, 400.74, 403.86, 400.50, 404.10, 400.74 +2015, 06, 42170, 2015.4548, 402.88, 400.33, 403.25, 400.73, 402.88, 400.33 +2015, 07, 42200, 2015.5370, 401.61, 400.81, 401.75, 400.98, 401.61, 400.81 +2015, 08, 42231, 2015.6219, 399.00, 400.52, 399.70, 401.25, 399.00, 400.52 +2015, 09, 42262, 2015.7068, 397.50, 400.95, 398.08, 401.55, 397.50, 400.95 +2015, 10, 42292, 2015.7890, 398.28, 401.86, 398.29, 401.85, 398.28, 401.86 +2015, 11, 42323, 2015.8740, 400.24, 402.48, 399.96, 402.17, 400.24, 402.48 +2015, 12, 42353, 2015.9562, 401.89, 402.81, 401.59, 402.49, 401.89, 402.81 +2016, 01, 42384, 2016.0410, 402.65, 402.60, 402.87, 402.81, 402.65, 402.60 +2016, 02, 42415, 2016.1257, 404.16, 403.41, 403.88, 403.12, 404.16, 403.41 +2016, 03, 42444, 2016.2049, 404.85, 403.30, 404.97, 403.39, 404.85, 403.30 +2016, 04, 42475, 2016.2896, 407.57, 404.77, 406.49, 403.67, 407.57, 404.77 +2016, 05, 42505, 2016.3716, 407.66, 404.28, 407.29, 403.92, 407.66, 404.28 +2016, 06, 42536, 2016.4563, 407.00, 404.48, 406.67, 404.16, 407.00, 404.48 +2016, 07, 42566, 2016.5383, 404.50, 403.73, 405.13, 404.39, 404.50, 403.73 +2016, 08, 42597, 2016.6230, 402.24, 403.79, 403.02, 404.61, 402.24, 403.79 +2016, 09, 42628, 2016.7077, 401.01, 404.49, 401.33, 404.83, 401.01, 404.49 +2016, 10, 42658, 2016.7896, 401.50, 405.08, 401.47, 405.04, 401.50, 405.08 +2016, 11, 42689, 2016.8743, 403.64, 405.88, 403.03, 405.25, 403.64, 405.88 +2016, 12, 42719, 2016.9563, 404.55, 405.47, 404.54, 405.44, 404.55, 405.47 +2017, 01, 42750, 2017.0411, 406.07, 406.02, 405.70, 405.63, 406.07, 406.02 +2017, 02, 42781, 2017.1260, 406.64, 405.89, 406.58, 405.81, 406.64, 405.89 +2017, 03, 42809, 2017.2027, 407.06, 405.52, 407.52, 405.97, 407.06, 405.52 +2017, 04, 42840, 2017.2877, 408.95, 406.17, 408.93, 406.14, 408.95, 406.17 +2017, 05, 42870, 2017.3699, 409.91, 406.53, 409.68, 406.30, 409.91, 406.53 +2017, 06, 42901, 2017.4548, 409.12, 406.56, 408.99, 406.46, 409.12, 406.56 +2017, 07, 42931, 2017.5370, 407.20, 406.39, 407.38, 406.61, 407.20, 406.39 +2017, 08, 42962, 2017.6219, 405.24, 406.77, 405.20, 406.77, 405.24, 406.77 +2017, 09, 42993, 2017.7068, 403.27, 406.75, 403.43, 406.92, 403.27, 406.75 +2017, 10, 43023, 2017.7890, 403.64, 407.24, 403.48, 407.07, 403.64, 407.24 +2017, 11, 43054, 2017.8740, 405.17, 407.43, 404.99, 407.22, 405.17, 407.43 +2017, 12, 43084, 2017.9562, 406.75, 407.68, 406.47, 407.37, 406.75, 407.68 +2018, 01, 43115, 2018.0411, 408.05, 408.00, 407.59, 407.52, 408.05, 408.00 +2018, 02, 43146, 2018.1260, 408.34, 407.59, 408.45, 407.68, 408.34, 407.59 +2018, 03, 43174, 2018.2027, 409.25, 407.72, 409.38, 407.83, 409.25, 407.72 +2018, 04, 43205, 2018.2877, 410.30, 407.52, 410.81, 408.01, 410.30, 407.52 +2018, 05, 43235, 2018.3699, 411.30, 407.91, 411.60, 408.20, 411.30, 407.91 +2018, 06, 43266, 2018.4548, 410.88, 408.31, 410.97, 408.43, 410.88, 408.31 +2018, 07, 43296, 2018.5370, 408.90, 408.09, 409.43, 408.66, 408.90, 408.09 +2018, 08, 43327, 2018.6219, 407.10, 408.63, 407.35, 408.92, 407.10, 408.63 +2018, 09, 43358, 2018.7068, 405.59, 409.08, 405.68, 409.19, 405.59, 409.08 +2018, 10, 43388, 2018.7890, 406.00, 409.61, 405.86, 409.46, 406.00, 409.61 +2018, 11, 43419, 2018.8740, 408.12, 410.38, 407.49, 409.73, 408.12, 410.38 +2018, 12, 43449, 2018.9562, 409.23, 410.15, 409.09, 409.99, 409.23, 410.15 +2019, 01, 43480, 2019.0411, 410.92, 410.87, 410.31, 410.25, 410.92, 410.87 +2019, 02, 43511, 2019.1260, 411.66, 410.90, 411.26, 410.49, 411.66, 410.90 +2019, 03, 43539, 2019.2027, 412.00, 410.46, 412.26, 410.71, 412.00, 410.46 +2019, 04, 43570, 2019.2877, 413.52, 410.73, 413.74, 410.93, 413.52, 410.73 +2019, 05, 43600, 2019.3699, 414.83, 411.43, 414.55, 411.15, 414.83, 411.43 +2019, 06, 43631, 2019.4548, 413.96, 411.39, 413.92, 411.37, 413.96, 411.39 +2019, 07, 43661, 2019.5370, 411.85, 411.04, 412.36, 411.58, 411.85, 411.04 +2019, 08, 43692, 2019.6219, 410.08, 411.62, 410.22, 411.80, 410.08, 411.62 +2019, 09, 43723, 2019.7068, 408.55, 412.05, 408.50, 412.02, 408.55, 412.05 +2019, 10, 43753, 2019.7890, 408.43, 412.06, 408.63, 412.23, 408.43, 412.06 +2019, 11, 43784, 2019.8740, 410.29, 412.55, 410.21, 412.45, 410.29, 412.55 +2019, 12, 43814, 2019.9562, 411.85, 412.78, 411.75, 412.66, 411.85, 412.78 +2020, 01, 43845, 2020.0410, 413.37, 413.32, 412.94, 412.87, 413.37, 413.32 +2020, 02, 43876, 2020.1257, 414.09, 413.33, 413.85, 413.08, 414.09, 413.33 +2020, 03, 43905, 2020.2049, 414.51, 412.94, 414.86, 413.27, 414.51, 412.94 +2020, 04, 43936, 2020.2896, 416.18, 413.35, 416.32, 413.47, 416.18, 413.35 +2020, 05, 43966, 2020.3716, 417.16, 413.75, 417.08, 413.67, 417.16, 413.75 +2020, 06, 43997, 2020.4563, 416.30, 413.75, 416.41, 413.88, 416.30, 413.75 +2020, 07, 44027, 2020.5383, 414.49, 413.71, 414.84, 414.09, 414.49, 413.71 +2020, 08, 44058, 2020.6230, 412.59, 414.16, 412.70, 414.31, 412.59, 414.16 +2020, 09, 44089, 2020.7077, 411.25, 414.77, 410.99, 414.53, 411.25, 414.77 +2020, 10, 44119, 2020.7896, 411.23, 414.85, 411.13, 414.75, 411.23, 414.85 +2020, 11, 44150, 2020.8743, 412.95, 415.22, 412.72, 414.97, 412.95, 415.22 +2020, 12, 44180, 2020.9563, 413.95, 414.88, 414.28, 415.19, 413.95, 414.88 +2021, 01, 44211, 2021.0411, -99.99, -99.99, -99.99, -99.99, -99.99, -99.99 +2021, 02, 44242, 2021.1260, -99.99, -99.99, -99.99, -99.99, -99.99, -99.99 +2021, 03, 44270, 2021.2027, -99.99, -99.99, -99.99, -99.99, -99.99, -99.99 +2021, 04, 44301, 2021.2877, -99.99, -99.99, -99.99, -99.99, -99.99, -99.99 +2021, 05, 44331, 2021.3699, -99.99, -99.99, -99.99, -99.99, -99.99, -99.99 +2021, 06, 44362, 2021.4548, -99.99, -99.99, -99.99, -99.99, -99.99, -99.99 +2021, 07, 44392, 2021.5370, -99.99, -99.99, -99.99, -99.99, -99.99, -99.99 +2021, 08, 44423, 2021.6219, -99.99, -99.99, -99.99, -99.99, -99.99, -99.99 +2021, 09, 44454, 2021.7068, -99.99, -99.99, -99.99, -99.99, -99.99, -99.99 +2021, 10, 44484, 2021.7890, -99.99, -99.99, -99.99, -99.99, -99.99, -99.99 +2021, 11, 44515, 2021.8740, -99.99, -99.99, -99.99, -99.99, -99.99, -99.99 +2021, 12, 44545, 2021.9562, -99.99, -99.99, -99.99, -99.99, -99.99, -99.99