No.20132 Rを用いた検出力の求め方  【赤岳】 2013/08/09(Fri) 13:06

いつもありがとうございます。
さて,今回,正規分布をN(μ,σ^2)としたとき,帰無仮説はN(0, 0.51^2)の分布をしているとして,その分布の右すそが1.96×0.51(=1.0)以上の値を取る確率を,対立仮説下のN(0.693, 0.54^2)の分布で求めたいのですが,Rでどのようにプログラムをかけばよいでしょうか。どなたかご教示ください。

No.20133 Re: Rを用いた検出力の求め方  【青木繁伸】 2013/08/09(Fri) 21:27

どの部分ですか?赤のハッチがかかっているところ?それとも黒のハッチがかかっている(クロスハッチのように見えるけど)ところ。
いやいや,それ以外だよということ?(記述が正確でないので,識別しがたい)


No.20135 Re: Rを用いた検出力の求め方  【赤岳】 2013/08/09(Fri) 22:46

青木先生,いつもありがとうございます。
求めたい部分は,まさに赤のハッチがかかっている部分です。
よろしくお願いします。

No.20136 Re: Rを用いた検出力の求め方  【青木繁伸】 2013/08/09(Fri) 23:00

赤の正規分布は N(0.693, 0.54^2) なのですから,求める確率は pnorm(1, mean=0.693, sd=0.54, lower.tail=FALSE) でしょう。
標 準正規分布 N(0, 1) で 1.96 より大きい値を取る確率は pnorm(1.96, lower.tail=FALSE) で求めているでしょう?それは pnorm(1.96, mean=0, sd=1, lower.tail=FALSE) と同じなんですよ。

No.20137 Re: Rを用いた検出力の求め方  【赤岳】 2013/08/09(Fri) 23:32

青木先生,よくわかりました。
標準正規分布での求め方はよく目にしていて,ご指摘のように省略形で書かれていたので応用ができず,悩んでいました。
ありがとうございました。

● 「統計学関連なんでもあり」の過去ログ--- 046 の目次へジャンプ
● 「統計学関連なんでもあり」の目次へジャンプ
● 直前のページへ戻る