2012-08-09 77 views
0

有沒有一種方法可以使用ccf在R中生成一個ccf圖,但只允許它使用負滯後?我試圖改變lag.max,但它仍然保持對稱的形狀。定製滯後acf/ccf在R

編輯:爲了澄清,我需要的是比被繪製的標準的默認(我覺得默認的是+ -34)

+0

從'ccf':默認的數滯後爲10 * log10(N/m),其中N是觀測數,m是系列數。最大滯後數比數據中的時間點數少一個。 – Seth 2012-08-10 16:37:50

+0

那麼如何改變默認的滯後數? – Alex 2012-08-10 20:56:56

+0

您可以像我在下面的回答中一樣,使用'lag.max = 200'更改滯後數。 – Seth 2012-08-10 20:58:36

回答

3

我也跟着在CCF一個例子較大的滯後?我創建然後ccf對象看其結構,我不使用plot=F繪製它尚未:

>xccf=ccf(mdeaths,fdeaths,ylab="wtf!",lag.max=200, plot=F) 
>str(xccf) 

str條命令表明,有一個有分量的呼叫滯後是某種形式的序列,看它與:

> xccf$lag 
, , 1 

       [,1] 
    [1,] -5.91666667 
    [2,] -5.83333333 
    [3,] -5.75000000 
    [4,] -5.66666667 
    [5,] -5.58333333 
...  ... ... 
[141,] 5.75000000 
[142,] 5.83333333 
[143,] 5.91666667 

的滯後的長度將是lag.max或觀察在較短的數據,無論哪個是更低的數目。對於這個示例數據,我將max.lag設置爲200,但由於沒有足夠的數據進行計算,因此只能得到143.如果只需要負數,那麼前71就是它們。所以我把這些延伸到方括號中,以便將它們分組。然後打電話給事情的陰謀。

圖(xccf [xccf $滯後[1:71]])

這給出了這樣的情節:

OOHH aahhhh!

+0

對不起,我應該是具體的,我需要滯後於標準範圍之外..所以我需要更多比說34滯後或任何它自動 – Alex 2012-08-09 22:58:29