タイプIIエラーの確率を見つけるにはどうすればよいですか?


12

タイプIIエラーはH1が真であることがわかっていますが、H0は拒否されません。

質問

標準偏差がわかっている正規分布を含むタイプIIエラーの確率を計算するにはどうすればよいですか?



1
私はこの質問を「対などの一般的なテストの力をどのように見つけるのですか?」と言い換えます。多くの場合、これはより頻繁に実行されるテストです。私はそのようなテストの力をどのように計算するのかわかりません。H0:μ=μ0H1:μ>μ0
確率

回答:


27

(タイプIエラーの確率)を指定することに加えて、完全に指定された仮説ペアが必要です。つまり、、、およびが既知である必要があります。(タイプIIエラーの確率)はです。片側を想定してい。R:αμ0μ1σβ1powerH1:μ1>μ0

> sigma <- 15    # theoretical standard deviation
> mu0   <- 100   # expected value under H0
> mu1   <- 130   # expected value under H1
> alpha <- 0.05  # probability of type I error

# critical value for a level alpha test
> crit <- qnorm(1-alpha, mu0, sigma)

# power: probability for values > critical value under H1
> (pow <- pnorm(crit, mu1, sigma, lower.tail=FALSE))
[1] 0.63876

# probability for type II error: 1 - power
> (beta <- 1-pow)
[1] 0.36124

編集:視覚化

ここに画像の説明を入力してください

xLims <- c(50, 180)
left  <- seq(xLims[1],   crit, length.out=100)
right <- seq(crit, xLims[2],   length.out=100)
yH0r  <- dnorm(right, mu0, sigma)
yH1l  <- dnorm(left,  mu1, sigma)
yH1r  <- dnorm(right, mu1, sigma)

curve(dnorm(x, mu0, sigma), xlim=xLims, lwd=2, col="red", xlab="x", ylab="density",
      main="Normal distribution under H0 and H1", ylim=c(0, 0.03), xaxs="i")
curve(dnorm(x, mu1, sigma), lwd=2, col="blue", add=TRUE)
polygon(c(right, rev(right)), c(yH0r, numeric(length(right))), border=NA,
        col=rgb(1, 0.3, 0.3, 0.6))
polygon(c(left,  rev(left)),  c(yH1l, numeric(length(left))),  border=NA,
        col=rgb(0.3, 0.3, 1, 0.6))
polygon(c(right, rev(right)), c(yH1r, numeric(length(right))), border=NA,
        density=5, lty=2, lwd=2, angle=45, col="darkgray")
abline(v=crit, lty=1, lwd=3, col="red")
text(crit+1,  0.03,  adj=0, label="critical value")
text(mu0-10,  0.025, adj=1, label="distribution under H0")
text(mu1+10,  0.025, adj=0, label="distribution under H1")
text(crit+8,  0.01,  adj=0, label="power", cex=1.3)
text(crit-12, 0.004,  expression(beta),  cex=1.3)
text(crit+5,  0.0015, expression(alpha), cex=1.3)

1
この回答にタイプミスはありますか?と呼ばれるものは実際にはあり、逆もまた同様です。どちらにしても、これは優れたグラフであり、Rコードの例です!powβ
jdods 2018年

1
@jdods確かに、lower.tail=FALSE行方不明がありました。どうもありがとうございました!
カラカル

@caracalは、〜layman用語で、ベータを考慮せずにp値(タイプ1エラーのリスク)を計算できる理由を説明できますが、タイプ2エラーのリスクを測定できるようにアルファを指定する必要がありますか?何かが足りないような気がします。すばらしい回答をありがとう。
Cystack

1
@Cystack p値、タイプ1エラー、タイプ2エラーの正確な意味は、コメントで伝えられる範囲を超えています。私はstats.stackexchange.com/q/46856/1909stats.stackexchange.com/q/129628/1909などの質問への回答を見ていきます。また、右上隅の[リンクされた]ボックスと[関連した]ボックスも参照してくださいより関連性の高いコンテンツ。
カラカル

1

カラカルの答えを補足するために、質問で暗示されているものを含む多くの一般的な設計のタイプIIエラー率または電力を計算するためのユーザーフレンドリーなGUIオプションを探している場合は、無料ソフトウェアのG Power 3をチェックしてください。

弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.