Buenas noches Andres,

Dale una mirada al siguiente ejemplo:

R> set.seed(123)
R> x <- rnorm(100)
R> d <- acf(x, type = "correlation", plot = F)
R> d

Autocorrelations of series ‘x’, by lag

     0      1      2      3      4      5      6      7      8      9     10
    11     12
 1.000 -0.026 -0.113  0.149 -0.094 -0.013  0.029  0.013 -0.025 -0.076 -0.027
 0.095 -0.184
    13     14     15     16     17     18     19     20
-0.148  0.160 -0.086 -0.027  0.140 -0.067 -0.062 -0.072
R> str(d)
List of 6
 $ acf   : num [1:21, 1, 1] 1 -0.0256 -0.1127 0.149 -0.0943 ...
 $ type  : chr "correlation"
 $ n.used: int 100
 $ lag   : num [1:21, 1, 1] 0 1 2 3 4 5 6 7 8 9 ...
 $ series: chr "x"
 $ snames: NULL
 - attr(*, "class")= chr "acf"
R> d$acf    # esto es un 'array'
, , 1

             [,1]
 [1,]  1.00000000
 [2,] -0.02559634
 [3,] -0.11265356
 [4,]  0.14898981
 [5,] -0.09427814
 [6,] -0.01325316
 [7,]  0.02944197
 [8,]  0.01314468
 [9,] -0.02524275
[10,] -0.07647498
[11,] -0.02749780
[12,]  0.09527532
[13,] -0.18418382
[14,] -0.14812903
[15,]  0.15993700
[16,] -0.08576017
[17,] -0.02679542
[18,]  0.14019189
[19,] -0.06662434
[20,] -0.06184345
[21,] -0.07206891

R> d$acf[,,1]
[1]  1.00000000 -0.02559634 -0.11265356  0.14898981 -0.09427814 -0.01325316
 0.02944197
 [8]  0.01314468 -0.02524275 -0.07647498 -0.02749780  0.09527532 -0.18418382
-0.14812903
[15]  0.15993700 -0.08576017 -0.02679542  0.14019189 -0.06662434 -0.06184345
-0.07206891

Saludos,
Jorge Ivan Velez


2010/4/26 Andres rodriguez trujillo <andreswk25@hotmail.com>

>
>
> Buen Dia
>
> Ando trabajando con el comando acf del paquete stats que por medio de ella
> obtengo la funcion de autocorrelacion y utilizo el ejemplo que esta sugerido
> en la ayuda del R la base de datos se llama "lh" y ps necesito cada una de
> las autocorrelaciones que me arroja el acf para poderlas multiplicar con un
> programa que ando realizando y pues lo que habia hecho era transformar la
> salida como un vector pues me facilitaria realizar las operaciones que deseo
> aunque lo consegui sucedia lo mismo antes de que lo llevara en forma de
> vector, por ultimo he mirando el comando interno del acf creo que lo arroja
> en forma de lista aun asi no entiendo mucho los comandos que utiliza esta
> funcion acf
>
>  d<-acf(lh,type="correlation",plot=F)
> > d
>
> Autocorrelations of series ‘lh’, by lag
>
>     0      1      2      3      4      5      6      7      8      9     10
>  1.000  0.576  0.182 -0.145 -0.175 -0.150 -0.021 -0.020 -0.004 -0.136
> -0.154
>    11     12     13     14     15     16
> -0.097  0.049  0.120  0.087  0.119  0.151
> > d[1]
>
> Autocorrelations of series ‘lh’, by lag
>
>    1
> 0.576
> > d[1]*2
> Error en d[1] * 2 : argumento no-numérico para operador binario
>
>
> ahora si lo transformo como un vector sucede esto
>
> d<-acf(lh,type="correlation",plot=F)
> > vector<-c(d)
> > vector
> $acf
> , , 1
>
>              [,1]
>  [1,]  1.000000000
>  [2,]  0.575524476
>  [3,]  0.181818182
>  [4,] -0.144755245
>  [5,] -0.174825175
>  [6,] -0.149650350
>  [7,] -0.020979021
>  [8,] -0.020279720
>  [9,] -0.004195804
> [10,] -0.135664336
> [11,] -0.153846154
> [12,] -0.097202797
> [13,]  0.048951049
> [14,]  0.119580420
> [15,]  0.086713287
> [16,]  0.118881119
> [17,]  0.151048951
>
>
> $type
> [1] "correlation"
>
> $n.used
> [1] 48
>
> $lag
> , , 1
>
>      [,1]
>  [1,]    0
>  [2,]    1
>  [3,]    2
>  [4,]    3
>  [5,]    4
>  [6,]    5
>  [7,]    6
>  [8,]    7
>  [9,]    8
> [10,]    9
> [11,]   10
> [12,]   11
> [13,]   12
> [14,]   13
> [15,]   14
> [16,]   15
> [17,]   16
>
>
> $series
> [1] "lh"
>
> $snames
> NULL
>
> > w<-vector[1]
> > w
> $acf
> , , 1
>
>              [,1]
>  [1,]  1.000000000
>  [2,]  0.575524476
>  [3,]  0.181818182
>  [4,] -0.144755245
>  [5,] -0.174825175
>  [6,] -0.149650350
>  [7,] -0.020979021
>  [8,] -0.020279720
>  [9,] -0.004195804
> [10,] -0.135664336
> [11,] -0.153846154
> [12,] -0.097202797
> [13,]  0.048951049
> [14,]  0.119580420
> [15,]  0.086713287
> [16,]  0.118881119
> [17,]  0.151048951
>
>
> > w*1
> Error en w * 1 : argumento no-numérico para operador binario
> >
>
> Estaria agradecido por la ayuda.
>
>
>
> _________________________________________________________________
> Explore the seven wonders of the world
>
> BRE
>        [[alternative HTML version deleted]]
>
>
> _______________________________________________
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
>

	[[alternative HTML version deleted]]

