均值模型的基本构件
p 阶自回归 AR(p) 用自身的过去解释当前值:xt=c+φ1xt−1+⋯+φpxt−p+εt。当系数使序列不会发散时它是平稳的,对 AR(1) 即 ∣φ∣<1,其均值为 μ=c/(1−∑φi)。q 阶移动平均 MA(q) 用过去的冲击解释当前值:xt=μ+εt+θ1εt−1+⋯+θqεt−q。它总是平稳的,记忆恰好是 q 步。两者结合得到 ARMA(p,q),往往比单纯的 AR 或 MA 用更少的参数就能描述一个序列。
μ=1−φc,Var(xt)=1−φ2σ2,x^T+h=μ+φh(xT−μ);ρMA(1)(1)=1+θ2θ,AIC=2k−2lnL,BIC=klnT−2lnL(1) 参数用最大似然估计。阶数通过比较信息准则来选择:AIC =2k−2lnL 与 BIC =klnT−2lnL,其中 k 是参数个数,L 是似然值;越小越好,对 k 的惩罚能防止加入只拟合噪声的项。当 T>7 时 BIC 的惩罚更重,因此倾向选择更小的模型。拟合之后残差应像白噪声,用 Ljung-Box 检验检查。平稳 AR(1) 的预测按几何速度回到均值,x^T+h=μ+φh(xT−μ),MA(q) 的预测从 q+1 步起等于均值。任何平稳 ARMA 的远期预测都只是均值,所以这类模型对接下来的几步有用,之后并不比平均值更好。
AR(1) 的预测设
xt=0.2+0.8xt−1+εt,冲击标准差为
0.6。均值
μ=0.2/(1−0.8)=1.0,方差
0.36/(1−0.64)=1.0。若今日值为
xT=2,一步预测为
1+0.8⋅1=1.8,两步预测为
1+0.64⋅1=1.64。与均值的差距每步缩小为
0.8 倍,十步之后预测为
1.11,几乎就是均值。
- AR(p) 依赖自身的过去值,MA(q) 依赖最近 q 个冲击,ARMA 两者都依赖。
- AIC 与 BIC 奖励拟合、惩罚参数;取值最小的模型胜出,BIC 更严格。
- 平稳 ARMA 的预测随着预测期变长收敛到均值。
实验:评分内容
- AR(1) 的均值15 分
- 两步预测15 分
- MA(1) 的自相关15 分
- 信息准则15 分
- 非平稳的系数15 分
- 配对术语15 分
- 为什么远期预测会到达均值10 分