252. AR

Autoregressive

AR(𝑝)𝜑(𝐵)𝑥𝑡=𝑐+𝜀𝑡𝜑(𝐵)=1𝜑1𝐵𝜑2𝐵2𝜑𝑝𝐵𝑝

i.e.,

𝑥𝑡=𝑐+𝜑1𝑥𝑡1++𝜑𝑝𝑥𝑡𝑝+𝜀𝑡

Stationary if all roots of 𝜑(𝐵)=0 lie outside the unit circle

Parameters: 𝜑1,,𝜑𝑝, 𝑐, 𝜎2
Orders: 𝑝

Example: AR(1)

Given

  • Order: 𝑝=1
  • Parameters: 𝜑1=0.5, 𝑐=0
  • Data:
𝑡12345678910111213141516
x𝑡121081114129131614111518161317

Step 1 — formula

Substitute 𝑝=1 into the AR(𝑝) recursion:

𝑥𝑡=𝑐+𝜑1𝑥𝑡1+𝜀𝑡

Forecast (set 𝜀𝑡=0):

𝑥̂𝑡=𝑐+𝜑1𝑥𝑡1

Innovation:

𝜀𝑡=𝑥𝑡𝑥̂𝑡

Step 2 — apply at 𝑡=2

Plug in 𝑐=0, 𝜑1=0.5, and 𝑥1=12:

𝑥̂2=0+0.512=6𝜀2=𝑥2𝑥̂2=106=4

Step 3 — iterate

𝑡𝑥𝑡𝑥̂𝑡=0.5𝑥𝑡1𝜀𝑡=𝑥𝑡𝑥̂𝑡
2100.512=6106=4
380.510=585=3
4110.58=4114=7
5140.511=5.5145.5=8.5
6120.514=7127=5
790.512=696=3
8130.59=4.5134.5=8.5
9160.513=6.5166.5=9.5
10140.516=8148=6
11110.514=7117=4
12150.511=5.5155.5=9.5
13180.515=7.5187.5=10.5
14160.518=9169=7
15130.516=8138=5
16170.513=6.5176.5=10.5
170.517=8.5

Residuals are large because AR(1) cannot capture the trend or seasonality in 𝑥 — that is the role of ARIMA / SARIMA / ETS, applied to the same data in their own examples.