打印重要的自相关值
Print significant auto-correlation value
如果我在 R (acf) 中进行自相关检验，我会得到一个很好的图表，水平线显示显着性的截止值。
acf 还在控制台中打印出各个滞后值，但是，在这里我看不出哪些是重要的。有没有不看图表的简单方法？
所以基本上为此我们需要知道截止值。通过写 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 是显着性水平。
如果我在 R (acf) 中进行自相关检验，我会得到一个很好的图表，水平线显示显着性的截止值。
acf 还在控制台中打印出各个滞后值，但是，在这里我看不出哪些是重要的。有没有不看图表的简单方法？
所以基本上为此我们需要知道截止值。通过写 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 是显着性水平。