National Lifetable

Data Loading

-작업 디렉토리를 지우고 시작하니, 상황에 따라 대처할 것.

rm(list=ls())
load("lifetable_kr.rda")
ls()

##  [1] "dx"           "dx.f"         "dx.m"         "dx.o"        
##  [5] "e0x"          "e0x.o"        "lifetable"    "lifetable.kr"
##  [9] "mux"          "mux.o"        "qx"           "qx.f"        
## [13] "qx.m"         "qx.o"

• 다음 코드를 이용하여 국민생존함수를 그려보자.

plot(lifetable.kr$age,lifetable.kr$lx.m/1000,xlab="",ylab="",col="blue",type="l")
lines(lifetable.kr$age,lifetable.kr$lx.f/1000,xlab="",ylab="",col="red")
x<-c(lifetable.kr$age,rev(lifetable.kr$age))
y<-c(lifetable.kr$lx.m/1000,rev(lifetable.kr$lx.f/1000))
polygon(x,y,border=NA,density=15)
legend(80,100,legend=c("남","여"),col=c("blue","red"),lty=1)
title(main="국민생존함수(2012)",xlab="연령(세)",ylab="생존률(%)")