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="mauigasfield.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$Maui/1000,xlim=c(1979,2010),type='n',las=1,axes=FALSE,lwd=1,xlab="",ylab="")
box(lwd=2)
abline(h=c(1:5),col='lightgray',lwd=1)
abline(v=c(1980,1985,1990,1995,2000,2005,2010),col='lightgray',lwd=1)
mtext(side=1,line=-1.8,"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="Maui gas field production",xlab="Year", ylab="billion cubic metres per year")
lines(e2$Year,e2$Maui/1000, col=4,lwd=2)
dev.off()