我有一个数据表,其中包含许多产品的消费量。我生成了每种产品的预测,现在想获得周期+1的均值和较高的80%。问题在于,根据所使用的方法,预测对象是具有不同结构的列表,因此我无法通过建立索引来检索值(我可以使用来命名data.table
)。
这是(虚拟)数据和代码:
# load required libraries
library(data.table)
library(xts)
library(forecast)
library(dplyr)
# create random data
set.seed(1)
a <- data.table(prod = sample(LETTERS[1:5], 20, TRUE), cons = sample(1:50, 20, TRUE), dt = sample(seq(as.Date("2016/06/01"), as.Date("2016/07/27"), by = "day"), 20, FALSE))
# create a time series of purchases
b <- a[, .(C=sum(cons)), by = .(dt, prod)][, x := .(list(xts(x = C, order.by = dt))), by = prod]
b <- b[, .SD[1,], by = prod]
# create a "reference" timeseries
dts <- xts(order.by = seq(as.Date("2016/06/01"), as.Date("2016/07/27"), by = "day"))
# merge reference and calculated timeseries, so zeros appear
b[, x2 := .(list(merge.xts(dts, x[[1]], fill = 0))), by = prod]
# calculate forecast for each extended timeseries
b[, fc := .(list(forecast(x2[[1]]))), by = prod]
现在,我要提取均值和较高水平。问题在于均值有时位于列表的插槽2中,有时不在列表中,因此我必须用名称来称呼它。在data.table
我这样做:
b[, mn := fc[[1]]$mean[1], by = prod]
b[, up := fc[[1]]$upper[1,1], by = prod]
但是,如果尝试在中执行相同的操作dplyr
,则会收到关闭错误:
b %>% mutate(mnD = .$fc[[1]]$mean[1])
## Error: invalid subscript type 'closure'
b %>% mutate(mnD = fc[[1]]$mean[1])
## Error: invalid subscript type 'closure'
我在做什么错,我该如何解决dplyr
?
与结合map_dbl()
从purrr
包将工作如下:
library(dplyr)
library(purrr)
b %>% as_data_frame() %>% mutate(mnD = map_dbl(fc, ~ .$mean[1]))
#> # A tibble: 5 x 7
#> prod dt C x x2 fc mnD
#> <chr> <date> <int> <list> <list> <list> <dbl>
#> 1 B 2016-07-17 47 <S3: xts> <S3: xts> <S3: forecast> 2.5241999
#> 2 C 2016-07-14 33 <S3: xts> <S3: xts> <S3: forecast> 1.1749266
#> 3 E 2016-06-30 7 <S3: xts> <S3: xts> <S3: forecast> 0.5952119
#> 4 D 2016-06-24 20 <S3: xts> <S3: xts> <S3: forecast> 3.3695962
#> 5 A 2016-07-04 18 <S3: xts> <S3: xts> <S3: forecast> 0.8421001
as_data_frame()
放在一边,不是必需的,而是添加它来以一种整洁的方式打印结果。没有它,列表列将打印所有数据。
本文收集自互联网,转载请注明来源。
如有侵权,请联系[email protected] 删除。
我来说两句