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

★ R ** はてな ★

　326　R ** はてな　　青木繁伸　　2002/01/22 (火) 18:14
　　329　Re: R ** はてな　　sb812109　　2002/01/23 (水) 01:28
　　　331　Re^2: R ** はてな　　青木繁伸　　2002/01/23 (水) 10:54

326.　R ** はてな　　青木繁伸　　2002/01/22 (火) 18:14

R の ?binom 関数で，落とし穴に入りそうなこと > p <- 0.5 # 母比率 > n <- 5 # 試行数 > f <-sapply(0:n, dbinom, n, p) # 密度関数 > names(f) <- 0:n > f 0 1 2 3 4 5 0.03125 0.15625 0.31250 0.31250 0.15625 0.03125 それで，Pr{x <= 1} は > pbinom(1, n, p) [1] 0.1875 で，Pr{x >=1} は > pbinom(1, n, p, lower.tail=FALSE) [1] 0.8125 と思いきや， > pbinom(0, n, p, lower.tail=FALSE) [1] 0.96875 でなくてはならない。
lower.tail logical; if TRUE (default), probabilities are P[X <= x], otherwise, P[X > x].
離散分布のときは，両方に等号がある方が良さそうに思う。
S でも同じですか。

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

329.　Re: R ** はてな　　sb812109　　2002/01/23 (水) 01:28

仕様なのでしょうか?　バグなのでしょうか?　間違えやすそうですね。

ちなみに，S-PLUS 2000 の pbinom では，lower.tail の指定は出来ません。

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

331.　Re^2: R ** はてな　　青木繁伸　　2002/01/23 (水) 10:54

> 仕様なのでしょうか?　バグなのでしょうか?　間違えやすそうですね。

バグではないです。
そのような仕様になっているのでしょう。

> ちなみに，S-PLUS 2000 の pbinom では，lower.tail の指定は出来ません。

そのようですね。
情報ありがとうございます。

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