绘制中间和最优 ARIMA(p,d,q) 的 AIC

Question

用优化算法找出最佳ARIMA(p,d,q)模型后，根据AIC的最小化，我也得到如下：

head():

  d p q       AIC
1 0 0 0 -2671.744
2 0 0 1 -2686.691
3 0 0 2 -2685.053
4 0 0 3 -2683.094
5 0 0 4 -2688.464
6 0 0 5 -2686.973

tail():

   d p q       AIC
31 0 5 0 -2685.982
32 0 5 1 -2684.926
33 0 5 2 -2683.003
34 0 5 3 -2683.654
35 0 5 4 -2686.069
36 0 5 5 -2698.871

在哪里，是的，d=0 无处不在。

我想知道是否有人知道如何绘制（3d 表面甚至 2d）不同的组合，以便可视化 AIC 方面的改进。以某种方式突出显示最大 AIC 可能也很酷。

Answer 1

不确定您到底想要什么类型的可视化，但热图似乎合适。这是ggplot2 中的一个。我将p 和q 制作成因子，然后使用geom_raster 制作瓷砖。

library(tidyverse)

d1 <- "dummy d p q       AIC
1 0 0 0 -2671.744
2 0 0 1 -2686.691
3 0 0 2 -2685.053
4 0 0 3 -2683.094
5 0 0 4 -2688.464
6 0 0 5 -2686.973" %>% read_table2()

d2 <- "dummy    d p q       AIC
31 0 5 0 -2685.982
32 0 5 1 -2684.926
33 0 5 2 -2683.003
34 0 5 3 -2683.654
35 0 5 4 -2686.069
36 0 5 5 -2698.871" %>% read_table2()

df <- bind_rows(d1, d2) %>% 
    select(-dummy) %>%
    mutate_at(vars(p, q), as.factor)

ggplot(df, aes(x = p, y = q, fill = AIC)) +
    geom_raster() +
    scale_fill_viridis_c(option = "A")

由reprex package (v0.2.0) 于 2018 年 4 月 30 日创建。