Infant Mortality in India¶
Context: I thought about this after reading Doing good better (highly recommend). Soon after I found out about kepler.gl which made this the perfect project to test it out.
Some (back of the envelope) calculations to show that we loose 8.3 lakh babies every year.
In [1]:
# India rates; all data from 2011 or so. no consistent data after that :(
infant_mortality_rate = 34
crude_birth_rate = 20.4
population = 1210854977
In [2]:
def lives_lost(infant_mortality_rate, crude_birth_rate, population):
# One must note this number does not include anybody dying over the age of 1.
number_of_births = crude_birth_rate * population / 1000
number_of_dead = infant_mortality_rate * number_of_births /1000
return number_of_dead
In [3]:
l = lives_lost(infant_mortality_rate, crude_birth_rate, population)
l
Out[3]:
In [4]:
lives_lost_per_min = l /(365 * 24 * 60)
lives_lost_per_min
Out[4]:
In [5]:
lives_created_per_min = population * crude_birth_rate /(1000 *365 *24*60)
lives_created_per_min
Out[5]:
In [6]:
life_expectancy = 68.35
qaly_per_life = 60 #something a little less than life expectancy, hacky
In [7]:
qaly_lost = l * qaly_per_life
qaly_lost
Out[7]:
Infant mortality, Population, Birth rate by State¶
After lots of painful data processing here.
In [8]:
import pandas as pd
In [9]:
birth_rate = pd.read_csv("birth_rate.csv",index_col="states")
infant_mortality = pd.read_csv("infant_mortality.csv", index_col="states")
population = pd.read_csv("population.csv", index_col="states")
In [10]:
birth_rate[:5]
Out[10]:
In [11]:
infant_mortality[:5]
Out[11]:
In [12]:
population[:5]
Out[12]:
In [13]:
assert len(birth_rate) == len(infant_mortality) == len(population)
In [14]:
lives_lost = pd.DataFrame(index=birth_rate.index)
number_of_births = pd.DataFrame(index=birth_rate.index)
In [15]:
for i in ["total","urban","rural"]:
number_of_births[i] = birth_rate[i] * population[i] / 1000
lives_lost[i] = infant_mortality[i] * number_of_births[i] /1000
In [16]:
lives_lost
Out[16]:
In [17]:
lives_lost.to_csv("lives_lost.csv",index=False)
number_of_births.to_csv("number_births.csv", index=False)
Analysis¶
In [18]:
import matplotlib.pyplot as plt
%matplotlib inline
In [19]:
lives_lost = lives_lost[:-1] # removing all india data
infant_mortality = infant_mortality[:-1]
In [20]:
lives_lost.plot.bar(figsize=(20,10), title="lives lost per state")
plt.show()
In [21]:
infant_mortality.plot.bar(figsize=(20,10), title="Infant mortality rate")
plt.show()
In [22]:
infant_mortality.idxmax()
Out[22]:
In [23]:
infant_mortality.idxmin()
Out[23]:
In [24]:
infant_mortality.hist(figsize=(20,10))
plt.show()
In [25]:
infant_mortality.mean()
Out[25]:
In [26]:
infant_mortality.median()
Out[26]:
In [27]:
infant_mortality.loc["Madhya Pradesh"]
Out[27]: