x <- read.delim("http://minato.sip21c.org/ldaR/japancensus2010tp.txt", encoding="CP932") if (!require(pyramid)) { install.packages("pyramid"); library(pyramid) } layout(matrix(1:48, 6, 8, byrow=TRUE)) par(mar=c(1, 1, 2, 1), cex=0.6) AreaList <- names(table(x$Area)) for (i in AreaList) { y <- subset(x, Area==i) pyramid(data.frame(Males=y$Males, Females=y$Females, Ages=y$Ages), Cadj=-0.01, Cgap=0.5, Csize=0.6, Llab="M", Rlab="F", Clab="", main=paste(i), AxisFM="d", Laxis=pretty(range(y$Males, y$Females), n=2)) }