ノンパラメトリック検定のサンプルサイズ計算

にほんブログ村 科学ブログ 数学へ

ノンパラメトリック検定のサンプルサイズ計算はどうやったらいいのか?

統計ソフトRではどうやるか調べた。

質問

ノンパラメトリック検定のサンプルサイズ計算をしたい。

どうやったらいいか?

回答

統計ソフトRの”samplesize"というパッケージで、 Mann-Whitney U検定のサンプルサイズ計算ができる。

samplesize パッケージのn.wilcox.ord()。

samplesizeパッケージをインストールして、 library()で呼び出す。

install.packages("samplesize")
library(samplesize)

検出力80%、有意水準5%、 両群同数で1:1(t=0.5)、 グループ1は33%、33%、34%と 三分割したグループが同数に対し、 グループ2は、66%、20%、14%だとすると、 サンプルサイズは36例ずつ必要と計算される。

> n.wilcox.ord(power=0.8, alpha=0.05, t=0.5, p=c(0.33, 0.33, 0.34), q=c(0.66, 0.20, 0.14))
$`total sample size`
[1] 72

$m
[1] 36

$n
[1] 36