library(reshape)
data(EduDat)
data(dictionary)
icountry <- "China"
ipop <- subset(EduDat,cc==getcode(icountry,dictionary) & scen2 != "BP")
ipop <- cast(ipop,yr~scen2,sum)
row.names(ipop) <- ipop[,1]
matplot(as.data.frame(ipop[,-1]),type="l",lwd=3,main=paste("Projected Population 15+ by Education Scenario, ",icountry,sep=""),cex.main=1.8,
xlab="Year",ylab="Total Population in 1000s",xaxt="n",cex=1.5,cex.axis=1.5,cex.lab=1.5)
axis(1,at=1:nrow(ipop),labels=row.names(ipop),cex.axis=1.5)
grid(lwd=2,lty=2)
legend(legend=names(ipop[,-1]),title="Education Scenario: ",
col=1:ncol(ipop[,-1]),lwd=3,lty=1:ncol(ipop[,-1]),
"topleft",inset=c(0.02,0.02),bty="o",bg="grey",cex=1.5)
Run the code above in your browser using DataLab