//Download "edfe2a.csv" from https://docs.google.com/spreadsheet/ccc?key=0AsLhMBA3NXL4dElpNnJHR0dqRG90LU9FVlg4UkZITmc
to your working R directory.
e2 <- read.csv("edfe2a.csv", skip= 5,header=TRUE, sep=",", na.strings=-99, dec=".", strip.white=TRUE)
str(e2) ## examine object ##
'data.frame': 41 obs. of 14 variables:
$ Year : int 1970 1971 1972 1973 1974 1975 1976 1977 1978 1979 ...
$ Maui : int NA NA NA NA NA NA NA NA NA 468 ...
$ Kapuni : int 158 301 414 488 515 566 1502 2347 2175 1077 ...
$ McKee : int NA NA NA NA NA NA NA NA NA NA ...
$ Mangahewa : int NA NA NA NA NA NA NA NA NA NA ...
$ Kaimiro..Ngatoro: int NA NA NA NA NA NA NA NA NA NA ...
$ Turangi : int NA NA NA NA NA NA NA NA NA NA ...
$ Kowhai : int NA NA NA NA NA NA NA NA NA NA ...
$ Pohokura : int NA NA NA NA NA NA NA NA NA NA ...
$ Tui : int NA NA NA NA NA NA NA NA NA NA ...
$ Maari : int NA NA NA NA NA NA NA NA NA NA ...
$ Kupe : int NA NA NA NA NA NA NA NA NA NA ...
$ Other : int NA NA NA NA NA NA NA NA NA NA ...
$ Total : int 158 301 414 488 515 566 1502 2347 2175 1545 ...
library("RSVGTipsDevice")
devSVGTips(file="nz-gas-production-m3.svg",toolTipMode=0,width=8,height=6)
par(mar=c(5,6,5,1)+0.1,cex.main=1.5,cex.lab=1.2)
plot(e2$Year,e2$Total/1000,type='n',las=1,axes=FALSE,lwd=1,xlab="",ylab="")
box(lwd=2)
abline(h=c(0.5,1,1.5,2,2.5),col='lightgray',lwd=1)
abline(v=c(1970,1980,1990,2000,2010),col='lightgray',lwd=1)
mtext(side=1,line=-2,"Data: Ministry of Economic Development \nEnergy Data File 2011 Table E.2.a")
axis(side=1, at = NULL, labels = TRUE, tick = TRUE)
axis(side=2, las=2,at = NULL, labels = TRUE, tick = TRUE)
title(col.main="black", cex.lab=1, main="Total New Zealand gas field production",xlab="Year", ylab="Billion cubic metres per year")
lines(e2$Year,e2$Total/1000, col=4,lwd=2)
dev.off()