1
我正在寻求两个特定的帮助点 1)如何创建一个列表给定我的数据库(all.df)的列表 2)如何矢量化这个列表上的功能如何创建列表并向其执行矢量化函数
我试图在使用Prophet库的客户/产品级别生成预测。我正在努力引导操作。 我目前运行一个for循环,我想避免并加速我的计算。
数据进行了分析
set.seed(1123)
df1 <- data.frame(
Date = seq(dmy("01/01/2017"), by = "day", length.out = 365*2),
Customer = "a",
Product = "xxx",
Revenue = sample(1:100, 365*2, replace=TRUE))
df2 <- data.frame(
Date = seq(dmy("01/01/2017"), by = "day", length.out = 365*2),
Customer = "a",
Product = "yyy",
Revenue = sample(25:200, 365*2, replace=TRUE))
df3 <- data.frame(
Date = seq(dmy("01/01/2017"), by = "day", length.out = 365*2),
Customer = "b",
Product = "xxx",
Revenue = sample(1:100, 365*2, replace=TRUE))
df4 <- data.frame(
Date = seq(dmy("01/01/2017"), by = "day", length.out = 365*2),
Customer = "b",
Product = "yyy",
Revenue = sample(25:200, 365*2, replace=TRUE))
all.df <- rbind(df1, df2, df3, df4)
这是我的预测功能
daily_forecast <- function(df, forecast.days = 365){
# fit actuals into prophet
m <- prophet(df,
yearly.seasonality = TRUE,
weekly.seasonality = TRUE,
changepoint.prior.scale = 0.55) # default value is 0.05
# create dummy data frame to hold prodictions
future <- make_future_dataframe(m, periods = forecast.days, freq = "day")
# run the prediction
forecast <- predict(m, future)
### Select the date and forecast from the model and then merge with actuals
daily_fcast <- forecast %>% select(ds, yhat) %>% dplyr::rename(Date = ds, fcast.daily = yhat)
actual.to.merge <- df %>% dplyr::rename(Date = ds, Actual.Revenue = y)
daily_fcast <- merge(actual.to.merge, daily_fcast, all = TRUE)
}
目前,我在同一时间通过一个客户/产品配合使用for循环
x <- df1 %>% select(-c(Customer, Product)) %>%
dplyr::rename(ds = Date, y = Revenue) %>%
daily_forecast()
我想不是,vectorise整个操作:
1 - 创建名单列表,即通过拆分all.df:
一)产品,然后
二)客户
2,然后请有上述
在列表1中创建的列表中daily_forecast功能图)我非常喜欢使用功能出purrr
。
谢谢!爱你的第二个解决方案,它非常优雅,并有自然流动。我希望我可以多次投票答复这个答案! – cephalopod