In statistica per un modello autoregressivo (casuale) valgono le seguenti relazioni, dette equazioni di Yule-Walker:

  • R y y = k = 1 N a k R y y ( n k ) R y x ( n ) {\displaystyle R_{yy}=-{\sum _{k=1}^{N}a_{k}R_{yy}(n-k) R_{yx}(n)}}
  • R y x ( n ) = { 0 p e r n > 0 σ n 2 p e r n = 0 {\displaystyle R_{yx}(n)=\left\{{\begin{matrix}0\,\,\,per\,\,\,n>0\\{\sigma _{n}^{2}}\,\,\,per\,\,\,n=0\,\end{matrix}}\right.}

In particolare, la matrice R {\displaystyle R} dei coefficienti delle equazioni di Yule-Walker è una matrice di Toeplitz; cioè è simmetrica (o hermitiana, nel caso di sequenze complesse) e tutti gli elementi appartenenti alla stessa diagonale, o subdiagonale, sono eguali tra loro. La matrice R [ N × N ] {\displaystyle R[N\times N]} è pertanto caratterizzata da N {\displaystyle N} numeri e può dunque essere rappresentata da:

R = [ r 0 r 1 . . r N r 1 r 0 . . r N 1 r 2 r 1 . . r N 2 . . . . . . . . r N r N 1 . . r 0 ] {\displaystyle R={\begin{bmatrix}r_{0}&r_{-1}&..&r_{-N}\\r_{1}&r_{0}&..&r_{-N 1}\\r_{2}&r_{1}&..&r_{-N 2}\\..&..&..&..\\r_{N}&r_{N-1}&..&r_{0}\\\end{bmatrix}}}

Nota: Per ricavare l'elemento m-esimo r m {\displaystyle r_{m}} si veda la procedura di derivazione sotto esposta.

Derivazione

Considerando un processo AR:

z t = i = 1 p ϕ i z t i α t . {\displaystyle z_{t}=\sum _{i=1}^{p}{\phi _{i}}\,z_{t-i} \alpha _{t}.}

Moltiplicando entrambi i membri per z t m {\displaystyle z_{t-m}} e, usando l'operatore del valore atteso, si ha:

E [ z t z t m ] = E [ i = 1 p ϕ i z t i z t m ] E [ α t z t m ] . {\displaystyle E[z_{t}z_{t-m}]=E\left[\sum _{i=1}^{p}{\phi _{i}}\,z_{t-i}z_{t-m}\right] E[\alpha _{t}z_{t-m}].}

Si ha che la funzione di autocovarianza è: E [ z t z t m ] = r m {\displaystyle E[z_{t}z_{t-m}]=r_{m}} . I valori della funzione del rumore bianco risultano indipendenti tra loro, e z t m {\displaystyle z_{t-m}} risulta indipendente da α t {\displaystyle {\alpha _{t}}} per m > 0. Se m 0 {\displaystyle m\neq 0\Rightarrow } E [ α t z t m ] = 0 {\displaystyle E[\alpha _{t}z_{t-m}]=0} . Per m = 0 {\displaystyle m=0} si ha:

E [ α t z t ] = E [ α t ( i = 1 p ϕ i z t i α t ) ] = i = 1 p ϕ i E [ α t z t i ] E [ α t 2 ] = 0 σ α 2 , {\displaystyle E[\alpha _{t}z_{t}]=E\left[\alpha _{t}(\sum _{i=1}^{p}{\phi _{i}}\,z_{t-i} \alpha _{t})\right]=\sum _{i=1}^{p}{\phi _{i}}\,E[\alpha _{t}\,z_{t-i}] E[\alpha _{t}^{2}]=0 \sigma _{\alpha }^{2},}

Pertanto, risulta:

r m = E [ i = 1 p ϕ i z t i z t m ] σ α 2 δ m . {\displaystyle r_{m}=E\left[\sum _{i=1}^{p}{\phi _{i}}\,z_{t-i}z_{t-m}\right] \sigma _{\alpha }^{2}\delta _{m}.}

Poiché:

E [ i = 1 p ϕ i z t i z t m ] = i = 1 p ϕ i E [ z t z t m i ] {\displaystyle E\left[\sum _{i=1}^{p}{\phi _{i}}\,z_{t-i}z_{t-m}\right]=\sum _{i=1}^{p}{\phi _{i}}\,E[z_{t}z_{t-m i}]} = i = 1 p ϕ i r m i , {\displaystyle \sum _{i=1}^{p}{\phi _{i}}\,r_{m-i},}

La risultante equazione di Yule-Walker è:

r m = i = 1 p ϕ i r m i σ α 2 δ m . {\displaystyle r_{m}=\sum _{i=1}^{p}{\phi _{i}}r_{m-i} \sigma _{\alpha }^{2}\delta _{m}.}

Bibliografia

  • G. U. Yule, On a method of investigating periodicities in disturbed series, with special reference to wolfer's sunspot numbers, Phil. Trans. Roy. Soc., 226-A:267–298, 1927.
  • Rob J Hyndman, Yule-Walker Type Estimates for Continuous Time Autoregressive Models, Dept. of Statistics, University of Melbourne, 1991.
  • Helmut Lütkepohl, Introduction to Multiple Time Series Analysis, ISBN 3540569405, Springer, 1993.
  • Jack HW Penm, Tim Brailsford, Richard Deane Terrell, The Adjustment of the Yule-Walker Relations in VAR Modeling: The Impact of the Euro on the Hong Kong, Canberra, A.C.T.: School of Finance and Applied Statistics, Australian National University, 2000.

Voci correlate

  • George Udny Yule
  • Schema di Yule

(PDF) Kernel autoregressive models using YuleWalker equations

arima yule walker equation Cross Validated

Solved In this question, we derive the autocovariances of an

2) Fluidodinamica le equazioni di Eulero CEMER

The modified YuleWalker method for multidimensional infinitevariance