x <- read.delim("./deaths.txt") library(pyramid) png("swedenkazakhstan.png", width=1200, height=400, pointsize=12) layout(t(1:3)) pyramid(data.frame(x$NSW,x$NK,x$AG), Llab="Sweden", Rlab="Kazakhstan", Clab="Age\n group", Cadj=0, AxisFM="d", AxisBM=",", Laxis=0:2*200000, Raxis=0:4*250000, main="Population structure of Sweden\n and Kazakhstan females in 1992") pyramid(data.frame(x$DSW,x$DK,x$AG), Llab="Sweden", Rlab="Kazakhstan", Clab="Age\n group", Cadj=0, AxisFM="d", AxisBM=",", main="Number of death by age in Sweden\n and Kazakhstan females in 1992") pyramid(data.frame(x$DSW/x$NSW,x$DK/x$NK,x$AG), Llab="Sweden", Rlab="Kazakhstan", Clab="Age\n group", Cadj=0, AxisFM="g", main="Age specific death rates of Sweden\n and Kazakhstan females in 1992") dev.off()