たこぜりー研究室

大腸菌DNAは4.64Mbp。酵母は13Mbpで、ヒトは3Gbp

スポンサーサイト

#
上記の広告は1ヶ月以上更新のないブログに表示されています。
新しい記事を書く事で広告が消せます。

自己相関係数

#271
自己相関係数 autocorrelation coeffient 。周期検定に使えそうだ。
そして、それ用の R の関数がこれ。

[Auto- and Cross- Covariance and -Correlation Function]
acf(x, lag.max = NULL, type = c("correlation", "covariance", "partial"), plot = TRUE, na.action = na.fail, demean = TRUE, ...)

で、プロットに関してはこれを参照。
Plot Autocovariance and Autocorrelation Function

プロットだから引数が多いが、このうち ci 引数と type 引数をいじってみる。ci は coverage probability for confidence interval 。既定では、青い線ででるみたい。

> plot(acf(c(...), 70), ci=.99, type="l")

そうすると、こんな感じ。

青梅のハエ

これを Excel でやろうとすると大変なのよ。その分理解しやすいんだけどね。

コメント


コメントの投稿


管理者にだけ表示を許可する

トラックバック

トラックバックURL:  http://takojelly.blog2.fc2.com/tb.php/271-619a42d0
上記広告は1ヶ月以上更新のないブログに表示されています。新しい記事を書くことで広告を消せます。