裏 RjpWiki

文字通り,RjpWiki の裏を行きます
R プログラム コンピュータ・サイエンス 統計学

棄却域の描画

2011年05月13日 | ブログラミング
たかが棄却域の描画に,複雑なこと,無駄なことをするプログラムが多い。
以下のように簡単に書けばよい。

rejection.region <- function(df, s0, s.max=4, length=100, col="#0000aa20")
{
x <- seq(-s0, s0, length=length)
plot(c(x, s0, -s0), c(dt(x, df), 0, 0), xlim=c(-s.max, s.max), type="l", xlab="t", ylab="f(t)")
x <- seq(s0, s.max, length=length)
y <- c(dt(x, df), 0, 0)
x <- c(x, s.max, s0)
polygon(x, y, col=col)
polygon(-x, y, col=col)
}
rejection.region(25, qt(0.975, 25), 4)
コメント