Distribució de Poisson

Distribució de Poisson

L'eix horitzontal és l'índex k . La funció només està definida en valors sencers de k . Les línies que connecten els punts són només guies per a l'ull i no indiquen continuïtat.
Funció de distribució de probabilitat L'eix horitzontal és l'índex k .
Tipus	distribution de Tweedie, Distribució binomial, natural exponential family with quadratic variance function (en) , Distribució de Conway–Maxwell–Poisson, Panjer distribution (en) , Distribució de Poisson composta, distribució univariant i distribució de probabilitat discreta
Epònim	Siméon Denis Poisson
Paràmetres	${\displaystyle \lambda \in (0,\infty )}$
Suport	k ∈ { 0, 1, 2, 3, ... }
FD	${\displaystyle {\frac {\Gamma (\lfloor k+1\rfloor ,\lambda )}{\lfloor k\rfloor !}}{\text{per }}k\geq 0}$ (on ${\displaystyle \Gamma (x,y)}$ és la Funció gamma incompleta)
Esperança matemàtica	${\displaystyle \lambda \,}$
Mediana	${\displaystyle {\text{normalment prop de}}\lfloor \lambda +1/3-002/\lambda \rfloor }$
Moda	${\displaystyle \lfloor \lambda \rfloor {\text{(amb }}\lambda -1{\text{ si }}\lambda {\text{ és un enter)}}}$
Variància	${\displaystyle \lambda \,}$
Coeficient de simetria	${\displaystyle \lambda ^{-1/2}\,}$
Curtosi	${\displaystyle 3+\lambda ^{-1}\,}$
Entropia	${\displaystyle \lambda [1\!-\!\ln(\lambda )]\!+\!i^{-\lambda }\sum _{k=0}^{\infty }{\frac {\lambda ^{k}\ln(k!)}{k!}}}$
FC	${\displaystyle \exp(\lambda (e^{it}-1))\,}$
EOM	Poisson_distribution
Mathworld	PoissonDistribution

En teoria de probabilitat i estadística, la distribució de Poisson o llei dels petits nombres o dels fenòmens rars és una distribució de probabilitat discreta que és un bon model per molts fenòmens naturals o socials. Una propietat especialment important és que una distribució binomial de paràmetres ${\displaystyle n}$ i ${\displaystyle p}$, ${\displaystyle B(n,p)}$, amb ${\displaystyle n}$ gran i ${\displaystyle p}$ petit es pot aproximar (en distribució) per una distribució de Poisson de paràmetre ${\displaystyle n\times p}$; moltes de les aplicacions de la distribució de Poisson es justifiquen en base aquesta propietat.

Aquesta distribució va ser introduïda per Siméon-Denis Poisson l'any 1837 en el seu treball Recherches sur la probabilité des jugements en matière criminelle et en matière civile^[1] ( Investigació sobre la probabilitat dels judicis en matèries criminals i civils ).

Les referències bàsiques d'aquesta pàgina són^[2][3]

Una variable aleatòria de Poisson ${\displaystyle X}$ de paràmetre ${\displaystyle \lambda >0}$ és una variable discreta que pot prendre qualsevol valor natural, 0, 1, 2, ... amb probabilitats$${\displaystyle P(X=k)=e^{-\lambda }{\frac {\lambda ^{k}}{k!}},\ k=0,1,2,\dots }$$S'escriu ${\displaystyle X\sim {\text{Poi}}(\lambda )}$ o ${\displaystyle X\sim {\text{Poiss}}(\lambda )}$ o alguna notació similar.

1. Esperança i variància. L'esperança i la variància d'una variable aleatòria amb distribució de Poisson de paràmetre ${\displaystyle \lambda }$ són iguals a ${\displaystyle \lambda }$:$${\displaystyle E[X]=\lambda \quad {\text{i}}\quad {\text{Var}}(X)=\lambda .}$$

Així, si la distribució de Poisson s'utilitza per modelar el comptatge d'un fenomen, aleshores ${\displaystyle \lambda }$ es pot interpretar com la freqüència esperada del fenomen.

Prova: $${\displaystyle E[X]=\sum _{k=0}^{\infty }k\,e^{-\lambda }{\frac {\lambda ^{k}}{k!}}=\lambda \,e^{-\lambda }\sum _{k=1}^{\infty }{\frac {\lambda ^{k-1}}{(k-1)!}}=\lambda \,e^{-\lambda }e^{\lambda }=\lambda ,}$$on hem utilitzat el desenvolupament es sèrie de la funció exponencial.

Per trobar el moment de segon ordre i la variància calcularem ${\displaystyle E[X(X-1)]}$:

$${\displaystyle E[X(X-1)]=\sum _{k=0}^{\infty }k\,(k-1)\,e^{-\lambda }{\frac {\lambda ^{k}}{k!}}=\lambda ^{2}\,e^{-\lambda }\sum _{k=2}^{\infty }{\frac {\lambda ^{k-2}}{(k-2)!}}=\lambda ^{2}.}$$D'altra banda, $${\displaystyle E[X(X-1)]=E[X^{2}]-E[X]=E[X^{2}]-\lambda ,}$$d'on$${\displaystyle E[X^{2}]=\lambda ^{2}+\lambda .}$$Llavors, $${\displaystyle {\text{Var}}(X)=E[X^{2}]-{\big (}E[X]{\big )}^{2}=\lambda .}$$

2. Moments d'ordre superior. Els moments d'ordre superior són polinomis de Touchard en ${\displaystyle \lambda }$, i els coeficients tenen una interpretació combinatòria. De fet, quan el valor esperat de la distribució de Poisson és 1, llavors segons la fórmula de Dobinski, el ${\displaystyle n}$-èsim moment iguala el nombre de particions de mida ${\displaystyle n}$.

3. Moda. La moda d'una variable aleatòria de distribució de Poisson amb un ${\displaystyle \lambda }$ no sencer és igual a ${\displaystyle \lfloor \lambda \rfloor }$, el més gran dels enters menors que ${\displaystyle \lambda }$ (els símbols ${\displaystyle \scriptstyle \lfloor \ \ \rfloor }$ representen la funció part sencera). Quan ${\displaystyle \lambda }$ és un enter positiu, les modes són ${\displaystyle \lambda }$ i ${\displaystyle \lambda -1}$.

4. Funció generatriu de moments. La Funció generatriu de la distribució de Poisson de paràmetre ${\displaystyle \lambda }$ és

$${\displaystyle \mathrm {E} \left[e^{tX}\right]=\sum _{k=0}^{\infty }e^{tk}e^{-\lambda }{\lambda ^{k} \over k!}=e^{-\lambda }\sum _{k=0}^{\infty }{\frac {{\big (}\lambda e^{t}{\big )}^{k}}{k!}}=e^{-\lambda }\,e^{\lambda e^{t}}=e^{\lambda (e^{t}-1)}.}$$

Anàlogament, la funció característica és $${\displaystyle E[e^{itX}]=e^{\lambda (e^{it}-1)}.}$$

5. Les variables aleatòries de Poisson són infinitament divisibles.

6. Divergència de Kullback-Leibler. La divergència Kullback-Leibler d'una variable aleatòria de Poisson de paràmetre ${\displaystyle \lambda _{0}}$ a una altra de paràmetre ${\displaystyle \lambda }$ és

$${\displaystyle D_{\mathrm {KL} }(\lambda ||\lambda _{0})=\lambda \left(1-{\frac {\lambda _{0}}{\lambda }}+{\frac {\lambda _{0}}{\lambda }}\log {\frac {\lambda _{0}}{\lambda }}\right).}$$

7. Sumes de variables aleatòries de Poisson. La suma de variables aleatòries de Poisson independents és una altra variable aleatòria de Poisson en què el paràmetre és la suma dels paràmetres de les originals. Dit d'una altra manera, si

$${\displaystyle X_{i}\sim {\text{Poi}}(\lambda _{i})\,,i=1,\dots ,N}$$

són N variables aleatòries de Poisson independents, llavors

$${\displaystyle \sum _{i=1}^{N}X_{i}\sim {\text{Poi}}\left(\sum _{i=1}^{N}\lambda _{i}\right)\,}$$.

La distribució de Poisson és un cas límit de la distribució binomial: una distribució binomial ${\displaystyle B(n,p)}$ amb ${\displaystyle n}$ gran, ${\displaystyle p}$ petita, i ${\displaystyle np}$ petita respecte a ${\displaystyle n}$, es pot aproximar raonablement bé per una distribució de Poisson de paràmetre ${\displaystyle \lambda =np}$. Històricament, Poisson va utilitzar aquest argument per introduir la seva distribució,^[4] i és el principal raonament per justificar la utilització d'aquesta distribució, tal com veurem en els exemples més avall.

Formalment, sigui ${\displaystyle (p_{n})_{n\geq 1}}$ una successió de nombres ${\displaystyle 0<p_{n}<1}$ tals que ${\displaystyle \lim _{n\to \infty }np_{n}=\lambda >0}$. Considerem una llei binomial de paràmetres ${\displaystyle n}$ i ${\displaystyle p_{n}}$, ${\displaystyle X_{n}\sim B(n,p_{n})}$, i sigui ${\displaystyle X\sim {\text{Poi}}(\lambda )}$ Aleshores$${\displaystyle \lim _{n\to \infty }P(X_{n}=k)=P(X=k),\ k=0,1,\dots }$$

En efecte, tenim que $${\displaystyle P(X_{n}=k)={n \choose k}p_{n}^{k}(1-p_{n})^{n-k},\quad k=0,\ldots ,n.}$$ Llavors, per a ${\displaystyle k=0,}$ $${\displaystyle \lim _{n\to \infty }P(X_{n}=0)=\lim _{n\to \infty }(1-p_{n})^{n}=\lim _{n\to \infty }{\Big (}(1-p_{n})^{1/p_{n}}{\Big )}^{np_{n}}=e^{-\lambda }=P(X=0),}$$on hem utilitzat un límit tipus nombre ${\displaystyle e}$ .

Per a qualsevol natural ${\displaystyle k\geq 1}$, tenim, per a ${\displaystyle n}$ prou gran,$${\displaystyle {\begin{array}{rl}P\{X_{n}=k\}&={\dfrac {n!}{k!\,(n-k)!}}p_{n}^{k}(1-p_{n})^{n-k}\\\\&={\dfrac {n(n-1)\cdots (n-k+1)}{k!\,n^{k}}}\,(np_{n})^{k}{\Big [}(1-p_{n})^{1 \over {p_{n}}}{\Big ]}^{np_{n}}(1-p_{n})^{-k},\end{array}}}$$ i passant al límit s'obté $${\displaystyle P\{X_{n}=k\}={\frac {\lambda ^{k}}{k!}}\,e^{-\lambda }=P(X=k).}$$Noteu que d'aquí es dedueix que la successió ${\displaystyle X_{1},\,X_{2},\dots }$ convergeix en distribució a ${\displaystyle X}$.

En particular, si partim de ${\displaystyle \lambda >0}$ i prenem ${\displaystyle p_{n}=\lambda /n}$ aleshores tenim que una distribució binomial ${\displaystyle B(n,p)}$ amb ${\displaystyle n}$ gran, ${\displaystyle n}$ petita i ${\displaystyle n\,p=\lambda }$ relativament petita respecte ${\displaystyle n}$, es pot aproximar per una distribució de Poisson ${\displaystyle {\text{Poi}}(\lambda )}$. Aquesta propietat fa que la distribució de Poisson també s'anomeni llei dels petits nombres; ^[5] o dels fenòmens rars, ja que el fenomen subjacent té probabilitat petita; Johnson and Kotz ^[3] atribueixen el nom <<llei del petits nombres>> a Bortkiewicz (o Borkiewitsch) (1868-1931) i el seu llibre Das Gesetz der kleinen Zahlen (La llei dels petits nombres). Un autor de pes, Feller,^[6] considera desafortunat aquest nom i dona nombrosos exemples per mostrar-ho. Cal tenir present que si bé la probabilitat ${\displaystyle p}$ és petita, el producte ${\displaystyle n\,p}$ pot ser gran, sempre que sigui relativament petit respecte de ${\displaystyle n}$ .

A conseqüència de teorema central del límit, per a valors grans de ${\displaystyle \lambda }$, una variable aleatòria de Poisson ${\displaystyle X}$ es pot aproximar per una variable normal de mitjana i variància ${\displaystyle \lambda }$ ja que $${\displaystyle \lim _{\lambda \to \infty }{\frac {X-\lambda }{\sqrt {\lambda }}}=N(0,1),\ {\text{en distribució}}.}$$Vegeu .^[7][8]

Suposem que per a cada valor ${\displaystyle t>0}$, que representa el temps, el nombre d'aparicions de cert fenomen aleatori segueix una distribució de Poisson de paràmetre ${\displaystyle \lambda t}$. Llavors, els temps que discorren entre dos aparicions successives segueix la distribució exponencial.

Si ${\displaystyle X}$ és una distribució de Poisson de paràmetre ${\displaystyle \lambda }$, aleshores ^[9] per a ${\displaystyle k=0,\,1,2,\dots }$,

$${\displaystyle P(X\leq k)=P(\chi _{2(k+1)}^{2}>2\lambda ),}$$

on ${\displaystyle \chi _{2(k+1)}^{2}}$ és una variable aleatòria amb distribució ${\displaystyle \chi ^{2}}$ amb ${\displaystyle 2(k+1)}$ graus de llibertat.

Aquesta fórmula és útil perquè permet calcular el valor de la funció de distribució d'una variable aleatòria de Poisson de manera ràpida.

En el llibre que hem citat abans de Bortkiewicz, s'estudia el nombre de soldats de l'exèrcit prussià morts per una guitza del seu cavall. A la taula^[10] hi ha el nombre de morts corresponent a 10 cossos de característiques similars ^[11] durant 20 anys

En total es van produir 122 morts. La següent taula és un resum de la Taula 1 en termes de freqüències absolutes i relatives:

La mitjana de morts per cos i any és $${\displaystyle {\overline {X}}={\frac {0\cdot 109+1\cdot 65+2\cdot 22+3\cdot 3+4\cdot 1}{200}}={\frac {122}{200}}=0,61.}$$

Suposem que cos de l'exèrcit tenia 1000 soldats (com veurem, aquest nombre, que desconeixem, no té cap paper, podríem posar ${\displaystyle n}$). Un model raonable per al nombre de soldats morts (per cos i any) és una distribució binomial ${\displaystyle B(1000,p)}$, on ${\displaystyle p}$ és la probabilitat que un soldat resulti mort. Aquesta probabilitat és molt petita i podem estimar-la de la següent manera: Tenim 10 cossos durant 10 anys, que representen ${\displaystyle 10\times 20\times 1000=200.000}$ soldats. Aleshores $${\displaystyle p={\frac {\text{nombre de soldats morts}}{\text{nombre total de soldats}}}={\frac {122}{200000}}=0,00061.}$$D'acord amb els comentaris anteriors, podem aproximar la distribució binomial per una distribució de Poisson de paràmetre $${\displaystyle \lambda =n\times p=1000\times 0,00061=0,61,}$$que és la mitjana que hem calculat abans. La Taula 3 dona les probabilitats d'una variable aleatòria de Poisson ${\displaystyle X}$ de paràmetre ${\displaystyle \lambda =0,61}$, ${\displaystyle X\sim {\text{Poi}}(0,61)}$:

Taula 3. Taula de probabilitats d'una distribució de Poisson Poi(0,61)

${\displaystyle X}$	0	1	2	3	${\displaystyle \geq 4}$
Probabilitat	0.543	0,331	0,101	0,020	0,003

En comparar la fila de les freqüències relatives de la Taula 2 i les probabilitats de la Taula 3 veiem que la concordança entre ambdues és realment notable.

Rutherford i Geiger ^[12] reporten que varen observar el nombre de partícules ${\displaystyle \alpha }$ emeses per una massa de material radioactiu durant 2.608 intervals de 7,5 segons cadascun. A la Taula 4 la ${\displaystyle k}$ denota el nombre de partícules emeses en un interval de 7,5 segons, la freqüència absoluta és el nombre de vegades que es va observar el corresponent valor ${\displaystyle k}$, després hi ha freqüència relativa i finalment la probabilitat que una variable de Poisson de paràmetre ${\displaystyle \lambda =3,87}$ prengui aquell valor:

Taula 4. Experiment de Rutherford et al.

${\displaystyle k}$	0	1	2	3	4	5	6	7	8	9	10	11	12	13	14	Total
Freq. absoluta	57	203	383	525	532	408	273	139	45	27	10	4	0	1	1	2608
Freq. relativa	0,022	0,078	0,147	0,201	0,204	0,156	0,105	0,053	0,017	0,010	0.004	0.002	0.000	0,000	0,000	1
Probabilitat Poisson	0,021	0,081	0,156	0,201	0,195	0,151	0,097	0,054	0,026	0,011	0.004	0.002	0.000	0,000	0,000	1

El paràmetre ${\displaystyle \lambda }$ s'ha calculat com la mitjana del nombre de partícules emeses per unitat de temps: $${\displaystyle \lambda ={\overline {X}}={\frac {0\times 57+1\times 203+\cdots }{2608}}={\frac {10097}{2608}}=3,87.}$$

Com veiem la concordança entre les dades observades i les probabilitats corresponents a la distribució de Poisson és molt bona. Breiman ^[13] argumenta perquè la distribució de Poisson apareix en aquesta situació: Suposem que durant un període llarg, la mitjana de partícules emeses per unitat de temps és ${\displaystyle \lambda }$ (en aquest cas, la unitat de temps seria 7,5 seg). Considerem ara el nombre de partícules emeses durant un interval de temps de longitud ${\displaystyle t}$ unitats; el resultat pot ser 0, 1, 2, etc. Per trobar una assignació raonable de probabilitats dividim l'interval ${\displaystyle [0,t]}$ en ${\displaystyle n}$ parts d'amplada ${\displaystyle t/n}$; prenem ${\displaystyle n}$ prou gran perquè els subintervals siguin molt petits; en cada subinterval podem observar 0, 1, 2, ... partícules. Suposem:

1. Les probabilitats d'observar 0, 1, ... partícules són les mateixes per a tots els subintervals.

2. La probabilitat d'observar dos o més partícules en un subinterval és menyspreable comparada amb les probabilitats d'observar-ne zero o una partícula.

3. La probabilitat d'observar o no una partícula en un subinterval és independent d'observar o no una partícula en un altre subinterval.

Durant la Segona Guerra Mundial els alemanys van bombardejar Londres i altres llocs amb bombes volants. Molta gent creia que les bombes tendien a agrupar-se en clústers; amb l'objectiu d'analitzar si aquesta suposició era veritat es va dividir el sud de Londres en 576 quadrats de 0'5 km²; durant el període considerat van caure 537 bombes, vegeu la Figura 1.

La taula 5 ^[14] mostra el nombre de quadrats on van caure 0 bombes, 1 bomba, etc. (freqüència absoluta), així com la freqüència relativa.

Taula 5. Bombardeig sobre Londres

${\displaystyle k}$	0	1	2	3	4	${\displaystyle \geq 5}$	Total
Freq. absoluta	229	211	93	35	7	1	576
Freq. relativa	0,397	0,366	0,161	0,061	0,012	0,002	1
Prob. Poisson	0,394	0,367	0,171	0,053	0,012	0,002	1
Freq. esperada Poisson	226,7	211,4	98,5	36,6	7,1	1,3	575,8

La mitjana de bombes per quadrat és $${\displaystyle {\overline {X}}={\frac {537}{576}}\approx 0,932.}$$ A la taula 5 hi ha una fila amb les probabilitats ${\displaystyle P(X=k)}$ d'una variable de Poisson ${\displaystyle X}$ de paràmetre ${\displaystyle \lambda =0,932}$ per ${\displaystyle k=0,1,\dots ,4}$ i el valor ${\displaystyle P(X\geq 5)}$. També hi ha una fila amb els productes$${\displaystyle 576\times P(X=k),\quad k=0,1,\dots ,4\quad {\text{i}}\quad 576\times P(X\geq 5),}$$els quals són les freqüències absolutes esperades per una variable de Poisson.

L'ajustament de les dades observades amb les donades per la distribució de Poisson és molt bona, la qual cosa suggereix que no va haver-hi clústers i que el lloc de caiguda de les bombes era completament aleatori. Chung^[15] argumenta el bon ajust amb la distribució de Poisson dient que, si el fenomen és purament aleatori, es pot estudiar com un problema de col·locació de 537 objectes en 576 caixes; la probabilitat d'observar exactament ${\displaystyle k}$ bombes (objectes) en un quadrat (una capsa) serà ${\displaystyle P(Y=k)}$, on ${\displaystyle Y}$ és una variable binomial de paràmetres ${\displaystyle n=537}$ i ${\displaystyle p=1/576}$ : ${\displaystyle Y\sim B(537,1/576)}$. Atès que ${\displaystyle n}$ és gran i ${\displaystyle p}$ petita, podem aproximar ${\displaystyle Y}$ per una variable de Poisson de paràmetre $${\displaystyle \lambda =n\,p=537\times {\frac {1}{576}}\approx 0,932,}$$que és el que hem fet anteriorment.

Feller ^[16] també analitza aquestes dades i comenta que <<la taula anterior indica que hi ha aleatorietat i homogeneïtat perfectes, [...] però que per a [l'ull] inexpert, l'aleatorietat apareix com regularitat [...]>>. Per profunditzar en aquest comentari de Feller és molt interessant el capítol 2 del llibre de J. A. Paulos ^[17]

Si el 2% dels llibres enquadernats en cert taller té enquadernació defectuosa, la probabilitat que 5 de 400 llibres enquadernats en aquest taller tinguin enquadernacions defectuoses es pot calcular usant la distribució de Poisson. En efecte, el nombre de llibres defectuosos es pot modelar per una distribució binomial de paràmetres ${\displaystyle n=400}$ i ${\displaystyle p=0,02}$, la qual el pot aproximar per una distribució de Poisson de paràmetre $${\displaystyle \lambda =400\times 0.02=8,}$$ que, d'altra banda, és la mitjana de llibres defectuosos. Per tant, la probabilitat desitjada és

$${\displaystyle P(X=5)={\frac {8^{5}e^{-8}}{5!}}=0,092.}$$

Aquest problema també podria resoldre's recorrent a una distribució binomial de paràmetres ${\displaystyle k=5}$, ${\displaystyle n=400}$ i ${\displaystyle p=0,02}$.

La distribució de Poisson s'aplica a diversos fenòmens discrets de la natura (és a dir, aquells fenòmens que ocorren 0, 1, 2, 3, ... vegades durant un període definit de temps o en una àrea determinada) quan la probabilitat d'ocurrència del fenomen és constant en el temps o l'espai. Exemples d'aquests esdeveniments que poden ser modelats per la distribució de Poisson inclouen:

• El nombre de cotxes que passen a través d'un cert punt en una ruta (prou distants dels semàfors) durant un període definit de temps.
• El nombre d'errors d'ortografia que un comet en escriure una única pàgina.
• El nombre de trucades telefòniques en una central telefònica per minut.
• El nombre de servidors web accedits per minut.
• El nombre d'animals morts trobats per unitat de longitud de ruta.
• El nombre de mutacions de determinada cadena d'ADN després de certa quantitat de radiació.
• El nombre de nuclis atòmics inestables que van decaure en un determinat període en una porció de substància radioactiva. La radioactivitat de la substància es debilitarà amb el temps, per tant el temps total de l'interval utilitzat en el model ha de ser significativament menor que la vida mitjana de la substància.
• El nombre d'estels en un determinat volum d'espai.
• El nombre d'estels fugaços per una unitat de temps.
• La distribució de receptors visuals a la retina de l'ull humà.
• La inventiva d'un inventor a través de la seva carrera.^[18]

A Wikimedia Commons hi ha contingut multimèdia relatiu a: Distribució de Poisson

• Procés de Poisson
• Regressió de Poisson