Rを使用してccfプロットを生成する方法はありますが、ccf
負のラグのみを使用できるようにしますか? 変更しlag.max
てみましたが、まだ対称的な形状を保持しています。
編集:明確にするために、プロットされる標準のデフォルトよりも大きいラグが必要です(デフォルトは+-34だと思います)
Rを使用してccfプロットを生成する方法はありますが、ccf
負のラグのみを使用できるようにしますか? 変更しlag.max
てみましたが、まだ対称的な形状を保持しています。
編集:明確にするために、プロットされる標準のデフォルトよりも大きいラグが必要です(デフォルトは+-34だと思います)
?ccf の例の 1 つに従いました。オブジェクトを作成して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 です。そのため、角括弧にラグを渡し、それらにサブセット化します。次に、プロットを呼び出します。
plot(xccf[xccf$lag[1:71],])
これにより、次のプロットが得られます。