Comment est calculé l'intervalle de confiance pour la fonction ACF?


12

Par exemple, dans R si vous appelez la acf()fonction, elle trace un corrélogramme par défaut et trace un intervalle de confiance à 95%. En regardant le code, si vous appelez plot(acf_object, ci.type="white"), vous voyez:

qnorm((1 + ci)/2)/sqrt(x$n.used)

comme limite supérieure pour le bruit blanc de type. Quelqu'un peut-il expliquer la théorie derrière cette méthode? Pourquoi obtenons-nous la qnorm de 1 + 0,95 puis divisons par 2 et après cela, divisons par le nombre d'observations?


1
FWIW, il ne s'agit pas vraiment de R.
gung - Reinstate Monica

Réponses:


9

Dans Analysis of Time Series de Chatfield (1980), il donne un certain nombre de méthodes d'estimation de la fonction d'autocovariance, y compris la méthode du jack-knife. Il note également qu'il peut être montré que la variance du coefficient d'autocorrélation au décalage k, , est normalement distribuée à la limite, et que Var ( ) ~ 1 / N (où N est le nombre d'observations). Ces deux observations sont à peu près au cœur du problème. Il ne donne pas de dérivation pour la première observation, mais fait référence à Kendall & Stuart, The Advanced Theory of Statistics (1966). rkrk

Maintenant, nous voulons α / 2 dans les deux queues, pour le test à deux queues, donc nous voulons le quantile 1 − α / 2.

Voyez ensuite que (1 + 1 − α) / 2 = 1 − α / 2 et multipliez par l'écart-type (c.-à-d. Racine carrée de la variance comme indiqué ci-dessus)


1
Bonne réponse, mais cela ne traite pas de la partie de la question de savoir pourquoi (0,95 + 1) / 2 (ou toute autre valeur ci prend); ce n'est pas vraiment suffisant pour une réponse séparée, je pense, donc je vais le mentionner ici: C'est simplement parce que nous voulons dans les deux queues, donc nous voulons le quantile . puis voyez que . (Robert: si vous voulez incorporer quelque chose dans ce sens dans votre réponse, allez-y)1 - α / 2 ( 1 + 1 - α ) / 2 = 1 - α / 2α/21α/2(1+1α)/2=1α/2
Glen_b -Reinstate Monica

Bonjour Glen_b, j'ai essayé de mettre à jour ma réponse en conséquence - j'ai également emprunté certains de vos mots, ce qui, je l'espère, va bien.
Robert de Graaf
En utilisant notre site, vous reconnaissez avoir lu et compris notre politique liée aux cookies et notre politique de confidentialité.
Licensed under cc by-sa 3.0 with attribution required.