8.11 Exercises
 Figure 8.24 shows the ACFs for 36 random numbers, 360 random numbers
and for 1,000 random numbers.
 Explain the differences among these figures. Do they all indicate the data are white noise?
 Why are the critical values at different distances from the mean of zero? Why are the autocorrelations different in each figure when they each refer to white noise?
 A classic example of a nonstationary series is the daily closing
IBM stock prices (data set
ibmclose
). Use R to plot the daily closing prices for IBM stock and the ACF and PACF. Explain how each plot shows the series is nonstationary and should be differenced.  For the following series, find an appropriate BoxCox transformation and order of differencing in order to obtain stationary data.

usnetelec

usgdp

mcopper

enplanements

visitors

 For the
enplanements
data, write down the differences you chose above using backshift operator notation.  Use R to simulate and plot some data from simple ARIMA models.
 Use the following R code to generate data from an AR(1) model with $\phi _{1} = 0.6$ and
$\sigma ^2=1$. The process starts with $y_0=0$.
R codey < ts(numeric(100))
e < rnorm(100)
for(i in 2:100)
y[i] < 0.6*y[i1] + e[i]  Produce a time plot for the series. How does the plot change as you change $\phi _1$?
 Write your own code to generate data from an MA(1) model with $\theta_{1} = 0.6$ and $\sigma^2=1$.
 Produce a time plot for the series. How does the plot change as you change $\theta _1$?
 Generate data from an ARMA(1,1) model with $\phi _{1}$ = 0.6 and $\theta _{1} = 0.6$ and $\sigma ^2=1$.
 Generate data from an AR(2) model with $\phi_{1} =0.8$ and $\phi_{2} = 0.3$ and $\sigma ^2=1$. (Note that these parameters will give a nonstationary series.)
 Graph the latter two series and compare them.
 Use the following R code to generate data from an AR(1) model with $\phi _{1} = 0.6$ and
$\sigma ^2=1$. The process starts with $y_0=0$.
 Consider the number of women murdered each year (per 100,000
standard population) in the United States (data set
wmurders
). By studying appropriate graphs of the series in R, find an appropriate ARIMA($p,d,q$) model for these data.
 Should you include a constant in the model? Explain.
 Write this model in terms of the backshift operator.
 Fit the model using R and examine the residuals. Is the model satisfactory?
 Forecast three times ahead. Check your forecasts by hand to make sure you know how they have been calculated.
 Create a plot of the series with forecasts and prediction intervals for the next three periods shown.
 Does
auto.arima
give the same model you have chosen? If not, which model do you think is better?
 Consider the quarterly number of international tourists to Australia
for the period 1999–2010. (Data set
austourists
.) Describe the time plot.
 What can you learn from the ACF graph?
 What can you learn from the PACF graph?
 Produce plots of the seasonally differenced data $(1  B^{4})Y_{t}$. What model do these graphs suggest?
 Does
auto.arima
give the same model that you chose? If not, which model do you think is better?  Write the model in terms of the backshift operator, and then without using the backshift operator.
 Consider the total net generation of electricity (in billion
kilowatt hours) by the U.S. electric industry (monthly for the
period 1985–1996). (Data set
usmelec
.) In general there are two peaks per year: in midsummer and midwinter. Examine the 12month moving average of this series to see what kind of trend is involved.
 Do the data need transforming? If so, find a suitable transformation.
 Are the data stationary? If not, find an appropriate differencing which yields stationary data.
 Identify a couple of ARIMA models that might be useful in describing the time series. Which of your models is the best according to their AIC values?
 Estimate the parameters of your best model and do diagnostic testing on the residuals. Do the residuals resemble white noise? If not, try to find another ARIMA model which fits better.
 Forecast the next 15 years of generation of electricity by the U.S. electric industry. Get the latest figures from http://data.is/zgRWCO to check on the accuracy of your forecasts.
 How many years of forecasts do you think are sufficiently accurate to be usable?
 For the
mcopper
data: if necessary, find a suitable BoxCox transformation for the data;
 fit a suitable ARIMA model to the transformed data using
auto.arima()
;  try some other plausible models by experimenting with the orders chosen;
 choose what you think is the best model and check the residual diagnostics;
 produce forecasts of your fitted model. Do the forecasts look reasonable?
 compare the results with what you would obtain using
ets()
(with no transformation).
 Choose one of the following seasonal time series:
condmilk
,hsales
,uselec
 Do the data need transforming? If so, find a suitable transformation.
 Are the data stationary? If not, find an appropriate differencing which yields stationary data.
 Identify a couple of ARIMA models that might be useful in describing the time series. Which of your models is the best according to their AIC values?
 Estimate the parameters of your best model and do diagnostic testing on the residuals. Do the residuals resemble white noise? If not, try to find another ARIMA model which fits better.
 Forecast the next 24 months of data using your preferred model.
 Compare the forecasts obtained using
ets()
.
 For the same time series you used in exercise Q10, try using a nonseasonal model applied to the seasonally adjusted data obtained from STL. The
stlf()
function will make the calculations easy (withmethod="arima"
). Compare the forecasts with those obtained in exercise Q10. Which do you think is the best approach?