{ "cells": [ { "cell_type": "code", "source": [ "import pandas as pd\n", "import numpy as np\n", "from pandas import read_csv\n", "import math\n", "import plotly.plotly as py\n", "import plotly\n", "import plotly.graph_objs as go\n", "\n", "pd.set_option('display.max_colwidth',-1)\n", "plotly.tools.set_credentials_file(username='****', api_key='***')" ], "outputs": [], "execution_count": 1, "metadata": { "collapsed": false, "outputHidden": false, "inputHidden": false } }, { "cell_type": "code", "source": [ "raw_data = read_csv(\"master.csv\")\n", "\n", "# Suicide rates are typically reported as number of deaths per 100,000 people; \n", "# that is, the number of suicides in a country in a year, divided by the population and multiplied by 100,000.\n", "# so suicides/100k pop make sense hence lets ignore suicide counts and population\n", "# Similarly gdp per capita is more informative than raw GDP\n", "\n\n", "countries = [\n", "{'timezones': ['Europe/Andorra'], 'code': 'AD', 'continent': 'Europe', 'name': 'Andorra', 'capital': 'Andorra la Vella'},\n", "{'timezones': ['Asia/Kabul'], 'code': 'AF', 'continent': 'Asia', 'name': 'Afghanistan', 'capital': 'Kabul'},\n", "{'timezones': ['America/Antigua'], 'code': 'AG', 'continent': 'North America', 'name': 'Antigua and Barbuda', 'capital': \"St. John's\"},\n", "{'timezones': ['Europe/Tirane'], 'code': 'AL', 'continent': 'Europe', 'name': 'Albania', 'capital': 'Tirana'},\n", "{'timezones': ['Asia/Yerevan'], 'code': 'AM', 'continent': 'Asia', 'name': 'Armenia', 'capital': 'Yerevan'},\n", "{'timezones': ['Africa/Luanda'], 'code': 'AO', 'continent': 'Africa', 'name': 'Angola', 'capital': 'Luanda'},\n", "{'timezones': ['America/Argentina/Buenos_Aires', 'America/Argentina/Cordoba', 'America/Argentina/Jujuy', 'America/Argentina/Tucuman', 'America/Argentina/Catamarca', 'America/Argentina/La_Rioja', 'America/Argentina/San_Juan', 'America/Argentina/Mendoza', 'America/Argentina/Rio_Gallegos', 'America/Argentina/Ushuaia'], 'code': 'AR', 'continent': 'South America', 'name': 'Argentina', 'capital': 'Buenos Aires'},\n", "{'timezones': ['Europe/Vienna'], 'code': 'AT', 'continent': 'Europe', 'name': 'Austria', 'capital': 'Vienna'},\n", "{'timezones': ['Australia/Lord_Howe', 'Australia/Hobart', 'Australia/Currie', 'Australia/Melbourne', 'Australia/Sydney', 'Australia/Broken_Hill', 'Australia/Brisbane', 'Australia/Lindeman', 'Australia/Adelaide', 'Australia/Darwin', 'Australia/Perth'], 'code': 'AU', 'continent': 'Oceania', 'name': 'Australia', 'capital': 'Canberra'},\n", "{'timezones': ['Asia/Baku'], 'code': 'AZ', 'continent': 'Asia', 'name': 'Azerbaijan', 'capital': 'Baku'},\n", "{'timezones': ['America/Barbados'], 'code': 'BB', 'continent': 'North America', 'name': 'Barbados', 'capital': 'Bridgetown'},\n", "{'timezones': ['Asia/Dhaka'], 'code': 'BD', 'continent': 'Asia', 'name': 'Bangladesh', 'capital': 'Dhaka'},\n", "{'timezones': ['Europe/Brussels'], 'code': 'BE', 'continent': 'Europe', 'name': 'Belgium', 'capital': 'Brussels'},\n", "{'timezones': ['Africa/Ouagadougou'], 'code': 'BF', 'continent': 'Africa', 'name': 'Burkina Faso', 'capital': 'Ouagadougou'},\n", "{'timezones': ['Europe/Sofia'], 'code': 'BG', 'continent': 'Europe', 'name': 'Bulgaria', 'capital': 'Sofia'},\n", "{'timezones': ['Asia/Bahrain'], 'code': 'BH', 'continent': 'Asia', 'name': 'Bahrain', 'capital': 'Manama'},\n", "{'timezones': ['Africa/Bujumbura'], 'code': 'BI', 'continent': 'Africa', 'name': 'Burundi', 'capital': 'Bujumbura'},\n", "{'timezones': ['Africa/Porto-Novo'], 'code': 'BJ', 'continent': 'Africa', 'name': 'Benin', 'capital': 'Porto-Novo'},\n", "{'timezones': ['Asia/Brunei'], 'code': 'BN', 'continent': 'Asia', 'name': 'Brunei Darussalam', 'capital': 'Bandar Seri Begawan'},\n", "{'timezones': ['America/La_Paz'], 'code': 'BO', 'continent': 'South America', 'name': 'Bolivia', 'capital': 'Sucre'},\n", "{'timezones': ['America/Noronha', 'America/Belem', 'America/Fortaleza', 'America/Recife', 'America/Araguaina', 'America/Maceio', 'America/Bahia', 'America/Sao_Paulo', 'America/Campo_Grande', 'America/Cuiaba', 'America/Porto_Velho', 'America/Boa_Vista', 'America/Manaus', 'America/Eirunepe', 'America/Rio_Branco'], 'code': 'BR', 'continent': 'South America', 'name': 'Brazil', 'capital': 'Bras\\xc3\\xadlia'},\n", "{'timezones': ['America/Nassau'], 'code': 'BS', 'continent': 'North America', 'name': 'Bahamas', 'capital': 'Nassau'},\n", "{'timezones': ['Asia/Thimphu'], 'code': 'BT', 'continent': 'Asia', 'name': 'Bhutan', 'capital': 'Thimphu'},\n", "{'timezones': ['Africa/Gaborone'], 'code': 'BW', 'continent': 'Africa', 'name': 'Botswana', 'capital': 'Gaborone'},\n", "{'timezones': ['Europe/Minsk'], 'code': 'BY', 'continent': 'Europe', 'name': 'Belarus', 'capital': 'Minsk'},\n", "{'timezones': ['America/Belize'], 'code': 'BZ', 'continent': 'North America', 'name': 'Belize', 'capital': 'Belmopan'},\n", "{'timezones': ['America/St_Johns', 'America/Halifax', 'America/Glace_Bay', 'America/Moncton', 'America/Goose_Bay', 'America/Blanc-Sablon', 'America/Montreal', 'America/Toronto', 'America/Nipigon', 'America/Thunder_Bay', 'America/Pangnirtung', 'America/Iqaluit', 'America/Atikokan', 'America/Rankin_Inlet', 'America/Winnipeg', 'America/Rainy_River', 'America/Cambridge_Bay', 'America/Regina', 'America/Swift_Current', 'America/Edmonton', 'America/Yellowknife', 'America/Inuvik', 'America/Dawson_Creek', 'America/Vancouver', 'America/Whitehorse', 'America/Dawson'], 'code': 'CA', 'continent': 'North America', 'name': 'Canada', 'capital': 'Ottawa'},\n", "{'timezones': ['Africa/Kinshasa', 'Africa/Lubumbashi'], 'code': 'CD', 'continent': 'Africa', 'name': 'Democratic Republic of the Congo', 'capital': 'Kinshasa'},\n", "{'timezones': ['Africa/Brazzaville'], 'code': 'CG', 'continent': 'Africa', 'name': 'Republic of the Congo', 'capital': 'Brazzaville'},\n", "{'timezones': ['Africa/Abidjan'], 'code': 'CI', 'continent': 'Africa', 'name': \"C\\xc3\\xb4te d'Ivoire\", 'capital': 'Yamoussoukro'},\n", "{'timezones': ['America/Santiago', 'Pacific/Easter'], 'code': 'CL', 'continent': 'South America', 'name': 'Chile', 'capital': 'Santiago'},\n", "{'timezones': ['Africa/Douala'], 'code': 'CM', 'continent': 'Africa', 'name': 'Cameroon', 'capital': 'Yaound\\xc3\\xa9'},\n", "{'timezones': ['Asia/Shanghai', 'Asia/Harbin', 'Asia/Chongqing', 'Asia/Urumqi', 'Asia/Kashgar'], 'code': 'CN', 'continent': 'Asia', 'name': \"People's Republic of China\", 'capital': 'Beijing'},\n", "{'timezones': ['America/Bogota'], 'code': 'CO', 'continent': 'South America', 'name': 'Colombia', 'capital': 'Bogot\\xc3\\xa1'},\n", "{'timezones': ['America/Costa_Rica'], 'code': 'CR', 'continent': 'North America', 'name': 'Costa Rica', 'capital': 'San Jos\\xc3\\xa9'},\n", "{'timezones': ['America/Havana'], 'code': 'CU', 'continent': 'North America', 'name': 'Cuba', 'capital': 'Havana'},\n", "{'timezones': ['Atlantic/Cape_Verde'], 'code': 'CV', 'continent': 'Africa', 'name': 'Cape Verde', 'capital': 'Praia'},\n", "{'timezones': ['Asia/Nicosia'], 'code': 'CY', 'continent': 'Asia', 'name': 'Cyprus', 'capital': 'Nicosia'},\n", "{'timezones': ['Europe/Prague'], 'code': 'CZ', 'continent': 'Europe', 'name': 'Czech Republic', 'capital': 'Prague'},\n", "{'timezones': ['Europe/Berlin'], 'code': 'DE', 'continent': 'Europe', 'name': 'Germany', 'capital': 'Berlin'},\n", "{'timezones': ['Africa/Djibouti'], 'code': 'DJ', 'continent': 'Africa', 'name': 'Djibouti', 'capital': 'Djibouti City'},\n", "{'timezones': ['Europe/Copenhagen'], 'code': 'DK', 'continent': 'Europe', 'name': 'Denmark', 'capital': 'Copenhagen'},\n", "{'timezones': ['America/Dominica'], 'code': 'DM', 'continent': 'North America', 'name': 'Dominica', 'capital': 'Roseau'},\n", "{'timezones': ['America/Santo_Domingo'], 'code': 'DO', 'continent': 'North America', 'name': 'Dominican Republic', 'capital': 'Santo Domingo'},\n", "{'timezones': ['America/Guayaquil', 'Pacific/Galapagos'], 'code': 'EC', 'continent': 'South America', 'name': 'Ecuador', 'capital': 'Quito'},\n", "{'timezones': ['Europe/Tallinn'], 'code': 'EE', 'continent': 'Europe', 'name': 'Estonia', 'capital': 'Tallinn'},\n", "{'timezones': ['Africa/Cairo'], 'code': 'EG', 'continent': 'Africa', 'name': 'Egypt', 'capital': 'Cairo'},\n", "{'timezones': ['Africa/Asmera'], 'code': 'ER', 'continent': 'Africa', 'name': 'Eritrea', 'capital': 'Asmara'},\n", "{'timezones': ['Africa/Addis_Ababa'], 'code': 'ET', 'continent': 'Africa', 'name': 'Ethiopia', 'capital': 'Addis Ababa'},\n", "{'timezones': ['Europe/Helsinki'], 'code': 'FI', 'continent': 'Europe', 'name': 'Finland', 'capital': 'Helsinki'},\n", "{'timezones': ['Pacific/Fiji'], 'code': 'FJ', 'continent': 'Oceania', 'name': 'Fiji', 'capital': 'Suva'},\n", "{'timezones': ['Europe/Paris'], 'code': 'FR', 'continent': 'Europe', 'name': 'France', 'capital': 'Paris'},\n", "{'timezones': ['Africa/Libreville'], 'code': 'GA', 'continent': 'Africa', 'name': 'Gabon', 'capital': 'Libreville'},\n", "{'timezones': ['Asia/Tbilisi'], 'code': 'GE', 'continent': 'Asia', 'name': 'Georgia', 'capital': 'Tbilisi'},\n", "{'timezones': ['Africa/Accra'], 'code': 'GH', 'continent': 'Africa', 'name': 'Ghana', 'capital': 'Accra'},\n", "{'timezones': ['Africa/Banjul'], 'code': 'GM', 'continent': 'Africa', 'name': 'The Gambia', 'capital': 'Banjul'},\n", "{'timezones': ['Africa/Conakry'], 'code': 'GN', 'continent': 'Africa', 'name': 'Guinea', 'capital': 'Conakry'},\n", "{'timezones': ['Europe/Athens'], 'code': 'GR', 'continent': 'Europe', 'name': 'Greece', 'capital': 'Athens'},\n", "{'timezones': ['America/Guatemala'], 'code': 'GT', 'continent': 'North America', 'name': 'Guatemala', 'capital': 'Guatemala City'},\n", "{'timezones': ['America/Guatemala'], 'code': 'GT', 'continent': 'North America', 'name': 'Haiti', 'capital': 'Port-au-Prince'},\n", "{'timezones': ['Africa/Bissau'], 'code': 'GW', 'continent': 'Africa', 'name': 'Guinea-Bissau', 'capital': 'Bissau'},\n", "{'timezones': ['America/Guyana'], 'code': 'GY', 'continent': 'South America', 'name': 'Guyana', 'capital': 'Georgetown'},\n", "{'timezones': ['America/Tegucigalpa'], 'code': 'HN', 'continent': 'North America', 'name': 'Honduras', 'capital': 'Tegucigalpa'},\n", "{'timezones': ['Europe/Budapest'], 'code': 'HU', 'continent': 'Europe', 'name': 'Hungary', 'capital': 'Budapest'},\n", "{'timezones': ['Asia/Jakarta', 'Asia/Pontianak', 'Asia/Makassar', 'Asia/Jayapura'], 'code': 'ID', 'continent': 'Asia', 'name': 'Indonesia', 'capital': 'Jakarta'},\n", "{'timezones': ['Europe/Dublin'], 'code': 'IE', 'continent': 'Europe', 'name': 'Republic of Ireland', 'capital': 'Dublin'},\n", "{'timezones': ['Asia/Jerusalem'], 'code': 'IL', 'continent': 'Asia', 'name': 'Israel', 'capital': 'Jerusalem'},\n", "{'timezones': ['Asia/Calcutta'], 'code': 'IN', 'continent': 'Asia', 'name': 'India', 'capital': 'New Delhi'},\n", "{'timezones': ['Asia/Baghdad'], 'code': 'IQ', 'continent': 'Asia', 'name': 'Iraq', 'capital': 'Baghdad'},\n", "{'timezones': ['Asia/Tehran'], 'code': 'IR', 'continent': 'Asia', 'name': 'Iran', 'capital': 'Tehran'},\n", "{'timezones': ['Atlantic/Reykjavik'], 'code': 'IS', 'continent': 'Europe', 'name': 'Iceland', 'capital': 'Reykjav\\xc3\\xadk'},\n", "{'timezones': ['Europe/Rome'], 'code': 'IT', 'continent': 'Europe', 'name': 'Italy', 'capital': 'Rome'},\n", "{'timezones': ['America/Jamaica'], 'code': 'JM', 'continent': 'North America', 'name': 'Jamaica', 'capital': 'Kingston'},\n", "{'timezones': ['Asia/Amman'], 'code': 'JO', 'continent': 'Asia', 'name': 'Jordan', 'capital': 'Amman'},\n", "{'timezones': ['Asia/Tokyo'], 'code': 'JP', 'continent': 'Asia', 'name': 'Japan', 'capital': 'Tokyo'},\n", "{'timezones': ['Africa/Nairobi'], 'code': 'KE', 'continent': 'Africa', 'name': 'Kenya', 'capital': 'Nairobi'},\n", "{'timezones': ['Asia/Bishkek'], 'code': 'KG', 'continent': 'Asia', 'name': 'Kyrgyzstan', 'capital': 'Bishkek'},\n", "{'timezones': ['Pacific/Tarawa', 'Pacific/Enderbury', 'Pacific/Kiritimati'], 'code': 'KI', 'continent': 'Oceania', 'name': 'Kiribati', 'capital': 'Tarawa'},\n", "{'timezones': ['Asia/Pyongyang'], 'code': 'KP', 'continent': 'Asia', 'name': 'North Korea', 'capital': 'Pyongyang'},\n", "{'timezones': ['Asia/Seoul'], 'code': 'KR', 'continent': 'Asia', 'name': 'South Korea', 'capital': 'Seoul'},\n", "{'timezones': ['Asia/Kuwait'], 'code': 'KW', 'continent': 'Asia', 'name': 'Kuwait', 'capital': 'Kuwait City'},\n", "{'timezones': ['Asia/Beirut'], 'code': 'LB', 'continent': 'Asia', 'name': 'Lebanon', 'capital': 'Beirut'},\n", "{'timezones': ['Europe/Vaduz'], 'code': 'LI', 'continent': 'Europe', 'name': 'Liechtenstein', 'capital': 'Vaduz'},\n", "{'timezones': ['Africa/Monrovia'], 'code': 'LR', 'continent': 'Africa', 'name': 'Liberia', 'capital': 'Monrovia'},\n", "{'timezones': ['Africa/Maseru'], 'code': 'LS', 'continent': 'Africa', 'name': 'Lesotho', 'capital': 'Maseru'},\n", "{'timezones': ['Europe/Vilnius'], 'code': 'LT', 'continent': 'Europe', 'name': 'Lithuania', 'capital': 'Vilnius'},\n", "{'timezones': ['Europe/Luxembourg'], 'code': 'LU', 'continent': 'Europe', 'name': 'Luxembourg', 'capital': 'Luxembourg City'},\n", "{'timezones': ['Europe/Riga'], 'code': 'LV', 'continent': 'Europe', 'name': 'Latvia', 'capital': 'Riga'},\n", "{'timezones': ['Africa/Tripoli'], 'code': 'LY', 'continent': 'Africa', 'name': 'Libya', 'capital': 'Tripoli'},\n", "{'timezones': ['Indian/Antananarivo'], 'code': 'MG', 'continent': 'Africa', 'name': 'Madagascar', 'capital': 'Antananarivo'},\n", "{'timezones': ['Pacific/Majuro', 'Pacific/Kwajalein'], 'code': 'MH', 'continent': 'Oceania', 'name': 'Marshall Islands', 'capital': 'Majuro'},\n", "{'timezones': ['Europe/Skopje'], 'code': 'MK', 'continent': 'Europe', 'name': 'Macedonia', 'capital': 'Skopje'},\n", "{'timezones': ['Africa/Bamako'], 'code': 'ML', 'continent': 'Africa', 'name': 'Mali', 'capital': 'Bamako'},\n", "{'timezones': ['Asia/Rangoon'], 'code': 'MM', 'continent': 'Asia', 'name': 'Myanmar', 'capital': 'Naypyidaw'},\n", "{'timezones': ['Asia/Ulaanbaatar', 'Asia/Hovd', 'Asia/Choibalsan'], 'code': 'MN', 'continent': 'Asia', 'name': 'Mongolia', 'capital': 'Ulaanbaatar'},\n", "{'timezones': ['Africa/Nouakchott'], 'code': 'MR', 'continent': 'Africa', 'name': 'Mauritania', 'capital': 'Nouakchott'},\n", "{'timezones': ['Europe/Malta'], 'code': 'MT', 'continent': 'Europe', 'name': 'Malta', 'capital': 'Valletta'},\n", "{'timezones': ['Indian/Mauritius'], 'code': 'MU', 'continent': 'Africa', 'name': 'Mauritius', 'capital': 'Port Louis'},\n", "{'timezones': ['Indian/Maldives'], 'code': 'MV', 'continent': 'Asia', 'name': 'Maldives', 'capital': 'Mal\\xc3\\xa9'},\n", "{'timezones': ['Africa/Blantyre'], 'code': 'MW', 'continent': 'Africa', 'name': 'Malawi', 'capital': 'Lilongwe'},\n", "{'timezones': ['America/Mexico_City', 'America/Cancun', 'America/Merida', 'America/Monterrey', 'America/Mazatlan', 'America/Chihuahua', 'America/Hermosillo', 'America/Tijuana'], 'code': 'MX', 'continent': 'North America', 'name': 'Mexico', 'capital': 'Mexico City'},\n", "{'timezones': ['Asia/Kuala_Lumpur', 'Asia/Kuching'], 'code': 'MY', 'continent': 'Asia', 'name': 'Malaysia', 'capital': 'Kuala Lumpur'},\n", "{'timezones': ['Africa/Maputo'], 'code': 'MZ', 'continent': 'Africa', 'name': 'Mozambique', 'capital': 'Maputo'},\n", "{'timezones': ['Africa/Windhoek'], 'code': 'NA', 'continent': 'Africa', 'name': 'Namibia', 'capital': 'Windhoek'},\n", "{'timezones': ['Africa/Niamey'], 'code': 'NE', 'continent': 'Africa', 'name': 'Niger', 'capital': 'Niamey'},\n", "{'timezones': ['Africa/Lagos'], 'code': 'NG', 'continent': 'Africa', 'name': 'Nigeria', 'capital': 'Abuja'},\n", "{'timezones': ['America/Managua'], 'code': 'NI', 'continent': 'North America', 'name': 'Nicaragua', 'capital': 'Managua'},\n", "{'timezones': ['Europe/Amsterdam'], 'code': 'NL', 'continent': 'Europe', 'name': 'Kingdom of the Netherlands', 'capital': 'Amsterdam'},\n", "{'timezones': ['Europe/Oslo'], 'code': 'NO', 'continent': 'Europe', 'name': 'Norway', 'capital': 'Oslo'},\n", "{'timezones': ['Asia/Katmandu'], 'code': 'NP', 'continent': 'Asia', 'name': 'Nepal', 'capital': 'Kathmandu'},\n", "{'timezones': ['Pacific/Nauru'], 'code': 'NR', 'continent': 'Oceania', 'name': 'Nauru', 'capital': 'Yaren'},\n", "{'timezones': ['Pacific/Auckland', 'Pacific/Chatham'], 'code': 'NZ', 'continent': 'Oceania', 'name': 'New Zealand', 'capital': 'Wellington'},\n", "{'timezones': ['Asia/Muscat'], 'code': 'OM', 'continent': 'Asia', 'name': 'Oman', 'capital': 'Muscat'},\n", "{'timezones': ['America/Panama'], 'code': 'PA', 'continent': 'North America', 'name': 'Panama', 'capital': 'Panama City'},\n", "{'timezones': ['America/Lima'], 'code': 'PE', 'continent': 'South America', 'name': 'Peru', 'capital': 'Lima'},\n", "{'timezones': ['Pacific/Port_Moresby'], 'code': 'PG', 'continent': 'Oceania', 'name': 'Papua New Guinea', 'capital': 'Port Moresby'},\n", "{'timezones': ['Asia/Manila'], 'code': 'PH', 'continent': 'Asia', 'name': 'Philippines', 'capital': 'Manila'},\n", "{'timezones': ['Asia/Karachi'], 'code': 'PK', 'continent': 'Asia', 'name': 'Pakistan', 'capital': 'Islamabad'},\n", "{'timezones': ['Europe/Warsaw'], 'code': 'PL', 'continent': 'Europe', 'name': 'Poland', 'capital': 'Warsaw'},\n", "{'timezones': ['Europe/Lisbon', 'Atlantic/Madeira', 'Atlantic/Azores'], 'code': 'PT', 'continent': 'Europe', 'name': 'Portugal', 'capital': 'Lisbon'},\n", "{'timezones': ['Pacific/Palau'], 'code': 'PW', 'continent': 'Oceania', 'name': 'Palau', 'capital': 'Ngerulmud'},\n", "{'timezones': ['America/Asuncion'], 'code': 'PY', 'continent': 'South America', 'name': 'Paraguay', 'capital': 'Asunci\\xc3\\xb3n'},\n", "{'timezones': ['Asia/Qatar'], 'code': 'QA', 'continent': 'Asia', 'name': 'Qatar', 'capital': 'Doha'},\n", "{'timezones': ['Europe/Bucharest'], 'code': 'RO', 'continent': 'Europe', 'name': 'Romania', 'capital': 'Bucharest'},\n", "{'timezones': ['Europe/Kaliningrad', 'Europe/Moscow', 'Europe/Volgograd', 'Europe/Samara', 'Asia/Yekaterinburg', 'Asia/Omsk', 'Asia/Novosibirsk', 'Asia/Krasnoyarsk', 'Asia/Irkutsk', 'Asia/Yakutsk', 'Asia/Vladivostok', 'Asia/Sakhalin', 'Asia/Magadan', 'Asia/Kamchatka', 'Asia/Anadyr'], 'code': 'RU', 'continent': 'Europe', 'name': 'Russia', 'capital': 'Moscow'},\n", "{'timezones': ['Africa/Kigali'], 'code': 'RW', 'continent': 'Africa', 'name': 'Rwanda', 'capital': 'Kigali'},\n", "{'timezones': ['Asia/Riyadh'], 'code': 'SA', 'continent': 'Asia', 'name': 'Saudi Arabia', 'capital': 'Riyadh'},\n", "{'timezones': ['Pacific/Guadalcanal'], 'code': 'SB', 'continent': 'Oceania', 'name': 'Solomon Islands', 'capital': 'Honiara'},\n", "{'timezones': ['Indian/Mahe'], 'code': 'SC', 'continent': 'Africa', 'name': 'Seychelles', 'capital': 'Victoria'},\n", "{'timezones': ['Africa/Khartoum'], 'code': 'SD', 'continent': 'Africa', 'name': 'Sudan', 'capital': 'Khartoum'},\n", "{'timezones': ['Europe/Stockholm'], 'code': 'SE', 'continent': 'Europe', 'name': 'Sweden', 'capital': 'Stockholm'},\n", "{'timezones': ['Asia/Singapore'], 'code': 'SG', 'continent': 'Asia', 'name': 'Singapore', 'capital': 'Singapore'},\n", "{'timezones': ['Europe/Ljubljana'], 'code': 'SI', 'continent': 'Europe', 'name': 'Slovenia', 'capital': 'Ljubljana'},\n", "{'timezones': ['Europe/Bratislava'], 'code': 'SK', 'continent': 'Europe', 'name': 'Slovakia', 'capital': 'Bratislava'},\n", "{'timezones': ['Africa/Freetown'], 'code': 'SL', 'continent': 'Africa', 'name': 'Sierra Leone', 'capital': 'Freetown'},\n", "{'timezones': ['Europe/San_Marino'], 'code': 'SM', 'continent': 'Europe', 'name': 'San Marino', 'capital': 'San Marino'},\n", "{'timezones': ['Africa/Dakar'], 'code': 'SN', 'continent': 'Africa', 'name': 'Senegal', 'capital': 'Dakar'},\n", "{'timezones': ['Africa/Mogadishu'], 'code': 'SO', 'continent': 'Africa', 'name': 'Somalia', 'capital': 'Mogadishu'},\n", "{'timezones': ['America/Paramaribo'], 'code': 'SR', 'continent': 'South America', 'name': 'Suriname', 'capital': 'Paramaribo'},\n", "{'timezones': ['Africa/Sao_Tome'], 'code': 'ST', 'continent': 'Africa', 'name': 'S\\xc3\\xa3o Tom\\xc3\\xa9 and Pr\\xc3\\xadncipe', 'capital': 'S\\xc3\\xa3o Tom\\xc3\\xa9'},\n", "{'timezones': ['Asia/Damascus'], 'code': 'SY', 'continent': 'Asia', 'name': 'Syria', 'capital': 'Damascus'},\n", "{'timezones': ['Africa/Lome'], 'code': 'TG', 'continent': 'Africa', 'name': 'Togo', 'capital': 'Lom\\xc3\\xa9'},\n", "{'timezones': ['Asia/Bangkok'], 'code': 'TH', 'continent': 'Asia', 'name': 'Thailand', 'capital': 'Bangkok'},\n", "{'timezones': ['Asia/Dushanbe'], 'code': 'TJ', 'continent': 'Asia', 'name': 'Tajikistan', 'capital': 'Dushanbe'},\n", "{'timezones': ['Asia/Ashgabat'], 'code': 'TM', 'continent': 'Asia', 'name': 'Turkmenistan', 'capital': 'Ashgabat'},\n", "{'timezones': ['Africa/Tunis'], 'code': 'TN', 'continent': 'Africa', 'name': 'Tunisia', 'capital': 'Tunis'},\n", "{'timezones': ['Pacific/Tongatapu'], 'code': 'TO', 'continent': 'Oceania', 'name': 'Tonga', 'capital': 'Nuku\\xca\\xbbalofa'},\n", "{'timezones': ['Europe/Istanbul'], 'code': 'TR', 'continent': 'Asia', 'name': 'Turkey', 'capital': 'Ankara'},\n", "{'timezones': ['America/Port_of_Spain'], 'code': 'TT', 'continent': 'North America', 'name': 'Trinidad and Tobago', 'capital': 'Port of Spain'},\n", "{'timezones': ['Pacific/Funafuti'], 'code': 'TV', 'continent': 'Oceania', 'name': 'Tuvalu', 'capital': 'Funafuti'},\n", "{'timezones': ['Africa/Dar_es_Salaam'], 'code': 'TZ', 'continent': 'Africa', 'name': 'Tanzania', 'capital': 'Dodoma'},\n", "{'timezones': ['Europe/Kiev', 'Europe/Uzhgorod', 'Europe/Zaporozhye', 'Europe/Simferopol'], 'code': 'UA', 'continent': 'Europe', 'name': 'Ukraine', 'capital': 'Kiev'},\n", "{'timezones': ['Africa/Kampala'], 'code': 'UG', 'continent': 'Africa', 'name': 'Uganda', 'capital': 'Kampala'},\n", "{'timezones': ['America/New_York', 'America/Detroit', 'America/Kentucky/Louisville', 'America/Kentucky/Monticello', 'America/Indiana/Indianapolis', 'America/Indiana/Marengo', 'America/Indiana/Knox', 'America/Indiana/Vevay', 'America/Chicago', 'America/Indiana/Vincennes', 'America/Indiana/Petersburg', 'America/Menominee', 'America/North_Dakota/Center', 'America/North_Dakota/New_Salem', 'America/Denver', 'America/Boise', 'America/Shiprock', 'America/Phoenix', 'America/Los_Angeles', 'America/Anchorage', 'America/Juneau', 'America/Yakutat', 'America/Nome', 'America/Adak', 'Pacific/Honolulu'], 'code': 'US', 'continent': 'North America', 'name': 'United States', 'capital': 'Washington, D.C.'},\n", "{'timezones': ['America/Montevideo'], 'code': 'UY', 'continent': 'South America', 'name': 'Uruguay', 'capital': 'Montevideo'},\n", "{'timezones': ['Asia/Samarkand', 'Asia/Tashkent'], 'code': 'UZ', 'continent': 'Asia', 'name': 'Uzbekistan', 'capital': 'Tashkent'},\n", "{'timezones': ['Europe/Vatican'], 'code': 'VA', 'continent': 'Europe', 'name': 'Vatican City', 'capital': 'Vatican City'},\n", "{'timezones': ['America/Caracas'], 'code': 'VE', 'continent': 'South America', 'name': 'Venezuela', 'capital': 'Caracas'},\n", "{'timezones': ['Asia/Saigon'], 'code': 'VN', 'continent': 'Asia', 'name': 'Vietnam', 'capital': 'Hanoi'},\n", "{'timezones': ['Pacific/Efate'], 'code': 'VU', 'continent': 'Oceania', 'name': 'Vanuatu', 'capital': 'Port Vila'},\n", "{'timezones': ['Asia/Aden'], 'code': 'YE', 'continent': 'Asia', 'name': 'Yemen', 'capital': \"Sana'a\"},\n", "{'timezones': ['Africa/Lusaka'], 'code': 'ZM', 'continent': 'Africa', 'name': 'Zambia', 'capital': 'Lusaka'},\n", "{'timezones': ['Africa/Harare'], 'code': 'ZW', 'continent': 'Africa', 'name': 'Zimbabwe', 'capital': 'Harare'},\n", "{'timezones': ['Africa/Algiers'], 'code': 'DZ', 'continent': 'Africa', 'name': 'Algeria', 'capital': 'Algiers'},\n", "{'timezones': ['Europe/Sarajevo'], 'code': 'BA', 'continent': 'Europe', 'name': 'Bosnia and Herzegovina', 'capital': 'Sarajevo'},\n", "{'timezones': ['Asia/Phnom_Penh'], 'code': 'KH', 'continent': 'Asia', 'name': 'Cambodia', 'capital': 'Phnom Penh'},\n", "{'timezones': ['Africa/Bangui'], 'code': 'CF', 'continent': 'Africa', 'name': 'Central African Republic', 'capital': 'Bangui'},\n", "{'timezones': ['Africa/Ndjamena'], 'code': 'TD', 'continent': 'Africa', 'name': 'Chad', 'capital': \"N'Djamena\"},\n", "{'timezones': ['Indian/Comoro'], 'code': 'KM', 'continent': 'Africa', 'name': 'Comoros', 'capital': 'Moroni'},\n", "{'timezones': ['Europe/Zagreb'], 'code': 'HR', 'continent': 'Europe', 'name': 'Croatia', 'capital': 'Zagreb'},\n", "{'timezones': ['Asia/Dili'], 'code': 'TL', 'continent': 'Asia', 'name': 'East Timor', 'capital': 'Dili'},\n", "{'timezones': ['America/El_Salvador'], 'code': 'SV', 'continent': 'North America', 'name': 'El Salvador', 'capital': 'San Salvador'},\n", "{'timezones': ['Africa/Malabo'], 'code': 'GQ', 'continent': 'Africa', 'name': 'Equatorial Guinea', 'capital': 'Malabo'},\n", "{'timezones': ['America/Grenada'], 'code': 'GD', 'continent': 'North America', 'name': 'Grenada', 'capital': \"St. George's\"},\n", "{'timezones': ['Asia/Almaty', 'Asia/Qyzylorda', 'Asia/Aqtobe', 'Asia/Aqtau', 'Asia/Oral'], 'code': 'KZ', 'continent': 'Asia', 'name': 'Kazakhstan', 'capital': 'Astana'},\n", "{'timezones': ['Asia/Vientiane'], 'code': 'LA', 'continent': 'Asia', 'name': 'Laos', 'capital': 'Vientiane'},\n", "{'timezones': ['Pacific/Truk', 'Pacific/Ponape', 'Pacific/Kosrae'], 'code': 'FM', 'continent': 'Oceania', 'name': 'Federated States of Micronesia', 'capital': 'Palikir'},\n", "{'timezones': ['Europe/Chisinau'], 'code': 'MD', 'continent': 'Europe', 'name': 'Moldova', 'capital': 'Chi\\xc5\\x9fin\\xc4\\x83u'},\n", "{'timezones': ['Europe/Monaco'], 'code': 'MC', 'continent': 'Europe', 'name': 'Monaco', 'capital': 'Monaco'},\n", "{'timezones': ['Europe/Podgorica'], 'code': 'ME', 'continent': 'Europe', 'name': 'Montenegro', 'capital': 'Podgorica'},\n", "{'timezones': ['Africa/Casablanca'], 'code': 'MA', 'continent': 'Africa', 'name': 'Morocco', 'capital': 'Rabat'},\n", "{'timezones': ['America/St_Kitts'], 'code': 'KN', 'continent': 'North America', 'name': 'Saint Kitts and Nevis', 'capital': 'Basseterre'},\n", "{'timezones': ['America/St_Lucia'], 'code': 'LC', 'continent': 'North America', 'name': 'Saint Lucia', 'capital': 'Castries'},\n", "{'timezones': ['America/St_Vincent'], 'code': 'VC', 'continent': 'North America', 'name': 'Saint Vincent and the Grenadines', 'capital': 'Kingstown'},\n", "{'timezones': ['Pacific/Apia'], 'code': 'WS', 'continent': 'Oceania', 'name': 'Samoa', 'capital': 'Apia'},\n", "{'timezones': ['Europe/Belgrade'], 'code': 'RS', 'continent': 'Europe', 'name': 'Serbia', 'capital': 'Belgrade'},\n", "{'timezones': ['Africa/Johannesburg'], 'code': 'ZA', 'continent': 'Africa', 'name': 'South Africa', 'capital': 'Pretoria'},\n", "{'timezones': ['Europe/Madrid', 'Africa/Ceuta', 'Atlantic/Canary'], 'code': 'ES', 'continent': 'Europe', 'name': 'Spain', 'capital': 'Madrid'},\n", "{'timezones': ['Asia/Colombo'], 'code': 'LK', 'continent': 'Asia', 'name': 'Sri Lanka', 'capital': 'Sri Jayewardenepura Kotte'},\n", "{'timezones': ['Africa/Mbabane'], 'code': 'SZ', 'continent': 'Africa', 'name': 'Swaziland', 'capital': 'Mbabane'},\n", "{'timezones': ['Europe/Zurich'], 'code': 'CH', 'continent': 'Europe', 'name': 'Switzerland', 'capital': 'Bern'},\n", "{'timezones': ['Asia/Dubai'], 'code': 'AE', 'continent': 'Asia', 'name': 'United Arab Emirates', 'capital': 'Abu Dhabi'},\n", "{'timezones': ['Europe/London'], 'code': 'GB', 'continent': 'Europe', 'name': 'United Kingdom', 'capital': 'London'},\n", "]\n", "\n", "def get_continent(row):\n", " cname = row['country']\n", " result = ''.join(country['continent'] for country in countries if country['name'] == cname)\n", " if result == '':\n", " result = 'Not Found'\n", " return result\n", "\n", "raw_data['continent'] = raw_data.apply(get_continent, axis=1)\n", "\n\n", "# Lets round gdp_per_capita to the nearest round number for EDA purposes\n", "raw_data['gdp_rounded'] = raw_data.apply(lambda row: math.ceil(row['gdp_per_capita ($)']/1000) * 1000, axis = 1)\n", "\n", "raw_data.sample(2)\n" ], "outputs": [ { "output_type": "execute_result", "execution_count": 2, "data": { "text/plain": [ " country year sex age suicides_no population \\\n", "4811 Bulgaria 2007 male 35-54 years 180 1061001 \n", "13472 Japan 1994 male 75+ years 1364 2476000 \n", "\n", " suicides/100k pop country-year HDI for year gdp_for_year ($) \\\n", "4811 16.97 Bulgaria2007 NaN 44,411,476,557 \n", "13472 55.09 Japan1994 NaN 4,907,039,384,470 \n", "\n", " gdp_per_capita ($) generation continent gdp_rounded \n", "4811 6073 Boomers Europe 7000 \n", "13472 41563 G.I. Generation Asia 42000 " ], "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", "
countryyearsexagesuicides_nopopulationsuicides/100k popcountry-yearHDI for yeargdp_for_year ($)gdp_per_capita ($)generationcontinentgdp_rounded
4811Bulgaria2007male35-54 years180106100116.97Bulgaria2007NaN44,411,476,5576073BoomersEurope7000
13472Japan1994male75+ years1364247600055.09Japan1994NaN4,907,039,384,47041563G.I. GenerationAsia42000
\n", "
" ] }, "metadata": {} } ], "execution_count": 2, "metadata": { "collapsed": false, "outputHidden": false, "inputHidden": false } }, { "cell_type": "code", "source": [ "countrywise_df = raw_data.groupby('country')['suicides/100k pop'].agg('sum').reset_index()\n", "yearwise_df = raw_data.groupby('year')['suicides/100k pop'].agg('sum').reset_index()\n", "genderwise_df = raw_data.groupby('sex')['suicides/100k pop'].agg('sum').reset_index()\n", "agewise_df = raw_data.groupby('age')['suicides/100k pop'].agg('sum').reset_index()\n", "gdpwise_df = raw_data.groupby('gdp_rounded')['suicides/100k pop'].agg('sum').reset_index()\n", "continent_df = raw_data.groupby('continent')['suicides/100k pop'].agg('sum').reset_index()\n", "continent_gdp_df = raw_data.groupby('continent')['gdp_rounded'].agg('sum').reset_index()" ], "outputs": [], "execution_count": 3, "metadata": { "collapsed": false, "outputHidden": false, "inputHidden": false } }, { "cell_type": "code", "source": [ "def show_plot(x, y, title, xlabel, ylabel, plottype='Bar'):\n", "\n", " if plottype != 'Bar':\n", " trace0 = go.Scatter(x=x,y=y)\n", " else: \n", " trace0 = go.Bar(x=x,y=y)\n", " \n", " data = [trace0]\n", "\n", " layout = go.Layout(\n", " title=title,\n", " xaxis=dict(\n", " title=xlabel,\n", " titlefont=dict(\n", " family='Courier New, monospace',\n", " size=18,\n", " color='#7f7f7f'\n", " )\n", " ),\n", " yaxis=dict(\n", " title=ylabel,\n", " titlefont=dict(\n", " family='Courier New, monospace',\n", " size=18,\n", " color='#7f7f7f'\n", " )\n", " )\n", " )\n", "\n", " fig = go.Figure(data=data, layout=layout)\n", " return fig\n" ], "outputs": [], "execution_count": 4, "metadata": { "collapsed": false, "outputHidden": false, "inputHidden": false } }, { "cell_type": "code", "source": [ "py.iplot(show_plot(countrywise_df['country'],countrywise_df['suicides/100k pop'], 'Suicide rates across countries','Country', 'Suicide Rate'))" ], "outputs": [ { "output_type": "stream", "name": "stderr", "text": [ "/anaconda3/lib/python3.7/site-packages/IPython/core/display.py:689: UserWarning:\n", "\n", "Consider using IPython.display.IFrame instead\n", "\n" ] }, { "output_type": "execute_result", "execution_count": 5, "data": { "text/plain": [ "" ], "text/html": [ "" ] }, "metadata": {} } ], "execution_count": 5, "metadata": { "collapsed": false, "outputHidden": false, "inputHidden": false } }, { "cell_type": "code", "source": [ "py.iplot(show_plot( continent_df['continent'],continent_df['suicides/100k pop'], 'Suicide rates across continents','Continent', 'Suicide Rate'))" ], "outputs": [ { "output_type": "execute_result", "execution_count": 6, "data": { "text/plain": [ "" ], "text/html": [ "" ] }, "metadata": {} } ], "execution_count": 6, "metadata": { "collapsed": false, "outputHidden": false, "inputHidden": false } }, { "cell_type": "code", "source": [ "py.iplot(show_plot(agewise_df['age'],agewise_df['suicides/100k pop'], 'Suicide rates by age group','Age group', 'Suicide Rate'))" ], "outputs": [ { "output_type": "execute_result", "execution_count": 7, "data": { "text/plain": [ "" ], "text/html": [ "" ] }, "metadata": {} } ], "execution_count": 7, "metadata": { "collapsed": false, "outputHidden": false, "inputHidden": false } }, { "cell_type": "code", "source": [ "py.iplot(show_plot( yearwise_df['year'],yearwise_df['suicides/100k pop'], 'Suicide rates over the years','Years', 'Suicide Rate', 'Scatter'))" ], "outputs": [ { "output_type": "execute_result", "execution_count": 8, "data": { "text/plain": [ "" ], "text/html": [ "" ] }, "metadata": {} } ], "execution_count": 8, "metadata": { "collapsed": false, "outputHidden": false, "inputHidden": false } }, { "cell_type": "code", "source": [ "py.iplot(show_plot(gdpwise_df['gdp_rounded'],gdpwise_df['suicides/100k pop'], 'Suicide rates by GDP-per-capita','GDP per capita', 'Suicide Rate', 'Scatter'))" ], "outputs": [ { "output_type": "execute_result", "execution_count": 9, "data": { "text/plain": [ "" ], "text/html": [ "" ] }, "metadata": {} } ], "execution_count": 9, "metadata": { "collapsed": false, "outputHidden": false, "inputHidden": false } }, { "cell_type": "markdown", "source": [ "while it looks like Suicide rate goes down with increase in GDP per capita, taking one of countries with top suicidal rates paints a different picture. Afterall increase in GDP percapita doesn't affect the suicide as above trend line shows." ], "metadata": {} }, { "cell_type": "code", "source": [ "Lithuania_df = raw_data[raw_data['country'] == 'Lithuania']\n", "Lithuania_df = Lithuania_df.groupby('gdp_rounded')['suicides/100k pop'].agg('sum').reset_index()\n", "py.iplot(show_plot(Lithuania_df['gdp_rounded'],Lithuania_df['suicides/100k pop'], 'Suicide rates of Lithuania by GDP-per-capita','GDP per capita', 'Suicide Rate', 'Scatter'))" ], "outputs": [ { "output_type": "execute_result", "execution_count": 10, "data": { "text/plain": [ "" ], "text/html": [ "" ] }, "metadata": {} } ], "execution_count": 10, "metadata": { "collapsed": false, "outputHidden": false, "inputHidden": false } }, { "cell_type": "code", "source": [ "data = [ dict(\n", " type = 'choropleth',\n", " locations = countrywise_df['country'],\n", " z = countrywise_df['suicides/100k pop'],\n", " text = countrywise_df['country'],\n", " locationmode = \"country names\",\n", " colorscale = [[0,\"rgb(5, 10, 172)\"],[0.35,\"rgb(40, 60, 190)\"],[0.5,\"rgb(70, 100, 245)\"],\\\n", " [0.6,\"rgb(90, 120, 245)\"],[0.7,\"rgb(106, 137, 247)\"],[1,\"rgb(220, 220, 220)\"]],\n", " autocolorscale = False,\n", " reversescale = True,\n", " marker = dict(\n", " line = dict (\n", " color = 'rgb(180,180,180)',\n", " width = 0.5\n", " ) ),\n", " colorbar = dict(\n", " autotick = False,\n", " tickprefix = '',\n", " title = 'Suicide rate'),\n", " ) ]\n", "\n", "layout = dict(\n", " title = 'Country wise heatmap of suicide rates',\n", " geo = dict(\n", " showframe = False,\n", " showcoastlines = False,\n", " projection = dict(\n", " type = 'Mercator'\n", " )\n", " )\n", ")\n", "\n", "fig = dict( data=data, layout=layout )\n", "py.iplot( fig, validate=False, filename='d3-world-map' )" ], "outputs": [ { "output_type": "execute_result", "execution_count": 11, "data": { "text/plain": [ "" ], "text/html": [ "" ] }, "metadata": {} } ], "execution_count": 11, "metadata": { "collapsed": false, "outputHidden": false, "inputHidden": false } }, { "cell_type": "markdown", "source": [ "Let's encode all the features appropriately and checkout the correlation." ], "metadata": { "collapsed": false, "outputHidden": false, "inputHidden": false } }, { "cell_type": "code", "source": [ "from sklearn_pandas import DataFrameMapper\n", "from sklearn.preprocessing import LabelEncoder\n", "from sklearn.preprocessing import StandardScaler\n", "\n", "mapper = DataFrameMapper([\n", " ('country', LabelEncoder()),\n", " ('year' , LabelEncoder()), \n", " ('sex', LabelEncoder()),\n", " ('age', LabelEncoder()),\n", " ('continent', LabelEncoder()),\n", " (['gdp_rounded'], StandardScaler()),\n", " (['suicides/100k pop'], StandardScaler())\n", "],df_out=True)\n", "\nout_df = mapper.fit_transform(raw_data.copy())" ], "outputs": [ { "output_type": "stream", "name": "stderr", "text": [ "/anaconda3/lib/python3.7/site-packages/sklearn/utils/validation.py:595: DataConversionWarning:\n", "\n", "Data with input dtype int64 was converted to float64 by StandardScaler.\n", "\n", "/anaconda3/lib/python3.7/site-packages/sklearn/utils/validation.py:595: DataConversionWarning:\n", "\n", "Data with input dtype int64 was converted to float64 by StandardScaler.\n", "\n" ] } ], "execution_count": 12, "metadata": { "collapsed": false, "outputHidden": false, "inputHidden": false } }, { "cell_type": "code", "source": [ "import seaborn as sns\n", "corr = out_df.corr()\n", "sns.heatmap(corr, \n", " xticklabels=corr.columns.values,\n", " yticklabels=corr.columns.values)" ], "outputs": [ { "output_type": "execute_result", "execution_count": 14, "data": { "text/plain": [ "" ] }, "metadata": {} }, { "output_type": "display_data", "data": { "text/plain": [ "
" ], "image/png": [ "iVBORw0KGgoAAAANSUhEUgAAAa0AAAFNCAYAAACzEjC1AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJzt3XucHmV9///XOwGkEAhyECHhVI0ociYIAUSONagICFoQqhx+pvVbDNbW/qBSPBWr1mo5fSlRQxBFqwiSIgiIBOQowUCAKBoplQgFQUSQctjd9/ePmd3cWe7dvZdsdmZu3k8e89h7Zq575nPvkv3sdc11kG0iIiKaYELVAURERHQqSSsiIhojSSsiIhojSSsiIhojSSsiIhojSSsiIhojSSsiIlYJSXMlPSrpniHOS9KZkpZKWixp55GumaQVERGryjxg5jDnDwKmldss4NyRLpikFRERq4TtG4DfDVPkEOBrLtwKrCdpk+GumaQVERFVmQI82LK/rDw2pNVWaTix0l547P5azrO13w4fqDqEtp71C1WH0DiLHvtV1SG0tf/G21cdQuNc9eCVWtlrjOZ3zhobveYvKZr1+s2xPWcUt2sX77D3T9KKiIjl+no7LlomqNEkqcGWAZu17E8FHhruDWkejIiI5dzX+bby5gPvK3sR7g48afvh4d6QmlZERCzXNybJCABJ3wT2ATaUtAz4OLA6gO1/B64A3gYsBZ4BjhvpmklaERExwGNTgyqv5aNGOG/gr0dzzSStiIhYbgxrWqtCklZERCw3hjWtVSFJKyIiluut97CRJK2IiFguzYMREdEUY9kRY1VI0oqIiOVS04qIiMaoeU0rM2KMMUkflrRW1XFERLwkvS90vlUgSWvsfRhom7QkTRznWCIiRqevr/OtAi/LpCXpfeUqmXdJulDSFpKuLY9dK2nzstw8SUe0vO/p8us+khZIuljSzyV9o5w7azawKXCdpOv63yPpU5JuA06VdGnL9Q6UdMm4fviIiOGM79yDo/aye6Yl6Y3Ax4A9bT8maX3gAoqFyC6QdDxwJnDoCJfaCXgjxYzEN5XXO1PSR4B9bT9WllsbuMf2aZIE/EzSRrZ/SzHP1vlj/iEjIl6qmnfEeDnWtPYDLu5PKrZ/B8wALirPXwjs1cF1fmJ7mYv+oXcCWw5Rrhf4bnkvl9c/RtJ65X2vHPwGSbMkLZS08Ctf+2bHHywiYmXZvR1vVXjZ1bQoFh0baZGz/vM9lIm9rCWt0VLmuZbXvQz9vXzWK/50zwf+E3gW+I7tnhfdvGWNmrouAhkRXSq9B2vnWuA9kjYAKJsHbwaOLM8fDdxYvn4A2KV8fQjllPojeApYZ6iTth+iaFI8FZg3utAjIlax3p7Otwq87Gpatu+VdDpwvaReYBEwG5gr6aNA/7MmgC8Dl0n6CUWy+2MHt5gDXCnpYdv7DlHmG8BGtpeszGeJiBhzo1i5uAovu6QFYPsCis4XrfZrU+4RYPeWQ6eUxxcAC1rKndjy+izgrJb9SW1C2IsiIUZE1EvNmwdflkmrSpLuoKix/W3VsUREvEjNew8maY0z27uMXCoioiKpaUVERGOkphUREU3hLAIZERGNkZpWREQ0Rp5pRUREY6SmFRERjZGaVkRENEZF0zN1KkkrIiKWS/NgrIz9dvhA1SG09aO76jkL1R7bH1t1CI2z0VqTqw6hrQ++8MqqQxjSKX1Lqw5h1UnSioiIxsgzrYiIaIya17RejutpRUTEUNzX+dYBSTMl3SdpqaST25zfXNJ1khZJWizpbcNdLzWtiIhYbgx7D0qaCJwDHAgsA26XNH/QWoKnAt+2fa6kbYArgC2HumaSVkRELDe2zYNvApbavh9A0rcoVoFvTVoG1i1fT6ZY2X1ISVoREbHc2CatKcCDLfvLgN0GlfkEcLWkDwFrAwcMd8E804qIiOXsjjdJsyQtbNlmDbqa2t1h0P5RwDzbU4G3ARdKGjI3paYVERHLjaKmZXsOMGeYIsuAzVr2p/Li5r8TgJnl9W6RtCawIfBouwumphUREcv19XW+jex2YJqkrSStARwJzB9U5tfA/gCS3gCsCfx2qAumphUREcuNYe9B2z2STgSuAiYCc23fK+lTwELb84G/Bb4s6W8omg6PtT24CXFAklZERCw3dL54iZfzFRTd2FuPndbyegmwZ6fXS9KKiIjlaj4jRpJWxSRNtN1bdRwREUDtk1Y6YoyCpE9LOqll/3RJsyV9VNLt5RQkn2w5/z1Jd0i6t7UrqKSnJX1K0m3AjHH+GBERQxvjaZzGWpLW6HwVeD9AOY7gSOARYBrFyO8dgV0k7V2WP972LsB0YLakDcrjawP32N7N9o3j+QEiIobjnt6OtyqkeXAUbD8g6XFJOwEbA4uAXYE/K18DTKJIYjdQJKrDyuOblccfB3qB7w51n7JWNgvgtZO35tVrT1kFnyYioo0sTdJ1vgIcC7wamEsxvuCfbZ/XWkjSPhTTkcyw/YykBRTjDwCeHe45VuuAvTdP2X9su/JERAynr96/ctI8OHqXUoze3pVi7MFVwPGSJgFImiLpVRQTPz5RJqzXA7tXFXBERMfGdnDxmEtNa5RsPy/pOuD3ZW3p6nIU9y2SAJ4GjgF+APyVpMXAfcCtVcUcEdGxmvceTNIapbIDxu7Au/uP2T4DOKNN8YPaXcP2pFUTXUTEShrjwcVjLUlrFMoFyi4HLrX9y6rjiYgYcxX1CuxUktYolNON/GnVcURErDLpPRgREY1R896DSVoRETHA6YgRERGNkZpWREQ0Rp5pRUREY6T3YERENEaaByMiojHSPBgREY2RmlasjGf9QtUhtLXH9sdWHUJbNy+eV3UIQ6rr92zqWhtWHUJbp7Os6hCGdB5bVB3CKpMu7xER0Rw9SVoREdEUeaYVERGNkWdaERHRFE7SioiIxkjSioiIxkjvwYiIaIz0HoyIiKaw6908OKHqACIiokb63PnWAUkzJd0naamkk4co8x5JSyTdK+mi4a6XmlZERCw3hh0xJE0EzgEOBJYBt0uab3tJS5lpwCnAnrafkPSq4a6ZpBUREQPGuMv7m4Cltu8HkPQt4BBgSUuZDwDn2H4CwPajw10wzYMREbHc2DYPTgEebNlfVh5r9TrgdZJuknSrpJnDXTA1rYiIGOCezmtakmYBs1oOzbE9p7VIu1sM2l8NmAbsA0wFfixpW9u/b3fPJK0xIGlt4NsU3/CJwKeBpcAXgUnAY8CxwG+BW4CP2l4g6Z+BPtsfqyLuiIgXGUXzYJmg5gxTZBmwWcv+VOChNmVutf0C8F+S7qNIYre3u2CaB8fGTOAh2zvY3hb4AXAWcITtXYC5wOm2eyiS17mSDizf98mKYo6IeLG+UWwjux2YJmkrSWsARwLzB5X5HrAvgKQNKZoL7x/qgqlpjY27gS9I+hxwOfAEsC1wjSQoal8PA9i+V9KFwH8CM2w/P/hirVXuLSa/lo3W2mRcPkRExFh2xLDdI+lE4CqK34Nzy9+BnwIW2p5fnvszSUuAXoqWqMeHumaS1hiw/QtJuwBvA/4ZuAa41/aMId6yHfB7YOMhrjdQ5d51073rPdIvIrrLGE+IYfsK4IpBx05reW3gI+U2ojQPjgFJmwLP2P468AVgN2AjSTPK86tLemP5+l3ABsDewJmS1qso7IiIF3GPO96qkJrW2NgO+BdJfcALwAeBHoqkNJni+/xvkh4BPgvsb/tBSWcDZwDvryjuiIgV1HwNyCStsWD7Kop22cH2bnPsdS3vO3OVBRUR8VIkaUVERFOkphUREc2RpBUREU2RmlZERDRGX0/VEQwvSSsiIpZzu+kC6yNJKyIiBqR5MCIiGsN9qWlFRERDpKYVERGN0debmlZERDREmgcjxtEe2x9bdQhDunnxvKpDaKvO37O6+v6a9fzV2W7euNFyzdeVqOd3PiIiKpGaVkRENEaSVkRENEaaByMiojH6euu9NnCSVkREDMg4rYiIaIy+zD0YERFN4SStiIhoivQejIiIxkjvwYiIaIze9B6MiIimyDOtiIhojDQPRkREY9S9y3u9Gy8bQNL3JN0h6V5Js8pjJ0j6haQFkr4s6ezy+EaSvivp9nLbs9roIyJWZKvjrQpJWivveNu7ANOB2ZKmAP8I7A4cCLy+pewZwJds7wocDnxlvIONiBhOb5863johaaak+yQtlXTyMOWOkGRJ04e7XpoHV95sSYeVrzcD/gK43vbvACR9B3hdef4AYBtp4Ie9rqR1bD/VesGyxjYLYIvJr2WjtTZZxR8hIqIwljUoSROBcyj+gF8G3C5pvu0lg8qtA8wGbhvpmqlprQRJ+1Akohm2dwAWAfcN85YJZdkdy23K4IQFYHuO7em2pydhRcR46rM63jrwJmCp7fttPw98CzikTblPA58Hnh3pgklaK2cy8ITtZyS9nqJJcC3gLZJeKWk1imbAflcDJ/bvSNpxXKONiBiBR7F1YArwYMv+svLYAEk7AZvZvryTCyZprZwfAKtJWkzxl8KtwG+Az1BUc38ILAGeLMvPBqZLWixpCfBX4x9yRMTQRlPTkjRL0sKWbdagy7Wrjg3kO0kTgC8Bf9tpfHmmtRJsPwccNPi4pIW255Q1rUspaljYfgz48/GNMiKic6N5pmV7DjBnmCLLKJ7195sKPNSyvw6wLbCgfNb/amC+pHfaXtjugklaq8YnJB0ArEmRsL5XcTwRER3pbVs5esluB6ZJ2oqiFepI4L39J20/CWzYvy9pAfB3QyUsSNJaJWz/XdUxRES8FH1jOCOG7R5JJwJXAROBubbvlfQpYKHt+aO9ZpJWREQM6Bvbmha2rwCuGHTstCHK7jPS9ZK0IiJigMc4aY21JK2IiBjQV3UAI0jSioiIAalpRUREY/RUHcAIkrQiImJAaloREdEYHU7eXpkkrYiIGDDWXd7HWpJWzS167FdVh9DWRmtNrjqEtqauteHIhSqyx/bHVh1CWzcvnld1CG2ttembqw5hSD3r91YdQlufG4NrjOHY4lUiSSsiIgb0KDWtiIhoiNS0IiKiMTK4OCIiGiO9ByMiojHSezAiIhojz7QiIqIxeupd0UrSioiI5VLTioiIxkhHjIiIaIx0eY+IiMaoe9KaUHUAdSNpS0nvbdmfLunMVXCfQyVtM9bXjYhYGVbnWxWStF5sS2AgadleaHv2KrjPoUCSVkTUSs8otip0XdKS9D5JiyXdJelCSVtIurY8dq2kzcty8ySdKelmSfdLOqK8xGeBN0u6U9LfSNpH0uXlez4haa6kBeV7Zrfc9xhJPynfd56kieXxpyWdXsZzq6SNJe0BvBP4l7L8a8b3uxQR0Z5HsVWhq5KWpDcCHwP2s70DcBJwNvA129sD3wBam/o2AfYC3kGRrABOBn5se0fbX2pzm9cDbwXeBHxc0uqS3gD8ObCn7R2BXuDosvzawK1lPDcAH7B9MzAf+Gh5n3quPxIRLzt96nyrQrd1xNgPuNj2YwC2fydpBvCu8vyFwOdbyn/Pdh+wRNLGHd7j+7afA56T9CiwMbA/sAtwu4pp/f8EeLQs/zxwefn6DuDAkW4gaRYwC0ATJzNhwtodhhYRsXLq3hGj25KWGLnW2nr+uUHv7UTre3opvocCLrB9SpvyL9j2oPLDB2jPAeYArLbGlLqP9YuILlL3pNVVzYPAtcB7JG0AIGl94GbgyPL80cCNI1zjKWCdl3DfIyS9qv++krZYBfeJiFiletX5VoWuSlq27wVOB66XdBfwRWA2cJykxcBfUDznGs5ioKfsOPE3Hd53CXAqcHV5n2sonpcN51vARyUtSkeMiKiLvlFsVei25kFsXwBcMOjwfm3KHTtof1L59QWKZ1StFpTnPjHoPdu2vP4P4D/a3GdSy+uLgYvL1zeRLu8RUTN1fx7RdUkrIiJeur6ap62uah6MiIiVM9bNg5JmSrpP0lJJJ7c5/xFJS1rG0g7bHyBJKyIiBozl4OJykoVzgIMoHocc1Wb6ukXA9HIs7cWsOCzpRZK0IiJiQI863zrwJmCp7fttP0/RAe2Q1gK2r7P9TLl7KzB1uAvmmVZERAwY42daU4AHW/aXAbsNU/4E4MrhLpikFRERA0aTslpn7ynNKSdHGCjS6S0kHQNMB94y3D2TtCIiYsBoxl+1zt4zhGXAZi37U4GHBheSdADFvLFvKafJG1KSVkREDBjj5sHbgWmStgJ+QzE70XtbC0jaCTgPmGn70RdfYkVJWhERMaB3DK9lu0fSicBVwERgru17JX0KWGh7PvAvwCTgO+WE47+2/c6hrpmkFRERA8Z6cLHtK4ArBh07reX1AaO5XpJWze2/8fZVh9DWB194ZdUhtHU6y6oOoXHW2vTNVYfQ1jMP/bjqEIa05/bHVR3CKlPv+TCStCIiokXdlyZJ0oqIiAGueV0rSSsiIgakphUREY3Rm5pWREQ0Rd2XJknSioiIAWkejIiIxkhHjIiIaIzUtCIiojFS04qIiMbocZJWREQ0RL1TVpJWRES0qHuX9wlVB9AJSVtKuqfqOIYjaZ6kI0ZRvvafKSJefjyK/6rwsqlpqVioRbbr3jkmIqIydf8FWYukJekfgaOBB4HHgDuA64C5wDPAjS1ljwUOA14BbAVcZPuTQ1x3S+DK8lozgEMl7QH8AyDg+7b//7Ls07Ynla+PAN5h+1hJ84A/ANOBVwN/b/viMgmeBewH/Fd5vf777gJ8kWJhs8eAY20/XB5/0WeKiKiL3pqnrcqbByVNBw4HdgLeRZEcAM4HZtue0eZtb6JIcjsC7y6vMZStga/Z3gl4AfgcRaLZEdhV0qEdhLkJsBfwDuCz5bHDymtvB3wA2KP8PKtTJLMjbPcnqdM7+EwDJM2StFDSwmVPP9hBeBERY6NvFFsVKk9aFMngMtv/a/sp4D+BtYH1bF9flrlw0Huusf247f8FLimvMZT/tn1r+XpXYIHt39ruAb4B7N1BjN+z3Wd7CbBxeWxv4Ju2e20/BPyoPL41sC1wjaQ7gVOBqZImj/CZBtieY3u67elTJ23WQXgREWPDdsdbFerQPKg2x/7I8D0vB58bruwfR7hXu2usOejcc0Nco919Bdw7uDYlab0R4oyIqFx6D47sRuBgSWtKmgS8vTz+pKT+GtTRg95zoKT1Jf0JcChwU4f3ug14i6QNJU0EjgL6az6PSHqDpAkUTX8juQE4UtJESZsA+5bH7wM2kjQDiuZCSW+0/fsRPlNEROXq3jxYeU3L9u2S5gN3Af8NLASeBI4D5kp6Brhq0NtupGheey1FR4yFHd7rYUmnUHTMEHCF7cvK0ycDl1N0BrmHohPFcC6leDZ2N/ALyuRn+/myI8eZZZPgasC/AfeO8JkiIiqXaZw68wXbn5C0FkUN5l9t/xTYoaXMJ1peP2r7xJEuavsBiudLrccuAi5qU/Zi4OI2x48dtD+p/GqgbQy276TNszLbdzD0Z4qIqFxvzUcF1SVpzZG0DcWzpAvKhBUREeOs3imrJknL9ntHUXYeMK/1mKQNgGvbFN/f9uMrFVxExMtImgfHQZmYdqw6joiIpqt778GuSFoRETE2qhp/1akkrYiIGJCaVkRENEbdew/WYXBxRETUhEexdULSTEn3SVoq6eQ2518h6T/K87eVE50PKUkrIiIG9OGOt5GUMw+dAxwEbAMcVQ5vanUC8ITt1wJfopjUfEhJWhERMWAskxbFihxLbd9v+3ngW8Ahg8ocAlxQvr4Y2L9c+qmtPNOKl+SUvqVVh9DWeWxRdQhD+v6a9fzn1rN+b9UhtLXn9sdVHcKQblp8ftUhrDJj3HtwCsXUeP2WAbsNVcZ2j6QngQ0o1iJ8kXr+K4qIiEqMZhFISbOAWS2H5tie01qkzdsGZ8VOygxI0oqIiAGjqWmVCWrOMEWWAa2LAk4FHhqizDJJqwGTgd8NdcE804qIiAFj/EzrdmCapK0krQEcCcwfVGY+8P7y9RHAjzxM5kxNKyIiBozlM63yGdWJFEsxTQTm2r5X0qeAhbbnA18FLpS0lKKGdeRw10zSioiIAWM9I4btK4ArBh07reX1s8C7O71eklZERAzILO8REdEYdZ/GKUkrIiIG9GWW94iIaIo0D0ZERGOkphUREY1R95rWmA4ulvTOdlPPDypz8xDH50k64iXed3VJd5Sv50p6VNI9g8qsL+kaSb8sv76yPC5JZ5bT4i+WtHN5fB9Jl7+UeCIimqrP7nirwpgmLdvzbX92hDJ7jOU9S3sB/clwHjCzTZmTgWttTwOuLfehmDJ/WrnNAs5dBfFFRDRCn3s73qowbNKStLak70u6S9I9kv68PP6ApA3L19MlLShfHyvp7PL1xpIuLd97l6Q9yuNPl18l6WxJSyR9H3hVy313kXS9pDskXSVpk/L47LL8Yknfagl1JnAlgO0baD9vVev09xcAh7Yc/5oLtwLr9d+vJZ5dJS2S9KeDjh8r6TJJPygXOft4y7mPlN+zeyR9uDy2paSfS7qg/AwXS1pruJ9BRMR4GuNpnMbcSM+0ZgIP2X47gKTJo7j2mcD1tg8rFwKbNOj8YcDWwHbAxsASYK6k1YGzgENs/7ZMlKcDx1PUjray/Zyk9VqutS/wyRHi2dj2wwC2H5bUnyTbTZ0/pX+nTLb98fy6zXXfBGwLPAPcXiZgA8dRTMEv4DZJ1wNPlJ/5BNs3SZoL/B/gCyPEHhExLsZ4aZIxN1Lz4N3AAZI+J+nNtp8cxbX3o2xqs93b5r17A98szz0E/Kg8vjVFErhG0p3AqRQzAwMsBr4h6RigB0DSpsDvbD8zithaDTct/hsoZjA+eIiEBXCN7cdt/y9wCUVT5V7Apbb/aPvp8viby/IP2r6pfP31suyKAUmzJC2UtHDZ0w8OPh0RscrUvaY1bNKy/QtgF4rk9c+S+ueL6ml575orcf92n1rAvbZ3LLftbP9Zee7tFEs37wLcUU5jfxDFZIwjeaSlmXET4NHy+HBT5z8MPAvsNIrPYNonwuHKr3jAnmN7uu3pUydtNvh0RMQqY7vjrQojPdPaFHjG9tcpmrB2Lk89QJE4AA4f4u3XAh8srzNR0rqDzt8AHFme24SiiQ/gPmAjSTPK964u6Y2SJgCb2b4O+HtgPYomx4HnWSNonf7+/cBlLcffVz5j2x14sr8ZEfg9RaL8jKR9hrjugWXPxD+heE52U/nZDpW0lqS1KZpCf1yW37z/swFHATd2EHtExLjodV/HWxVGah7cDvhJ2Uz3MeCfyuOfBM6Q9GNgqC4kJwH7SrobuAN446DzlwK/pKjFnQtcD2D7eYo1VT4n6S7gTmAPimntv15ebxHwJeApYJrtn/dfVNI3gVuArSUtk3RCeeqzFAnml8CB5T4Usw/fDywFvkzxjGmA7UeAg4FzJA1eJhqKpHNhGed3bS+0/VOKXow/AW4DvmJ7UVn+Z8D7JS0G1ie9FSOiRupe01LdH7oNR9JewDG2/6qi+x8LTLd9YofltwQut71tp/d462YH1fIH9OvnHq86hLbOm7BF1SEM6ftr1nMs/9XPDvW4tlqrT6jn9wvgpsXnVx1CW6tv+KfDPZroyEaTt+74d85vn7xvpe83WvX9v6IDtm8kzWsREWOm7hWZRietqtmeR9EM2Gn5Byh6RkZE1FLmHoyIiMZITSsiIhoji0BGRERjpHkwIiIao+5LkyRpRUTEgNS0IiKiMdIRIyIiGqMvHTEiIqIpUtOKiIjGqHfKavjcgzF6kmbZnlN1HIMlrtGpa1xQ39gSV3cYaZb36D6zqg5gCIlrdOoaF9Q3tsTVBZK0IiKiMZK0IiKiMZK0Xn7q2naeuEanrnFBfWNLXF0gHTEiIqIxUtOKiIjGSNKKiIjGSNKKcSdpgqQ9qo6jqSStXXUMrSRt1cmxKEhaU9JHJF0i6buS/kbSmlXH1RRJWl2u/Efxdkm1+Vnb7gP+teo4hiLp05JWa9lfV9L5VcZUxrGHpCXAz8r9HST934rDAvhum2MXj3sULSTdLWnxUFuVsQFfA94InAWcDbwBuLDSiBok0zh1v3OB44AzJX0HmGf75xXHBHC1pMOBS1y/3kCrAbdJOg54NcUvl7OqDQmALwFvBeYD2L5L0t5VBSPp9RS/fCdLelfLqXWBqmsO7yi//nX5tT8pHA08M/7hrGBr2zu07F8n6a7KommYJK0uZ/uHwA8lTQaOAq6R9CDwZeDrtl+oKLSPAGsDPZKeBVSE63UrimeA7VMkXQvcBjwB7G17acVhAWD7QUmth3qrigXYmiI5rAcc3HL8KeADlURUsv3fAJL2tL1ny6mTJd0EfKqayABYJGl327cCSNoNuKnCeBolSetlQNIGwDHAXwCLgG8AewHvB/apIibb61Rx306UtZczKH6xbQecLel42w9VGxkPls8CLWkNYDZlU2EVbF8GXCZphu1bqopjBGtL2sv2jVA0sVL8sVSl3YD3Sfp1ub858DNJd1P84bZ9daHVX8ZpdTlJlwCvp2gemWf74ZZzC21PrzC2VwLTaGlKsn1DVfH0k/QT4FjbS8r9dwGfsf36iuPakCKZHkBRM70aOMn24xXHtRFFzWpLWv4Qtn18VTH1k7QLMBeYTDGB+ZPA8bZ/WmFMWwx3vr+WGO0laXWxsvPFqbarbAppS9L/B5wETAXuBHYHbrG9X6WBAZIm2u4ddGyDqpNDXUm6GfgxcActzZW223XQqISkdSl+3z1ZdSxQdKIB3lzu/th2nml1KEmry0m6xfaMquMYrGwK2RW41faO5UP9T9r+84pDQ9LGwGeAKbZnStoGmGH7qxXHdWabw08CC8umukpIutP2jlXdfzgtP8tNbR9Uh5+lpJMoaqaXlIcOA+bYrkNnn9qrTTfoWGWulnS4Bj29r4FnbT8LIOkVZY/GrSuOqd884Cpgk3L/F8CHK4tmuTWBHYFfltv2wPrACZL+rcK4Lpf0tgrvP5x5FD/LTcv9OvwsTwB2s32a7dMoWhkq7bjSJOmI0f3q2ktvmaT1gO9R9Gh8Aqi6o0O/DW1/W9IpALZ7JFXZS6/fa4H9bPcASDqX4rnWgcDdFcZ1EvAPkp4Hnqc+/49BPX+WYsXrAc1yAAAM70lEQVRen73lsehAklaXq2svPduHlS8/Iek6igflP6gwpFZ/LHtcGkDS7hTNcFWbQvEHSH8sa1M0e/VKeq6qoOr6/1ipjj/L8ynGAV5KkawOASptem6SJK0uJ+la2/uPdKwKkvYCptk+v+yBNgX4r4rDgqJ2Oh94TTmmZyPgiGpDAuDzwJ2SFlD8stsb+Ew5rdMPqwqqbHo+GtjK9qclbQZsYvsnVcXUonY/S9tfLH+Ge5WHjrO9qMKQGiVJq0uVc5mtBWxYdi3vb35Yl+Xt+5WR9HFgOsVzrPOB1YGvA3sO975x8hrgIGAz4HCKcTWV/1ux/VVJV1KMt/s5RdPgMtt/BD5aYWj/F+gD9gM+DTwNnEPR0aZStn8q6S0U/58JuK/CAfWDieL7lqbBUaj8H2KsMn9J8cB5U4quyP3/MP5A8QulaocBOwE/BbD9kKS6NDP9o+3vlMn+AIp5Es+lSF6VGWqYAEWyqNJutneWtAjA9hPl4OfKDJpWqtXrJGH7kiHOr3KSTgPeTTFno4DzJX3H9j9VFVOTJGl1KdtnAGdI+lBNu9I+b9uS+p81VD1LQav+h+RvB/7d9mWSPlFhPP1OYvkwgX37hwlUHBPAC5Imsvy50UYUNYgq9U8r9SpgD+BH5f6+wAKWdzevwlHATi29Zz9L8cdbklYHkrS6nO2zyqlrtmTF2Qq+VllQhW9LOg9YT9IHgOMp5kOsg9+UsR0AfE7SK6jH8JBnbT8raWCYgKQ6DBM4E7gUeJWk0ymeGZ1aZUC2jwOQdDmwTf9MMJI2ofqWhgcohi88W+6/AvhVZdE0TJJWl5N0IcUzmjtZXoMwxfIIVXqOovPAHyieN5xm+5pqQxrwHmAm8AXbvy9/0VX5zKhfLYcJ2P6GpDuA/Smauw61XdmciINs2Tp1GfAI8Lqqgik9B9wr6RqKf4sHAjf2Dx63PbvK4OouM2J0OUk/o/hLs1Y/aEn/BBxJ0SwyF7iqbjHWWdm5YDLwA9vP1yCeicDGrFib//XQ7xgfks6mmN/ymxQJ4khgqe0PVRjT+4c7b/uC8YqliZK0upyKNbRmD/prsxbKrtJ/RrHe13Tg28BXbaeppEEkfQj4OEUtpn+gbG1mKy87ZfTP83eD7UurjCdWTpoHu9+GwJJy5vKBAai231ldSAMxWNL/AP8D9ACvBC6WdI3tv682uhiFkygWNqzlhMJlT8EqO17EGEpNq8uVzUgvYvv68Y6llaTZFOt5PQZ8Bfie7RfKmel/afs1VcYXnStnNDmwf3qpOilrWZ+j6EUo6jXFVLwEqWl1uaqT0zA2BN41eO0g232S3jHEe6Ke7gcWSPo+K9bmv1hdSAM+Dxxco44hSNrS9gODju1q+/aKQmqUOnTjjVVI0lOS/lBuz0rqlfSHquMqZ7huu9hdnX7BREd+DVwDrAGs07LVwSM1/P/pEklT+nfK1pC5FcbTKGkefJmRdCjwJtv/UHUsEauapDOAV1MME2itBVY5I8auFFNfHQzsTLHe18G2H6wqpiZJ0noZknSr7d2rjiOaTdK/2f6wpP+knA2jVR06+0g6v81h2z5+3INpIWkGcB7FAOO32/5tlfE0SZ5pdblBc7BNoOhanr9UYixcWH79QqVRDKN/Zow6aJPc16JYJuWr5XyIlSf5JkjS6n4Ht7zuoZhC5pBqQoluYvuO8uWO5VyXA8ol5SvvBFTWtNrVAquoadU2uTdJmgcjYqVI+qntnQcdW2R7p6piaonj8JbdNSlWF3ioyqmSJG1MsXacy1geqSqWJkrS6nKSpgJnUaxTZeBG4CTbyyoNLBpP0lHAeykWM/xxy6l1gF7bB1QS2DDKcYA/tD3uy7lI2hH4d4rpt35THp4K/B74YBaC7EyaB7vf+cBFFOv3ABxTHjuwsoiiW9wMPEwx5u5fW44/BSyuJKKRTQM2r+je84C/tH1b60FJu5fndqggpsZJTavLSbrT9o4jHYvoRpKeomhhUPn1f4BTbH+3glh+aXvaEOeW2n7teMfURKlpdb/HJB1DMcs1FAvQ1XKOuGimOk+VZLsug5wBrixnDfka0D8mazPgfcAPKouqYVLT6nKSNgfOBmZQ/KV5M8Ws75UvGxHdQdJSajZVUitJ7wT2LncX2L68wlgOoui9O4UiuS8D5tu+oqqYmiZJq8tJugD4sO0nyv31KRY3rHRwZXQPSTfZ3rPqONopl7LfFfhGeegoYKHtU6qLKlZGklaXa9f1uC7dkaM71HGqpH6SFlOMI+sr9ycCi6pY60vSZOAUiprWq8rDjwKXAZ+1/fvxjqmJMmFu95sg6ZX9O2VNK88yYyytCzxDsaDnweVWp5n612t5PbmyKIpFTp8A9rW9ge0NgH0purx/p8K4GiU1rS4n6X0Uf91dTPFM6z3A6bYvHPaNEV2gHEv2WeA6imdIe1P0HvxWBbHcZ3vr0Z6LFSVpvQxI2gbYj+If7bW2l1QcUnSRug5glySKwbs9FM+1BNxm+38qiudq4IfABf2zYJSzYxxLsYhm7QZj11GSVkSsFEnXUAxg76+9HwMcbbvyAeyS7rC9S9VxAJTN9Cez4jOtR4D5wOds/66q2JokSSsiVkqdB7BLOgeYl1WBu0c6YkTEynpM0jGSJpbbMdRnAPu+wC2SfiVpsaS7yx6FtSKpNkuo1F1qWhGxUuo8gF3SFu2O2/7v8vwr+8cwVknSr21XNSdioyRpRcRKafIA9nbLqqzCew1VwxPwOtuvGI84mi7jdSJiZW3fWlux/TtJTRm8rnG818bAWynGag2O4eZxjKPRkrQiYmVNaG1ma9gA9vFsarocmGT7zsEnJC0YxzgarSn/Y0VEff0rcLOkFQawVxtS/dg+YZhz7x3PWJosz7QiYqU1dQD7eM7DKWkhcBNwJcVs88+Ox327TZJWRHQ1STsDe1HUAm+y/dOWc+uP16BeSauVccyk6Ir/OHAVcKXtX4xHDN0gSSsiupak04B3A/0zzh8KfMf2P1UXVUHSJsBBFElsGnCL7f9TbVT1l6QVEV1L0s+Anfqb4iT9CfBT22+oNrKCpAnAJOBpYIbtmyoOqfYyI0ZEdLMHgDVb9l8B/KqaUAqSLpK0rqS1gSXAfcDfJmF1JkkrIrrZc8C9kuZJOh+4B3ha0pmSzqwopm1s/4GiqfIKYHPgLyqKpXHS5T0iutml5dZvQUVxtFpd0uoUSets2y9IynOaDiVpRUTXsn1B1TG0cR5Fs+VdwA3l/Ih/qDSiBklHjIjoOpLuZpjZLmxvP47hjEjSarZ7qo6jCVLTiohu9I7y61+XX/sXqDwaeGb8w1muXK34M8Cmtg8qB2bPAL5aZVxNkZpWRHQtSTfZ3nOkY+Mc05XA+cDHbO9QDjpeZHu7qmJqkvQejIhutrakvfp3JO0JrF1hPAAb2v420AdQNgv2VhtSc6R5MCK62fHA+ZImUzzjehKoepXgP0raoIwHSbtTxBUdSPNgRHQtSR8pX04C/sjyxHVHuyVCximmnYGzgG0pxo1tBBxhe6hFIqNFklZEdC1JFwHTgfkUM9C/HbgdeD3FHISfryiu1YCty5jus/1CFXE0UZJWRHQtSVcBh9t+utyfBFwMHEZR29pmHGN513DnbV8y3Pko5JlWRHSzzYHnW/ZfALaw/b+SnhvnWA4uv74K2AP4Ubm/L8VMHUlaHUjSiohudhFwq6TLyv2DgW+2TFY7bmwfByDpcor5Bx8u9zcBzhnPWJoszYMR0dUk7UKx+KKAG20vrDiee2xv27I/AVjceiyGlqQVETGOJJ1NsejjNyl6Mx4JLLX9oUoDa4gkrYiIcVZ2ynhzuXuD7UuHKx/LJWlFRERjpCNGRMQ4kHSj7b0kPcWKM9ALsO11KwqtUVLTioiIxsiEuRER40jS7pLWadmfJGm3KmNqktS0IiLGkaRFwM4uf/mWXd4X2t652siaITWtiIjxJbfUFmz3kf4FHUvSiogYX/dLmi1p9XI7Cbi/6qCaIkkrImJ8/RXF3IO/AZYBuwGzKo2oQfJMKyIiGiPtqBER40DS39v+vKSzWHGcFgC2Z1cQVuMkaUVEjI+flV8rnbC36dI8GBERjZGaVkTEOJJ0He2bB/erIJzGSdKKiBhff9fyek3gcKCnolgaJ82DEREVk3S97bdUHUcTpKYVETGOJK3fsjsB2AV4dUXhNE6SVkTE+LqD4pmWKJoF/ws4odKIGiTNgxER0RiZxikiYhxJenf/0iSSTpV0iaTM8N6hJK2IiPH1j7afkrQX8FbgAuDcimNqjCStiIjx1Vt+fTtwru3LgDUqjKdRkrQiIsbXbySdB7wHuELSK8jv4o6lI0ZExDiStBYwE7jb9i8lbQJsZ/vqikNrhCStiIhojFRJIyKiMZK0IiKiMZK0IiKiMZK0IiKiMZK0IiKiMf4fh4nCIk3nl+EAAAAASUVORK5CYII=\n" ] }, "metadata": { "needs_background": "light" } } ], "execution_count": 14, "metadata": { "collapsed": false, "outputHidden": false, "inputHidden": false } }, { "cell_type": "markdown", "source": [ "## Observations & Inferences\n", "\n", "#### The case of unaccounted data\n", "Two largest populations China and India are missing in this data set. For instance In 2014 World population was ~ 7 Billion thse 2 countires alone accounted for ~ 2.5 Billion people. Appears this sample is already having some selection bias (an inadvertant bias ofcourse)\n", "\n", "#### Why GDP is an almost useless data point ? The case of Inequal Income Distribution\n", "While GDP is one indicator of a country's or geography's economic prosperity, it does not account of inequality in income or wealth distribution. (Thomas Picketty has captured this in his book and this **TED talk** [Capital](https://www.ted.com/talks/thomas_piketty_new_thoughts_on_capital_in_the_twenty_first_century). Why this is important ? Because the GDP figures can be influenced by few extremely well doing private enterprises but that doesnt mean people living below the poverty line are lifted out of poverty. So suicides stemming out of poverty needs a different metric. Say, we have a direct measure of poverty and if we know what % of the population are below the poverty line, we can see how it correlates with suicide rate. How much does the richest person make and the poorest person make can give an idea on how deep the schisms are. The [2019 World Economic forum at Davos summit](https://www.theguardian.com/business/2019/jan/21/world-26-richest-people-own-as-much-as-poorest-50-per-cent-oxfam-report) figured top 26 richest people own as much as poorest 50%.\n", "\n", "#### The case of age of the country \n", "Age seems to be good indicator of suicidal rate but overall age of the country is important when determining which country has higher suicidal rate. For instance India's average age is younger than Japan's. If our data suggests people in age group 74+ are highly suicidal then most people in Japan would be dead now. So this needs to be taken with a grain of salt.\n", "\n\n", "## Conclusions\n", "\n", "While Age and Sex seems to be a good indicators of the suicide rate purely by the virtue of correlation, still to know which country or geography more has higher suicidal rates and what are the causations are hard to draw from this dataset. Without which charting a course of remedy will be shooting in the dark at best. Some other factors that can be considered are prevalency of mental illness and gender dysphoria - because transgender suicide rates are [very high realtively speaking](https://www.youtube.com/watch?v=Kx_7biZoNaY). Access to (or Lack of) healthcare, Happines Index can be a good one to look into. Finally the case misreported or unreported suicides should be factored in as well to see what percent of our sample is a suicide really.\n", "\n\n" ], "metadata": {} } ], "metadata": { "kernel_info": { "name": "python3" }, "language_info": { "name": "python", "version": "3.7.1", "mimetype": "text/x-python", "codemirror_mode": { "name": "ipython", "version": 3 }, "pygments_lexer": "ipython3", "nbconvert_exporter": "python", "file_extension": ".py" }, "kernelspec": { "name": "python3", "language": "python", "display_name": "Python 3" }, "nteract": { "version": "0.12.3" } }, "nbformat": 4, "nbformat_minor": 4 }