Links

6/15の日記

Posted on June 15, 2015

統計学入門

p.255の検出力の算出の答えが合わず。scipy使って自分でやった解答はこれ。

In [10]: 1 - norm(scale=sqrt(25. / 3), loc=6).cdf(norm(scale=sqrt(25. / 3)).ppf(0.95))
Out[10]: 0.66771319889291236

Rのpowerを使うとこう。

> power.t.test(n=6, delta=6, sd=5, type="two.sample", alternative="one.sided")

     Two-sample t test power calculation 

              n = 6
          delta = 6
             sd = 5
      sig.level = 0.05
          power = 0.6142154
    alternative = one.sided

NOTE: n is number in *each* group

テキストの答えは \(0.96\) でほど遠い。そして知恵袋にも質問があったが、こちらは\(\mu_1 - \mu_2 = 5\)となっていて間違えているように見える。書籍では \(6\) になっている。追加の質問もあって、結局よくわからないことになっている。

ということで、誤植であるという結論に至ったのだけどどうなんだろ。