Calcular la correlación de rodadura utilizando rollapply
Tengo objeto zoológico con más de 10000 filas.
> head(tt)
A B
2007-01-04 0.005945924 0.0021167475
2007-01-05 -0.004201991 -0.0080020024
2007-01-08 0.001740897 0.0045804104
2007-01-09 0.000000000 -0.0008163931
2007-01-10 -0.004503531 0.0032615812
2007-01-11 -0.005841138 0.0043863282
He intentado variaciones de la siguiente línea, pero en vano.
rollapply(tt, 21, function(x) cor(x[,1],x[,2]))
Cada entrada dio una correlación de 1, parece que está tomando el 1 de la diagonal de la matriz de correlación.
2013-11-25 1 1
2013-11-26 1 1
2013-11-27 1 1
2013-11-29 1 1
2013-12-02 1 1
2013-12-03 1 1
Lo que realmente quiero es -0.4649, como el siguiente
> cor(tt)
A B
A 1.0000000 -0.4649881
B -0.4649881 1.0000000