tfr <- wpp.indicator("fert")
# Histogram of TFR
print(qplot(value, data=tfr) + facet_wrap(~ Year))
mxM01 <- wpp.indicator("mortagesex", sex="M", age="0")
# Plot map
plot(gvisGeoMap(tfr, locationvar='charcode', numvar='value'))
Run the code above in your browser using DataLab