block4.py

japanew = pd.read_csv(R"./Resources/japanew.csv") #import sample dataset
japanew["Japan_MX"] = japanew.iloc[:,1]/japanew.iloc[:,2] #create Japan death rate
japanew["EngW_MX"] = japanew.iloc[:,3]/japanew.iloc[:,4] #create England Wales death rate

#create log death rate
japanew["Japan_Lg_MX"] = np.log(japanew["Japan_MX"])
japanew["EngW_Lg_MX"] = np.log(japanew["EngW_MX"])

#age structure
japanew["Japan_AgeStr"] = japanew.iloc[:,2]/japanew.iloc[:,2].sum() 
japanew["EngW_AgeStr"] = japanew.iloc[:,4]/japanew.iloc[:,4].sum() 
japanew["MeanAgeStr"] = (japanew["Japan_AgeStr"]+japanew["EngW_AgeStr"])/2

#calculate standardized death rate
japanew["Japan_Std_DR"] = japanew["MeanAgeStr"]*japanew["Japan_MX"]
japanew["EngW_Std_DR"] = japanew["MeanAgeStr"]*japanew["EngW_MX"]

jap_std_DR = japanew["Japan_Std_DR"].sum()
ew_std_DR = japanew["EngW_Std_DR"].sum()

print("Japan Standardized Death Rate: %.5f" %jap_std_DR)
print("England/Wales Standardized Death Rate: %.5f" %ew_std_DR)