R 如何从容器函数中提取conf并添加到ggplot2中

发布于03月05日

对于此数据:

big = runif(60);mangum = runif(61:120)
df = data.frame(big, mangum)

我们可以绘制:

ccf(df$big, df$mangum) ## we see 2 horizontal dashed lines

现在我想通过ggplot2绘制，如何将两条水平虚线(它们确定哪些相关值是重要的)添加到ggplot2图中

    df_x <- eval(substitute(big), df)
    df_y <- eval(substitute(mangum), df)
    ccf.object <- ccf(df_x, df_y, plot = FALSE)
    output_table <- cbind(lag = ccf.object$lag, 
                x.corr = ccf.object$acf) %>%
    as_tibble() %>%      mutate(cat = ifelse(x.corr > 0, "green", "red"))
    output_table %>%    ggplot(aes(x = lag, y = x.corr)) +
    geom_bar(stat = "identity", aes(fill = cat))

推荐答案

如果查看stats:::plot.acf内部，您会看到这些行的默认公式基本上是c(-1, 1) * qnorm((1 + 0.95)/2)/sqrt(n.used)，其中n.used是从ccf的输出中获得的.因此，要获得与ccf的输出匹配的ggplot，您可以执行以下操作:

output_table %>%    
  ggplot(aes(x = lag, y = x.corr)) +
  geom_col(aes(fill = paste0(cat, 3)), width = 0.4) +
  geom_hline(yintercept = c(-1, 1) * qnorm((1 + 0.95)/2)/sqrt(
    ccf(df$big, df$mangum, plot = FALSE)$n.used), 
    linetype = 2, color = 'blue') +
  geom_hline(yintercept = 0) +
  scale_fill_identity() +
  theme_classic(base_size = 14) +
  theme(axis.line = element_blank(),
        panel.border = element_rect(fill = NA, color = 'black'))