返回信息流下面一段拟合带异常值的ARIMA模型代码：
model=arimax(log(airmiles),order=c(0,1,1),seasonal=list(order=c(0,1,1),period=12),
xtransf=data.frame(I911=1*(seq(airmiles)==69)),transfer=list(c(1,1)),
xreg=data.frame(Dec96=1*(seq(airmiles)==12),Jan97=1*(seq(airmiles)==13),Dec02=1*(seq(airmiles)==84)),
method='ML')
为神马会报错啊：
错误于optim(init[mask], armafn, method = "BFGS", hessian = TRUE, control = optim.control, :
non-finite finite-difference value [6]
但是如果把xtransf改成
xtransf=data.frame(I911=1*(seq(airmiles)==69),I911=1*(seq(airmiles)==69)),transfer=list(c(0,0),c(1,0))
就完全正常！但是对应着模型的差别只不过是一个很简单的分式通分求和啊！肿么会这样！
以上程序需要加载TSA程序包。内容来自《时间序列分析与应用，R语言》（Time Series Analysis and Applications in R, 2nd ed.）。
这是一条镜像帖。来源：北邮人论坛 / matlab / #9033同步于 2012/5/14
Matlab机器人发帖
有木有会R的童鞋啊~
moiraphael
2012/5/14镜像同步0 回复
订阅后，新回复会通过你的通知中心匿名送达。
0 条回复
暂无回复 · 你可以订阅本帖等待新回复。