We will consider two examples of decomposition. The first an example where the time series is additive and the second where the time series is multiplicative.

Example from OTexts Forecasting: principles and practice fpp

What is white noise? The term white noise means the Irregular part of the time series is Normally distributed with mean zero and it has constand variance.

x.error <- as.ts( rnorm(100) )
plot(x.error)

So when we decompose a time series the Irregular part should look like white noise.

Additive decomposition

Consider the beer production data.

library(fpp)
class(ausbeer)
[1] "ts"
data("ausbeer")

Plot the data.

plot(ausbeer, type="o", xlab="Year")

plot(decompose( ausbeer, "additive" ), xlab="Year")

Example from ts4

Mutliplicative decomposition

Consider the Johnson & Johnson data.

library(astsa)
class(jj)
[1] "ts"
data(jj)

Plot the data

plot(jj, type="o", xlab="Quarter", ylab="Quarterly Earnings per Share")

We see the seasonal component increasing in time. So we should consider a multiplicative decomposition.

plot(decompose( jj, "multiplicative" ), xlab="Quarter")

What if we use an additive model.

plot(decompose( jj, "additive" ), xlab="Quarter")

The Irregular part does not look like white noise.

Using piping from dplyr

jj %>% decompose("multiplicative") %>% plot()

jj %>% decompose() %>% plot()

Using a different R function that also does decomposition. stl()

ausbeer %>% stl("periodic") %>% plot()

jj %>% stl("periodic") %>% plot()

log(jj) %>% stl("periodic") %>% plot()

LS0tCnRpdGxlOiAiRGVjb21wb3NpdGlvbiIKb3V0cHV0OiBodG1sX25vdGVib29rCi0tLQoKV2Ugd2lsbCBjb25zaWRlciB0d28gZXhhbXBsZXMgb2YgZGVjb21wb3NpdGlvbi4gIFRoZSBmaXJzdCBhbiBleGFtcGxlIHdoZXJlIHRoZSB0aW1lIHNlcmllcyBpcyBhZGRpdGl2ZSBhbmQgdGhlIHNlY29uZCB3aGVyZSB0aGUgdGltZSBzZXJpZXMgaXMgbXVsdGlwbGljYXRpdmUuCgojIEV4YW1wbGUgZnJvbSBbT1RleHRzXShodHRwczovL3d3dy5vdGV4dHMub3JnLykgRm9yZWNhc3Rpbmc6IHByaW5jaXBsZXMgYW5kIHByYWN0aWNlIFtmcHBdKGh0dHBzOi8vd3d3Lm90ZXh0cy5vcmcvZnBwKQoKIyBXaGF0IGlzIHdoaXRlIG5vaXNlPyAgVGhlIHRlcm0gd2hpdGUgbm9pc2UgbWVhbnMgdGhlIElycmVndWxhciBwYXJ0IG9mIHRoZSB0aW1lIHNlcmllcyBpcyBOb3JtYWxseSBkaXN0cmlidXRlZCB3aXRoIG1lYW4gemVybyBhbmQgaXQgaGFzIGNvbnN0YW5kIHZhcmlhbmNlLgoKYGBge3J9CnguZXJyb3IgPC0gYXMudHMoIHJub3JtKDEwMCkgKQpwbG90KHguZXJyb3IpCmBgYAoKU28gd2hlbiB3ZSBkZWNvbXBvc2UgYSB0aW1lIHNlcmllcyB0aGUgSXJyZWd1bGFyIHBhcnQgc2hvdWxkIGxvb2sgbGlrZSB3aGl0ZSBub2lzZS4KCiMgQWRkaXRpdmUgZGVjb21wb3NpdGlvbgoKQ29uc2lkZXIgdGhlIGJlZXIgcHJvZHVjdGlvbiBkYXRhLgoKYGBge3J9CmxpYnJhcnkoZnBwKQoKY2xhc3MoYXVzYmVlcikKCmRhdGEoImF1c2JlZXIiKQpgYGAKClBsb3QgdGhlIGRhdGEuCgpgYGB7cn0KcGxvdChhdXNiZWVyLCB0eXBlPSJvIiwgeGxhYj0iWWVhciIpCmBgYAoKYGBge3J9CnBsb3QoZGVjb21wb3NlKCBhdXNiZWVyLCAiYWRkaXRpdmUiICksIHhsYWI9IlllYXIiKQpgYGAKCgojIEV4YW1wbGUgZnJvbSBbdHM0XShodHRwOi8vd3d3LnN0YXQucGl0dC5lZHUvc3RvZmZlci90c2E0L3RzYUVaLnBkZikKCiMgTXV0bGlwbGljYXRpdmUgZGVjb21wb3NpdGlvbgoKQ29uc2lkZXIgdGhlIEpvaG5zb24gJiBKb2huc29uIGRhdGEuCgpgYGB7cn0KbGlicmFyeShhc3RzYSkKCmNsYXNzKGpqKQoKZGF0YShqaikKYGBgCgpQbG90IHRoZSBkYXRhCgpgYGB7cn0KcGxvdChqaiwgdHlwZT0ibyIsIHhsYWI9IlF1YXJ0ZXIiLCB5bGFiPSJRdWFydGVybHkgRWFybmluZ3MgcGVyIFNoYXJlIikKYGBgCgpXZSBzZWUgdGhlIHNlYXNvbmFsIGNvbXBvbmVudCBpbmNyZWFzaW5nIGluIHRpbWUuICBTbyB3ZSBzaG91bGQgY29uc2lkZXIgYSBtdWx0aXBsaWNhdGl2ZSBkZWNvbXBvc2l0aW9uLgoKYGBge3J9CnBsb3QoZGVjb21wb3NlKCBqaiwgIm11bHRpcGxpY2F0aXZlIiApLCB4bGFiPSJRdWFydGVyIikKYGBgCgpXaGF0IGlmIHdlIHVzZSBhbiBhZGRpdGl2ZSBtb2RlbC4KCmBgYHtyfQpwbG90KGRlY29tcG9zZSggamosICJhZGRpdGl2ZSIgKSwgeGxhYj0iUXVhcnRlciIpCmBgYAoKVGhlIElycmVndWxhciBwYXJ0IGRvZXMgbm90IGxvb2sgbGlrZSB3aGl0ZSBub2lzZS4KCgojIFVzaW5nIHBpcGluZyBmcm9tIGRwbHlyCgpgYGB7cn0KamogJT4lIGRlY29tcG9zZSgibXVsdGlwbGljYXRpdmUiKSAlPiUgcGxvdCgpCmBgYAoKYGBge3J9CmpqICU+JSBkZWNvbXBvc2UoKSAlPiUgcGxvdCgpCmBgYAoKVXNpbmcgYSBkaWZmZXJlbnQgUiBmdW5jdGlvbiB0aGF0IGFsc28gZG9lcyBkZWNvbXBvc2l0aW9uLiAgc3RsKCkKCmBgYHtyfQphdXNiZWVyICU+JSBzdGwoInBlcmlvZGljIikgJT4lIHBsb3QoKQpgYGAKCmBgYHtyfQpqaiAlPiUgc3RsKCJwZXJpb2RpYyIpICU+JSBwbG90KCkKYGBgCgoKYGBge3J9CmxvZyhqaikgJT4lIHN0bCgicGVyaW9kaWMiKSAlPiUgcGxvdCgpCmBgYAoKCgo=