R auto.arima 包版本 3.22 在solve.default(res$hessian * n.used) 中出错答案

【问题标题】：R auto.arima package version 3.22 Error in solve.default(res$hessian * n.used)R auto.arima 包版本 3.22 在solve.default(res$hessian * n.used) 中出错
【发布时间】：2012-07-08 18:56:56
【问题描述】：

我正在使用这样的 auto.arima()：

F12.Xreg <- fourier(NEW.JMB,12)
F4.Xreg<- fourier(NEW.JMB,4)
fb<-cbind(F12.Xreg, F4.Xreg)
fit <-auto.arima(NEW.JMB, D=0, max.P=0, max.Q=0, xreg=fb)

我不知道如何避免错误：

Error in solve.default(res$hessian * n.used) : 
 system is computationally singular: reciprocal condition number = 6.64873e-30
Error in if (diffs == 1 & constant) { : argument is of length zero
In addition: Warning message:
In auto.arima(NEW.JMB, D = 0, max.P = 0, max.Q = 0, xreg = fb) :
  Unable to calculate AIC offset

谢谢

【问题讨论】：

标签： r forecasting

【解决方案1】：

两个傅立叶矩阵具有相同的列，因此设计矩阵不是满秩的。您可以通过一次调用 fourier() 来创建傅立叶项。

【讨论】：

嗨 Rob，我对两个傅立叶矩阵fourier(x, 7) 和fourier(x,365) 有同样的问题。当我像上面那样将它们 cbind 在一起时，我可以看到奇异值都是非零的，所以 xreg 仍然是满秩的。知道为什么仍然会出现错误吗？
啊抱歉，实际上是因为我在 xreg 中包含了另一个外部变量。奇怪的是奇异值仍然非零