LoginSignup
3

More than 3 years have passed since last update.

簡単なポアソン分布まとめ(Rでの出力付き)

Posted at
tags: 緑本

ポアソン分布

1. ポアソン分布とは

二項分布で, nが大きくpが小さい時, ポアソン分布に従う.

2. ポアソン分布の確率の定義

ポアソン分布の確率は, 以下のように定義される.
$$
p(y | \lambda) = \frac{\lambda^{y}\exp(-\lambda)}{y!}
$$
平均が$\lambda$である時, 確率変数が$y$という値になる確率.

3. ポアソン小数の法則

$np → \lambda$となるように, $n → \infty$, $p → 0$となる極限においては, 以下が成り立つ.
$$
{}_n \mathrm{C} _y\, p^y(1 - p)^{n - y} → \frac{\lambda^{y}\exp(-\lambda)}{y!}
$$

4. ポアソン分布の期待値と確率変数

ポアソン分布を, $Po(\lambda)$で表す.
確率変数$Y$がポアソン分布$Po(\lambda)$に従うならば, 期待値と分散は等しくなる.
$$
E(Y) = \lambda, V(Y) = \lambda\
(\lambda\fallingdotseq np, 1 - p \fallingdotseq 1 \Longrightarrow E(y) = np → \lambda, V(y) = np(1 - p) → \lambda)
$$

5. ポアソン分布を使う場合

主に, 以下のような場合ポアソン分布を使う.

  1. カウントデータである
  2. 平均と分散が概ね等しい
  3. 標本が大きく, 確率が小さい時

6. ポアソン分布の例

プロシア騎兵連隊において馬に蹴られて亡くなった兵士の数

$\lambda = 0.61, E(Y) = 0.61, V(Y) = 0.61, \sqrt{V(Y)} = 0.77$

死亡数 0 1 2 3 4 5
観測数 109 65 22 3 1 0
理論値 108.7 66.3 20.2 4.1 0.6 0.1

7. ポアソン分布の出力

7-1. ポアソン分布の確率の出力

ポアソン分布の確率の出力は, dpois(x, lambda = "")を用いて行う. 6. での例を用いると, 以下のように出力できる.

y <- c(0, 1, 2, 3, 4, 5)
prob <- dpois(y, lambda = 0.61)
prob %>% round(digits = 2)
> prob %>% round(digits = 2)
[1] 0.54 0.33 0.10 0.02 0.00 0.00

7-2. ポアソン分布のグラフの出力

plot()を用いて, グラフを出力させると, 以下のようになる.

y <- c(0, 1, 2, 3, 4, 5)
prob <- dpois(y, lambda = 0.61)
plot(y, prob, type="o")

8. 参考文献

・久保拓弥(2012)『データ解析のための統計モデリング入門』 岩波書店.
・東京大学教養学部統計学教室 (編集)(1991)『統計学入門 (基礎統計学Ⅰ)』 東京大学出版.

Register as a new user and use Qiita more conveniently

  1. You get articles that match your needs
  2. You can efficiently read back useful information
  3. You can use dark theme
What you can do with signing up
3