254. ARMA

Autoregressive moving average

ARMA(𝑝,𝑞)𝜑(𝐵)𝑥𝑡=𝑐+𝜃(𝐵)𝜀𝑡𝜑(𝐵)=1𝜑1𝐵𝜑𝑝𝐵𝑝


𝜃(𝐵)=1+𝜃1𝐵++𝜃𝑞𝐵𝑞


Combines AR and MA structure. Requires stationary input series.

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

Example: ARMA(1,1)

Given

  • Orders: 𝑝=1, 𝑞=1
  • Parameters: 𝜑1=0.5, 𝜃1=0.5, 𝑐=0
  • Initial conditions: 𝑥0=0, 𝜀0=0
  • Data:
𝑡12345678910111213141516
x𝑡121081114129131614111518161317

Step 1 — formula

Substitute 𝑝=1, 𝑞=1 into the ARMA(𝑝,𝑞) recursion:

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

Forecast (set 𝜀𝑡=0):

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

Innovation:

𝜀𝑡=𝑥𝑡𝑥̂𝑡

Step 2 — apply at 𝑡=1

Plug in 𝑐=0, 𝜑1=0.5, 𝜃1=0.5, 𝑥0=0, 𝜀0=0:

𝑥̂1=0+0.50+0.50=0𝜀1=𝑥1𝑥̂1=120=12

Step 3 — iterate

Each row uses the previous observation 𝑥𝑡1 and the previous innovation 𝜀𝑡1. Values rounded to 4 decimal places.

𝑡𝑥𝑡𝑥̂𝑡=0.5𝑥𝑡1+0.5𝜀𝑡1𝜀𝑡=𝑥𝑡𝑥̂𝑡
1120.50+0.50=0120=12
2100.512+0.512=6+6=121012=2
380.510+0.5(2)=51=484=4
4110.58+0.54=4+2=6116=5
5140.511+0.55=5.5+2.5=8148=6
6120.514+0.56=7+3=101210=2
790.512+0.52=6+1=797=2
8130.59+0.52=4.5+1=5.5135.5=7.5
9160.513+0.57.5=6.5+3.75=10.251610.25=5.75
10140.516+0.55.75=8+2.875=10.8751410.875=3.125
11110.514+0.53.125=7+1.5625=8.5625118.5625=2.4375
12150.511+0.52.4375=5.5+1.2188=6.7188156.7188=8.2813
13180.515+0.58.2813=7.5+4.1406=11.64061811.6406=6.3594
14160.518+0.56.3594=9+3.1797=12.17971612.1797=3.8203
15130.516+0.53.8203=8+1.9102=9.9102139.9102=3.0898
16170.513+0.53.0898=6.5+1.5449=8.0449178.0449=8.9551
170.517+0.58.9551=8.5+4.4775=12.9775