### Average Air Temperature Racife Brazil Splus

# to read racife.dat 

x <- scan("E:\\Stat207\\Midterm\\Racife.dat")

# you can make x a regular time series object

x <- rts(x, start=c(1953,1), freq=12, units="months")

# plot the data

par(mfrow=c(1,1)) 
ts.plot(x, xlab="Month", ylab="Air Temperature")

# Classical Decomposition Model

x.sabl <- sabl(x)
sablplot(x.sabl, "Temperature")