我正在尝试将 tbats 模型应用于时间序列数据。此数据框包含 到 范围内的0值1。在应用模型之前，我用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()在这个数据集中工作得很好。