[R-es] Intervalos de confianza en una autocorreloción

[Rubén Gómez Antolí]

Hola:

El 18/05/13 23:38, Carlos Ortega escribió:
> Hola,
>
> Aquí tienes un ejemplo de cómo se calculan (busca la primera entrada de
> "confidence interval"):

Muchas gracias Carlos, ya casi lo tenía con la pista que me dio Jorge 
(gracias también).

>
> http://zoonek2.free.fr/UNIX/48_R/15.html

Ya te quedaría agradecido si me dijeras que búsqueda has hecho porque os 
prometo que he invertido más de una hora buscando sin éxito.

La función ha quedado finalmente así:

niv.conf.acf<-function (x,ci= 0.95, ci.type = c("white","ma"),
   ...)
{
   # Niveles de confianza de la función de autocorrelación.

   ci.type <- match.arg(ci.type)
   if ((nser <- ncol(x$lag)) < 1L)
       stop("x$lag must have at least 1 column")
   if (ci.type != "ma" && ci > 0 && x$type != "covariance") {
       clim <- qnorm((1+ci)/2)/sqrt(x$n.used)
       return(clim*c(-1,1))
   } else {
       stop("Can't use ci.type=\"ma\" in this function")
   }
}

Sin embargo, observo que, utilizando un objeto acf, este no «transmite» 
ni el argumento «ci» ni «ci.type».

Y, otra duda, esta ya más de «teoría» (algo en lo que estoy flojo). 
Estoy intentando replicar la tabla de la página 26 de este artículo:

http://hoeger.com.ve/ama/pdf/taller-cambiosclimaticos-02.pdf

Entiendo que por valores tabulados se refiere a los intervalos de 
confianza, pero quizás estoy equivocado, ya que no me coincide el 
resultado con los que se reflejan la tabla:

dput(colotti.sep)
c(23.9, 24.4, 24.1, 24.5, 24.4, 24.6, 24.4, 24.9, 24.9, 24.2,
24.7, 26, 25.1, 25.5, 25.3, 25.5, 24.6, 25.7, 25.4, 24.8, 25,
25.3, 25.4, 25.3, 24.9, 24.7, 24.6, 25, 25.3, 25.2, 25.8, 25.3,
25.4, 25.5, 24.6)
niv.conf.acf(acf(colotti.sep,ci=.95,plot=F))
[1] -0.3312944  0.3312944

De nuevo gracias por todo.

Salud y Revolución.

Lobo.
-- 
Libertad es poder elegir en cualquier momento. Ahora yo elijo GNU/Linux,
para no atar mis manos con las cadenas del soft propietario.
Porque la libertad no es tu derecho, es tu responsabilidad.
http://www.mucharuina.com
---------
Desde El Ejido, en Almería, usuario registrado Linux #294013
http://www.counter.li.org