Rでは、時系列のACF相関ログラムを実行でき、信頼区間バンドは薄い青でプロットされます。しかし、私がACFオブジェクトの構造を引っ張ると、私はこれらの値を見つけることができません。誰も信頼区間バンドの値を抽出する方法を知っていますか？ACFからの信頼区間値の抽出
List of 6
$ acf : num [1:27, 1, 1] 1 0.06453 -0.06354 0.00213 -0.01324 ...
$ type : chr "correlation"
$ n.used: int 501
$ lag : num [1:27, 1, 1] 0 1 2 3 4 5 6 7 8 9 ...
$ series: chr "tser[i:(i + 500)]"
$ snames: NULL
- attr(*, "class")= chr "acf"
私は、機能を見て持っていたと私は信頼区間を抽出する簡単な方法を見ることはできません。領域はplot.acf関数で計算されます。この機能を表示するには、この機能で
getS3method("plot", "acf")
を使用し、変数climがあり、これはあなたが後にしているものです。最も簡単な方法はplot.acfをmyplot.acfにコピーすることですが、climの値を返します。
ありがとうございます。あまりにも簡単な作業ではありませんが、他の方法が見つからない場合は、コードを修正する必要があります。どのように正確に修正するかについての簡単な投稿は高く評価されます。 – pat
@pat 'return'文をコメントアウトし、代わりに' clim'を返します – csgillespie
これは答えとして受け入れられるべきです。関数をプルするコードを以下に示します： 'dump（" plot.acf "、file =" function_dump_plot_acf.txt "、envir = asNamespace（" stats "））'。また、 'invisible（）'を 'return（clim）'に置き換える必要があります。 – user1700890
私はこの質問が超古いものだと知っていますが、信頼区間の値が必要な人は、信頼水準のz値を使用した観測数のsqrtで割った値にすぎません。 plot.acf関数でこれがここで計算される：with.ci論理値（ユーザが信頼区間をプロットしたいかどうかとciが所望の信頼レベルである場合、例えば0.95、0.9、等示す情報である
clim0 <- if (with.ci)
qnorm((1 + ci)/2)/sqrt(x$n.used)
...）
EDIT：これは訂正がある場合ではない場合は、ホワイトノイズ遅れの値であると仮定し、あなたがいくつかを読むことができる
clim <- clim0 * sqrt(cumsum(c(1, 2 * x$acf[-1, i, j]^2)))
を適用することができます場合、これは信頼区間でありますその詳細についてhere
どのようにACFをプロットしますか？どのパッケージが訴えていますか？ – agstudy
統計パッケージに組み込まれていると思います。 http://stat.ethz.ch/R-manual/R-patched/library/stats/html/acf.html – pat