改变 R 中 svars 包生成的 ggplot 中的美学

Question

我正在使用svars 包生成一些 IRF 图。这些图是使用ggplot2 渲染的，但是我需要一些帮助来改变一些美学。

有什么方法可以更改阴影置信带的填充和 alpha 以及实线的颜色？我知道在ggplot2 中，您可以将fill 和alpha 参数传递给geom_ribbon（以及col 到geom_line），只是不确定如何在此包源的plot 函数中执行相同操作代码。

# Load Dataset and packages
library(tidyverse)
library(svars)
data(USA)

# Create SVAR Model
var.model <- vars::VAR(USA, lag.max = 10, ic = "AIC" )
svar.model <- id.chol(var.model)

# Wild Bootstrap

cores <- parallel::detectCores() - 1
boot.svar <- wild.boot(svar.model, n.ahead = 30, nboot = 500, nc = cores)

# Plot the IRFs

plot(boot.svar)

我还在查看历史分解图的命令（见下文）。有什么办法可以省略前两个方面，只在同一方面绘制底部三行？

hist.decomp <- hd(svar.model, series = 1)
plot(hist.decomp)

Answer 1

在调用plot 后重置aes_params 可以轻松实现您想要的第一个结果。为了你的第二个目标。可能有一种方法可以操纵ggplot 对象。相反，我下面的方法从头开始构建情节。基本上，我从vars:::plot.hd 复制并粘贴了数据整理代码，并为所需系列过滤了准备好的数据集：

# Plot the IRFs

p <- plot(boot.svar)

p$layers[[1]]$aes_params$fill <- "pink"
p$layers[[1]]$aes_params$alpha <- .5
p$layers[[2]]$aes_params$colour <- "green"
p

# Helper to convert to long dataframe. Source: svars:::plot.hd
hd2PlotData <- function(x) {
  PlotData <- as.data.frame(x$hidec)
  if (inherits(x$hidec, "ts")) {
    tsStructure = attr(x$hidec, which = "tsp")
    PlotData$Index <- seq(from = tsStructure[1], to = tsStructure[2], 
                          by = 1/tsStructure[3])
    PlotData$Index <- as.Date(yearmon(PlotData$Index))
  }
  else {
    PlotData$Index <- 1:nrow(PlotData)
    PlotData$V1 <- NULL
  }
  dat <- reshape2::melt(PlotData, id = "Index")
  dat
}

hist.decomp <- hd(svar.model, series = 1)

dat <- hd2PlotData(hist.decomp)

dat %>% 
  filter(grepl("^Cum", variable)) %>% 
  ggplot(aes(x = Index, y = value, color = variable)) + 
  geom_line() + 
  xlab("Time") + 
  theme_bw()

编辑 更改构面标签的一种方法是通过自定义labeller 函数。有关通过数据更改构面标签的不同方法，请参阅here：

myvec <- LETTERS[1:9]

mylabel <- function(labels, multi_line = TRUE) {
  data.frame(variable = labels) 
}
p + facet_wrap(~variable, labeller = my_labeller(my_labels))