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

Question

我正在使用这样的 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

谢谢

Answer 1

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