library(demography) USER <- readline("username> ") PWD <- readline("password> ") # RStudioでは上のコードはうまく動かないかもしれない # 素のRでやるか,USER <- "hogehoge@foo.bar"のように # 直接ユーザとパスワードは付値してしまうと良い。 jpn <- hmd.mx(country="JPN", username=USER, password=PWD) rus <- hmd.mx(country="RUS", username=USER, password=PWD) ukr <- hmd.mx(country="UKR", username=USER, password=PWD) jpnm2013 <- jpn$rate$male[, "2013"] rusm2013 <- rus$rate$male[, "2013"] ukrm2013 <- ukr$rate$male[, "2013"] matplot(0:80, cbind(jpnm2013[1:81], rusm2013[1:81], ukrm2013[1:81]), type="b", lty=1:3, pch=1:3, col=1:3, xlab="Age (0-80)", ylab="Age-specific mortality rate (ASMR)") legend("topleft", lty=1:3, pch=1:3, col=1:3, legend=c("Japan", "Russia", "Ukraine")) matplot(0:80, cbind(rusm2013[1:81]-jpnm2013[1:81], ukrm2013[1:81]-jpnm2013[1:81]), type="b", lty=1:2, pch=1:2, col=1:2, xlab="Age (0-80)", ylab="Rate differences with Japan's ASMR") legend("topleft", lty=1:2, pch=1:2, col=1:2, legend=c("Russia", "Ukraine")) matplot(0:110, cbind(rusm2013/jpnm2013, ukrm2013/jpnm2013), type="b", lty=1:2, pch=1:2, col=1:2, xlab="Age", ylab="Rate ratios to Japan's ASMR") legend("topleft", lty=1:2, pch=1:2, col=1:2, legend=c("Russia", "Ukraine"))