title: Loi Poisson binomiale
text: # Loi Poisson binomiale
  
  En théorie des probabilités et en statistique, la loi Poisson binomiale est une loi de probabilité discrète de la somme d'épreuves de Bernoulli indépendantes.
  
  En d'autres termes, c'est la loi de probabilité du nombre de succès (nombre de pile) d'une suite de $n$ lancers de pile ou face dont les probabilités de succès (d'obtenir pile) sont $p_{1},p_{2},\dots ,p_{n}$. La loi binomiale ordinaire est un cas spécial de la loi Poisson binomiale lorsque toutes les probabilités sont les mêmes : $p_{1}=p_{2}=\dots =p_{n}$ .
  
  ## Espérance et variance
  
  Puisque la loi Poisson binomiale est une somme de n variables aléatoires indépendantes de loi de Bernoulli, son espérance et sa variance sont simplement les sommes des espérances et variances des lois de Bernoulli :
  
      $\mu =\sum \limits _{i=1}^{n}p_{i}$
      $\sigma ^{2}=\sum \limits _{i=1}^{n}(1-p_{i})p_{i}.$
  
  ## Fonction de masse
  
  La probabilité d'obtenir $k$ succès sur un total de n essais peut être écrit comme la somme :
  
      $\mathbb {P} (K=k)=\sum \limits _{A\in F_{k}}{\prod \limits _{i\in A}{{p}_{i}}\prod \limits _{j\in {{A}^{c}}}{(1-{{p}_{j}})}}$
  
  où $F_{k}$ est l'ensemble de tous les sous-ensembles de $\{1,2,\dots ,n\}$ contenant $k$ éléments. Par exemple si n=3, alors $F_{2}=\left\{\{1,2\},\{1,3\},\{2,3\}\right\}$. $A^{c}$ est le complémentaire de $A$.
  
  L'ensemble $F_{k}$ contient (n/k)=n!/(n−k)!k! {\textstyle {n \choose k}={\frac {n!}{(n-k)!k!}}} éléments, ainsi les calculs deviennent très grands en pratique, par exemple pour n=30, $F_{15}$ contient un nombre de l'ordre de $10^{20}$ éléments. Il existe cependant des méthodes efficaces pour calculer $\mathbb {P} (K=k)$.
  
  On peut utiliser une formule $itérative^,$ :
  
      $\mathbb {P} (K=k)=\left\{{\begin{aligned}&\prod _{i=1}^{n}(1-p_{i})&\mathrm {si} \ k=0\\&{\frac {1}{k}}\sum _{i=1}^{k}(-1)^{i-1}\mathbb {P} (K=k-i)T(i)&\mathrm {si} \ k>0\end{aligned}}\right.$
  
  où $T(i)=\sum _{j=1}^{n}\left({\frac {p_{j}}{1-p_{j}}}\right)^{i}$ .
  
  Une autre possibilité est d'utiliser la transformée de Fourier discrète :
  
      $\mathbb {P} (K=k)={\frac {1}{n+1}}\sum \limits _{l=0}^{n}{{\omega ^{lk}}\prod \limits _{m=1}^{n}{\left(1+({\omega ^{l}}-1){p_{m}}\right)}}$
  
  où $\omega =\exp \left(-{\frac {2{\rm {i}}\pi }{n+1}}\right)$ avec i l'unité imaginaire.
  
  D'autres méthodes sont décrites dans les ouvrages de Chen.
  
  