2015-10-27 45 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

相關問題