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

★ Mann Whitney U検定の検出力 ★

　193　Mann Whitney U検定の検出力　　パリ　　2003/01/27 (月) 10:02
　　199　Re: Mann Whitney U検定の検出力　　　　　2003/01/27 (月) 20:20
　　194　Re: Mann Whitney U検定の検出力　　青木繁伸　　2003/01/27 (月) 13:49
　　　197　Re^2: Mann Whitney U検定の検出力　　パリ　　2003/01/27 (月) 17:44
　　　　200　Re^3: Mann Whitney U検定の検出力　　　　　2003/01/27 (月) 20:48
　　　　198　Re^3: Mann Whitney U検定の検出力　　　　　2003/01/27 (月) 20:11
　　　　　208　Re^4: Mann Whitney U検定の検出力　　パリ　　2003/01/28 (火) 10:13

193.　Mann Whitney U検定の検出力　　パリ　　2003/01/27 (月) 10:02

検出力についてどなたか教えてください。

牛の血中動態を調べるため，A，B群各7頭ずつに薬を投与しました。
Aは成牛，Bは仔牛で，投与した薬剤は同じです。

2群の最高血中濃度Cmaxの平均値はU検定で検定しました。
A(n=7):Cmax1.9,1.8,2.0,1.4,2.4,1.7,1.7 Average1.8, SD 0.3　
B(n=7):Cmax1.2,1.1,0.92,1.5,1.0,1.3,1.8 Average1.3, SD 0.3
この検出力と最小検出差の計算方法を教えてください。

　　　　　[このページのトップへ]

199.　Re: Mann Whitney U検定の検出力　　　　　2003/01/27 (月) 20:20

> A(n=7):Cmax1.9,1.8,2.0,1.4,2.4,1.7,1.7 Average1.8, SD 0.3　
> B(n=7):Cmax1.2,1.1,0.92,1.5,1.0,1.3,1.8 Average1.3, SD 0.3

有効桁数をけちりすぎていますね。
> x
[1] 1.9 1.8 2.0 1.4 2.4 1.7 1.7
> y
[1] 1.20 1.10 0.92 1.50 1.00 1.30 1.80
> mean(x)
[1] 1.842857
> sd(x)
[1] 0.3101459
> mean(y)
[1] 1.26
> sd(y)
[1] 0.3063767

　　　　　[このページのトップへ]

194.　Re: Mann Whitney U検定の検出力　　青木繁伸　　2003/01/27 (月) 13:49

> この検出力と最小検出差の計算方法を教えてください。

G*Power で計算できます。
以下を参照

http://www.psycho.uni-duesseldorf.de/aap/projects/gpower/reference/reference_manual_06.html
http://www.psycho.uni-duesseldorf.de/aap/projects/gpower/how_to_use_gpower.html

　　　　　[このページのトップへ]

197.　Re^2: Mann Whitney U検定の検出力　　パリ　　2003/01/27 (月) 17:44

> > この検出力と最小検出差の計算方法を教えてください。
>
> G*Power で計算できます。

ありがとうございました。
なんとか計算してみましたが，effct sizeが16.6098となりエラーが出てしまいます。計算は不可能なのでしょうか?
検出力と最小検出差はどのような検定においても求められるものなのですか?
その必要性があまりないので，参考書にあまり詳しく述べられないのか教えてください。

　　　　　[このページのトップへ]

200.　Re^3: Mann Whitney U検定の検出力　　　　　2003/01/27 (月) 20:48

> 検出力と最小検出差はどのような検定においても求められるものなのですか?
> その必要性があまりないので，参考書にあまり詳しく述べられないのか教えてください。

必要性があるかないかは，利用者個々が決めるべきものでしょう。
必要なら必要，不要なら不要。

　　　　　[このページのトップへ]

198.　Re^3: Mann Whitney U検定の検出力　　　　　2003/01/27 (月) 20:11

> なんとか計算してみましたが，effct sizeが16.6098となりエラーが出てしまいます。計算は不可能なのでしょうか?

うむ。エラー出ましたか?
t 検定でパワーを求めて，それを U test のパワーに変換するのかな。効率が0.95ということで。
私，やってみたら，effect size = 1.6667 になりまして，エラーにはならなかったけど。
G*Power 使わなくても，effect size=(1.8-1.3)/0.3=1.6666666なんですが（あなたの場合，両群ともSD=0.3でわかりやすかったわ!）。

　　　　　[このページのトップへ]

208.　Re^4: Mann Whitney U検定の検出力　　パリ　　2003/01/28 (火) 10:13

> 私，やってみたら，effect size = 1.6667 になりまして，エラーにはならなかったけど。
> G*Power 使わなくても，effect size=(1.8-1.3)/0.3=1.6666666なんですが（あなたの場合，両群ともSD=0.3でわかりやすかったわ!）。

ご指摘の通りです。計算ミスしてました。以後気をつけます。
今回の検出力は決して高くないんですが，どうしても知りたいっていう上の命令で。
私にはあまり必要性がないように思っています。

　　　　　[このページのトップへ]