思路:
问题1:如何计算模型的方差(以AR(1)模型为例)
yt=ϕ1yt−1+εty_{t}=\phi _{1}y_{t-1}+\varepsilon _{t}yt=ϕ1yt−1+εt
则Green函数为
Gj=ϕ1j,j=0,1,...G_{j}=\phi_1^j,j=0,1,...Gj=ϕ1j,j=0,1,...
故AR(1)模型的方差为
Var(yt)=∑j=0∞Gj2Var(εt−j)=∑j=0∞(ϕ1j)2σϵ2)=σϵ21−ϕ1jVar(y_{t})=\sum_{j=0}^{\infty }G_{j}^2Var(\varepsilon _{t-j})=\sum_{j=0}^{\infty }(\phi_1^j)^2\sigma_\epsilon ^{2})=\frac{\sigma_\epsilon ^{2}}{1-\phi_1^j}Var(yt)=∑j=0∞Gj2Var(εt−j)=∑j=0∞(ϕ1j)2σϵ2)=1−ϕ1jσϵ2
注:
AR(p)模型
yt=ϕ1yt−1+ϕ2yt−2+...+ϕpyt−p+εty_{t}=\phi _{1}y_{t-1}+\phi _{2}y_{t-2}+...+\phi _{p}y_{t-p}+\varepsilon _{t}yt=ϕ1yt−1+ϕ2yt−2+...+ϕpyt−p+εt
用滞后算子可表示为:
ϕ(B)yt=εt\phi(B)y_{t}=\varepsilon _{t}ϕ(B)yt=εt
其中ϕ(B)=1−ϕ1B−ϕ2B2−...−ϕpBp\phi(B)=1-\phi _{1}B-\phi _{2}B^{2}-...-\phi _{p}B^{p}ϕ(B)=1−ϕ1B−ϕ2B2−...−ϕpBp
有因式分解可知,ϕ(B)=∏i=1p(1−λiB)\phi(B)=\prod_{i=1}^{p}(1-\lambda _{i}B)ϕ(B)=∏i=1p(1−λiB)
所以,yt=εt∏i=1p(1−λiB)=∑i=1pεtki1−λiBy_{t}=\frac{\varepsilon _{t}}{\prod_{i=1}^{p}(1-\lambda _{i}B)}=\sum_{i=1}^{p}\frac{\varepsilon _{t}k_{i}}{1-\lambda _{i}B }yt=∏i=1p(1−λiB)εt=∑i=1p1−λiBεtki
式中,ki(i=1,2,3,4,...,p)为常数k_{i}(i=1,2,3,4,...,p)为常数ki(i=1,2,3,4,...,p)为常数
这里需要引入Green函数的定义
设λ1,λ2,λ3...λp\lambda _{1},\lambda _{2},\lambda _{3}...\lambda _{p}λ1,λ2,λ3...λp为平稳AR(p)模型的特征根,则平稳AR(p)模型可以写成:
yt=εtϕ(B)=∑i=1pεtki1−λiB=∑i=1p∑j=0∞ki(λiB)jεt=∑j=0∞∑i=1pkiλijBjεt=∑j=0∞∑i=1pkiλijεt−j=∑j=0∞Gjεt−jy_{t}=\frac{\varepsilon _{t}}{\phi (B)}=\sum_{i=1}^{p}\frac{\varepsilon _{t}k_{i}}{1-\lambda _{i}B }=\sum_{i=1}^{p}\sum_{j=0}^{\infty }k_{i}(\lambda _{i}B)^{j}\varepsilon _{t}=\sum_{j=0}^{\infty }\sum_{i=1}^{p}k_{i}\lambda _{i}^{j}B^{j}\varepsilon _{t}=\sum_{j=0}^{\infty }\sum_{i=1}^{p}k_{i}\lambda _{i}^{j}\varepsilon _{t-j}=\sum_{j=0}^{\infty }G_{j}\varepsilon _{t-j}yt=ϕ(B)εt=∑i=1p1−λiBεtki=∑i=1p∑j=0∞ki(λiB)jεt=∑j=0∞∑i=1pkiλijBjεt=∑j=0∞∑i=1pkiλijεt−j=∑j=0∞Gjεt−j
式中,G0=1,Gj=∑i=1pkiλij(j=1,2....)G_{0}=1,G_{j}=\sum_{i=1}^{p}k_{i}\lambda _{i}^{j}(j=1,2....)G0=1,Gj=∑i=1pkiλij(j=1,2....)
本文详细介绍了AR(1)模型的方差计算方法,通过Green函数阐述了其计算过程,展示了从AR(p)模型到Green函数的转换,并给出了特征根和系数的关系,为理解时间序列分析中的方差计算提供了清晰的理论基础。
595

被折叠的 条评论
为什么被折叠?



