たかが棄却域の描画に,複雑なこと,無駄なことをするプログラムが多い。
以下のように簡単に書けばよい。
![](https://blogimg.goo.ne.jp/user_image/2f/0d/05dc3b9998c8232929fe51db9b413baf.png)
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)
以下のように簡単に書けばよい。
![](https://blogimg.goo.ne.jp/user_image/2f/0d/05dc3b9998c8232929fe51db9b413baf.png)
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)