############################
# Author: Joshua M. Tebbs  #
# Date: 20 Dec 2009        #
# Update: 25 Jul 2011      #
# Purpose: STAT 520 R code #
# Modified by Hanson 2015  #
# CHAPTER 1                #
############################

library("TSA") # loads Chan's R package (once it is installed as shown in class)

# Example 1.1.
# Page 2
globaltemps=ts(read.table("http://people.stat.sc.edu/hansont/stat520/globaltemps.txt"),start=1856)
plot(globaltemps,ylab="Global temperature deviations",xlab="Year",type="o")

# Example 1.2.
# Page 3
data(milk)
plot(milk,ylab="Amount of milk produced",xlab="Year",type="o")

# Example 1.3.
# Page 4
data(CREF)
plot(CREF,ylab="CREF stock values",type="o")

# Example 1.4.
# Page 5
homeruns=ts(read.table("http://people.stat.sc.edu/hansont/stat520/homeruns.txt"),start=1909)
plot(homeruns,ylab="Number of homeruns",xlab="Year",type="o")

# Example 1.5.
# Page 6
earthquake=ts(read.table("http://people.stat.sc.edu/hansont/stat520/earthquake.txt"),start=1900)
plot(earthquake,ylab="Number of earthquakes (7.0 or greater)",xlab="Year",type="o")

# Example 1.6.
# Page 7
enrollment=ts(read.table("http://people.stat.sc.edu/hansont/stat520/enrollment.txt"),start=1954)
plot(enrollment,ylab="USC Columbia fall enrollment",xlab="Year",type="o")

# Example 1.7.
# Page 8
data(star)
plot(star,ylab="Star brightness",type="o")

# Example 1.8.
# Page 9
data(airmiles)
plot(airmiles,ylab="Airline miles",xlab="Year",type="o")

# Example 1.9.
# Page 10 
sp500=ts(read.table("http://people.stat.sc.edu/hansont/stat520/sp500.txt"))
plot(sp500,ylab="SP500 Index",xlab="Time",type="o")

# Example 1.10.
# Page 11
ventilation=ts(read.table("http://people.stat.sc.edu/hansont/stat520/ventilation.txt"))
plot(ventilation,ylab="Ventilation (L/min)",xlab="Observation time",type="o")

# Example 1.11.
# Page 12
exchangerate=ts(read.table("http://people.stat.sc.edu/hansont/stat520/exchangerate.txt"),start=1980,frequency=52)
plot(exchangerate,ylab="British pounds",xlab="Year",type="o")

# Example 1.12.
# Page 13
data(oil.price)
plot(oil.price,ylab="Oil prices",xlab="Year",type="o")

# Example 1.13.
# Page 14
data(larain)
plot(larain,ylab="LA rainfall amounts",xlab="Year",type="o")

# Example 1.14.
# Page 15
brick=ts(read.table("http://people.stat.sc.edu/hansont/stat520/brick.txt"),start=1956,freq=4)
plot(brick,ylab="Australian clay brick production (in millions)",xlab="Time",type="o")

# Example 1.15.
# Page 16
supremecourt=ts(read.table("http://people.stat.sc.edu/hansont/stat520/supremecourt.txt"),start=1926)
plot(supremecourt,ylab="Percentage granted review",xlab="Time",type="o")

# Example 1.8.
# Add special plotting symbols for months
# Page 19
data(airmiles)
plot(airmiles,ylab="Airline miles",xlab="Year",type='l')
points(y=airmiles,x=time(airmiles),pch=as.vector(season(airmiles)),cex=1)