# save the data as a comma separated values format; "icpswitch.csv"
12907 13880 13390 12491 13471 12923 12424 12625 12722 11062 12699 14624
12564 15151 16038 16156 15862 14450 16246 14533 15175 12607 13509 14211
13906 15768 21020 18618 17511 14173 15193 13672 13988 14188 15843 16283
15497 19363 17693 18925 16187 15469 15586 21955 23884 17032 21004 23513
22248 21316 23892 24621 22430 23147 23706 26438 25688 21795 26539 34892
30573 25986 27783 51831 30005 30267 30457 29109 27350 23214 30821 32605
30698 32949 43920 42306 41741 34999 31752 34532 32435 26761 35791 35998
29072 32045 32526 34503 35612 NA NA NA NA NA NA NA
# read data from "icpswitch.csv" as a numeric vector
i2<-scan(file="icpswitch.csv",skip=0,sep=",",nlines=8)
#convert data into a 'time series', first month is April 2005
i2ts<-ts(i2, frequency = 12, start = c(2005, 4))
# load 'RSVGTipsDevice' module
library("RSVGTipsDevice")
# create plot type and dimensions
devSVGTips(file="icp-switch.svg",toolTipMode=0,width=8,height=6)
plot(i2ts,ylim=c(0,50500),xlim=c(2005,2012.5),type="n",lwd=1,col=2,axes="FALSE",cex.main=1.6, main="Consumer switches of supplier",las=1,xlab="Date",ylab="")
points(i2ts,pch=16,col=4)
axis(side=1,las=1,at=c(2005:2012),labels=c(2005:2012))
axis(side=2, las=2)
grid()
mtext(side=1,line=-1.5,cex=1.2,"Source: Electricity Authority")
mtext(side=2,line=-1.2,cex=1,"Consumers per month")
box(lwd=2)
dev.off()