时间序列与R语言应用(part5)--移动平均MA模型及其可逆性
學習筆記
參考書目:《計量經(jīng)濟學》、《計量經(jīng)濟學模型與R語言應用》
文章目錄
- 一階移動平均過程
- qqq階移動平均過程
- 可逆性
一階移動平均過程
在研究qqq階移動平均過程之前,我們先以一階移動平均過程開刀,即MA(1)MA(1)MA(1)序列, 一階移動平均模型表達式為:
Yt=et?θet?1(1)Y_t=e_t-\theta e_{t-1}\tag{1} Yt?=et??θet?1?(1)
其中ete_tet?和et?1e_{t-1}et?1?都是白噪聲,顯然,E(Yt)=0,Var(Yt)=σe2(1+θ2)E(Y_t)=0, Var(Y_t)=\sigma_e^2(1+\theta^2)E(Yt?)=0,Var(Yt?)=σe2?(1+θ2),此時有:
Cov(Yt,Yt?1)=Cov(et?θet?1,et?1?θet?2)=?θCov(et?1,et?1)=?θσe2(2)Cov(Y_t,Y_{t-1})=Cov(e_t-\theta e_{t-1}, e_{t-1}-\theta e_{t-2})=-\theta Cov(e_{t-1}, e_{t-1})=-\theta \sigma_e^2\tag{2} Cov(Yt?,Yt?1?)=Cov(et??θet?1?,et?1??θet?2?)=?θCov(et?1?,et?1?)=?θσe2?(2)
和:
Cov(Yt,Yt?2)=Cov(et?θet?1,et?2?θet?3)=0(3)Cov(Y_t,Y_{t-2})=Cov(e_t-\theta e_{t-1}, e_{t-2}-\theta e_{t-3})=0\tag{3} Cov(Yt?,Yt?2?)=Cov(et??θet?1?,et?2??θet?3?)=0(3)
對于k≥2,Cov(Yt,Yt?k)=0k\geq2, Cov(Y_t,Y_{t-k})=0k≥2,Cov(Yt?,Yt?k?)=0,即過程大于1階滯后時,不存在自相關。
通過上面的信息,我們可以得到該MA(1)MA(1)MA(1)過程的1階自相關函數(shù)為:
ρ1=Cov(Yt,Yt?1)Var(Yt)=?θ1+θ2(4)\rho_1=\frac{Cov(Y_t,Y_{t-1})}{Var(Y_t)}=\frac{-\theta}{1+\theta^2}\tag{4} ρ1?=Var(Yt?)Cov(Yt?,Yt?1?)?=1+θ2?θ?(4)
顯然,若存在:
Yt=et?ηet?1(5)Y_t=e_t-\eta e_{t-1}\tag{5} Yt?=et??ηet?1?(5)
其中η=1/θ\eta=1/\thetaη=1/θ,則有1階自相關函數(shù):
ρ1=?η1+η2=?1/θ1+1/θ2=?θ1+θ2(6)\rho_1=\frac{-\eta}{1+\eta^2}=\frac{-1/\theta}{1+{1/\theta}^2}=\frac{-\theta}{1+\theta^2}\tag{6} ρ1?=1+η2?η?=1+1/θ2?1/θ?=1+θ2?θ?(6)
通過這個結果,我們看到θ\thetaθ被1/θ1/\theta1/θ代替時,自相關函數(shù)完全相同。
qqq階移動平均過程
對于q階移動平均過程:
Yt=et?θ1et?1?θ2et?2?...?θqet?q(7)Y_t=e_t-\theta_1 e_{t-1}-\theta_2 e_{t-2}-...-\theta_q e_{t-q}\tag{7} Yt?=et??θ1?et?1??θ2?et?2??...?θq?et?q?(7)
其中,et,et?1,et?2,...,et?qe_t, e_{t-1}, e_{t-2}, ..., e_{t-q}et?,et?1?,et?2?,...,et?q?,都是白噪聲,于是有:
E(Yt)=0γ0=Var(Yt)=(1+θ12+...+θq2)σe2γ1=Cov(Yt,Yt?1)=(?θ1+θ1θ2+θ2θ3+...+θq?1θq)σe2......γq?1=Cov(Yt,Yt?q+1)=(?θq?1+θ1θq)σe2γq=Cov(Yt,Yt?q)=?θqσe2E(Y_t)=0 \\\gamma_0=Var(Y_t)=(1+\theta_1^2+...+\theta_q^2)\sigma_e^2 \\\gamma_1=Cov(Y_t,Y_{t-1})=(-\theta_1+\theta_1 \theta_2 +\theta_2 \theta_3 +...+ \theta_{q-1} \theta_q)\sigma_e^2 \\......\\\gamma_{q-1}=Cov(Y_t,Y_{t-q+1})=(-\theta_{q-1}+\theta_1 \theta_q)\sigma_e^2 \\\gamma_{q}=Cov(Y_t,Y_{t-q})=-\theta_q \sigma_e^2 E(Yt?)=0γ0?=Var(Yt?)=(1+θ12?+...+θq2?)σe2?γ1?=Cov(Yt?,Yt?1?)=(?θ1?+θ1?θ2?+θ2?θ3?+...+θq?1?θq?)σe2?......γq?1?=Cov(Yt?,Yt?q+1?)=(?θq?1?+θ1?θq?)σe2?γq?=Cov(Yt?,Yt?q?)=?θq?σe2?
對于k>q,Cov(Yt,Yt?k)=0k>q, Cov(Y_t,Y_{t-k})=0k>q,Cov(Yt?,Yt?k?)=0,即過程大于q階滯后時,不存在自相關。則根據(jù)平穩(wěn)性條件,有限階的移動平均模型總是平穩(wěn)的。
可逆性
根據(jù)上一個Blog的證明,自回歸AR過程也可以被認為是一個無窮階移動平均過程.但是由于某些原因,自回歸表達式更便利。那么,移動平均模型可以被重新表示為自回歸模型嗎?
考慮MA(1)MA(1)MA(1)模型:
Yt=et?θet?1(8)Y_t=e_t-\theta e_{t-1}\tag{8} Yt?=et??θet?1?(8)
先把方程改寫成et=Yt+θet?1e_t=Y_t+\theta e_{t-1}et?=Yt?+θet?1?,再遞推可得:
et=Yt+θ(Yt?1+θet?2)=Yt+θYt?1+θ2et?2e_t=Y_t+\theta(Y_{t-1}+\theta e_{t-2})=Y_t + \theta Y_{t-1} + \theta^2 e_{t-2} et?=Yt?+θ(Yt?1?+θet?2?)=Yt?+θYt?1?+θ2et?2?
若∣θ∣<1|\theta|<1∣θ∣<1,我們可以對過去值無限重復以上的替代過程,得到表達式:
et=Yt+θYt?1+θ2Yt?2+...e_t=Y_t + \theta Y_{t-1} + \theta^2 Y_{t-2} + ... et?=Yt?+θYt?1?+θ2Yt?2?+...
或:
Yt=(?θYt?1?θ2Yt?2?...)+etY_t=(-\theta Y_{t-1} - \theta^2 Y_{t-2} -...) + e_t Yt?=(?θYt?1??θ2Yt?2??...)+et?
若∣θ∣<1|\theta|<1∣θ∣<1,我們可以看到MA(1)MA(1)MA(1)過程可以逆轉換成一個無窮階的自回歸模型,當且僅當∣θ∣<1|\theta|<1∣θ∣<1,我們稱MA(1)MA(1)MA(1)可逆。
對于一般的MA(q)MA(q)MA(q)模型,定義MA特征多項式為:
θ(x)=1?θ1x?θ2x2?...?θqxq\theta(x)=1-\theta_1 x -\theta_2 x^2 -...-\theta_q x^q θ(x)=1?θ1?x?θ2?x2?...?θq?xq
和相應的MA特征方程:
1?θ1x?θ2x2?...?θqxq=01-\theta_1 x -\theta_2 x^2 -...-\theta_q x^q=0 1?θ1?x?θ2?x2?...?θq?xq=0
可以證明MA(q)MA(q)MA(q)模型可逆,當且僅當MA特征方程根的模大于1.
在MA(1)MA(1)MA(1)過程中,我們看到θ\thetaθ被1/θ1/\theta1/θ代替時,會得到完全一樣的自相關函數(shù)ρ1\rho_1ρ1?。比如:Yt=et+2et?1Y_t=e_t+2 e_{t-1}Yt?=et?+2et?1?和Yt=et+0.5et?1Y_t=e_t+0.5 e_{t-1}Yt?=et?+0.5et?1?有相同的自相關函數(shù),但是只有第2個以特征根為-2的MA過程是可逆的。
后記:本次Blog暫無R語言應用,本系列未完待續(xù)…困dog我了
總結
以上是生活随笔為你收集整理的时间序列与R语言应用(part5)--移动平均MA模型及其可逆性的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 磊科 NW705P 无线路由器恢复出厂设
- 下一篇: 《弗洛伊德及其后继者》读书笔记(part