Как вручную вычислить ci.type = 'ma' в ccf()

Question

Я пытался вычислить доверительный интервал для ci.type = 'ma' вручную. Я прочитал родственный пост (cross-correlation using ccf in R), который был полезным, но я получаю зацикливаться на деталях

getAnywhere(plot.acf) 

До сих пор у меня есть:

x <- ccf(x=x,y=y,ci.type='ma') 

clim0 <- qnorm((1+0.95)/2)/sqrt(x$n.used) # upper CI when ci.type=NULL 
clim <- clim0 + * sqrt(cumsum(c(1, 2 * x$acf[-1, 
               i, i]^2))) # from getAnywhere(plot.acf) 

Но я не понимаю, как я определен в функции подробности. Кто-нибудь делал это раньше?

Answer 1

вызова getAnywhere(plot.acf), я итеративно с 1 по длине вашей серии:

for (i in 1L:nser) { 
       clim <- clim0 * sqrt(cumsum(c(1, 2 * x$acf[-1, 
        i, i]^2))) 
       ylim <- range(c(-clim, clim, ylim)) 
      } 

однако, в функции ACF вы можете указать «ма» для ХИ с ci.type='ma' см https://stat.ethz.ch/R-manual/R-devel/library/stats/html/acf.html