R中ACF图中的虚线


9

我正在阅读Cowpertwait和Metcalfe的《 R入门时间序列》这本书。在第36页上,其行位于:1/n±2/n。我在这里阅读过R论坛,其行位于。 ±1.96/n

我运行了以下代码:

b = c(3,1,4,1)

acf(b)

并且我发现这些行看起来好像是。那么,显然这本书是错的吗?还是我误读了所写的内容?作者在谈论的内容略有不同吗?±1.96/4

*请注意,我对1.96对2的次要细节差异不感兴趣。我假设这只是作者使用2 sd与实际1.96 sd的经验法则。

编辑:我运行了此模拟:

acf1 = 0
acf2 = 0
acf3 = 0
for(i in 1:5000){
  resids= runif(1000)
  residsacf = c(acf(resids,plot= FALSE))
  acf1[i] = residsacf$acf[2,,1]
  acf2[i] = residsacf$acf[3,,1]
  acf3[i] = residsacf$acf[4,,1]
}
meanacf1 = mean(acf1)
meanacf2 = mean(acf2)
meanacf3 = mean(acf3)
meanacf1
meanacf2
meanacf3

我似乎总是对所有3 获得接近值。 1/n

进一步编辑:我看到趋势为1/n(k1)/n2


1
真, 1n±2n?集中于1n
mpiktas 2011年

在《恩德斯的应用经济时间序列》(第2版,第67-68页)中,2/ñ来自Box and Jenkins(1976),《时间序列预测,分析和控制》。恩德斯使用以下估计var[Rs
v一个[R[Rs=Ť-1个1个+2Ĵ=1个s-1个[RĴ2
恩德斯的用途 Ť作为系列的长度。
杰森·摩根

通常的极限是白噪声零假设下的临界值,在这种情况下,Enders中的方差表达式会崩溃为 1个/Ť
Rob Hyndman

时间序列分析中的 Shumway和Stoffer 及其应用:与R一起使用±2/ñ也一样 在此处查看其ACF代码。
詹森·摩根

Answers:


7

样本自相关为负偏差,并且第一个样本自相关系数具有均值 -1个/ñ 哪里 ñ是观察数。但是梅特卡夫(Metcalfe)和考珀特韦特(Cowpertwait)说所有自相关系数都具有均值是不正确的,并且说R在-1个/ñ±1.96/ñ

渐近平均数为0,这就是R在绘制直线处使用的 ±1.96/ñ


感谢您的回应Rob。在滞后1处ACF的期望为-1 / n时,我是否理解正确?如果是这样,那么在第一次滞后的情况下,虚线不应该居中吗?另外,由于他们似乎写的不是错字。您认为它们的含义有所不同还是完全错误?我去了他们的网站,但没有看到它被列为勘误表。
亚当

1
对于任何合理的样本量, 1个/ñ2/ñ所以没关系。我已经与安德鲁·梅特卡夫(Andrew Metcalfe)通讯,他承认有关R的错误。我想他们还没有更新勘误表。
罗布·海恩德曼

从技术上讲,R的缺陷和作者的假设R是否正确?
亚当

有两个问题。首先,均值-1 / n仅适用于第一个自相关函数,但作者说,它适用于所有相关函数。那是他们的错误,而不是R的错误。其次,R使用渐近结果(如我所见过的所有其他软件包一样)而不是小样本结果。因此,R没错,它只是使用可以改进的近似值。
罗布·海恩德曼

这类似于在分母中使用n而不是n-1计算样本方差吗?
亚当
By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.