R 从ACF相关图中提取置信区间值

在R中，我们可以运行时间序列的ACF相关图，置信区间带将以浅蓝色绘制。但是当我拉动ACF对象的结构时，我找不到这些值。有人知道如何提取置信区间带的值吗

e、 g

---

我已经查看了函数，但找不到提取置信区间的简单方法。该区域在

plot.acf

函数中计算。要查看此函数，请使用

getS3method("plot", "acf")

---

在该函数中，有一个变量

clim

，这就是您要查找的变量。最简单的方法是将

plot.acf

复制到

myplot.acf

，但返回

clim

值。

我知道这个问题非常古老，但如果有人想要置信区间值，它只是置信水平的z值除以所用观察数的sqrt。在

plot.acf

函数中，计算如下：

clim0 <- if (with.ci) 
    qnorm((1 + ci)/2)/sqrt(x$n.used)

---

您可以阅读更多关于这方面的内容

如何绘制acf？你在起诉哪个软件包？我想它是内置在统计软件包中的。谢谢你。这看起来不太容易，但如果我找不到其他方法。。。我将不得不尝试修改代码。任何关于如何准确修改的快速帖子都将不胜感激。@pat只需注释掉

return

语句，而不是return

clim

这应该被接受为答案。下面是拉动函数的代码：

dump（“plot.acf”，file=“function\u dump\u plot\u acf.txt”，envir=asNamespace（“stats”）

。还需要将

不可见（）

替换为

返回（clim）

如果有人编辑此答案以清楚显示用户1700890建议的代码，这将非常有用，包括在哪里替换

不可见（）

以及

返回（clim）

什么？这个问题可能非常古老，但这个答案也非常有用！因此，置信区间就是$\pm\frac{z{a/2}}{\sqrt{n}}}$，其中$a$是显著性水平，$n$是数据的大小。

clim0 <- if (with.ci) 
    qnorm((1 + ci)/2)/sqrt(x$n.used)

clim <- clim0 * sqrt(cumsum(c(1, 2 * x$acf[-1, i, j]^2)))