有沒有一種方法可以使用ccf
在R中生成一個ccf圖,但只允許它使用負滯後?我試圖改變lag.max
,但它仍然保持對稱的形狀。定製滯後acf/ccf在R
編輯:爲了澄清,我需要的是比被繪製的標準的默認(我覺得默認的是+ -34)
有沒有一種方法可以使用ccf
在R中生成一個ccf圖,但只允許它使用負滯後?我試圖改變lag.max
,但它仍然保持對稱的形狀。定製滯後acf/ccf在R
編輯:爲了澄清,我需要的是比被繪製的標準的默認(我覺得默認的是+ -34)
我也跟着在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]])
這給出了這樣的情節:
對不起,我應該是具體的,我需要滯後於標準範圍之外..所以我需要更多比說34滯後或任何它自動 – Alex 2012-08-09 22:58:29
從'ccf':默認的數滯後爲10 * log10(N/m),其中N是觀測數,m是系列數。最大滯後數比數據中的時間點數少一個。 – Seth 2012-08-10 16:37:50
那麼如何改變默認的滯後數? – Alex 2012-08-10 20:56:56
您可以像我在下面的回答中一樣,使用'lag.max = 200'更改滯後數。 – Seth 2012-08-10 20:58:36