打印显着的自相关值

Question

如果我在R（acf）中进行自相关测试，我得到一个很好的图形，水平线显示重要性的截止值。

acf也打印出控制台中的各个滞后值，但是，在这里我无法看到哪些是重要的。没有查看图表，有没有一种简单的方法可以做到这一点？

Answer 1

所以基本上我们需要知道截止值。通过编写acf和stats:::plot.acf，您可以看到它对于不同的参数值可能会有所不同，但是对于默认值，您应该使用以下内容：

set.seed(123)
x <- arima.sim(list(ar = 0.5), 100)
r <- acf(x, plot = FALSE)$acf
which(abs(r)[-1] >= qnorm(1 - 0.05 / 2) / sqrt(length(x)))
# [1]  1  2  3  9 10 12 13

其中0.05是这种情况下的显着性水平。