x <- read.delim("./japancensus2010tp.txt") require(pyramid) png("./japanpyramidsallpref2010.png", width=1200, height=900) 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="’j«",Rlab="—«",Clab="",main=paste(i),AxisFM="fg") } dev.off()