Що таке автокореляція для випадкової прогулянки?


11

Здається, це дійсно високо, але це для мене контрсуєтивне. Може хтось, будь ласка, пояснить? Я дуже розгублений цим питанням і буду вдячний за детальне, проникливе пояснення. Заздалегідь дякую!

Відповіді:


11

(Я написав це як відповідь на інший пост, який був позначений як дублікат цього, коли я складав його; я зрозумів, що опублікую його тут, а не викидаю. Схоже, це говорить зовсім схоже на те, що було відповідь, але достатньо різний, що хтось може щось із цього отримати.)

Випадкова прогулянка має виглядyt=i=1tϵi

Зауважте, щоyt=yt1+ϵt

Отже .Cov(yt,yt1)=Cov(yt1+ϵt,yt1)=Var(yt1)

Також зауважте, щоσt2=Var(yt)=tσϵ2

Отже .corr(yt,yt1)=σt12σt1σt=σt1σt=t1t=11t112t

Що означає, що ви повинні побачити співвідношення майже 1, оскільки як тільки починає великі , і це майже те саме - відносна різниця між ними, як правило, досить мала.tytyt1

Ви можете переконатись у тому, що vs .ytyt1

введіть тут опис зображення

Тепер ми можемо бачити це дещо інтуїтивно - уявіть, що знизився до (як ми бачимо, це було в моєму моделюванні випадкової прогулянки зі стандартним нормальним рівнем шуму). Тоді буде досить близьким до ; це може бути або може бути але майже напевно, що він знаходиться в межах кількох одиниць . Отже, коли серія рухається вгору і вниз, сюжет vs майже завжди буде знаходитися в досить вузькому діапазоні лінії ... все ж, коли зростає, точки будуть охоплювати більше і більший розтяг уздовж цього - 20 y t - 20 - 22 - 18,5 - 20 y t y t - 1 y = x t y = x yt120yt202218.520ytyt1y=xty=xлінія (спред по лінії збільшується з , але вертикальний спред залишається приблизно постійним); кореляція повинна підходити 1.t


9

(x0,x1,x2,,xn)(x0,x1,,xn1)(x1,x2,,xn).

xi+1xixix0n(xi,xi+1)y=x±1y=x±11(n/2)2=n/4R2

R211n/4=14n.

n=1000R214/n

Малюнок


Ось Rкод, який створив зображення.

set.seed(17)
n <- 1e3
x <- cumsum((runif(n) <= 1/2)*2-1)          # Binomial random walk at x_0=0
rho <- format(cor(x[-1], x[-n]), digits=3)  # Lag-1 correlation

par(mfrow=c(1,2))
plot(x, type="l", col="#e0e0e0", main="Sample Path")
points(x, pch=16, cex=0.75,  col=hsv(1:n/n, .8, .8, .2))
plot(x[-n], x[-1], asp=1, pch=16, col=hsv(1:n/n, .8, .8, .2),
     main="Lag-1 Scatterplot",
     xlab="Current value", ylab="Next value")
mtext(bquote(rho == .(rho)))
Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.