if (require(wpp2019)==FALSE) { install.packages("wpp2019", dep=TRUE) library(wpp2019) } if (require(pyramid)==FALSE) { install.packages("pyramid", dep=TRUE) library(pyramid) } data(popM) data(popF) COUNTRIES <- c("Brazil", "China", "France", "Japan") YEARS <- c("1950", "1985", "2020") windows(width=2048/72, height=1109/72) layout(matrix(1:(4*3), 3, 4)) for (COUNTRY in COUNTRIES) { for (YEAR in YEARS) { MP <- subset(popM, name==COUNTRY)[, c("age", YEAR)] FP <- subset(popF, name==COUNTRY)[, c("age", YEAR)] POP <- data.frame(M=MP[, YEAR], F=FP[, YEAR], A=MP$age) pyramid(POP, Laxis=pretty(c(0,max(c(POP$M, POP$F)))), Cadj=-0.04, Cgap=0.2, AxisFM="d", Llab="Males (x1000)", Rlab="Females (x1000)", Lcol="#bfe4ff", Rcol="#ffcabf", main=paste("Population of",COUNTRY,"in",YEAR,"(wpp2019)",sep=" ")) } }