Time Series - Moving Average


Advertisements

For a stationary time series, a moving average model sees the value of a variable at time ‘t’ as a linear function of residual errors from ‘q’ time steps preceding it. The residual error is calculated by comparing the value at the time ‘t’ to moving average of the values preceding.

Mathematically it can be written as −

y_t=c+ϵ_t+θ_1 ϵ_(t-1)+θ_2 ϵ_(t-2)+...+θ_q ϵ_(t-q)

Where‘q’ is the moving-average trend parameter

ϵ_t is white noise, and

ϵ_(t-1), ϵ_(t-2)… ϵ_(t-q) are the error terms at previous time periods.

Value of ‘q’ can be calibrated using various methods. One way of finding the apt value of ‘q’ is plotting the partial auto-correlation plot.

A partial auto-correlation plot shows the relation of a variable with itself at prior time steps with indirect correlations removed, unlike auto-correlation plot which shows direct as well as indirect correlations, let’s see how it looks like for ‘temperature’ variable of our data.

// Code Snippet 10: Showing PACP

A partial auto-correlation is read in the same way as a correlogram.

In [143]:

from statsmodels.graphics.tsaplots import plot_pacf

plot_pacf(train, lags = 100)
plt.show()
Code Snippet 10
Advertisements