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?

Nick Nikolaev
la source
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)

Robert de Graaf
la source
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