{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "Cette semaine, après avoir travaillé sur [l'âge moyen des députés français](http://flothesof.github.io/age-moyen-deputes-francais.html), je vous propose de nous intéresser à celui des sénateurs.\n", "\n", "Premier constat : la tâche ne va pas être aussi facile que dans le cas des députés. En effet, je n'ai pas réussi à trouver une base de données avec les dates de naissance des sénateurs ou de fichier PDF comme dans le cas des députés avec toutes les données d'un coup. Il va donc falloir travailler un peu plus afin d'extraire les données du [site web du sénat et son index alphabétique](http://www.senat.fr/senateurs/senatl.html). \n", "\n", "Afin de faire ceci, nous allons utiliser les outils Python classique pour ce genre de tâche : `requests` et `beautiful soup`." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Comment extraire la date de naissance d'un sénateur de sa notice bibliographique ?" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Tout part de la notice bibliographique du sénateur. Par exemple, celle de [Michel Berson, sénateur de l'Essonne](http://www.senat.fr/senateur/berson_michel11068s.html).\n", "\n", "On peut télécherger le contenu de cette page à l'aide du code suivant :" ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": true }, "outputs": [], "source": [ "import requests" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" } ], "source": [ "r = requests.get(\"http://www.senat.fr/senateur/berson_michel11068s.html\")\n", "r" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "La réponse 200 signifie que le téléchargement de la page s'est bien passé. \n", "\n", "Pour extraire la date de naissance, il faut regarder le code HTML de la page et utiliser `beautiful soup`. On construit notre soupe :" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "collapsed": true }, "outputs": [], "source": [ "from bs4 import BeautifulSoup\n", "soup = BeautifulSoup(r.text, 'html.parser')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "La date de naissance se trouve dans un tag 'dd' lui-même contenu dans un tag 'dl'. On peut trouver ces éléments ainsi :" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "[
\n", " m.berson@senat.fr \n", "
,
\n", " \n", " \n", "
,
Né le 21 avril 1945
,
Cadre dans le secteur bancaire
,
\"Place\n", "
,
\"Circonscription\"\n", " \n", "
]" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "items = soup.find('dl').find_all('dd')\n", "items" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "La date de naissance qui nous intéresse est dans le deuxième item :" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "['Né le 21 avril 1945']" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" } ], "source": [ "[item.text for item in items if item.text.startswith('Né')]" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "On peut, sur le même principe, également écrire une fonction qui extrait le nom du sénateur de la page :" ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "'M. Michel Berson'" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" } ], "source": [ "soup.find('h1','title-01').text" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Maintenant que nous savons faire ça, nous pouvons écrire un fonction qui extrait la date de naissance de l'url de la page d'un sénateur, de même que son nom :" ] }, { "cell_type": "code", "execution_count": 7, "metadata": { "collapsed": false }, "outputs": [], "source": [ "def url2data(url):\n", " \"\"\"Télécharge la page de l'url et extrait la date de naissance du sénateur.\"\"\"\n", " r = requests.get(url)\n", " if r.ok:\n", " soup = BeautifulSoup(r.text, 'html.parser')\n", " birthdate = [item.text for item in soup.find('dl').find_all('dd') if item.text.startswith('Né')][0]\n", " name = soup.find('h1','title-01').text\n", " return (name, birthdate)\n", " else:\n", " raise(\"Something went wrong with the web request\")" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "On peut tester la fonction :" ] }, { "cell_type": "code", "execution_count": 8, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "('M. Michel Berson', 'Né le 21 avril 1945')" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" } ], "source": [ "url2data(\"http://www.senat.fr/senateur/berson_michel11068s.html\")" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Pour avoir toutes les dates de naissances, il nous faut maintenant appliquer cette procédure à toutes les urls que l'on trouve sur l'index alphabétique du sénat." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Répéter l'opération d'extraction pour chaque sénateur" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "A partir de la page d'index biographique, nous pouvons, de la même manière que précédemment, extraire les liens vers les différentes pages de sénateurs :" ] }, { "cell_type": "code", "execution_count": 11, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 11, "metadata": {}, "output_type": "execute_result" } ], "source": [ "r = requests.get(\"http://www.senat.fr/senateurs/senatl.html\")\n", "r" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "On fait une petite sousoupe :" ] }, { "cell_type": "code", "execution_count": 12, "metadata": { "collapsed": false }, "outputs": [], "source": [ "soup = BeautifulSoup(r.text, 'html.parser')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Et maintenant on cherche des liens qui sont dans des balises du type ul li a. Tant qu'à faire, on ne garde que ceux qui pointen vers une page de sénateur :" ] }, { "cell_type": "code", "execution_count": 13, "metadata": { "collapsed": false }, "outputs": [], "source": [ "urls = [link.get('href') for link in soup.select('ul li a') if link.get('href').startswith('/senateur/')]" ] }, { "cell_type": "code", "execution_count": 14, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "348" ] }, "execution_count": 14, "metadata": {}, "output_type": "execute_result" } ], "source": [ "len(urls)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Nous dénombrons bien 348 sénateurs. Passons à l'extraction des dates de naissance." ] }, { "cell_type": "code", "execution_count": 15, "metadata": { "collapsed": false }, "outputs": [], "source": [ "birthdates = {}\n", "for url in urls:\n", " name, bd = url2data(\"http://www.senat.fr\" + url)\n", " birthdates[name] = bd" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Regardons le résultat :" ] }, { "cell_type": "code", "execution_count": 16, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "{'M. Abdourahamane Soilihi': 'Né le 4 octobre 1959',\n", " 'M. Alain Anziani': 'Né le 30 mai 1951',\n", " 'M. Alain Bertrand': 'Né le 23 février 1951',\n", " 'M. Alain Chatillon': 'Né le 15 mars 1943',\n", " 'M. Alain Dufaut': 'Né le 2 janvier 1944',\n", " 'M. Alain Duran': 'Né le 2 janvier 1958',\n", " 'M. Alain Fouché': 'Né le 4 décembre 1942',\n", " 'M. Alain Gournac': 'Né le 13 septembre 1943',\n", " 'M. Alain Houpert': 'Né le 13 août 1957',\n", " 'M. Alain Joyandet': 'Né le 15 janvier 1954',\n", " 'M. Alain Marc': 'Né le 29 janvier 1957',\n", " 'M. Alain Milon': 'Né le 16 septembre 1947',\n", " 'M. Alain Néri': 'Né le 1er mai 1942',\n", " 'M. Alain Richard': 'Né le 29 août 1945',\n", " 'M. Alain Vasselle': 'Né le 27 juin 1947',\n", " 'M. Albéric de Montgolfier': 'Né le 6 juillet 1964',\n", " 'M. Alex Türk': 'Né le 25 janvier 1950',\n", " 'M. André Gattolin': 'Né le 24 juin 1960',\n", " 'M. André Reichardt': 'Né le 5 décembre 1949',\n", " 'M. André Trillard': 'Né le 24 octobre 1947',\n", " 'M. Antoine Karam': 'Né le 21 février 1950',\n", " 'M. Antoine Lefèvre': 'Né le 18 février 1966',\n", " 'M. Benoît Huré': 'Né le 5 juin 1953',\n", " 'M. Bernard Cazeau': 'Né le 27 avril 1939',\n", " 'M. Bernard Delcros': 'Né le 12 février 1953',\n", " 'M. Bernard Fournier': 'Né le 13 septembre 1946',\n", " 'M. Bernard Lalande': 'Né le 6 avril 1954',\n", " 'M. Bernard Saugey': 'Né le 3 mars 1943',\n", " 'M. Bruno Gilles': 'Né le 26 décembre 1960',\n", " 'M. Bruno Retailleau': 'Né le 20 novembre 1960',\n", " 'M. Bruno Sido': 'Né le 19 février 1951',\n", " 'M. Charles Guené': 'Né le 6 avril 1952',\n", " 'M. Charles Revet': 'Né le 9 novembre 1937',\n", " 'M. Christian Cambon': 'Né le 8 mars 1948',\n", " 'M. Christian Favier': 'Né le 20 février 1951',\n", " 'M. Christian Manable': 'Né le 19 juin 1948',\n", " 'M. Christian Namy': 'Né le 19 octobre 1938',\n", " 'M. Christophe Béchu': 'Né le 11 juin 1974',\n", " 'M. Christophe-André Frassa': 'Né le 4 février 1968',\n", " 'M. Claude Bérit-Débat': 'Né le 19 février 1946',\n", " 'M. Claude Haut': 'Né le 22 décembre 1944',\n", " 'M. Claude Kern': 'Né le 13 mars 1959',\n", " 'M. Claude Malhuret': 'Né le 8 mars 1950',\n", " 'M. Claude Nougein': 'Né le 4 décembre 1946',\n", " 'M. Claude Raynal': 'Né le 10 octobre 1957',\n", " 'M. Cyril Pellevat': 'Né le 18 février 1981',\n", " 'M. Cédric Perrin': 'Né le 20 janvier 1974',\n", " 'M. Daniel Chasseing': 'Né le 10 avril 1945',\n", " 'M. Daniel Dubois': 'Né le 5 février 1952',\n", " 'M. Daniel Gremillet': 'Né le 31 juillet 1953',\n", " 'M. Daniel Laurent': 'Né le 4 février 1949',\n", " 'M. Daniel Percheron': 'Né le 31 août 1942',\n", " 'M. Daniel Raoul': 'Né le 28 juillet 1941',\n", " 'M. Daniel Reiner': 'Né le 17 janvier 1941',\n", " 'M. David Assouline': 'Né le 16 juin 1959',\n", " 'M. David Rachline': 'Né le 2 décembre 1987',\n", " 'M. Didier Guillaume': 'Né le 11 mai 1959',\n", " 'M. Didier Mandelli': 'Né le 7 août 1964',\n", " 'M. Didier Marie': 'Né le 19 mai 1960',\n", " 'M. Didier Robert': 'Né le 26 avril 1964',\n", " 'M. Dominique Bailly': 'Né le 2 janvier 1960',\n", " 'M. Dominique Watrin': 'Né le 14 juin 1953',\n", " 'M. Dominique de Legge': 'Né le 18 février 1952',\n", " 'M. Francis Delattre': 'Né le 11 septembre 1946',\n", " 'M. Franck Montaugé': 'Né le 14 septembre 1961',\n", " 'M. François Baroin': 'Né le 21 juin 1965',\n", " 'M. François Bonhomme': 'Né le 5 mai 1969',\n", " 'M. François Calvet': 'Né le 1er avril 1953',\n", " 'M. François Commeinhes': 'Né le 10 novembre 1949',\n", " 'M. François Fortassin': 'Né le 2 août 1939',\n", " 'M. François Grosdidier': 'Né le 25 février 1961',\n", " 'M. François Marc': 'Né le 19 mars 1950',\n", " 'M. François Patriat': 'Né le 21 mars 1943',\n", " 'M. François Pillet': 'Né le 13 mai 1950',\n", " 'M. François Zocchetto': 'Né le 14 décembre 1958',\n", " 'M. François-Noël Buffet': 'Né le 28 août 1963',\n", " 'M. Félix Desplan': 'Né le 22 février 1943',\n", " 'M. Gaëtan Gorce': 'Né le 2 décembre 1958',\n", " 'M. Georges Labazée': 'Né le 16 juin 1943',\n", " 'M. Georges Patient': 'Né le 1er avril 1949',\n", " 'M. Gilbert Barbier': 'Né le 3 mars 1940',\n", " 'M. Gilbert Bouchet': 'Né le 8 janvier 1947',\n", " 'M. Gilbert Roger': 'Né le 9 novembre 1953',\n", " 'M. Guillaume Arnell': 'Né le 20 janvier 1958',\n", " 'M. Guy-Dominique Kennel': 'Né le 14 avril 1952',\n", " 'M. Gérard Bailly': 'Né le 28 janvier 1940',\n", " 'M. Gérard Collomb': 'Né le 20 juin 1947',\n", " 'M. Gérard Cornu': 'Né le 6 février 1952',\n", " 'M. Gérard César': 'Né le 19 décembre 1934',\n", " 'M. Gérard Dériot': 'Né le 1er novembre 1944',\n", " 'M. Gérard Larcher': 'Né le 14 septembre 1949',\n", " 'M. Gérard Longuet': 'Né le 24 février 1946',\n", " 'M. Gérard Miquel': 'Né le 17 juin 1946',\n", " 'M. Gérard Roche': 'Né le 29 décembre 1942',\n", " 'M. Henri Cabanel': 'Né le 9 mars 1959',\n", " 'M. Henri Tandonnet': 'Né le 14 décembre 1949',\n", " 'M. Henri de Raincourt': 'Né le 17 novembre 1948',\n", " 'M. Hervé Marseille': 'Né le 20 août 1954',\n", " 'M. Hervé Maurey': 'Né le 7 mai 1961',\n", " 'M. Hervé Poher': 'Né le 19 mai 1952',\n", " 'M. Hilarion Vendegou': 'Né le 4 septembre 1941',\n", " 'M. Hubert Falco': 'Né le 15 mai 1947',\n", " 'M. Hugues Portelli': 'Né le 22 décembre 1947',\n", " 'M. Jackie Pierre': 'Né le 27 mai 1946',\n", " 'M. Jacques Bigot': 'Né le 31 juillet 1952',\n", " 'M. Jacques Chiron': 'Né le 3 octobre 1949',\n", " 'M. Jacques Cornano': 'Né le 18 novembre 1956',\n", " 'M. Jacques Gautier': 'Né le 18 septembre 1946',\n", " 'M. Jacques Genest': 'Né le 6 août 1950',\n", " 'M. Jacques Gillot': 'Né le 4 mars 1948',\n", " 'M. Jacques Grosperrin': 'Né le 17 octobre 1955',\n", " 'M. Jacques Legendre': 'Né le 2 décembre 1941',\n", " 'M. Jacques Mézard': 'Né le 3 décembre 1947',\n", " 'M. Jacques-Bernard Magner': 'Né le 1er juin 1952',\n", " 'M. Jean Bizet': 'Né le 30 août 1947',\n", " 'M. Jean Desessard': 'Né le 6 septembre 1952',\n", " 'M. Jean Louis Masson': 'Né le 25 mars 1947',\n", " 'M. Jean Pierre Vogel': 'Né le 7 août 1956',\n", " 'M. Jean-Baptiste Lemoyne': 'Né le 15 septembre 1977',\n", " 'M. Jean-Claude Boulard': 'Né le 28 mars 1943',\n", " 'M. Jean-Claude Carle': 'Né le 9 juin 1948',\n", " 'M. Jean-Claude Frécon': 'Né le 3 septembre 1944',\n", " 'M. Jean-Claude Gaudin': 'Né le 8 octobre 1939',\n", " 'M. Jean-Claude Lenoir': 'Né le 27 décembre 1944',\n", " 'M. Jean-Claude Leroy': 'Né le 3 juin 1952',\n", " 'M. Jean-Claude Luche': 'Né le 7 août 1952',\n", " 'M. Jean-Claude Requier': 'Né le 4 octobre 1947',\n", " 'M. Jean-François Husson': 'Né le 17 avril 1961',\n", " 'M. Jean-François Longeot': 'Né le 27 décembre 1956',\n", " 'M. Jean-François Mayet': 'Né le 30 janvier 1940',\n", " 'M. Jean-François Rapin': 'Né le 11 mai 1966',\n", " 'M. Jean-Jacques Filleul': 'Né le 22 juin 1943',\n", " 'M. Jean-Jacques Lasserre': 'Né le 11 mars 1944',\n", " 'M. Jean-Jacques Lozach': 'Né le 8 février 1954',\n", " 'M. Jean-Jacques Panunzi': 'Né le 5 avril 1956',\n", " 'M. Jean-Louis Carrère': 'Né le 4 décembre 1944',\n", " 'M. Jean-Louis Tourenne': 'Né le 25 août 1944',\n", " 'M. Jean-Léonce Dupont': 'Né le 31 janvier 1955',\n", " 'M. Jean-Marc Gabouty': 'Né le 17 mai 1949',\n", " 'M. Jean-Marie Bockel': 'Né le 22 juin 1950',\n", " 'M. Jean-Marie Morisset': 'Né le 18 août 1947',\n", " 'M. Jean-Marie Vanlerenberghe': 'Né le 29 mars 1939',\n", " 'M. Jean-Noël Cardoux': 'Né le 18 décembre 1946',\n", " 'M. Jean-Noël Guérini': 'Né le 1er janvier 1951',\n", " 'M. Jean-Paul Emorine': 'Né le 20 mars 1944',\n", " 'M. Jean-Paul Fournier': 'Né le 16 octobre 1945',\n", " 'M. Jean-Pierre Bosino': 'Né le 15 janvier 1959',\n", " 'M. Jean-Pierre Caffet': 'Né le 8 novembre 1951',\n", " 'M. Jean-Pierre Cantegrit': 'Né le 2 juillet 1933',\n", " 'M. Jean-Pierre Godefroy': 'Né le 23 septembre 1944',\n", " 'M. Jean-Pierre Grand': 'Né le 18 novembre 1950',\n", " 'M. Jean-Pierre Leleux': 'Né le 8 mai 1947',\n", " 'M. Jean-Pierre Masseret': 'Né le 23 août 1944',\n", " 'M. Jean-Pierre Raffarin': 'Né le 3 août 1948',\n", " 'M. Jean-Pierre Sueur': 'Né le 28 février 1947',\n", " 'M. Jean-Pierre Vial': 'Né le 17 février 1951',\n", " 'M. Jean-Vincent Placé': 'Né le 12 mars 1968',\n", " 'M. Jean-Yves Leconte': 'Né le 31 octobre 1966',\n", " 'M. Jean-Yves Roux': 'Né le 18 janvier 1970',\n", " 'M. Jeanny Lorgeoux': 'Né le 2 janvier 1950',\n", " 'M. Joseph Castelli': 'Né le 4 juin 1948',\n", " 'M. Joël Guerriau': 'Né le 9 novembre 1957',\n", " 'M. Joël Labbé': 'Né le 18 octobre 1952',\n", " 'M. Jérôme Bignon': 'Né le 1er janvier 1949',\n", " 'M. Jérôme Durain': 'Né le 2 juin 1969',\n", " 'M. Ladislas Poniatowski': 'Né le 10 novembre 1946',\n", " 'M. Louis Duvernois': 'Né le 17 mai 1941',\n", " 'M. Louis Nègre': 'Né le 8 février 1947',\n", " 'M. Louis Pinton': 'Né le 24 octobre 1948',\n", " 'M. Louis-Jean de Nicolaÿ': 'Né le 18 septembre 1949',\n", " 'M. Loïc Hervé': 'Né le 8 juin 1980',\n", " 'M. Luc Carvounas': 'Né le 8 juin 1971',\n", " 'M. Marc Daunis': 'Né le 22 avril 1955',\n", " 'M. Marc Laménie': 'Né le 11 juillet 1956',\n", " 'M. Martial Bourquin': 'Né le 23 juillet 1952',\n", " 'M. Mathieu Darnaud': 'Né le 11 juillet 1975',\n", " 'M. Maurice Antiste': 'Né le 22 juin 1953',\n", " 'M. Maurice Vincent': 'Né le 20 septembre 1955',\n", " 'M. Michel Amiel': 'Né le 3 juillet 1954',\n", " 'M. Michel Berson': 'Né le 21 avril 1945',\n", " 'M. Michel Billout': 'Né le 19 février 1958',\n", " 'M. Michel Boutant': 'Né le 23 novembre 1956',\n", " 'M. Michel Bouvard': 'Né le 17 mars 1955',\n", " 'M. Michel Canevet': 'Né le 14 décembre 1960',\n", " 'M. Michel Delebarre': 'Né le 27 avril 1946',\n", " 'M. Michel Fontaine': 'Né le 6 mai 1952',\n", " 'M. Michel Forissier': 'Né le 12 septembre 1943',\n", " 'M. Michel Houel': 'Né le 8 novembre 1942',\n", " 'M. Michel Le Scouarnec': 'Né le 12 juillet 1949',\n", " 'M. Michel Magras': 'Né le 6 janvier 1954',\n", " 'M. Michel Mercier': 'Né le 7 mars 1947',\n", " 'M. Michel Raison': 'Né le 5 novembre 1949',\n", " 'M. Michel Savin': 'Né le 25 novembre 1958',\n", " 'M. Michel Vaspart': 'Né le 22 février 1948',\n", " 'M. Michel Vergoz': 'Né le 10 janvier 1950',\n", " 'M. Nuihau Laurey': 'Né le 29 décembre 1964',\n", " 'M. Olivier Cadic': 'Né le 22 avril 1962',\n", " 'M. Olivier Cigolotti': 'Né le 16 août 1963',\n", " 'M. Pascal Allizard': 'Né le 19 décembre 1962',\n", " 'M. Patrick Abate': 'Né le 28 décembre 1956',\n", " 'M. Patrick Chaize': 'Né le 22 mars 1963',\n", " 'M. Patrick Masclet': 'Né le 31 mai 1952',\n", " 'M. Paul Vergès': 'Né le 5 mars 1925',\n", " 'M. Philippe Adnot': 'Né le 25 août 1945',\n", " 'M. Philippe Bas': 'Né le 20 juillet 1958',\n", " 'M. Philippe Bonnecarrère': 'Né le 12 juillet 1955',\n", " 'M. Philippe Dallier': 'Né le 8 décembre 1962',\n", " 'M. Philippe Dominati': 'Né le 12 avril 1954',\n", " 'M. Philippe Esnol': 'Né le 2 avril 1954',\n", " 'M. Philippe Kaltenbach': 'Né le 9 janvier 1966',\n", " 'M. Philippe Leroy': 'Né le 3 février 1940',\n", " 'M. Philippe Madrelle': 'Né le 21 avril 1937',\n", " 'M. Philippe Mouiller': 'Né le 20 septembre 1969',\n", " 'M. Philippe Nachbar': 'Né le 26 septembre 1950',\n", " 'M. Philippe Paul': 'Né le 25 janvier 1965',\n", " 'M. Pierre Camani': 'Né le 19 février 1952',\n", " 'M. Pierre Charon': 'Né le 1er mars 1951',\n", " 'M. Pierre Frogier': 'Né le 16 novembre 1950',\n", " 'M. Pierre Laurent': 'Né le 1er juillet 1957',\n", " 'M. Pierre Médevielle': 'Né le 17 avril 1960',\n", " 'M. Pierre-Yves Collombat': 'Né le 18 juillet 1945',\n", " 'M. Rachel Mazuir': 'Né le 12 février 1940',\n", " 'M. Raymond Vall': 'Né le 13 janvier 1942',\n", " 'M. René Danesi': 'Né le 30 août 1947',\n", " 'M. René Vandierendonck': 'Né le 1er avril 1951',\n", " 'M. René-Paul Savary': 'Né le 3 janvier 1953',\n", " 'M. Richard Yung': 'Né le 22 septembre 1947',\n", " 'M. Robert Hue': 'Né le 19 octobre 1946',\n", " 'M. Robert Laufoaulu': 'Né le 7 juillet 1947',\n", " 'M. Robert Navarro': 'Né le 1er mai 1952',\n", " 'M. Robert del Picchia': 'Né le 10 novembre 1942',\n", " 'M. Roger Karoutchi': 'Né le 26 août 1951',\n", " 'M. Roger Madec': 'Né le 27 octobre 1950',\n", " 'M. Roland Courteau': 'Né le 24 février 1943',\n", " 'M. Ronan Dantec': 'Né le 5 août 1963',\n", " 'M. Rémy Pointereau': 'Né le 30 mars 1953',\n", " 'M. Serge Dassault': 'Né le 4 avril 1925',\n", " 'M. Serge Larcher': 'Né le 17 octobre 1945',\n", " 'M. Simon Sutour': 'Né le 18 août 1952',\n", " 'M. Stéphane Ravier': 'Né le 4 août 1969',\n", " 'M. Thani Mohamed Soilihi': 'Né le 20 juin 1972',\n", " 'M. Thierry Carcenac': 'Né le 19 décembre 1950',\n", " 'M. Thierry Foucaud': 'Né le 18 janvier 1954',\n", " 'M. Vincent Capo-Canellas': 'Né le 4 mai 1967',\n", " 'M. Vincent Delahaye': 'Né le 23 août 1959',\n", " 'M. Vincent Eblé': 'Né le 15 octobre 1957',\n", " 'M. Xavier Pintat': 'Né le 15 mars 1954',\n", " 'M. Yannick Botrel': 'Né le 9 décembre 1951',\n", " 'M. Yannick Vaugrenard': 'Né le 25 juin 1950',\n", " 'M. Yves Daudigny': 'Né le 23 février 1947',\n", " 'M. Yves Détraigne': 'Né le 21 décembre 1954',\n", " 'M. Yves Pozzo di Borgo': 'Né le 3 mai 1948',\n", " 'M. Yves Rome': 'Né le 25 avril 1950',\n", " 'M. Yvon Collin': 'Né le 10 avril 1944',\n", " 'M. Éric Bocquet': 'Né le 8 novembre 1957',\n", " 'M. Éric Doligé': 'Né le 25 mai 1943',\n", " 'M. Éric Jeansannetas': 'Né le 21 octobre 1962',\n", " 'Mme Agnès Canayer': 'Née le 21 septembre 1965',\n", " 'Mme Aline Archimbaud': 'Née le 2 novembre 1948',\n", " 'Mme Anne Chain-Larché': 'Née le 11 avril 1958',\n", " 'Mme Anne Emery-Dumas': 'Née le 30 août 1959',\n", " 'Mme Anne-Catherine Loisier': 'Née le 13 mars 1969',\n", " 'Mme Annick Billon': 'Née le 3 août 1967',\n", " 'Mme Annie David': 'Née le 17 janvier 1963',\n", " 'Mme Annie Guillemot': 'Née le 27 janvier 1956',\n", " 'Mme Bariza Khiari': 'Née le 3 septembre 1946',\n", " 'Mme Brigitte Gonthier-Maurin': 'Née le 23 avril 1956',\n", " 'Mme Brigitte Micouleau': 'Née le 10 septembre 1951',\n", " 'Mme Caroline Cayeux': 'Née le 1er novembre 1948',\n", " 'Mme Catherine Deroche': 'Née le 24 février 1953',\n", " 'Mme Catherine Di Folco': 'Née le 30 novembre 1960',\n", " 'Mme Catherine Génisson': 'Née le 22 avril 1949',\n", " 'Mme Catherine Morin-Desailly': 'Née le 6 juillet 1960',\n", " 'Mme Catherine Procaccia': 'Née le 13 octobre 1949',\n", " 'Mme Catherine Tasca': 'Née le 13 décembre 1941',\n", " 'Mme Catherine Troendlé': 'Née le 20 février 1961',\n", " 'Mme Chantal Deseyne': 'Née le 29 janvier 1957',\n", " 'Mme Chantal Jouanno': 'Née le 12 juillet 1969',\n", " 'Mme Christiane Hummel': 'Née le 8 mai 1942',\n", " 'Mme Christiane Kammermann': 'Née le 10 juillet 1932',\n", " 'Mme Christine Prunaud': 'Née le 4 septembre 1950',\n", " 'Mme Claire-Lise Campion': 'Née le 27 juillet 1951',\n", " 'Mme Claudine Lepage': 'Née le 10 août 1949',\n", " 'Mme Colette Giudicelli': 'Née le 24 novembre 1943',\n", " 'Mme Colette Mélot': 'Née le 20 avril 1947',\n", " 'Mme Corinne Bouchoux': 'Née le 5 janvier 1964',\n", " 'Mme Corinne Féret': 'Née le 17 septembre 1961',\n", " 'Mme Corinne Imbert': 'Née le 20 novembre 1958',\n", " 'Mme Cécile Cukierman': 'Née le 26 avril 1976',\n", " 'Mme Danielle Michel': 'Née le 3 avril 1947',\n", " 'Mme Delphine Bataille': 'Née le 16 septembre 1969',\n", " 'Mme Dominique Estrosi Sassone': 'Née le 14 novembre 1958',\n", " 'Mme Dominique Gillot': 'Née le 11 juillet 1949',\n", " 'Mme Esther Benbassa': 'Née le 27 mars 1950',\n", " 'Mme Evelyne Yonnet': 'Née le 9 octobre 1954',\n", " 'Mme Fabienne Keller': 'Née le 20 octobre 1959',\n", " 'Mme Françoise Cartron': 'Née le 27 mars 1949',\n", " 'Mme Françoise Férat': 'Née le 5 mars 1949',\n", " 'Mme Françoise Gatel': 'Née le 14 mars 1953',\n", " 'Mme Françoise Laborde': 'Née le 8 juillet 1958',\n", " 'Mme Frédérique Espagnac': 'Née le 19 avril 1972',\n", " 'Mme Gisèle Jourda': 'Née le 24 mars 1955',\n", " 'Mme Hermeline Malherbe': 'Née le 18 février 1969',\n", " 'Mme Hélène Conway-Mouret': 'Née le 13 septembre 1960',\n", " 'Mme Isabelle Debré': 'Née le 23 mai 1957',\n", " 'Mme Jacky Deromedi': 'Née le 15 juin 1944',\n", " 'Mme Jacqueline Gourault': 'Née le 20 novembre 1950',\n", " 'Mme Josette Durrieu': 'Née le 20 mars 1937',\n", " 'Mme Joëlle Garriaud-Maylam': 'Née le 20 mars 1955',\n", " 'Mme Karine Claireaux': 'Née le 15 novembre 1963',\n", " 'Mme Lana Tetuanui': 'Née le 23 février 1970',\n", " 'Mme Laurence Cohen': 'Née le 15 janvier 1953',\n", " 'Mme Leila Aïchi': 'Née le 14 mai 1970',\n", " 'Mme Marie Mercier': 'Née le 16 juin 1958',\n", " 'Mme Marie-Annick Duchêne': 'Née le 23 août 1940',\n", " 'Mme Marie-Christine Blandin': 'Née le 22 septembre 1952',\n", " 'Mme Marie-France Beaufils': 'Née le 22 novembre 1946',\n", " 'Mme Marie-Françoise Perol-Dumont': 'Née le 26 mai 1952',\n", " 'Mme Marie-Hélène Des Esgaulx': 'Née le 26 mai 1950',\n", " 'Mme Marie-Noëlle Lienemann': 'Née le 12 juillet 1951',\n", " 'Mme Marie-Pierre Monier': 'Née le 28 février 1958',\n", " 'Mme Maryvonne Blondin': 'Née le 20 janvier 1947',\n", " 'Mme Michelle Demessine': 'Née le 18 juin 1947',\n", " 'Mme Michelle Meunier': 'Née le 24 janvier 1956',\n", " 'Mme Michèle André': 'Née le 6 février 1947',\n", " 'Mme Mireille Jouve': 'Née le 30 décembre 1960',\n", " 'Mme Nathalie Goulet': 'Née le 24 mai 1958',\n", " 'Mme Nelly Tocqueville': 'Née le 2 mai 1950',\n", " 'Mme Nicole Bonnefoy': 'Née le 11 août 1958',\n", " 'Mme Nicole Bricq': 'Née le 10 juin 1947',\n", " 'Mme Nicole Duranton': 'Née le 13 octobre 1958',\n", " 'Mme Odette Herviaux': 'Née le 6 janvier 1948',\n", " 'Mme Pascale Gruny': 'Née le 18 février 1960',\n", " 'Mme Patricia Morhet-Richaud': 'Née le 28 juillet 1961',\n", " 'Mme Patricia Schillinger': 'Née le 18 janvier 1963',\n", " 'Mme Samia Ghali': 'Née le 10 juin 1968',\n", " 'Mme Sophie Joissains': 'Née le 25 octobre 1969',\n", " 'Mme Sophie Primas': 'Née le 6 juin 1962',\n", " 'Mme Stéphanie Riocreux': 'Née le 20 mai 1966',\n", " 'Mme Sylvie Goy-Chavent': 'Née le 23 mai 1963',\n", " 'Mme Sylvie Robert': 'Née le 28 août 1963',\n", " 'Mme Valérie Létard': 'Née le 13 octobre 1962',\n", " 'Mme Vivette Lopez': 'Née le 27 décembre 1954',\n", " 'Mme Éliane Assassi': 'Née le 2 octobre 1958',\n", " 'Mme Éliane Giraud': 'Née le 12 mars 1952',\n", " 'Mme Élisabeth Doineau': 'Née le 7 avril 1961',\n", " 'Mme Élisabeth Lamure': 'Née le 20 novembre 1947',\n", " 'Mme Évelyne Didier': 'Née le 14 mars 1948'}" ] }, "execution_count": 16, "metadata": {}, "output_type": "execute_result" } ], "source": [ "birthdates" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "On peut facilement construire une dataframe pandas avec ces données :" ] }, { "cell_type": "code", "execution_count": 17, "metadata": { "collapsed": true }, "outputs": [], "source": [ "import pandas as pd" ] }, { "cell_type": "code", "execution_count": 18, "metadata": { "collapsed": false }, "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", "
NomDate de naissance
0Mme Catherine GénissonNée le 22 avril 1949
1Mme Anne Chain-LarchéNée le 11 avril 1958
2M. Georges PatientNé le 1er avril 1949
3Mme Corinne BouchouxNée le 5 janvier 1964
4M. Martial BourquinNé le 23 juillet 1952
\n", "
" ], "text/plain": [ " Nom Date de naissance\n", "0 Mme Catherine Génisson Née le 22 avril 1949\n", "1 Mme Anne Chain-Larché Née le 11 avril 1958\n", "2 M. Georges Patient Né le 1er avril 1949\n", "3 Mme Corinne Bouchoux Née le 5 janvier 1964\n", "4 M. Martial Bourquin Né le 23 juillet 1952" ] }, "execution_count": 18, "metadata": {}, "output_type": "execute_result" } ], "source": [ "df = pd.DataFrame(data=[(key, birthdates[key]) for key in birthdates], \n", " columns=('Nom', 'Date de naissance'))\n", "df.head()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "On peut, comme dans mon billet précédent, calculer des dates de naissance numériques à l'aide de ces données :" ] }, { "cell_type": "code", "execution_count": 19, "metadata": { "collapsed": true }, "outputs": [], "source": [ "def convert2date(date_str):\n", " d = {'janvier': 1, 'février': 2, 'mars': 3, 'avril': 4,\n", " 'mai': 5, 'juin': 6, 'juillet': 7, 'août': 8, \n", " 'septembre': 9, 'octobre': 10, 'novembre': 11, 'décembre':12}\n", " splits = date_str.split(' ')[2:]\n", " \n", " if splits[0].startswith('1er'):\n", " splits.pop(0)\n", " splits.insert(0, '1')\n", " \n", " return pd.Timestamp(splits[2] + '/' + str(d[splits[1]]) + '/' + splits[0])" ] }, { "cell_type": "code", "execution_count": 20, "metadata": { "collapsed": false }, "outputs": [], "source": [ "df['Date de naissance2'] = df['Date de naissance'].map(convert2date)" ] }, { "cell_type": "code", "execution_count": 21, "metadata": { "collapsed": false }, "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", "
NomDate de naissanceDate de naissance2
0Mme Catherine GénissonNée le 22 avril 19491949-04-22
1Mme Anne Chain-LarchéNée le 11 avril 19581958-04-11
2M. Georges PatientNé le 1er avril 19491949-04-01
3Mme Corinne BouchouxNée le 5 janvier 19641964-01-05
4M. Martial BourquinNé le 23 juillet 19521952-07-23
\n", "
" ], "text/plain": [ " Nom Date de naissance Date de naissance2\n", "0 Mme Catherine Génisson Née le 22 avril 1949 1949-04-22\n", "1 Mme Anne Chain-Larché Née le 11 avril 1958 1958-04-11\n", "2 M. Georges Patient Né le 1er avril 1949 1949-04-01\n", "3 Mme Corinne Bouchoux Née le 5 janvier 1964 1964-01-05\n", "4 M. Martial Bourquin Né le 23 juillet 1952 1952-07-23" ] }, "execution_count": 21, "metadata": {}, "output_type": "execute_result" } ], "source": [ "df.head()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Enfin, on peut calculer l'âge de chaque sénateur :" ] }, { "cell_type": "code", "execution_count": 22, "metadata": { "collapsed": true }, "outputs": [], "source": [ "import datetime" ] }, { "cell_type": "code", "execution_count": 23, "metadata": { "collapsed": false }, "outputs": [], "source": [ "df['Âge'] = pd.Timestamp(datetime.date.today()) - df['Date de naissance2']\n", "df['Âge (années)'] = [x.days / 365.25 for x in df['Âge']]" ] }, { "cell_type": "code", "execution_count": 24, "metadata": { "collapsed": false }, "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", "
NomDate de naissanceDate de naissance2ÂgeÂge (années)
0Mme Catherine GénissonNée le 22 avril 19491949-04-2224379 days66.746064
1Mme Anne Chain-LarchéNée le 11 avril 19581958-04-1121103 days57.776865
2M. Georges PatientNé le 1er avril 19491949-04-0124400 days66.803559
3Mme Corinne BouchouxNée le 5 janvier 19641964-01-0519008 days52.041068
4M. Martial BourquinNé le 23 juillet 19521952-07-2323191 days63.493498
\n", "
" ], "text/plain": [ " Nom Date de naissance Date de naissance2 \\\n", "0 Mme Catherine Génisson Née le 22 avril 1949 1949-04-22 \n", "1 Mme Anne Chain-Larché Née le 11 avril 1958 1958-04-11 \n", "2 M. Georges Patient Né le 1er avril 1949 1949-04-01 \n", "3 Mme Corinne Bouchoux Née le 5 janvier 1964 1964-01-05 \n", "4 M. Martial Bourquin Né le 23 juillet 1952 1952-07-23 \n", "\n", " Âge Âge (années) \n", "0 24379 days 66.746064 \n", "1 21103 days 57.776865 \n", "2 24400 days 66.803559 \n", "3 19008 days 52.041068 \n", "4 23191 days 63.493498 " ] }, "execution_count": 24, "metadata": {}, "output_type": "execute_result" } ], "source": [ "df.head()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Analyse des données recueillies " ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "On peut maintenant tracer de jolis histogrammes :" ] }, { "cell_type": "code", "execution_count": 25, "metadata": { "collapsed": true }, "outputs": [], "source": [ "%matplotlib inline\n", "import matplotlib.pyplot as plt\n", "plt.style.use('ggplot')" ] }, { "cell_type": "code", "execution_count": 26, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 26, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYEAAAEiCAYAAAAf0W5kAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XtcVHX+P/DXzHBzZGQYLiqii4qXYMVa0TU1AS8V6n4l\n+0aPJJO+lZtoa1hmpqSWbdlXBRS1tkwta9O+G6St22PbZFovtYKXJJDyhlKmIMjFC5eZef/+8OH8\nHLmIMMMMndfzHz2Hc87nNXNmeM2cc2ZQiYiAiIgUSe3sAERE5DwsASIiBWMJEBEpGEuAiEjBWAJE\nRArGEiAiUjCWAFETTCYTXn/9dRQUFDg7SodTWFiI1157DVeuXHF2FLoFlsANYmJiMGPGjGaXWbp0\nKfr3799OiehmX3/9NdRqNc6ePevwsV577TXs3r0b06dPh9lsdvh4baFWq/HRRx85OwYAoKKiApMn\nT0bPnj2h1WqdHYdugSVwm+bNm4dvv/22xcv369cPr7zyigMTKY9KpXL4GEeOHME///lPfPbZZ3jg\ngQfw6quvOnzMX4tp06YhKSkJjz32mEO2/9RTT2HMmDEO2bYSuTk7QEej1Wpd+tWNyWSCmxt3a1tF\nRERg7969AICXXnrJyWk6lh07djg7gl3V19fD3d3d2TEchu8EGrFs2TJ0794dfn5+mD59us1xzaVL\nl6Jfv37W6Z9//hn//d//jYCAAHTq1AmhoaFYuXIlgGuHl06cOIGlS5dCrVZDo9HgzJkzAIBvv/0W\nUVFR0Gq1MBgMSEhIQGlpqU2OtLQ09OzZE97e3pg0aRI++ugjm0MhmzZtgru7O4xGI373u9/By8sL\nX331FYqKivDggw+iR48e6Ny5MyIiIrBlyxabbcfExODJJ59ESkoKunbtCl9fX6SkpEBE8Morr6Bb\nt24IDAzEokWLbNbr3bs3Xn75ZSQlJcHX1xddu3bFunXrUFdXhz/96U8wGAwIDg7G2rVrbda7fPky\n5syZg+DgYHTu3BlDhgxBZmbmLffFmjVr0LNnT3Tu3BmxsbHW++9GBw4cwH333QedTofAwEA8+OCD\nNss1t4+aMmPGDISGhkKr1aJv375YuHAh6uvrb2v/OCpbdnY2Bg8ejE6dOuHOO++E0WhssExJSQkS\nExMRGBiILl264J577sHu3butPzeZTJg7dy569uwJLy8vBAUFYerUqc2O++677yIsLAydOnWCn58f\noqOjb+u2Xn/ubN++HXfccQe8vb0RExOD48ePW5epqKjAtGnT8Jvf/AZarRYDBw7EqlWrbLaxYcMG\n62FBjUaD999/H8CtH2OnT5+GWq3Gvn37bG7Xze/W1Wo11qxZg4SEBOj1eus7mj//+c/o27cvvLy8\nEBgYiNjYWNTW1jZ7n3UIQlbR0dGi1+tl7ty58sMPP8iXX34pBoNBXn75ZesyS5YskX79+lmn//CH\nP8j48ePlyJEjcvr0aTEajfLxxx+LiEh5ebn07t1b5s2bJ+fPn5fz58+LxWKRc+fOSZcuXeTRRx+V\n/Px82bt3r0REREhUVJR1u3/729/Ezc1N1qxZI8ePH5fNmzdLUFCQqNVq+fnnn0VEZNOmTaJWq+X3\nv/+9GI1GOXXqlFy4cEHy8vJk7dq1kpeXJydPnpSMjAxxd3cXo9HY4La++OKLcuzYMdm4caOoVCqJ\njY2V+fPny7Fjx2Tz5s2iUqnkiy++sK4XEhIivr6+kpqaKidOnJDXXntNVCqVTJgwwTrv9ddfF7Va\nLUePHrUZLyYmRvbt2yenTp2Sd955Rzw9PWXXrl1N7o+srCxxc3OTtLQ0OXbsmLz33nvStWtXm/sg\nPz9fvL29ZenSpfLjjz/K999/L/Hx8dK/f3+pra295T5qjMVikUWLFklOTo6cPn1aduzYIUFBQbJk\nyZLb2j+OyHb27Fnp3LmzPPHEE3L06FH517/+JREREaJWq+XDDz8UEZGrV69KWFiYPPTQQ3Lw4EE5\nceKE/PnPfxYvLy8pLCwUEZGVK1dKz5495d///rcUFxdLbm6upKenNznugQMHxM3NTbZs2SJnzpyR\n77//XjZs2HBbt3XJkiXSuXNniY2NlUOHDsmRI0dkyJAhMnr0aOs4586dk+XLl8vhw4elqKhIPvzw\nQ9HpdLJp0yYREbl06ZIkJCTIyJEjpaSkRM6fPy81NTUicuvHWFFRkajVatm7d6/NbQsNDZWlS5da\np1Uqlfj7+8vatWvl5MmTcvz4cfn000+lS5cu8ve//12Ki4vlu+++k/T0dOvYHRlL4AbR0dFy5513\n2sybOXOmjBgxwjp9cwkMHjzY5gF0s5sfYCIiixYtkp49e0p9fb113nfffScqlUp2794tIiIjR46U\nxx57zGa9F198sdESuPlB3ZjJkyfLjBkzbG7rXXfdZbNMeHi4RERE2MwbPHiwzJs3zzodEhIiDzzw\ngHXaYrFIly5d5L/+679s5vn6+sratWtFRCQ7O1s6deokVVVVNtv+n//5H5tt3WzUqFHy6KOP2sx7\n/vnnbe6DxMREeeSRR2yWqampEa1WK5999pn1NjS3j1oiNTVV+vfvb51uyf5xRLaFCxdKSEiImM1m\n67zPP/9cVCqVtQQ2btwoPXv2tFlGRGTMmDGSnJwsIiJz5syRsWPHtnjczMxM0ev1Ul1d3ejPW3Jb\nlyxZIu7u7lJWVmZdZuvWraLRaKxF0Zg5c+bIvffea51+8sknJSYmxmaZljzGioqKRKVStagEnnrq\nKZtlUlNTZcCAAWIymZrM2VHx4PFNBg8ebDMdFBSEf/7zn00u/+yzz+KPf/wjdu7ciejoaEycOBH3\n3HNPs2MUFBRg+PDhNsfuIyIi4OPjg/z8fIwaNQoFBQVISEiwWe/uu+9udHuRkZE201evXsXSpUvx\n+eef45dffkFdXR3q6uoQExPT7G3t1q0bunfv3mBeSUlJk+upVCoEBAQgIiLCZl5gYKB1vdzcXNTW\n1iIoKMhmO/X19c1eaVVQUNDgEMWoUaNsDg/k5OTgxIkT0Ol0NsvV1tbi2LFjAFq3j9555x1s2LAB\nRUVFuHz5MkwmE+SGL9xtyf5xRLajR49i2LBhUKv//5HcUaNG2SyTm5uLX375BT4+Pjbz6+rqrOez\nHn/8cYwfPx6hoaEYP348xo8fjz/84Q9NHvseP348evfujZCQEIwfPx5jxozBlClT4Ofn1+LbClx7\nPhkMBptpEUFJSQmCg4MhIli+fDm2bt2Kn376CTU1Naivr0dISEiT98n129yax1hThg4dajMdHx+P\n1atXo1evXrj33nsxduxYxMXFwdvb+7a37WpYAjfx8PCwmVapVLBYLE0un5iYiNjYWHzxxRfIzs5G\nbGwspkyZYj1O2RYtuQpGo9E0yPz8889jx44dSE1NRf/+/dG5c2fMnTsXVVVVNsvd/IRXqVSNzrv5\n9t/uehaLBXq9Hrm5uTa/SIGG9/ftslgsmDZtGhYsWNBg29d/Qd3uPvrkk08we/ZsvPnmmxg9ejS6\ndOmCbdu2NTg/cqv944hsLWGxWBAWFoasrKwG414vgcGDB6OoqAhffvklsrOz8eyzzyIlJQX/+c9/\nGv3F1rlzZxw4cAB79+7Fv/71L7z11lt44YUXsGvXLtx1110tuq1A48+v65kBYMWKFVi+fDnS0tJw\n5513QqfTYdWqVdi5c+ctb/OtHmPXi/Pmn998ruf67b1RUFAQfvjhB2RnZ2PXrl1YtmwZ5s+fj/37\n96NHjx7NZnN1LAE76Nq1K6ZPn47p06cjNjYWU6dOxbp16+Dt7Q0PD48G15iHh4dj06ZNNlfyfPfd\nd6isrMSgQYMAAGFhYfjmm2/w9NNPW9f75ptvWpRn9+7dSEhIwIMPPgjg2oP+xx9/RLdu3exxc29b\nZGQkKioqcPXqVYSFhbV4vbCwMOzbtw8zZ860ztuzZ0+DbR85cgS9e/dudlvN7aOb7d69G7/73e8w\nZ84c67xTp041yHar/eOIbGFhYdiyZQtExPoLtLH75IMPPoBOp4O/v3+T42q1WkyePBmTJ0/GggUL\n0L17d3z99deYOHFio8urVCqMGjUKo0aNwpIlSxAWFoaPPvoId911V4tv663s3r0b999/P6ZPn26d\n9+OPP9os09hzqiWPsYCAAACwOZldUlKCn3/+uUXZ3N3dce+99+Lee+/FK6+8gq5duyIrKwuzZs1q\n0fquilcHtdEzzzyDf/zjHzh58iTy8/Pxt7/9Db169bI+gXv37o29e/eiuLgYZWVlAIDZs2ejqqoK\niYmJyM/Px549e/DYY48hKioKI0aMAAA899xz+Pjjj5GRkYETJ07g/fffxwcffADg1q9ABwwYgM8+\n+ww5OTkoKCjAjBkz2uXDVU0ZM2YMxo0bhylTpuCzzz7DqVOncPDgQWRkZGDDhg1Nrvfcc89h69at\nWL16NY4fP46NGzc2uMrppZdewtGjR/Hoo48iJycHRUVF1le2RUVFAG69j242YMAA5OXlYfv27Th5\n8iTS09MbXMnUkv3jiGwzZ85EaWkpnnrqKRQWFuKrr77CokWLbB4TCQkJ6N27NyZOnIgvv/wSp0+f\nxv79+/HGG29g+/btAK694v7oo49QUFCAoqIibNiwAW5ubk0eOtm+fTvS0tJw8OBBFBcXIzMzEz/9\n9BPCw8NbfFubcuMr8wEDBsBoNMJoNOLYsWNISUnB/v37bZbv3bs3CgsLUVBQgLKyMtTV1WHMmDEY\nO3Zss48xLy8vjBw5Em+++SaOHDmCAwcOYPr06fDy8mo2HwC89957ePfdd3HkyBGcOXMGW7ZswaVL\nl27rRY3LcsqZCBcVExPT4ITQsmXLpHfv3tbpm08Mz5o1SwYMGCBarVb8/f1l0qRJUlBQYP15bm6u\nDBkyRDp16iRqtVpOnz4tIiL/+c9/JCoqSrRarfj6+sqjjz4qpaWlNmOnpaVJcHCwaLVauf/+++Uv\nf/mLqNVqKS8vF5FrJ4bd3d0b3I7i4mK5//77xdvb23pVy80n0xq7rePGjZPHH3/cZt79998v06ZN\ns0737t1bXnvtNZtl+vXr1+Dk5h133CEpKSnW6ZqaGlmwYIH06dNHPD09pXv37hIbGyvZ2dkN8t9o\n9erV1vtg/Pjx8v7779ucfBUR+f777yUuLk4MBoNotVrp16+f/PGPf5SLFy+KyK330c3q6+vl6aef\nFj8/P/Hx8ZGEhARZu3atqNVqm+VutX8ckU1EZNeuXRIRESFeXl4yaNAgyc7Otrk6SOTalWlJSUkS\nHBwsnp6eEhwcLFOmTJHDhw+LiMjbb78tkZGR4uPjIzqdToYNGyY7duxocsx///vfMmbMGAkMDJRO\nnTpJ//795c0337RZ5la39ebnjojInj17bJ4XlZWV8vDDD4uPj4/4+/vL7Nmz5eWXX7Z5DpaXl8vE\niRPFx8dH1Gq1bN68WURa9hg7duyYREdHi7e3t/Tv318yMzMbPH5vvi9FRD799FMZMWKEGAwG6dy5\nswwaNEg2btzY7H7qKFQizf95yfXr1+PgwYPw8fHBihUrAACXLl1CWloaSktLERgYiOTkZOuxxszM\nTGRnZ0Oj0SAxMbHByUdqvVdeeQUZGRkNTtSSa+D+oY7oloeDYmJisHDhQpt5WVlZGDRoENLT0xEe\nHm59q/zTTz/hm2++QWpqKhYsWIB33323wUmYpuTn57cifvtrr5wmkwnLly9HXl4eTp48iXfffRcr\nVqzAU0891aL1eX/a180527p/HKUj3J8dISOgnJy3LIGBAwc2OFOem5uLqKgoAEB0dDRycnKs80eM\nGAGNRoPAwEB0797d5tOAzVHKHd5SKpUKRqMR48aNw29/+1ukpqZi0aJFLf4OG96f9nVzzrbuH0fp\nCPdnR8gIKCdnq64OqqyshF6vBwDo9XpUVlYCAMrLy21OLBkMBpSXl7cpoFJpNBr84x//cHYMagL3\nD/1a2OXqoPb4VkciIrK/Vr0T0Ov1qKiosP57/ZOJBoMBFy5csC5XVlZm8+nAG+Xn59u8jYmPj29N\nlHbHnPbFnPbVEXJ2hIxAx8q5bds263R4eLj10t2WaFEJyLXvGLJODxkyBEajEXFxcTAajdavLYiM\njMTq1asxadIklJeX49y5cwgNDW10m40Fdea17C2l0+lQXV3t7Bi3xJz25eicmhNHUffGfIdt/zqP\nF5fD3PcOh4/THO5z+woKCmpTYd2yBNLT01FQUIDq6mrMnDkT8fHxiIuLQ2pqKrKzsxEQEIDk5GQA\nQHBwMO6++24kJyfDzc0NTz75JA8VERG5sFuWwI0fnb9RSkpKo/MfeOABPPDAA21LRURE7YJfG0FE\npGAsASIiBWMJEBEpGEuAiEjBWAJERArGEiAiUjCWABGRgrEEiIgUjCVARKRgLAEiIgVjCRARKRhL\ngIhIwVgCREQKxhIgIlIwlgARkYKxBIiIFIwlQESkYCwBIiIFYwkQESkYS4CISMFYAkRECsYSICJS\nMJYAEZGCsQSIiBSMJUBEpGBuzg5A5Mo0Fy8A5aWo1bhBYzY5bByVqd5h2yZqDkuAqDnlpah7Y77D\nh/Gcs9jhYxA1hoeDiIgUjCVARKRgLAEiIgVjCRARKRhLgIhIwVgCREQKxhIgIlIwlgARkYK16cNi\nn3/+ObKzs6FSqdCrVy8kJSWhpqYGaWlpKC0tRWBgIJKTk6HVau2Vl4iI7KjV7wTKy8vxxRdfYPny\n5VixYgXMZjP27NmDrKwsDBo0COnp6QgPD0dmZqY98xIRkR216XCQxWJBTU0NzGYz6urqYDAYkJub\ni6ioKABAdHQ0cnJy7BKUiIjsr9WHgwwGAyZNmoSkpCR4enoiIiICERERqKyshF6vBwDo9XpUVlba\nLSwREdlXq98JXL58Gbm5uVi3bh3efvtt1NbWYvfu3Q2WU6lUbQpIRESO0+p3Anl5eQgMDIS3tzcA\nYNiwYfjhhx+g1+tRUVFh/dfHx6fR9fPz85Gfn2+djo+Ph06na22cduPh4cGcduTqOWs17fNFu+31\nYkmjcYPWyfe3q+/z6zpKTgDYtm2b9f/h4eEIDw9v8bqtfoT7+/vj2LFjqKurg7u7O/Ly8tC3b194\neXnBaDQiLi4ORqMRkZGRja7fWNDq6urWxmk3Op2OOe3I1XM68m8I3EhE2mUcs9nk9Pvb1ff5dR0p\nZ3x8fKvXb3UJhIaGYvjw4Zg/fz40Gg1CQkIwbtw41NTUIDU1FdnZ2QgICEBycnKrwxERkWO16b3u\nQw89hIceeshmnre3N1JSUtoUioiI2gc/MUxEpGAsASIiBWMJEBEpGEuAiEjBWAJERArGEiAiUjCW\nABGRgrEEiIgUjCVARKRgLAEiIgVjCRARKRhLgIhIwVgCREQKxhIgIlIwlgARkYKxBIiIFIwlQESk\nYCwBIiIFYwkQESkYS4CISMFYAkRECsYSICJSMJYAEZGCsQSIiBSMJUBEpGAsASIiBWMJEBEpGEuA\niEjBWAJERArGEiAiUjCWABGRgrEEiIgUzM3ZAYhul+biBaC8tF3GUpnq22UcImdhCVDHU16Kujfm\nt8tQnnMWt8s4RM7Cw0FERArGEiAiUrA2HQ66cuUK3nrrLRQXF0OlUmHmzJno3r070tLSUFpaisDA\nQCQnJ0Or1dorLxER2VGbSmDjxo246667MHfuXJjNZtTW1uLTTz/FoEGDMHnyZGRlZSEzMxMJCQn2\nyktERHbU6sNBV65cQWFhIWJiYgAAGo0GWq0Wubm5iIqKAgBER0cjJyfHPkmJiMjuWv1OoKSkBDqd\nDuvWrcPp06fRp08fJCYmorKyEnq9HgCg1+tRWVlpt7BERGRfrX4nYLFYcOrUKdx3331Yvnw5PD09\nkZWV1WA5lUrVpoBEROQ4rX4nYDAY4Ofnh759+wIAhg8fjqysLOj1elRUVFj/9fHxaXT9/Px85Ofn\nW6fj4+Oh0+laG6fdeHh4MKcdtSZnrab9Pt7SXi9i2mscjcYNWic/Ln7Nj01n2bZtm/X/4eHhCA8P\nb/G6rX426fV6+Pn54ezZswgKCkJeXh6Cg4MRHBwMo9GIuLg4GI1GREZGNrp+Y0Grq6tbG6fd6HQ6\n5rSj1uTUmE0OStOQiPyqxjGbTU5/XPyaH5vOoNPpEB8f3+r12/SS6vHHH8eaNWtgMpnQtWtXJCUl\nwWKxIDU1FdnZ2QgICEBycnJbhiAiIgdqUwmEhITg9ddfbzA/JSWlLZslIqJ2wk8MExEpGEuAiEjB\nWAJERArGEiAiUjCWABGRgrEEiIgUjCVARKRgLAEiIgVjCRARKRhLgIhIwVgCREQKxhIgIlIwlgAR\nkYKxBIiIFIwlQESkYCwBIiIFYwkQESkYS4CISMFYAkRECsYSICJSMJYAEZGCsQSIiBSMJUBEpGAs\nASIiBWMJEBEpGEuAiEjBWAJERArGEiAiUjCWABGRgrEEiIgUjCVARKRgLAEiIgVjCRARKRhLgIhI\nwVgCREQKxhIgIlIwt7ZuwGKxYMGCBTAYDJg/fz4uXbqEtLQ0lJaWIjAwEMnJydBqtfbISkREdtbm\ndwI7d+5Ejx49rNNZWVkYNGgQ0tPTER4ejszMzLYOQUREDtKmEigrK8OhQ4cwduxY67zc3FxERUUB\nAKKjo5GTk9O2hERE5DBtKoHNmzdj2rRpUKlU1nmVlZXQ6/UAAL1ej8rKyrYlJCIih2l1CRw8eBA+\nPj4ICQmBiDS53I0FQURErqXVJ4YLCwuRm5uLQ4cOoa6uDlevXsWaNWug1+tRUVFh/dfHx6fR9fPz\n85Gfn2+djo+Ph06na22cduPh4cGcdtSanLWaNl/P0GLt9SKmvcbRaNygdfLj4tf82HSWbdu2Wf8f\nHh6O8PDwFq/b6mfT1KlTMXXqVABAQUEBduzYgWeeeQZbtmyB0WhEXFwcjEYjIiMjG12/saDV1dWt\njdNudDodc9pRa3JqzCYHpWmouXe5HXEcs9nk9MfFr/mx6Qw6nQ7x8fGtXt/unxOIi4tDXl4e5syZ\ng++//x5xcXH2HoKIiOzELu+rw8LCEBYWBgDw9vZGSkqKPTZLREQO1n4HV4nI6VRubtCcOOr4gQwB\nMPv6O34cajOWAJGSVFehLn2pw4fxeHE5wBLoEPjdQURECsYSICJSMJYAEZGCsQSIiBSMJUBEpGAs\nASIiBWMJEBEpGEuAiEjBWAJERArGEiAiUjCWABGRgrEEiIgUjCVARKRgLAEiIgVjCRARKRhLgIhI\nwVgCREQKxhIgIlIw/nlJshvNxQtAeeltrVOrcYPGbLqtdVSm+ttanoiaxhIg+ykvRd0b8x0+jOec\nxQ4fg0gpeDiIiEjBWAJERArGEiAiUjCWABGRgrEEiIgUjCVARKRgLAEiIgVjCRARKRhLgIhIwVgC\nREQKxhIgIlIwlgARkYKxBIiIFKzV3yJaVlaGjIwMVFZWQqVSYezYsZgwYQIuXbqEtLQ0lJaWIjAw\nEMnJydBqtfbMTEREdtLqEtBoNJg+fTpCQkJQU1OD+fPnY/DgwcjOzsagQYMwefJkZGVlITMzEwkJ\nCfbMTEREdtLqw0F6vR4hISEAAC8vL/To0QNlZWXIzc1FVFQUACA6Oho5OTl2CUpERPZnl3MCJSUl\nOH36NPr374/Kykro9XoA14qisrLSHkMQEZEDtLkEampqsGrVKiQmJsLLy6vBz1UqVVuHICIiB2nT\nn5c0m81YuXIlRo8ejaFDhwK49uq/oqLC+q+Pj0+j6+bn5yM/P986HR8fD51O15Y47cLDw4M5m1Cr\naZ+/VtqeLyzaa6xf2zgajRu0TTz++Byyv23btln/Hx4ejvDw8Bav26Zn7fr16xEcHIwJEyZY5w0Z\nMgRGoxFxcXEwGo2IjIxsdN3GglZXV7clTrvQ6XTM2YTb/YPxrSUi7TJOe471axvHbDY1+fjjc8i+\ndDod4uPjW71+q0ugsLAQu3fvRq9evfDCCy9ApVLhkUceQVxcHFJTU5GdnY2AgAAkJye3OhwRETlW\nq0tg4MCB2Lp1a6M/S0lJaXUgIiJqP/zEMBGRgrEEiIgUjCVARKRg7XNNHzmV5uIFoLzU4eOoTPUO\nH4OI7IsloATlpah7Y77Dh/Gcs9jhYxCRffFwEBGRgrEEiIgUjCVARKRgLAEiIgVjCRARKRhLgIhI\nwVgCREQKxhIgIlIwlgARkYKxBIiIFIwlQESkYCwBIiIFYwkQESkYS4CISMFYAkRECsYSICJSMJYA\nEZGCsQSIiBSMJUBEpGAsASIiBWMJEBEpGEuAiEjBWAJERArGEiAiUjCWABGRgrk5O4CSaS5eAMpL\nHbLtWo0bNGYTAEBlqnfIGERNUbm5QXPiaKM/u/GxaReGAJh9/e23PYVhCThTeSnq3pjv8GE85yx2\n+BhENqqrUJe+tF2G8nhxOcASaDUeDiIiUjCWABGRgrEEiIgUjCVARKRgDjsxfPjwYWzatAkigpiY\nGMTFxTlqKCJSsOauRGqLBlcx/UqvQnJICVgsFmzYsAEvv/wyfH19sWDBAgwdOhQ9evRwxHBEpGTt\ndCXSr/UqJIeUwPHjx9G9e3cEBAQAAEaOHImcnJwOUQJuFy9Aqiub/Hmdmxs0prZf46zS+UDavBUi\norZxSAmUl5fDz8/POm0wGHD8+HFHDGV38lMR6la/0uTP6+w0jsefXga0ne20NSKi1uGHxW6i7hoE\n96l/bPrnKhUs0vbX8OquQbA0846DiKg9qETs8BvtJj/++CM++eQTLFy4EACQlZUFADYnh/Pz85Gf\nn2+djo+Pt3cMIiJF2LZtm/X/4eHhCA8Pb/nK4gBms1lmz54tJSUlUl9fL88//7wUFxc3u87WrVsd\nEcXumNO+mNO+OkLOjpBRRDk5HXI4SK1W44knnsCyZcsgIhgzZgyCg4MdMRQREbWBw84J3HnnnUhP\nT3fU5omIyA40S5YsWeLsENcFBgY6O0KLMKd9Mad9dYScHSEjoIycDjkxTEREHQO/O4iISMFYAkRE\nCuaUD4uVlZUhIyMDlZWVUKlUGDt2LCZMmIBLly4hLS0NpaWlCAwMRHJyMrRarTMiAgDq6+uxePFi\nmEwmmM3fh+U8AAAI+klEQVRmDB8+HA899JDL5QSufV/TggULYDAYMH/+fJfMOGvWLGi1WqhUKmg0\nGrz++usumfPKlSt46623UFxcDJVKhZkzZ6J79+4ulfPs2bNIS0uDSqWCiOD8+fN4+OGHMXr0aJfK\nCQCff/45srOzoVKp0KtXLyQlJaGmpsblcu7cuRNfffUVALjU76T169fj4MGD8PHxwYoVKwCg2VyZ\nmZnIzs6GRqNBYmIiBg8e3PwAdrhM9bZdvHhRTp06JSIiV69elT/96U/y008/yQcffCBZWVkiIpKZ\nmSlbtmxxRjwbNTU1InLtsw8vvfSSHDt2zCVz7tixQ9LT0+WNN94QEXHJjLNmzZLq6mqbea6YMyMj\nQ3bt2iUiIiaTSS5fvuySOa8zm80yY8YMKS0tdbmcZWVlMmvWLKmvrxcRkVWrVkl2drbL5Txz5ow8\n99xzUldXJ2azWV599VX55ZdfXCLn0aNH5dSpU/Lcc89Z5zWVq7i4WObNmycmk0nOnz8vs2fPFovF\n0uz2nXI4SK/XIyQkBADg5eWFHj16oKysDLm5uYiKigIAREdHIycnxxnxbHh6egK49q7AbDYDgMvl\nLCsrw6FDhzB27FjrPFfLCAAiArnpOgRXy3nlyhUUFhYiJiYGAKDRaKDVal0u543y8vLQtWtX+Pv7\nu2ROi8WCmpoamM1m1NXVwWAwuFzOn3/+GaGhoXB3d4darcYdd9yB/fv348CBA07POXDgQHTubPs9\nY03df7m5uRgxYgQ0Gg0CAwPRvXv3W35vm9O/O6ikpASnT59G//79UVlZCb1eD+BaUVRWOv+7dSwW\nC1588UWcP38e9913H0JDQ10u5+bNmzFt2jRcuXLFOs/VMgKASqXCsmXLoFarMW7cOIwdO9blcpaU\nlECn02HdunU4ffo0+vTpg8TERJfLeaN9+/Zh1KhRAFxvvxsMBkyaNAlJSUnw9PREREQEIiIiXC5n\nz5498fHHH+PSpUtwd3fHoUOH0KdPH1RUVLhUzuuauv/Ky8vRv39/63IGgwHl5eXNbsupJVBTU4NV\nq1YhMTERXl5eDX6uUqmckMqWWq3Gm2++iStXrmDFihUoLi5usIwzc14/VhgSEmLzXUw3c4X78tVX\nX4Wvry+qqqqwbNkyBAUFNVjG2TktFgtOnTqFJ554An379sWmTZus3311I2fnvM5kMiE3NxcJCQmN\n/tzZOS9fvozc3FysW7cOWq0Wq1atwu7duxss5+ycPXr0wOTJk7Fs2TJ4eXkhJCQEanXDAyXOztmU\ntuRyWgmYzWasXLkSo0ePxtChQwFca7TrzVtRUQEfHx9nxWtAq9UiLCwMhw8fdqmchYWFyM3NxaFD\nh1BXV4erV69izZo1LpXxOl9fXwBAly5dMHToUBw/ftzlchoMBvj5+aFv374AgOHDhyMrK8vlcl53\n+PBh9OnTB126dAHges+hvLw8BAYGwtvbGwAwbNgw/PDDDy6XEwBiYmKshwH/+te/ws/PzyVzAk3v\nZ4PBgAsXLliXKysrg8FgaHZbTrtEdP369QgODsaECROs84YMGQKj0QgAMBqNiIyMdFK6a6qqqqyH\nWOrq6pCXl4cePXq4VM6pU6di/fr1yMjIwLPPPovf/va3eOaZZ1wqIwDU1taipqYGwLV3gEeOHEGv\nXr1cLqder4efnx/Onj0L4NovseDgYJfLed2ePXswcuRI67Sr5fT398exY8dQV1cHEXHp+7OqqgoA\ncOHCBezfvx+jRo1ymZw3n09rKldkZCT27dsHk8mEkpISnDt3DqGhoc1u2ymfGC4sLMTixYvRq1cv\nqFQqqFQqPPLIIwgNDUVqaiouXLiAgIAAJCcnNzgh0p7OnDmDtWvXwmKxQEQwYsQITJkyBZcuXXKp\nnNcVFBRgx44d1ktEXSljSUkJ/vd//xcqlQpmsxn33HMP4uLiXC4nABQVFeHtt9+GyWRC165dkZSU\nBIvF4nI5a2trkZSUhIyMDHTq1AkAXPL+/OSTT7Bv3z5oNBqEhITg6aefRk1NjcvlXLx4MS5dugSN\nRoPp06cjPDzcJe7P9PR0FBQUoLq6Gj4+PoiPj8fQoUObzJWZmYldu3bBzc2tRZeI8msjiIgUjJ8Y\nJiJSMJYAEZGCsQSIiBSMJUBEpGAsASIiBWMJEBEpGEuAFO/MmTNISEhAXl6es6MQtTuWACnehx9+\niPnz5+P//u//GnzLKdGvndO/RZTImaqqqhAdHY2IiAgAsP6RDiKl4CeGSdGysrLw1VdfoaqqCv7+\n/nj44YcxbNgwWCwWbNmyBV9//TW0Wi0mTpyIjRs34q9//SvUajWuXLmC999/H4cOHYJarUZUVBQe\nfvhhl/2WSaKm8J0AKVq3bt3w6quvQq/X45tvvkFGRgZWr16N/fv347vvvsOKFSvg6emJlStX2qy3\ndu1a+Pr6Ys2aNaipqcHy5cvh7++PcePGOemWELUOzwmQog0fPtz6xznuvvtudOvWDcePH8e3336L\n2NhY+Pr6QqvVIi4uzrpORUUFDh8+jOnTp8PDwwNdunTBhAkTsHfvXmfdDKJW4zsBUrSvv/4af//7\n31FaWgrg2tdcV1VV4eLFi/D397cu5+fnZ/3/hQsXYDKZMGPGDOs8EbFZnqijYAmQYl24cAF/+ctf\nsHjxYuuf5HvhhRcAXPsDOGVlZTbLXufv7w8PDw+89957PAdAHR4PB5Fi1dTUQKVSQafTwWKxIDs7\n2/rnQ4cPH46dO3eivLwcly9fxvbt263r6fV6REREYPPmzbh69SpEBOfPn0dBQYGzbgpRq/GdAClW\ncHAwJk2ahIULF0KtVmP06NEYOHAgAGDcuHE4d+4c5s2bB61Wi9jYWBQUFFj/7uzs2bPx4YcfYu7c\nuaipqUFgYCAmT57szJtD1Cq8RJSoBQ4fPox33nkHa9eudXYUIrvi4SCiRtTV1eHQoUOwWCwoLy/H\nJ598gt///vfOjkVkd3wnQNSIuro6LF68GGfPnoWHhweGDBmCxMREeHl5OTsakV2xBIiIFIyHg4iI\nFIwlQESkYCwBIiIFYwkQESkYS4CISMFYAkRECvb/AMH57yO9u+AhAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "df['Âge (années)'].hist()\n", "plt.xlabel('âge')\n", "plt.title('histogramme des âges des sénateurs')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "On peut demander à pandas de nous calculer les données statistiques usuelles pour ces âges :" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Afin de comparer ces âges à la population française, nous réutilisons le code de mon post précédent :" ] }, { "cell_type": "code", "execution_count": 27, "metadata": { "collapsed": true }, "outputs": [], "source": [ "df_pop = pd.read_excel('http://www.insee.fr/fr/ppp/bases-de-donnees/donnees-detaillees/bilan-demo/fichiers-xls/pop-1janvier-fm.xls',\n", " skiprows=7,\n", " header=1,\n", " skipfooter=3)" ] }, { "cell_type": "code", "execution_count": 28, "metadata": { "collapsed": true }, "outputs": [], "source": [ "decimation_factor = 1000\n", "french_pop_ages = []\n", "for ind, value in enumerate(df_pop['Ensemble']):\n", " french_pop_ages += [ind for _ in range(int(value / decimation_factor))]" ] }, { "cell_type": "code", "execution_count": 29, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 29, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX4AAAEECAYAAAAvY19bAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XtclFX+B/DPACICozKECjO5tOGNIbwAYeKqEGr5cpVW\nm1K3X7pkeb+UVqbi0ote4nopvEBWGFq2SRd13faVS6ltoluYUDqGinkDBYFJBSfkMs/vD9cHRi4z\nwAxzeT7v18vXi4fnPM+c+QpfzpxznnNkgiAIICIiyXCxdQWIiKhjMfETEUkMEz8RkcQw8RMRSQwT\nPxGRxDDxExFJjJupAmlpaTh+/Di6deuGdevWNVlm27ZtyMvLQ+fOnTF37lwEBgZaup5ERGQhJlv8\n0dHRWL58ebPnc3NzUVJSgo0bN+L555/Hu+++a/aLa7Vas8s6O8aiHmNRj7Gox1jUa28sTCb+/v37\nw8vLq9nzOTk5GDlyJACgT58+0Ov1uH79ulkvzv/IeoxFPcaiHmNRj7GoZ/XEb4pOp4Ovr694rFAo\noNPp2ntbIiKyEg7uEhFJjMnBXVMUCgXKy8vF4/LycigUiibLarVao48oGo2mvS/vNBiLeoxFPcai\nHmNRT6PRIDMzUzxWq9VQq9VmX29W4hcEAc2t5RYeHo79+/dj2LBhOHPmDLy8vNC9e/cmyzZVuStX\nrphdWWcml8tRUVFh62rYBcaiHmNRr7lY1C2cCugrO6QO/zc8EZVuXZo85+3ugp1P9u2QegQEBLTr\nD6HJxJ+SkoJTp06hoqICs2fPhkajQW1tLWQyGWJjYzFkyBDk5uZi/vz58PDwwOzZs9tcGSKiptQt\nnIrrzSV3T2+4vvsPi73WtE/OoLLa0OQ5b3cX7O2g5G5NMlsvy8wW/x1s2dVjLOpJKRYtttw9vdF9\n274OicXEnfnYO62/1V+nPQICAtp1fbv7+ImILEJfadGWe0tMteqdHRM/EXUYU616SzLZZWPnrXpr\nssvE7+3tDZlMZutqdChXV1fI5XJbV8Mu2HMsBEFAZWXHDCQ6pQ5s1VdWGySd3Ftil4lfJpNJpl+T\nHIu9/kEiag27TPxEROZqrktHCn31bcXET0QOjV06rcfET0R2T+qzcCyNid8K0tLSkJ+fj5SUFFtX\nhajDWWPmTmW1AQdmRXDsz0KY+C3s4MGDOHnyJLZs2WKxe6pUKmRnZ+N3v/udxe5JZDUdOHOH2oaJ\n38Kio6MRHR1t0Xtac2prXV0dXF1drXZ/InOxO6fjMJptsGXLFoSFhaFfv34YOXIksrOzIQgCNm/e\njKioKISEhGD27Nm4ceMGAKCwsBAqlQqffPIJHn74YYSGhmLjxo3i/fLy8jBhwgQEBwcjLCwMK1as\nQG1tLQBg0qRJEAQBsbGx6NevH/bt2wcAyMrKwpgxYxAcHIy4uDj8/PPP4v1UKhUuXrwoHi9evBhr\n164FABw9ehTh4eFITU3F4MGD8eKLL0Kn0+HZZ59FcHAw1Go1Jk2aZPUYEt3r7iBtU/86avEzqWCL\nv5XOnTuHjIwMfPnll/Dz80NRURHq6uqQnp6Of//73/j888+hUCiwcuVKvPbaa0ZdPjk5OTh8+DAK\nCgowfvx4jBs3DkFBQXB1dUViYiIGDRqEK1eu4M9//jO2b9+O+Ph4fPbZZ1CpVPj666/Ru3dvAMDJ\nkyexZMkS7NixA6Ghofjss88wY8YMfPvtt+jUqZPJTwilpaW4ceMGvv/+exgMBmzYsAEBAQE4efIk\nBEHA8ePHrRpDIrIth038dTMntPsebemHdHV1RU1NDfLz8+Hj4wOlUgkA+PDDD/HGG2+gZ8+eAO60\nsiMjI7Fp0yYAd7prXnrpJbi7uyM4OBjBwcE4deoUgoKC8NBDD4n3VyqVmDZtGv773/8iPj5e/H7D\ntfR27tyJZ555BgMHDgQATJ48GRs3bsTx48cRGRnZ7BLaDd/DkiVL0KlTJwBAp06dcO3aNVy6dAmB\ngYGIiIhodVyIzMHuHPvgsInfVoNHgYGBSExMxIYNG3DmzBmMGjUKCQkJKCwsxHPPPQcXlzs/vIIg\nwM3NDaWlpeK1fn5+4tceHh7Q6/UAgF9++QWJiYn46aefUFVVhdraWoSGhjZbh8LCQnz66ad4//33\nxdeqqalBcXGxWe9BoVCISR8A5syZg3Xr1mHq1KmQyWSYOnUq5s6da35QiMzEOff2wWETvy1NnDgR\nEydOxK1bt/Dyyy/jjTfegFKpxPr16xEeHt6ofGFhYYv3W7ZsGR566CG8/fbb6NKlC9577z3861//\narZ8QEAAFixYgPnz5zd5vkuXLvjtt9/E49LSUqNlXO/tCvL09ERCQgISEhJw5swZPPnkkxg0aBCi\noqJarDdJV0cutkaWx8TfSufOnUNxcTEiIiLQqVMneHh4wGAw4JlnnkFycjJSUlKgVCpRXl6OH374\nAWPGjAGAFrtfbt26BW9vb3Tp0gUFBQXYsWMH7rvvPvF8jx49cPHiRXE657Rp0/Dcc89h+PDhGDx4\nMPR6PY4ePYpHHnkEnp6eCAkJwe7du/HKK6/gm2++wdGjR8VuoaZ89dVXCAoKQmBgILy8vODm5iZ+\nciFqUgtTNqd9cgaVO/ObPMfuHPvAxN9K1dXVWL16NQoKCuDm5obw8HD87W9/g5+fHwRBwJQpU3Dt\n2jX4+vpiwoQJYuK/t5Xd8HjlypV4+eWXkZaWhpCQEEycOBHZ2dni+RdffBGLFi3C7du3sWbNGowf\nPx5r167FihUrcOHCBXh4eCAiIgKPPPIIACAxMRGLFi3C9u3bMXbsWDz++OMtvqfz589jxYoV0Ol0\n6NatG5599lnxXkStxe4c+2eXO3BJadchciwd/bNpr78LdTMnNNvit9YOVvYaC1vgDlxE1OH+L+qv\n7M5xYEz8RNRqlZ082Z3jwPinmYhIYtjiJ6JmNbvJSY3eBrUhS2HiJ6JmNTdDp27mBGA6V+B0VOzq\nISKSGCZ+IiKJYVcPkcRx4TTpYeJ3IJMnT8bkyZPx9NNPt/raoqIixMTEID8/3+Ibu+Tk5GDx4sUo\nLS3Fpk2bxKeVbe3VV1+Fv78/Fi5caOuq2DU+aSs9TPxOaujQoVi3bh2GDx8O4M5yz6dPn7bKa61b\ntw7x8fGYMWOGVe7fVsnJybauApFdYuKndissLESfPn2aPS8IglW3jyTraXYVTq7A6dDYgdcGQ4cO\nxebNmxEdHQ21Wo2XXnoJ1dXV4vmdO3eKWzD+5S9/QUlJiXhOpVJh27ZtGDZsGEJDQ5GUlCSe27Bh\ng9FSy3e3bDQYGve/Xrx4ERqNBiEhIQgNDcX8+fPFdUwWLFiAoqIiTJ8+Hf369cPbb7/d6F4lJSWY\nMWMG1Go1hg8fjo8++sioHrNmzcLChQvRr18/PProozhx4kSTsYiKisLly5fx7LPPol+/fqipqcHk\nyZOxZs0axMXFISgoCJcuXcKuXbswatQo9OvXD1FRUfjwww/Fe9zdDnLr1q0YOHAgwsLCsGvXLvF8\nVVUVEhMTERkZieDgYPzpT3/C7du3AQAvvPACBg8ejODgYEyePBlnzpwRr2u45WRL20uWlJRg5syZ\nCA0NxbBhw7Bt27Ym36sk/W8Vzkb/Uj4yfS3ZLSb+NtqzZw/+/ve/48iRIzh37hxSUlIAAIcPH0Zy\ncjLeeecd5ObmQqlUYs6cOUbXfvnll+K//fv34+OPPxbPtbSKZ0OCIGD+/PnIy8vDoUOHcPXqVaxf\nvx4AsHHjRiiVSmzfvh2nT5/GrFmzGt1r9uzZUCqVyMvLw9atW5GcnIwjR46I57OysvDEE08gPz8f\nsbGxeO2115qsR3Z2NgICArBjxw6cPn1a3ODl888/x7p163DmzBkolUr4+fnhgw8+wOnTp7Fhwwb8\n9a9/xcmTJ8X7lJaW4tatWzh+/DjWrl2L5cuX4+bNmwCA119/HSdPnsS+ffug1WqxfPlycdnomJgY\nHDlyBD/++CNCQkIwb968Juu5detWcXvJn376Ca+++qoYx+nTpyMkJAS5ubnYtWsX0tPT8Z///KfJ\n+ziqaZ+cwcSd+U3+4wCu9DhsV8/EZhaIao32DGjNmDEDvXr1AnCnhb1y5UosXboUe/bswZQpU6BW\nqwHc2WQlODgYRUVF4jaN8+bNQ9euXdG1a1c899xz2LNnT6sHbAMDAxEYGAjgzo5aM2fOxJtvvmlU\nprmFV4uKivDDDz9g586d6NSpE9RqNaZMmYJPP/0Uw4YNAwA8/PDDGDVqFIA7g8rp6ekt1ufe19Jo\nNAgKCgIAuLi4ICYmRjwXGRmJkSNH4vvvv0dISAiAO9s/Llq0SCzr5eWFc+fOYdCgQdi1axe++OIL\n9OjRAwAQFhYm3uupp54Sv168eDHee+89VFZWwtvbuCuiue0l8/LyoNPpxAHg+++/H1OmTMHevXsx\nYsSIFt+zI+EALjXksInf1j/E/v7+4tcqlUrszikpKTHaQ9fT0xM+Pj64evWqmPibu7Y1ysrKkJCQ\ngO+++w56vR51dXXo3r27Wddeu3YN3bt3R5cuXYzq0bA7p+E2kV26dMHt27dhMBjM3qDl3mVjDxw4\ngDfffBO//PILBEFAVVUVBgwYIJ738fExuneXLl1w69Yt6HQ6VFdXi5vQNGQwGJCcnIwvvvgCOp0O\nMpkMMpkMOp2uUeKfPXs21q9f32h7ycLCQhQXF4t/qAVBgMFgQGRkpFnvk8gROWzit7WG+wgUFhaK\nm6z37NkTRUVF4jm9Xo9ff/3VKNlfuXJFHAwtKioSr/X09DTaMrGlPwjJyclwcXHBwYMH0bVrV+zf\nvx8rVqwQz7c0mNqzZ09cv34der0enp6eYj3ufoKxhIavX11djeeffx6bNm3C2LFj4eLigvj4eJOb\nwgN3Ps107twZFy5cMPpDAQC7d+9GVlYWMjMzoVQqcfPmTQQHBzd5Xy8vrya3lwwICEDv3r3x7bff\ntv9NEzkIdu61UUZGBq5evYpff/0VmzZtwsSJEwHc2Y93165dOHXqFG7fvo3k5GQMGTJEbO0DQFpa\nGm7cuIGioiKkp6eL1wYHB+O7775DUVERbt68iS1btjT7+pWVlfD09IS3tzeuXr2KtLQ0o/N+fn64\ndOmS0ffuJsSAgACEh4dj9erVuH37Nk6dOoWPP/7YaMDzXu3Zr6empgY1NTVQKBRwcXHBgQMH8M03\n35h1rUwmw1NPPYXExESUlJTAYDDghx9+QHV1NSorK+Hu7o5u3bpBr9dj9erVzf7B++qrr3DhwgUA\nMNpecvDgwfD29kZqaiqqqqpQV1eH06dP48cff2zz+yWyd0z8bfTEE09g6tSpGD58OB544AEsWLAA\nAPCHP/wBS5cuxcyZMxEWFoZLly4hNTXV6Nq72yE+9thjGD16tNi/P2LECEyYMAGjR4/GuHHjMHr0\naKPrGia1F198ESdOnMCAAQMwffp0jBs3zqjsvHnz8NZbb0GtVmPr1q2Nrt+yZQsuXbqEIUOGYObM\nmVi6dGmLm6u39AnC1IC0l5cXXn/9dbzwwgtQq9XYu3cvxo4d2+z97r3HypUr0b9/f4wZMwa9e/fG\n6tWrIQgCnnzySSiVSoSFhSEmJqbJje7vOn/+PJ5++mn07dsXcXFx4vaSLi4u2L59O7RaLR555BGE\nhoZi6dKlDrnTEwdwyVxmbb2Yl5eHjIwMCIKA6OhoxMXFGZ3X6/XYtGkTysrKYDAY8Mc//lEcGDTF\nEbdevPfhqNZQqVTIzs5uss+aTJs6dSp27txps+cC7HnrRWtsedjSFosdzd7zQkdq79aLJpsBBoMB\n6enpWL58OdavX4/s7GyjPmwA2L9/P+6//36sXbsWq1atwo4dO1BXV9euihE1VFVVhdraWtTW1jbq\nwiKi1jGZ+AsKCuDv7w8/Pz+4ubkhKioKOTk5RmVkMpk4KFlVVQW5XA5XV1fr1NgOtKe1ySdY26ag\noAADBgyAXq9vd2uHSOpMzurR6XTw9fUVjxUKBQoKCozKPPbYY1izZg1eeOEFVFVVYdGiRZavqR05\nevRom6+9fPmyBWsiHSEhITh79qytq0HkFCwynTMvLw8PPPAAVq1aheLiYiQlJWHdunXw8PAwKqfV\naqHVasVjjUYDuVze6H7O/GmBHJurq2uTP7PW4u7u3qrXa0vdbsRPgHCr6b5zmZe8Q99vS1obC2eX\nmZkpfq1Wq8VnUcxhMvErFAqUlZWJxzqdDgqFwqjMoUOHxAHfXr16oUePHigqKsKDDz5oVK6pyjU1\nWMP/XLJXdXV1Nh3cNbV2flvqJtyqaHEA114GVDm4W08ul0Oj0bT5epOJPygoCMXFxSgtLYWPjw+y\ns7MbrW9+33334cSJE+jfvz+uX7+Oq1evig8lEZHlcOkFsgSTif/uU5ZJSUkQBAExMTFQqVTIysqC\nTCZDbGwsJk2ahNTUVCxZsgQAMG3atEaPzLeGIAiSa/W7urpyJtT/2HMs2vMgG5G9MGsevzU1NY9f\nivgxth5jUe/eWDj7XP2W8OeintXn8RMRkXPhIm1Edoabn5O1MfET2ZmGA7js3iBrYPOBiEhimPiJ\niCSGiZ+ISGKY+ImIJIaDu0Q20tzsHc7cIWtj4ieyES6/QLbCpgURkcQw8RMRSQwTPxGRxLCPn8iK\nuPwC2SMmfiIr4gAu2SM2OYiIJIaJn4hIYpj4iYgkhomfiEhiOLhL1E6cuUOOhomfqJ04c4ccDZsj\nREQSw8RPRCQxTPxERBLDxE9EJDEc3CUyA2fukDNh4icyg6PP3KlbOBXQVzZ90tO7YytDNsfETyQF\n+kq4vvsPW9eC7AQ/oxIRSQwTPxGRxDDxExFJDPv4if6HM3dIKpj4if7H0WfuEJmLzRgiIolh4ici\nkhgmfiIiiTGrjz8vLw8ZGRkQBAHR0dGIi4trVEar1WL79u2oq6tD165dsWrVKotXloiI2s9k4jcY\nDEhPT0dCQgJ8fHywbNkyREREQKlUimX0ej3S09OxYsUKKBQK3Lx506qVJiKitjPZ1VNQUAB/f3/4\n+fnBzc0NUVFRyMnJMSpz+PBhREZGQqFQAAC6du1qndoSEVG7mWzx63Q6+Pr6iscKhQIFBQVGZa5c\nuYK6ujokJiaiqqoKjz/+OEaMGGH52hIRUbtZZB6/wWDA+fPnkZCQgNu3b2PFihXo27cvevXqZYnb\nE1kMH9IiMiPxKxQKlJWVicc6nU7s0mlYRi6Xw93dHe7u7hgwYAAuXLjQKPFrtVpotVrxWKPRQC6X\nt/c9OAV3d3fG4n+sGYvKagMOzIqwyr2twVKxuA44/M8Xf0eMZWZmil+r1Wqo1WqzrzWZ+IOCglBc\nXIzS0lL4+PggOzsbCxcuNCoTERGBbdu2wWAwoKamBmfPnsX48eMb3aupylVUVJhdWWcml8sZi/+x\ndiwcKc6WjIUjve+m8Heknlwuh0ajafP1JhO/i4sL4uPjkZSUBEEQEBMTA5VKhaysLMhkMsTGxkKp\nVGLgwIFYsmQJXFxcEBsbC5VK1eZKERGR9ZjVxz9o0CCkpKQYfW/06NFGxxMmTMCECRMsVzMiIrIK\njmYREUkMV+ckp9Tc7B3O3CFi4icnxSWWiZrH5g8RkcQw8RMRSQy7eoicRN3CqYC+sumTnt4dWxmy\na0z8RM5CXwnXd/9h61qQA2DiJ4fFdXeI2oaJnxwWZ+4QtQ2bRUREEsPET0QkMUz8REQSw8RPRCQx\nHNwlu8aZO0SWx8RPdo0zd4gsj00mIiKJYeInIpIYJn4iIolh4icikhgO7pLNceYOUcdi4iebu3fm\njlwuR0VFhQ1rROTc2JwiIpIYJn4iIolh4icikhgmfiIiieHgLnUIztwhsh9M/NQhuOYOkf1gU4uI\nSGKY+ImIJIaJn4hIYtjHTxbDAVzrq1s4FdBXNn3S07tjK0MOi4mfLIYDuB1AXwnXd/9h61qQg2Mz\njIhIYpj4iYgkhomfiEhimPiJiCTGrMHdvLw8ZGRkQBAEREdHIy4urslyBQUFWLlyJRYtWoTIyEiL\nVpTsR3Ozdzhzh8gxmEz8BoMB6enpSEhIgI+PD5YtW4aIiAgolcpG5T766CMMHDjQapUl+8DZO0SO\nzWQTraCgAP7+/vDz84ObmxuioqKQk5PTqNyXX36JoUOHomvXrlapKBERWYbJxK/T6eDr6yseKxQK\n6HS6RmVycnIwZswYy9eQiIgsyiKdshkZGZg2bZp4LAiCJW5LRERWYLKPX6FQoKysTDzW6XRQKBRG\nZX755Re89dZbEAQBFRUVyM3NhZubG8LDw43KabVaaLVa8Vij0UAul7f3PTgFd3d3h4qFNevqaLGw\npntjcR3Wjb0948+FsczMTPFrtVoNtVpt9rUmE39QUBCKi4tRWloKHx8fZGdnY+HChUZlNm/eLH6d\nmpqKsLCwRkm/ucpVVFSYXVlnJpfL7SoWptbdsWZd7S0WttRULKQaG/5c1JPL5dBoNG2+3mTid3Fx\nQXx8PJKSkiAIAmJiYqBSqZCVlQWZTIbY2Ng2vzjZL87cIXJeZs3jHzRoEFJSUoy+N3r06CbLzpkz\np/21IiIiq+ETN0REEsNlmYnsTMM196/fe5Jr7pMFMPFLGDdOsVMN1tzngCZZAxO/hHEAl0ia2Kwj\nIpIYJn4iIolh4icikhj28Ts5DuAS0b2Y+J0cB3CJ6F5s8hERSQwTPxGRxLCrh8hGGj6ha4RP55KV\nMfET2UqDJ3SJOhITvxPgzB0iag0mfifAmTtE1BpsDhIRSQwTPxGRxDDxExFJDPv4HQQHcInIUpj4\nHQQHcInIUthUJCKSGCZ+IiKJYeInIpIY9vHbmeYGcTmAS0SWwsRvZziIS0TWxmYkEZHEMPETEUkM\nEz8RkcSwj98G+BQuEdkSE78NNDWAK5fLUVFRYaMakbU0u8sWwJ22yGaY+ImsibtskR1ivwIRkcQw\n8RMRSQy7eqyEA7hEZK+Y+K2ET+ASkb1i05OISGLMavHn5eUhIyMDgiAgOjoacXFxRucPHz6MvXv3\nAgA8PDwwc+ZM9O7d2/K1JSKidjPZ4jcYDEhPT8fy5cuxfv16ZGdno6ioyKhMjx49kJiYiLVr12LS\npEnYunWr1SpMRETtY7LFX1BQAH9/f/j5+QEAoqKikJOTA6VSKZbp27ev+HWfPn2g0+msUFX7wwFc\nInJEJhO/TqeDr6+veKxQKFBQUNBs+a+//hqDBg2yTO3sHAdwicgRWXRWz8mTJ3Ho0CG8/vrrTZ7X\narXQarXisUajgVwut2QVOpyl6u/u7u7wsbAUZ4rFdbTvZ8SZYtFejIWxzMxM8Wu1Wg21Wm32tSYT\nv0KhQFlZmXis0+mgUCgalbt48SLeeecdvPbaa/D2bnoNkqYq5+jr01iq/lyrp56zxaI978XZYtEe\njEU9uVwOjUbT5utNJv6goCAUFxejtLQUPj4+yM7OxsKFC43KlJWVYf369Zg3bx569erV5soQOSIu\nxEaOxmTid3FxQXx8PJKSkiAIAmJiYqBSqZCVlQWZTIbY2Fh8+umnqKysRHp6OgRBgKurK1avXt0R\n9SeyPS7ERg5GJgiCYMsKXLlyxZYvb5KpmTs7n+zb5LnW4sfYeo4Wi7qZE6yW+B0tFtbEWNQLCAho\n1/VcssEEztwhImfDyeZERBLDxE9EJDFM/EREEsPET0QkMRzcBdfcISJpYeIHZ+4QkbSwOUtEJDFs\n8ROZgcsykDNh4icyB5dlICciqcTf3CAuB3CJSEoklfg5iEtExMFdIiLJYeInIpIYJn4iIolh4ici\nkhinG9zl8gtERC1zusTPmTvUVnxIi6TC6RI/UZvxIS2SCPZ9EBFJDBM/EZHEOGRXDwdwiYjaziET\nPwdwqa04gEvkoImfqM04gEvEPn4iIqlhi5+cUrNdOuzOIWLiJyfFLh2iZtlt4ufMHSIi67DbxM+Z\nO0RE1sGmMxGRxDDxExFJDBM/EZHE2LyPf+LO/Ca/zwFcIiLrsHni5wAutRWXXyBqG5snfqI241x9\nojYxK/Hn5eUhIyMDgiAgOjoacXFxjcps27YNeXl56Ny5M+bOnYvAwEBL15UkiK16IsszmfgNBgPS\n09ORkJAAHx8fLFu2DBEREVAqlWKZ3NxclJSUYOPGjTh79izeffddvPHGG1atODmPe5P79YYnPb3Z\nqieyMJOJv6CgAP7+/vDz8wMAREVFIScnxyjx5+TkYOTIkQCAPn36QK/X4/r16+jevbuVqk2OxlTL\nvWFyl8vlqKio6KCaEUmPycSv0+ng6+srHisUChQUFJgso9PpmPglpjXJnYhsh4O7Tq7FZGxpTO5E\nDsFk4lcoFCgrKxOPdTodFApFozLl5eXicXl5eaMyAKDVaqHVasVjjUaDgICANlXcGcnlcsvf9JND\nlr9nB7BKLBwUY1GPsaiXmZkpfq1Wq6FWq82+1uRTUkFBQSguLkZpaSlqa2uRnZ2N8PBwozLh4eH4\n5ptvAABnzpyBl5dXk908arUaGo1G/New4lLHWNRjLOoxFvUYi3qZmZlGubQ1SR8wo8Xv4uKC+Ph4\nJCUlQRAExMTEQKVSISsrCzKZDLGxsRgyZAhyc3Mxf/58eHh4YPbs2W1+Q0REZF1m9fEPGjQIKSkp\nRt8bPXq00XF8fLzlakVERFZj0wVxWvvxxJkxFvUYi3qMRT3Gol57YyETBEGwUF2IiMgBcAlMIiKJ\nYeInIpIYmz3AZc7Cb86qvLwcmzdvxo0bNyCTyfDoo49i3LhxqKysxFtvvYXS0lL06NEDixcvhqen\np62ra3UGgwHLli2DQqHAK6+8Itk46PV6vP3227h8+TJkMhlmz54Nf39/Scbin//8Jw4ePAiZTIbe\nvXtjzpw5qKqqkkws0tLScPz4cXTr1g3r1q0DgBZ/L3bv3o2DBw/C1dUV06dPx8CBA1t+AcEG6urq\nhHnz5gnXrl0TampqhCVLlgiFhYW2qIpN/Prrr8L58+cFQRCE3377TViwYIFQWFgofPDBB8KePXsE\nQRCE3bt3Cx9++KENa9lx9u3bJ6SkpAjJycmCIAiSjcPmzZuFAwcOCIIgCLW1tcKtW7ckGYvy8nJh\n7ty5Qk3uLl5YAAADaElEQVRNjSAIgrBhwwbh4MGDkorFzz//LJw/f1546aWXxO819/4vX74sLF26\nVKitrRVKSkqEefPmCQaDocX726Srp+HCb25ubuLCb1LRvXt3cdlqDw8PKJVKlJeX49ixY+Jid6NG\njZJETMrLy5Gbm4tHH31U/J4U46DX65Gfn4/o6GgAgKurKzw9PSUZC+DOp8CqqirU1dWhuroaCoVC\nUrHo378/vLy8jL7X3Ps/duwYhg0bBldXV/To0QP+/v6N1lO7l026esxZ+E0qrl27hosXL6Jv3764\nceOG+MRz9+7dcePGDRvXzvq2b9+OZ555Bnq9XvyeFONw7do1yOVypKam4uLFi/j973+P6dOnSzIW\nCoUC48ePx5w5c9C5c2eEhoYiNDRUkrFoqLn3r9Pp0LdvX7Hc3UUyW8LBXRuqqqrChg0bMH36dHh4\neDQ6L5PJbFCrjnO3DzMwMBBCC7OKnT0OwJ0W7vnz5zF27FisWbMGnTt3xp49exqVk0Isbt26hWPH\njiE1NRVbt27F7du38e233zYqJ4VYtKQ9798mLX5zFn5zdnV1dVi/fj1GjBiBiIgIAHf+it/dx+D6\n9evo1q2bjWtpXfn5+Th27Bhyc3NRXV2N3377DZs2bZJcHIA7vxO+vr548MEHAQBDhw7Fnj17JBmL\nEydOoEePHvD2vrPD2sMPP4zTp09LMhYNNff+782nzS2S2ZBNWvzmLPzm7NLS0qBSqTBu3Djxe2Fh\nYTh06BAA4NChQ04fk6lTpyItLQ2bN2/GokWLEBISgvnz50suDsCdX2pfX19cuXIFwJ3kp1KpJBmL\n++67D2fPnkV1dTUEQZBsLARBMPok3Nz7Dw8Px5EjR1BbW4tr166huLgYQUFBLd7bZk/u5uXl4f33\n3xcXfpPSdM78/HysWrUKvXv3hkwmg0wmw5QpUxAUFIQ333wTZWVl8PPzw+LFixsN8DirU6dOYd++\nfeJ0TinG4cKFC9i6dStqa2vRs2dPzJkzBwaDQZKx+OSTT3DkyBG4uroiMDAQs2bNQlVVlWRikZKS\nglOnTqGiogLdunWDRqNBREREs+9/9+7dOHDgANzc3MyazsklG4iIJIaDu0REEsPET0QkMUz8REQS\nw8RPRCQxTPxERBLDxE9EJDFM/EREEsPET0QkMf8P9GL9OhZZbHQAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "import numpy as np\n", "import statsmodels.api as sm # recommended import according to the docs\n", "\n", "age_grid = np.linspace(0, 100)\n", "ecdf_mp = sm.distributions.ECDF(df['Âge (années)'].values)\n", "ecdf_french_pop = sm.distributions.ECDF(french_pop_ages)\n", "y_mp = ecdf_mp(age_grid)\n", "y_french_pop = ecdf_french_pop(age_grid)\n", "plt.step(age_grid, y_mp, label='sénateurs')\n", "plt.step(age_grid, y_french_pop, label='population française')\n", "plt.legend(loc='upper left')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Finalement, on peut passer à la comparaison des statistiques par quartiles :" ] }, { "cell_type": "code", "execution_count": 30, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "count 348.000000\n", "mean 62.791947\n", "std 9.050571\n", "min 28.134155\n", "25% 57.147844\n", "50% 63.954825\n", "75% 68.805613\n", "max 90.877481\n", "Name: Âge (années), dtype: float64" ] }, "execution_count": 30, "metadata": {}, "output_type": "execute_result" } ], "source": [ "df['Âge (années)'].describe()" ] }, { "cell_type": "code", "execution_count": 31, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "count 64461.000000\n", "mean 40.705931\n", "std 24.080480\n", "min 0.000000\n", "25% 20.000000\n", "50% 41.000000\n", "75% 60.000000\n", "max 105.000000\n", "dtype: float64" ] }, "execution_count": 31, "metadata": {}, "output_type": "execute_result" } ], "source": [ "pd.Series(french_pop_ages).describe()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Finalement, on peut produire le tableau récapitulatif suivant :" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "| population | âge médian | âge minimal | âge maximal | 25% de la population a moins de | 50% a moins de | 50% a plus de | 25% a plus de |\n", "|------------------------|------------|-------------|-------------|---------------------------------|----------------|---------------|---------------|\n", "| les Français | 41 ans | 0 ans | 105 ans | 20 ans | 41 ans | 41 ans | 59 ans |\n", "| les sénateurs français | 64 ans | 28 ans | 91 ans | 57 ans | 64 ans | 64 ans | 69 ans |" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Conclusions " ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Cette nouvelle analyse consolide l'idée selon laquelle les français sont représentés par des élus dans une tranche d'âge bien éloignée d'une répartition comparable à la population française.\n", "\n", "Evidemment, il est difficile de tirer des conclusions générales à partir de ce chiffre, mais il est frappant à quel point les sénateurs sont vieux par rapport disons, l'âge de la retraite. D'une certaine manière, la moitié des sénateurs ont plus que l'âge de la retraite." ] } ], "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.5.1" } }, "nbformat": 4, "nbformat_minor": 0 }