Russia is sadly notorious for its ridiculously high adult male mortality. According to Human Mortality Database data (2010), the probability for a Russian men to survive from 20 to 60 was just 0.64 ¹. For women the probability is 0.87. This huge gender disproportion in mortality results in a peculiar sex ratio profile (see my old DemoTrends post and the previous blog post).

Now let’s compare age-specific mortality rates of Russian men to that of the Japanese. For years and years Japan performs best in reducing mortality. It became standard to compare mortality schedules of other countries to the Japanese one ².

First, I need to get HMD data for both Russian and Japanese males. Again, I am using the amazing R package HMDPFDplus of Tim Riffe to download HMD data with just a couple of lines of R code.

# load required packages
library(tidyverse) # version 1.0.0
library(HMDHFDplus) # version 1.1.8

# load life tables for men, RUS and JPN
rus <- readHMDweb('RUS', "mltper_1x1", ik_user_hmd, ik_pass_hmd)
jpn <- readHMDweb('JPN', "mltper_1x1", ik_user_hmd, ik_pass_hmd)

Please note, the arguments ik_user_hmd and ik_pass_hmd are my login credidantials at the website of Human Mortality Database, which are stored locally at my computer. In order to access the data, one needs to create an account at http://www.mortality.org/ and provide his own credidantials to the readHMDweb() function.

Next, I select the most recent year for comparison, 2014, and compute the rate ratio of age specific mortality rates.

# compare mortality rates for 2014
ru <- rus %>% filter(Year == 2014) %>% transmute(age = Age, rus = mx)
jp <- jpn %>% filter(Year == 2014) %>% transmute(age = Age, jpn = mx)
df <- left_join(jp, ru, 'age') %>% mutate(ru_rate = rus / jpn)

Finally, I plot the resulting rate ratio of male mortality in Russia and Japan.

# get nice font
library(extrafont)
myfont <- "Roboto Condensed"

# plot
gg <- ggplot(df, aes(age, ru_rate)) + 
        geom_hline(yintercept = 1, color = 'red') +
        geom_line(aes(group=1)) + 
        scale_y_continuous('mortality rate ratio',
                           breaks = 0:10, labels = 0:10, limits = c(0, 10)) +
        annotate('text',x=c(0, 55), y = c(1.75,5), 
                 color = c('red','black'), hjust = 0, vjust = 1, size = 7,
                 label = c('Japan','Russia'), family = myfont) +
        ggtitle('Compare age-specific mortality of males',
                subtitle = "Russia and Japan, 2014, HMD")+
        theme_bw(base_size = 15, base_family = myfont)

In the middle ages, male mortality in Russian is up to 10 times higher than in Japan!

This post is based on my earlier twit and gist.