0

我正在尝试将 tbats 模型应用于时间序列数据。此数据框包含 到 范围内的01。在应用模型之前,我用NA's 替换了所有 0 值(因为我想获取数据集的日志并且 0 值的日志将给出 -inf),然后获取所有值的日志以准备数据帧TBATS

完成上述所有更改后,数据框如下所示:

     y_ds0           timestamp
1    -2.9957323 2017-08-07 01:01:00
2    -3.2188758 2017-08-07 01:01:10
3    -3.2188758 2017-08-07 01:01:20
4    -3.5065579 2017-08-07 01:01:30
5    -3.5065579 2017-08-07 01:01:40
6    -3.9120230 2017-08-07 01:01:50
7    -3.9120230 2017-08-07 01:02:00
8    -4.6051702 2017-08-07 01:02:10
9    -4.6051702 2017-08-07 01:02:20
10   -4.6051702 2017-08-07 01:02:30
11   -4.6051702 2017-08-07 01:02:40
12   -4.6051702 2017-08-07 01:02:50
13           NA 2017-08-07 01:03:00
14           NA 2017-08-07 01:03:10
15           NA 2017-08-07 01:03:20
16           NA 2017-08-07 01:03:30
17           NA 2017-08-07 01:03:40
18           NA 2017-08-07 01:03:50
19           NA 2017-08-07 01:04:00
20           NA 2017-08-07 01:04:10
21           NA 2017-08-07 01:04:20
22           NA 2017-08-07 01:04:30
23           NA 2017-08-07 01:04:40
24           NA 2017-08-07 01:04:50
25           NA 2017-08-07 01:05:00
26   -2.5257286 2017-08-07 01:05:10
27   -2.6592600 2017-08-07 01:05:20
28   -2.8134107 2017-08-07 01:05:30
29   -2.9957323 2017-08-07 01:05:40
30   -3.2188758 2017-08-07 01:05:50
31   -3.5065579 2017-08-07 01:06:00
32   -3.5065579 2017-08-07 01:06:10
33   -3.9120230 2017-08-07 01:06:20
34   -3.9120230 2017-08-07 01:06:30
35   -3.9120230 2017-08-07 01:06:40
36   -4.6051702 2017-08-07 01:06:50
37   -4.6051702 2017-08-07 01:07:00

它在数据中有 1000 多行。

在上面的数据框中,我TBATS使用以下代码应用了模型:

tk_ts_dataframe <- tk_ts(ts_df ,
            start = 1,
            freq = 6,
            silent = TRUE)

fit <- tbats(tk_ts_dataframe)

forecast <- forecast(fit, h = 1)

我收到如下错误:

属性错误(best.model$errors)<-属性(origy):暗淡[产品1435]与对象的长度不匹配[236]

请帮我解决这个问题。auto.arima()在这个数据集中工作得很好。

4

0 回答 0