Я изучаю книгу «Вводные временные ряды с R» Каупертвейта и Меткалфа. На странице 36 написано, что строки: . Я прочитал здесь R форум, что строки в .
Я запустил следующий код:
b = c(3,1,4,1)
acf(b)
и я вижу, что строки выглядят как . Итак, очевидно, что книга не так? Или я неправильно читаю то, что было написано? Авторы говорят о чем-то немного другом?
* Обратите внимание, меня не интересует несоответствие 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.
Дальнейшее редактирование: я вижу тенденцию
r
time-series
Адам
источник
источник
Ответы:
Автокорреляция выборки имеет отрицательное смещение, и коэффициент автокорреляции первого образца имеет среднее значение где - количество наблюдений. Но Меткалф и Каупертвейт неверны, говоря, что все коэффициенты автокорреляции имеют это среднее значение, и они также неверны, говоря, что R строит линии при .n - 1 / n ± 1,96 / √- 1 / н N - 1 / n ± 1,96 / n--√
Асимптотически среднее значение равно 0, и именно это R использует при построении линий в .± 1,96 / н--√
источник