目的

ポアソン分布に従うことが期待される標本から，ポアソン定数（λ）の信頼区間を求める。

使用法

poisson.conf(x, conf=0.95)

引数

x      標本（観察値のベクトル）
conf   信頼率（信頼度）。パーセント表現ではないので注意。

ソース

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

poisson.conf <- function(x, conf=0.95)
{
        N2 <- length(x)*2
        df <- 2*sum(x)+2
        alpha2 <- (1-conf)/2
        return(qchisq(c(alpha2, 1-alpha2), df)/N2)
}


使用例

poisson.conf(c(2,1,2,3,4,3))
set.seed(1)
poisson.conf(rpois(50, lambda=5))

出力結果例

> poisson.conf(c(2,1,2,3,4,3))
[1] 1.524230 4.123370
> set.seed(1)
> poisson.conf(rpois(50, lambda=5)) # λ=5 のポアソン乱数 50 個から，λを推定する
[1] 4.643444 5.914431