Use the help function to explore what the series
autoplot()to plot some of the series in these data sets.
- What is the time interval of each series?
filter()to find what days corresponded to the peak closing price for each of the four stocks in
Download the file
tute1.csvfrom the book website, open it in Excel (or some other spreadsheet application), and review its contents. You should find four columns of information. Columns B through D each contain a quarterly series, labelled Sales, AdBudget and GDP. Sales contains the quarterly sales for a small company over the period 1981-2005. AdBudget is the advertising budget and GDP is the gross domestic product. All series have been adjusted for inflation.
You can read the data into R with the following script:
<- readr::read_csv("tute1.csv") tute1 View(tute1)
Convert the data to time series
<- tute1 %>% mytimeseries mutate(Quarter = yearmonth(Quarter)) %>% as_tsibble(index = Quarter)
Construct time series plots of each of the three series
%>% mytimeseries pivot_longer(-Quarter) %>% ggplot(aes(x = Quarter, y = value, colour = name)) + geom_line() + facet_grid(name ~ ., scales = "free_y")
Check what happens when you don’t include
USgaspackage contains data on the demand for natural gas in the US.
- Install the
- Create a tsibble from
us_totalwith year as the index and state as the key.
- Plot the annual natural gas consumption by state for the New England area (comprising the states of Maine, Vermont, New Hampshire, Massachusetts, Connecticut and Rhode Island).
- Install the
tourism.xlsxfrom the book website and read it into R using
- Create a tsibble which is identical to the
tourismtsibble from the
- Find what combination of
Purposehad the maximum number of overnight trips on average.
- Create a new tsibble which combines the Purposes and Regions, and just has total trips by State.
Create time plots of the following four time series:
help()) to find out about the data in each series.
- For the last plot, modify the axis labels and title.
aus_arrivalsdata set comprises quarterly international arrivals to Australia from Japan, New Zealand, UK and the US.
gg_subseries()to compare the differences between the arrivals from these four countries.
- Can you identify any unusual observations?
Monthly Australian retail data is provided in
aus_retail. Select one of the time series as follows (but choose your own seed value):
set.seed(12345678) <- aus_retail %>% myseries filter(`Series ID` == sample(aus_retail$`Series ID`,1))
Explore your chosen retail time series using the following functions:
ACF() %>% autoplot()
Can you spot any seasonality, cyclicity and trend? What do you learn about the series?
Use the following graphics functions:
ACF()and explore features from the following time series: “Total Private”
- Can you spot any seasonality, cyclicity and trend?
- What do you learn about the series?
- What can you say about the seasonal patterns?
- Can you identify any unusual years?
The following time plots and ACF plots correspond to four different time series. Your task is to match each time plot in the first row with one of the ACF plots in the second row.
aus_livestockdata contains the monthly total number of pigs slaughtered in Victoria, Australia, from Jul 1972 to Dec 2018. Use
filter()to extract pig slaughters in Victoria between 1990 and 1995. Use
ACFfor this data. How do they differ from white noise? If a longer period of data is used, what difference does it make to the ACF?
Use the following code to compute the daily changes in Google closing stock prices.
<- gafa_stock %>% dgoog filter(Symbol == "GOOG", year(Date) >= 2018) %>% mutate(trading_day = row_number()) %>% update_tsibble(index = trading_day, regular = TRUE) %>% mutate(diff = difference(Close))
Why was it necessary to re-index the tsibble?
Plot these differences and their ACF.
Do the changes in the stock prices look like white noise?