「統計学関連なんでもあり」の過去ログ--- 046

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

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