2015-10-27 47 views
0

我想检查一下“鸟类”在不同时滞下的“狼”之间是否存在相关性。获取相关值很容易,但我该如何解决滞后问题(I需要检查1:4滞后的校正值)?我寻找的输出是一个包含滞后值和相关相关值的数据表。如何获得两个滞后变量之间的相关性

df <- read.table(text = " day birds wolfs  
          0 2  21   
          1 8   4 
          2 2   5 
          3 2   4 
          4 3   6 
          5 1   12 
          6 7   10 
          7 1   9 
          8 2   12 header = TRUE) 

输出(不是真正的结果): 滞后CorValue

0  0.9 
1  0.8 
2  0.7 
3  0.9 

回答

1

如果你这样做:

corLag<-ccf(df$birds,df$wolfs,lag.max=max(df$day)) 

它会返回此:

系列“X的

自相关',滞后

-8  -7  -6  -5  -4  -3  -2  -1  0  1  2  3  4  5  6  7  8 
-0.028 0.123 -0.045 -0.019 0.145 -0.176 -0.082 -0.126 -0.296 0.757 -0.134 -0.180 0.070 -0.272 0.549 -0.170 -0.117 

第一行是滞后,第二行是相关值。你可以检查cor(df$birds,df$wolfs)确实等于-0.296

+0

你好@etienne,有没有办法限制滞后数?例如我有1k时间点,但我想检查滞后1:10的相关性?我试过:corLag $ ccf [1:10]但没有成功。 – mql4beginner

+1

是的:'lag.max = 10'而不是'max(df $ day)'。它返回的滞后从-10到+10 – etienne

+1

然后你可以使用'corLag [0:10]'有相关性和相关性从1到10的差距 – etienne

相关问题