stat545a-2013-hw02_lee-woo

Locale change and library embeddings

Sys.setenv(lang = "EN")
library("lattice")

This article uses the following data: link to the data

dat <- read.delim("http://www.stat.ubc.ca/~jenny/notOcto/STAT545A/examples/gapminder/data/gapminderDataFiveYear.txt")

The basic features of data are:

• Number of observations: 1704 <- calculated automatically by inline R code
• Number of columns: 6 <- calculated automatically by inline R code
• The data looks like:

tail(dat, 3)

##       country year      pop continent lifeExp gdpPercap
## 1702 Zimbabwe 1997 11404948    Africa   46.81     792.4
## 1703 Zimbabwe 2002 11926563    Africa   39.99     672.0
## 1704 Zimbabwe 2007 12311143    Africa   43.49     469.7

A basic summary of the data:

summary(dat)

##         country          year           pop              continent  
##  Afghanistan:  12   Min.   :1952   Min.   :6.00e+04   Africa  :624  
##  Albania    :  12   1st Qu.:1966   1st Qu.:2.79e+06   Americas:300  
##  Algeria    :  12   Median :1980   Median :7.02e+06   Asia    :396  
##  Angola     :  12   Mean   :1980   Mean   :2.96e+07   Europe  :360  
##  Argentina  :  12   3rd Qu.:1993   3rd Qu.:1.96e+07   Oceania : 24  
##  Australia  :  12   Max.   :2007   Max.   :1.32e+09                 
##  (Other)    :1632                                                   
##     lifeExp       gdpPercap     
##  Min.   :23.6   Min.   :   241  
##  1st Qu.:48.2   1st Qu.:  1202  
##  Median :60.7   Median :  3532  
##  Mean   :59.5   Mean   :  7215  
##  3rd Qu.:70.8   3rd Qu.:  9325  
##  Max.   :82.6   Max.   :113523  
## 

A density plot of GDP per capita (wikipedia link: what is GDP per capita?):

densityplot(dat$gdpPercap)

A scatterplot of African countries with X:GDP per Capita and Y:life expectancy:

xyplot(lifeExp ~ gdpPercap, data = dat, subset = continent == "Africa")