dianwang

1个月预测

echart <- function(x1, x2, data){
  p <- data %>% 
    group_by(`业扩类型`) %>% 
    e_charts(`年月`, timeline = TRUE) %>% 
    e_line(`lstm`) %>% 
    e_line(`prophet`) %>% 
    e_line(`ingarch`) %>%
    e_legend(right = 0, show = T) %>%  # 设置图例位置
    e_tooltip(trigger = "axis") %>%    # 设置鼠标悬浮提示
    e_title(paste(x1, x2))
  p
}

详细版

ii <- 0
for(i in unique(b$一级行业)){
  for(j in unique(b$`申请/接电`)){
    x1 <- paste(j,"1 month")
    pic_temp <- paste0('pic.', ii)
    temp <- echart(i,x1,filter(b,`一级行业` == i&`申请/接电` == j))
    assign(pic_temp, temp)
    ii <- ii+1
  }
}
e_arrange(pic.0, pic.1, pic.2, pic.3,pic.4,pic.5,pic.6,pic.7,pic.8,pic.9,pic.10,pic.11,pic.12,pic.13,pic.14,pic.15,pic.16,pic.17,pic.18,pic.19,pic.20,pic.21,pic.22,pic.23, cols = 1, rows = 24)

均值版

不区分业扩类型及行业

 filter(b, prophet != Inf & lstm != Inf & ingarch != Inf)%>%
  group_by(年月) %>% reframe(
  `年月` = `年月`[1],
  `prophet` = mean(`prophet`,na.rm = TRUE),
  `lstm` = mean(`lstm`,na.rm = TRUE),
  `ingarch` = mean(`ingarch`,na.rm = TRUE)
) %>% 
  e_charts(`年月`) %>% 
  e_line(`prophet`) %>% 
  e_line(`lstm`) %>% 
  e_line(`ingarch`) %>%
  #e_line(sGARCH) %>%
  e_legend(right = 0, show = T) %>%  # 设置图例位置
  e_tooltip(trigger = "axis") %>%    # 设置鼠标悬浮提示
  e_title("一个月中各模型的预测均值")

3个月预测误差

详细版

ii <- 0
for(i in unique(a$一级行业)){
  for(j in unique(a$`申请/接电`)){
    x1 <- paste(j,"3 month")
    pic_temp <- paste0('pic.', ii)
    temp <- echart(i,x1,filter(a,`一级行业` == i&`申请/接电` == j))
    assign(pic_temp, temp)
    ii <- ii+1
  }
}
e_arrange(pic.0, pic.1, pic.2, pic.3,pic.4,pic.5,pic.6,pic.7,pic.8,pic.9,pic.10,pic.11,pic.12,pic.13,pic.14,pic.15,pic.16,pic.17,pic.18,pic.19,pic.20,pic.21,pic.22,pic.23, cols = 1, rows = 24)

均值版

不区分业扩类型及行业

 filter(a, prophet != Inf & lstm != Inf & ingarch != Inf)%>%
  group_by(年月) %>% reframe(
  `年月` = `年月`[1],
  `prophet` = mean(`prophet`,na.rm = TRUE),
  `lstm` = mean(`lstm`,na.rm = TRUE),
  `ingarch` = mean(`ingarch`,na.rm = TRUE)
) %>% 
  e_charts(`年月`) %>% 
  e_line(`prophet`) %>% 
  e_line(`lstm`) %>% 
  e_line(`ingarch`) %>%
  #e_line(sGARCH) %>%
  e_legend(right = 0, show = T) %>%  # 设置图例位置
  e_tooltip(trigger = "axis") %>%    # 设置鼠标悬浮提示
  e_title("三个月中各模型的预测均值")