Juan Buzzetti — Oct 18, 2012, 10:34 AM
# R Assignment - Part 2 - Quantitative Research Methods
# Juan Buzzetti
# ID 9485651
# import dataset
wced<-read.csv("/Users/Buzz/Desktop/HENV665Q/Lab7/wced.csv",header=T)
names(wced)
[1] "City" "Year" "Population"
[4] "Emission_total" "Emissions_percapita" "Region"
attach(wced)
View(wced)
# total emissions for each region
plot(Region,Emission_total, xlab="Region", ylab="Total Emissions - Millions T.CO2", main="CO2 Emissions by region")
# total emissions by population
plot(Population ,Emission_total, xlab="Population", ylab="Per Capita - T.CO2 per person", main="CO2 Emissions by population")
# total emissions by population (different colours for regions)
plot(Population ,Emission_total, pch=1, col=c("blue","red","green","yellow","black")[Region], xlab="Population", ylab="CO2 Emissions", main="CO2 Emissions by population")
# legend
legend(14300000,53, c("Africa", "Asia", "Europe", "N. America", "S. America"), pch=c(1), cex=1, col=c("blue","red","green","yellow","black"))
# data subset excluding S. American and African cities
wced_new<-subset(wced, Region!="S. America" & Region!="Africa")
attach(wced_new)
The following object(s) are masked from 'wced':
City, Emission_total, Emissions_percapita, Population, Region,
Year
View(wced_new)