Linee tratteggiate nel diagramma ACF in R


9

Sto sfogliando il libro "Introductory Time Series with R" di Cowpertwait e Metcalfe. A pagina 36 Si dice che le linee sono: . Ho letto quiR forumche le linee sono a±1,96/-1/n±2/n . ±1.96/n

Ho eseguito il seguente codice:

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

acf(b)

e vedo che le linee sembrano apparire a . Quindi, ovviamente, il libro è sbagliato? Oppure sto leggendo male ciò che è stato scritto? Gli autori parlano di qualcosa di leggermente diverso?±1.96/4

* Nota, non sono interessato alla discrepanza minore tra 1,96 e 2. Presumo che questo fosse solo l'autore che utilizzava la regola empirica di 2 sd contro l'attuale 1,96 sd.

Modifica: ho eseguito questa simulazione:

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

Mi sembra sempre di ottenere valori vicini a per tutti e 3. 1/n

Ulteriore modifica: sto vedendo una tendenza di 1/n-(K-1)/n2


1
Davvero, ? Centrato su-1-1n±2n ? -1n
mpiktas,

In Enders ' Applied Economic Time Series (2a edizione, pp 67-68) spiega che il proviene da Box e Jenkins (1976),Previsione, analisi e controllo delle serie storiche. Enders utilizzava la seguente stima divar(rs):var(rs)=T - 1 ( 1+2 s - 1 j = 1 r 2 j ) . Enders usaTcome lunghezza della serie. 2/Nvun'r(rS)
vun'r(rS)=T-1(1+2Σj=1S-1rj2).
T
Jason Morgan,

1/T

±2/N

Risposte:


7

-1/nn-1/n±1.96/n

±1.96/n


Grazie per la risposta Rob. Sono corretto nel comprendere che l'aspettativa dell'ACF in ritardo 1 è -1 / n? In tal caso, le linee tratteggiate non dovrebbero essere centrate lì per il primo ritardo? Inoltre, poiché sembra che ciò che hanno scritto non fosse un errore di battitura. Pensi che significhino qualcosa di diverso o semplicemente sbagliano? Sono andato sul loro sito Web e non lo vedo elencato come errata.
Adam,

1
1/n2/n

Tecnicamente, il difetto di R non è corretto e l'ipotesi degli autori R era corretta?
Adam,

Ci sono due problemi Innanzitutto, la media di -1 / n si applica solo alla prima funzione di autocorrelazione, ma gli autori affermano che si applica a tutte le funzioni di correlazione. Questo è il loro errore, non quello di R. In secondo luogo, R utilizza il risultato asintotico (come ogni altro pacchetto software che ho visto) anziché il piccolo risultato di esempio. Quindi R non è sbagliato, utilizza solo un'approssimazione che può essere migliorata.
Rob Hyndman,

è analogo al calcolo della varianza del campione usando n nel denominatore anziché n-1?
Adam,
Utilizzando il nostro sito, riconosci di aver letto e compreso le nostre Informativa sui cookie e Informativa sulla privacy.
Licensed under cc by-sa 3.0 with attribution required.