尽管尝试为半小时时间序列数据创建ARIMA模型,但auto.arima不会自动包括季节性,尽管明显(每天)季节性。当在调用auto.arima时指定D = 1时(按照here的指示),ARIMA模型返回以下错误:
Error in if (length(dx) > ncol(dxreg)) lm(dx ~ dxreg - 1, na.action = na.omit) else list(rank = 0L) :
argument is of length zero
在下面提供了一个最小的可复制示例
structure(c(860.110666666667, 0, 0, 0, 0, 0, 0, 0, 1795.942,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1233.68133333333, 0, 0, 0, 0, 2460.29066666667, 1910.628, 3026.91,
1313.402, 587.111666666667, 1389.328, 0, 425.473666666667, 2676.30433333333,
359.064, 598.866666666667, 0, 0, 0, 0, 0, 0, 826.639333333333,
0, 0, 0, 0, 0, 525.482333333333, 600.823333333333, 0, 1293.36533333333,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 383.521333333333,
0, 0, 0, 0, 0, 0, 739.561, 3241.267, 2820.05166666667, 3262.343,
1779.022, 1406.778, 631.806666666667, 3261.663, 25.2043333333333,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1120.446, 931.430666666667, 329.675,
167.672, 2957.59733333333, 0, 0, 0, 0, 0, 0, 486.451666666667,
1210.21233333333, 0, 1192.211, 0, 0, 0, 0, 0, 0, 1899.427, 0,
0, 0, 0, 0, 0, 717.466333333333, 0, 0, 0, 2421.93233333333, 3231.43766666667,
1160.77066666667, 2315.53966666667, 3252.93833333333, 3210.40633333333,
1805.097, 1445.225), .Tsp = c(1505352600, 1505352600.00816, 17520
), class = "ts")
aa <- forecast::auto.arima(dfts, D = 1)
答案 0 :(得分:1)
您定义了一个长度为144的时间序列,其频率等于17520,然后要求一个季节性差异-即相隔17520个时间段的观测值之间的差异-无法做到。
该错误消息不是很有帮助-我会解决该问题。