文档介绍：如果您需要使用本文档，请点击下载按钮下载!
如果您需要使用本文档，请点击下载按钮下载!
如果您需要使用本文档，请点击下载按钮下载!
时间序列回归模型
干预分析
概念及模型
Box和Tiao引入的干预分析提供了对于干预影响时间序列的效果进行评估的一个框架，假设干预是可以通过时间序列的均值函数或者趋势而对过程施加影响，干预可以自然产生也可以人为施加的，如国家的宏观调控等。
其模型可以如下表示：
其中mt代表均值的变化，Nt是ARIMA过程。
干预的分类
阶梯响应干预
脉冲响应干预
如果您需要使用本文档，请点击下载按钮下载!
如果您需要使用本文档，请点击下载按钮下载!
如果您需要使用本文档，请点击下载按钮下载!
干预的实例分析
模型初探
对数化航空客运里程的干预模型的估计
> data(airmiles)
> acf((diff(diff(window(log(airmiles),end=c(2001,8)),12))),=48)#用window得到在911事件以前的未爱干预的时间序列子集
对暂用的模型进行诊断
>fitmode<-arima(airmiles,order=c(0,1,1),seasonal=list(order=c(0,1,0)))
> tsdiag(fitmode)
如果您需要使用本文档，请点击下载按钮下载!
如果您需要使用本文档，请点击下载按钮下载!
如果您需要使用本文档，请点击下载按钮下载!
从诊断图可以看出存在三个异常点，acf在12阶存在高度相关因此在季节中加入MA（1）系数。
拟合带有干预信息的模型
函数：
arimax(x, order = c(0, 0, 0), seasonal = list(order = c(0, 0, 0), period = NA),
xreg = NULL, = TRUE, = TRUE, fixed = NULL,
init = NULL, method = c("CSS-ML", "ML", "CSS"), , = list(),
kappa = 1e+06, io = NULL, xtransf, transfer = NULL)
arimax函数扩展了arima函数，可以处理时间序列中干扰分析及异常值。假设干扰影响过程的均值，相对未受干扰的无价值函数的偏离用一些协变量的ARMA滤波器的输出这种来表示，偏差被称作传递函数。。
=arimax(log(airmiles),order=c(0,1,1),seasonal=list(order=c(0,1,1),
period=12),xtransf=(I911=1*(seq(airmiles)==69),
I911=1*(seq(airmiles)==69)),
transfer=list(c(0,0),c(1,0)),xreg=(Dec96=1*(seq(airmiles)==12),
Jan97=1*(seq(airmiles)==13),Dec02=1*(seq(airmiles)==84)),method='ML')
>
Call:
arimax(x = log(airmiles), order = c(0, 1, 1), seasonal = list(order = c(0, 1,
1), period = 12), xreg = (Dec96 = 1 * (seq(airmiles) == 12), Jan97 = 1 *
如果您需要使用本文档，请点击下载按钮下载!
如果您需要使用本文档，请点击下载按钮下载!
如果您需要使用本文档，请点击下载按钮下载!
(seq(airmiles) == 13), Dec02 = 1 * (seq(airmiles) == 84)), method = "ML",
xtransf = (I911 = 1 * (seq(airmiles) == 69), I911 = 1 * (seq(airmiles) ==
69)), transfer = list(c(0, 0), c(1, 0)))
Coeff