布瓦松分布

機率統計

教學錄影

數學符號

數學基礎

排列組合

機率統計簡介

機率

機率公理

隨機變數

連續測度

單一分布

條件機率

聯合分布

貝氏定理

動差生成函數

特徵函數

機率法則匯總

離散分布

二項分布

多項分布

負二項分布

幾何分布

超幾何分布

布瓦松分布

連續分布

均勻分布

常態分布

Gamma 分布

指數分布

卡方分布

柯西分布

Weibull 分布

T 分布

F 分布

Beta 分布

多維分布

統計

抽樣

敘述統計

推論統計

中央極限定理

估計方法

單組樣本估計

兩組樣本估計

檢定方法

單組樣本檢定

兩組樣本檢定

平均値的推論

變異數的推論

無母數推論

迴歸分析

變異數分析

實驗設計

因子實驗

品質管制

時間序列

數據分類

統計定理匯總

統計情況分類

計算統計

蒙地卡羅法

最大似然法則

假說與學習

EM 算法

簡單貝氏分類

貝氏網路

隨機過程

馬可夫鏈

蒙地卡羅馬可夫

資源

範例

投影片

教學錄影

練習題

考題解答

訊息

相關網站

參考文獻

最新修改

簡體版

English

布瓦松分布 (Poisson distribution)

  • 意義:在單位時間內,事件出現平均 λ 次的機率分布。
  • 課本公式:
(1)
\begin{align} f(x) = \frac{e^{-k} k^x}{x!} \end{align}
  • R 的公式:p(x) = λ^x e^{-λ}/x!
(2)
\begin{align} pois(λ) = \frac{λ^x e^{-λ}}{x!} \end{align}
  • 變數意義:k = λ;

特性

布瓦松分布可以與泰勒展開式中的 Maclaurin 級數對映起來,所謂的Maclaurin級數就是泰勒展開式在 0 點的展開式。

If the Taylor series is centered at zero, then that series is also called a Maclaurin series, named after the Scottish mathematician Colin Maclaurin, who made extensive use of this special case of Taylor series in the 18th century.

(3)
\begin{equation} e^x = 1+x+x^2/2!+x^3/3!+ ... + x^k/k! + .... \end{equation}

布瓦松分配的公式來源

布瓦松分配可視為二項分配的極限形式,當 binom(n, p) 當中 n 趨近於無限大,而 p 非常小的時候,就會趨近布瓦松分配。

關鍵公式:

(4)
\begin{align} \lim_{n\to\infty}\left(1-{\lambda \over n}\right)^n=e^{-\lambda} \end{align}

證明過程:

(5)
\begin{eqnarray} \lim_{n\to\infty} P(X_n=k) &=& \lim_{n\to\infty}{n \choose k} p^k (1-p)^{n-k} \\ &=&\lim_{n\to\infty}{n! \over (n-k)!k!} \left({\lambda \over n}\right)^k \left(1-{\lambda\over n}\right)^{n-k}\\ &=&\lim_{n\to\infty} \underbrace{\left[\frac{n!}{n^k\left(n-k\right)!}\right]}_{A_n} \left(\frac{\lambda^k}{k!}\right) \underbrace{\left(1-\frac{\lambda}{n}\right)^n}_{\to\exp\left(-\lambda\right)} \underbrace{\left(1-\frac{\lambda}{n}\right)^{-k}}_{\to 1} \\ &=& \left[ \lim_{n\to\infty} A_n \right] \left(\frac{\lambda^k}{k!}\right)\exp\left(-\lambda\right) \\ &\to& \left(\frac{\lambda^k}{k!}\right)\exp\left(-\lambda\right) \end{eqnarray}

其中的 An 趨近於 1 ,證明如下:

(6)
\begin{eqnarray} A_n &=& \frac{n!}{n^k\left(n-k\right)!}\\ &=& \frac{n\cdot (n-1)\cdots \big(n-(k-1)\big)}{n^k}\\ &=& 1\cdot(1-\tfrac{1}{n})\cdots(1-\tfrac{k-1}{n})\\ &\to & 1\cdot 1\cdots 1 = 1 \end{eqnarray}

期望值與變異數

(7)
\begin{eqnarray} 1. && E(X) = k = λ \\ 2. && Var(X) = k = λ \end{eqnarray}

動差生成函數

(8)
\begin{equation} m_x(t) = e^{k (e^t-1) } = e^{λ (e^t-1) } \end{equation}

R 程式範例一

lambda=5.0; k=seq(0,20); 
plot(k, dpois(k, lambda), type='h', main='dpois(lambda=4.0)', xlab='k')
dpoisPlot.jpg

R 程式範例二

require(graphics)

-log(dpois(0:7, lambda=1) * gamma(1+ 0:7)) # == 1
Ni <- rpois(50, lambda = 4); table(factor(Ni, 0:max(Ni)))

1 - ppois(10*(15:25), lambda=100)  # becomes 0 (cancellation)
    ppois(10*(15:25), lambda=100, lower.tail=FALSE)  # no cancellation

par(mfrow = c(2, 1))
x <- seq(-0.01, 5, 0.01)
plot(x, ppois(x, 1), type="s", ylab="F(x)", main="Poisson(1) CDF")
plot(x, pbinom(x, 100, 0.01),type="s", ylab="F(x)",
     main="Binomial(100, 0.01) CDF")

執行結果:

> require(graphics)
> 
> -log(dpois(0:7, lambda=1) * gamma(1+ 0:7)) # == 1
[1] 1 1 1 1 1 1 1 1
> Ni <- rpois(50, lambda = 4); table(factor(Ni, 0:max(Ni)))

 0  1  2  3  4  5  6  7  8 
 1  3  6  8 11 11  4  3  3 
> 
> 1 - ppois(10*(15:25), lambda=100)  # becomes 0 (cancellation)
 [1] 1.233094e-06 1.261664e-08 7.085799e-11 2.252643e-13 4.440892e-16
 [6] 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00
[11] 0.000000e+00
>     ppois(10*(15:25), lambda=100, lower.tail=FALSE)  # no cancellation
 [1] 1.233094e-06 1.261664e-08 7.085800e-11 2.253110e-13 4.174239e-16
 [6] 4.626179e-19 3.142097e-22 1.337219e-25 3.639328e-29 6.453883e-33
[11] 7.587807e-37
> 
> par(mfrow = c(2, 1))
> x <- seq(-0.01, 5, 0.01)
> plot(x, ppois(x, 1), type="s", ylab="F(x)", main="Poisson(1) CDF")
> plot(x, pbinom(x, 100, 0.01),type="s", ylab="F(x)",
+      main="Binomial(100, 0.01) CDF")
>
poisson.jpg

參考文獻

  1. Wikipedia:卜瓦松分佈
  2. Wikipedia:Poisson_distribution

Facebook

Unless otherwise stated, the content of this page is licensed under Creative Commons Attribution-NonCommercial-ShareAlike 3.0 License