Loading [Contrib]/a11y/accessibility-menu.js

Pages

2011-03-26

Ukraine's population pyramid with R

require(plotrix)
data = read.csv("UNPop_ukraine_age_sex_1985-2010.csv")
data = data[-379, ]
data1 = subset(data, select=c("Sex", "Year", "Age", "Value"))
xy.pop1985 = subset(data1[data1$Sex=="Male" & data1$Year==1985,], select="Value")$Value
xx.pop1985 = subset(data1[data1$Sex=="Female" & data1$Year==1985,], select="Value")$Value
agelabels = unique(subset(data1, select="Age")$Age)
xx.pop1985ratio = xx.pop1985 / sum(xx.pop1985) * 100
xy.pop1985ratio = xy.pop1985 / sum(xy.pop1985) * 100
par(mar=pyramid.plot(xy.pop1985ratio,xx.pop1985ratio,labels=agelabels,main="Ukraine population pyramid 1985",gap=1.5,show.values=TRUE))
UNPop_ukraine_age_sex_1985-2010.csv
(via World Population Prospects: The 2008 Revision Population Database)