目的

負の超幾何分布の確率を求める

使用法

NegativeGeometric(x, N, n, r)
NegativeGeometric2(x, N, n, r)

引数

x	確率変数
N	母集団の大きさ
n	母集団中で特性を持つものの個数
r	標本中で特性を持つものの個数

ソース

# ガンマ関数の自然対数を使って階乗を計算
# n! = exp(lgamma(n+1))

インストールは，以下の 1 行をコピーし，R コンソールにペーストする
source("http://aoki2.si.gunma-u.ac.jp/R/src/NegativeGeometric.R", encoding="euc-jp")

# 負の超幾何分布

# ガンマ関数の自然対数を使って階乗を計算
# n! = exp(lgamma(n+1))
NegativeGeometric <- function(x, N, n, r)
{
        exp(lgamma(x)-lgamma(r)-lgamma(x-r+1)-lgamma(N+1)+lgamma(n+1)+lgamma(N-n+1)+lgamma(N-x+1)-lgamma(n-r+1)-lgamma(N-x-n+r+1))
}

# 二項係数 choose(n, k) の自然対数を利用
# nCk = exp(lchoose(n, k))

NegativeGeometric2 <- function(x, N, n, r)
{
        exp(lchoose(x-1, r-1)+lchoose(N-x, n-r)-lchoose(N, n))
}



使用例

> NegativeGeometric(4, 5000, 400, 3)
[1] 0.001404224

> NegativeGeometric2(4, 5000, 400, 3)
[1] 0.001404224

 解説ページ