# Distribución binomial

En teoría de la probabilidad y estadística, la distribución binomial o distribución binómica es una distribución de probabilidad discreta que cuenta el número de éxitos en una secuencia de $n$ ensayos de Bernoulli independientes entre sí con una probabilidad fija $p$ de ocurrencia de éxito entre los ensayos. Un experimento de Bernoulli se caracteriza por ser dicotómico, esto es, solo dos resultados son posibles, a uno de estos se le denomina “éxito” y tiene una probabilidad de ocurrencia $p$ y al otro se le denomina “fracaso” y tiene una probabilidad $q=1-p$.
La distribución binomial se utiliza con frecuencia para modelizar el número de aciertos en una muestra de tamaño n extraída con reemplazo de una población de tamaño N. Si el muestreo se realiza sin reemplazo, las extracciones no son independientes, por lo que la distribución resultante es una distribución hipergeométrica, no una distribución binomial. Sin embargo, para N mucho mayores que n, la distribución binomial sigue siendo una buena aproximación, y se utiliza ampliamente.
Más matemáticamente, la distribución binomial es una distribución discreta de probabilidad descrita por dos parámetros: n el número de experimentos realizados, y p la probabilidad de éxito. Para cada experimento llamado ensayo Bernoulli, utilizamos una variable aleatoria que toma el valor 1 cuando se consigue un éxito y el valor 0 en caso contrario. La variable aleatoria, suma de todas estas variables aleatorias, cuenta el número de éxitos y sigue una distribución binomial. Es posible entonces obtener la probabilidad de k éxitos en una repetición de n experimentos:
    $\mathbb {P} (X=k)={n \choose k}\,p^{k}(1-p)^{n-k}.$
En esta fórmula interviene el coeficiente binomial ${n \choose k}$ del que se deriva el nombre de la ley.
La importancia de esta ley es en primer lugar histórica ya que fue objeto de estudio del Teorema de De Moivre-Laplace, resultado del xviii fundador de la teoremas de convergencia. Una distribución binomial también se puede utilizar para modelar situaciones simples de acertar o fallar, un juego de lanzar una moneda, por ejemplo. El cálculo de su función de probabilidad se vuelve rápidamente tedioso cuando n} es grande, es posible entonces utilizar aproximaciones por otras leyes de probabilidad|leyes de probabilidad como la distribución de Poisson o la distribución normal y utilizar tablas de valores.
La ley binomial se utiliza en diversos campos de estudio, especialmente a través de pruebas estadísticas que permiten interpretar datos y tomar decisiones en situaciones que dependen del azar. Debido a su sencilla definición, es una de las leyes de probabilidad que se estudian en los cursos introductorios de teoría de la probabilidad.

## Explicación intuitiva
### Primeros ejemplos
#### Cara o cruz

Consideremos n lanzamientos sucesivos de una moneda. Entonces el número de veces que la moneda aparece en la cara de la cruz sigue la distribución binomial donde el número de experimentos realizados es n y la probabilidad de éxito es $p={\frac {1}{2}}$

#### Lanzar un dado

Consideremos n lanzamientos sucesivos de un dado de 6 caras. Entonces el número de veces que sale un 1 sigue la distribución binomial donde el número de experimentos realizados es n y la probabilidad de éxito es $p={\frac {1}{6}}$.

### Definición intuitiva

Una ley de Bernoulli describe el comportamiento de un experimento aleatorio que tiene dos posibles resultados tradicionalmente llamados éxito y fracaso. Un experimento de este tipo se denomina ensayo de Bernoulli. Por ejemplo, en un lanzamiento de cara o cruz, cara puede considerarse un éxito y cruz un fracaso. En este modelo, la probabilidad de éxito es un valor fijo, es decir, permanece constante en cada renovación del experimento aleatorio. Esta probabilidad de éxito se denomina "p".
Consideremos la situación en la que un experimento aleatorio de este tipo (dos resultados posibles y una probabilidad fija) se repite un cierto número de veces de forma independiente; denotemos n este número de veces. Esta repetición independiente de ensayos de Bernoulli se denomina Esquema de Bernoulli o simplemente Ensayos de Bernoulli. Una distribución binomial describe el número de veces que aparece el éxito en los n experimentos realizados. Dado que el número de éxitos obtenidos es un valor aleatorio, una distribución binomial se describe dando las probabilidades de que el éxito se produzca precisamente k veces en los n ensayos.

### Árbol de probabilidades

Una forma visual de representar una secuencia de experimentos es utilizar un árbol de probabilidad. Cada suceso está representado por dos ramas: una para el éxito y otra para el fracaso. En cada extremo, se añaden dos ramas (éxito y fracaso) para el siguiente ensayo. Esto se repite hasta que se alcanza el número total de sucesos. En cada extremo final, se puede contar el número de éxitos. Se multiplica el número de veces que k tiene éxito por la probabilidad de que k tenga éxito para obtener la probabilidad correspondiente de la distribución binomial.

## Definición
### Notación

Si una variable aleatoria discreta $X$ tiene una distribución binomial con parámetros $n\in \mathbb {N}$ y $p$ con $0<p<1$ entonces escribiremos $X\sim \operatorname {Bin} (n,p)$.

### Función de Probabilidad

Si $X\sim \operatorname {Bin} (n,p)$ entonces su función de probabilidad está dada por
    $\operatorname {P} [X=x]={n \choose x}p^{x}(1-p)^{n-x}$
para $x=0,1,2,\dots ,n$, siendo
    $\!{n \choose x}={\frac {n!}{x!(n-x)!}}\,\!$
el coeficiente binomial y se lee “las combinaciones de $n\,\!$ en $x\,\!$“.
En ocasiones, para calcular las probabilidades binomiales se utiliza la siguiente fórmula recursiva para calcular $\operatorname {P} [X=x+1]$ en términos de $\operatorname {P} [X=x]$
    ${\begin{aligned}\operatorname {P} [X=x+1]&={\binom {n}{x+1}}p^{x+1}(1-p)^{n-(x+1)}\\&={\frac {n!}{(n-x-1)!(x+1)!}}p^{x+1}(1-p)^{n-x-1}\\&={\frac {n!(n-x)}{(n-x)!x!(x+1)}}p^{x}(1-p)^{n-x}{\frac {p}{1-p}}\\&=\left({\frac {n-x}{x+1}}\right)\left({\frac {p}{1-p}}\right){\frac {n!}{x!(n-x)!}}p^{x}(1-p)^{n-x}\\&=\left({\frac {n-x}{x+1}}\right)\left({\frac {p}{1-p}}\right){\binom {n}{x}}p^{x}(1-p)^{n-x}\\&=\left({\frac {n-x}{x+1}}\right)\left({\frac {p}{1-p}}\right)\operatorname {P} [X=x]\end{aligned}}$

### Función de Distribución Acumulada

La función de distribución acumulada de una variable aleatoria $X\sim \operatorname {Bin} (n,p)$ está dada por
    $F_{X}(x)=\operatorname {P} [X\leq x]=\sum _{k=0}^{\lfloor x\rfloor }{n \choose k}p^{k}(1-p)^{n-k}$ si $0\leq x\leq n$
También puede ser expresada en términos de la función beta incompleta como
    ${\begin{aligned}F_{X}(x)&=\operatorname {P} [X\leq x]\\&=I_{1-p}(n-x,x+1)\\&=(n-x){n \choose x}\int _{0}^{1-p}t^{n-x-1}(1-t)^{x}dt\end{aligned}}$
que es equivalente a la función de distribución acumulada de la distribución F.
La distribución binomial es la base del test binomial de significación estadística.

## Experimento binomial

Existen muchas situaciones en las que se presenta una experiencia binomial. Cada uno de los experimentos es independiente de los restantes (la probabilidad del resultado de un experimento no depende del resultado del resto). El resultado de cada experimento ha de admitir sólo dos categorías (a las que se denomina éxito y fracaso). El valor de ambas posibilidades ha de ser constante en todos los experimentos, y se denotan como $p$ y $q$ respectivamente o como $p$ y $1-p$ de forma alternativa.
Se designa por $X$ a la variable que mide el número de éxitos que se han producido en los $n$ experimentos.
Cuando se dan estas circunstancias, se dice que la variable $X$ sigue una distribución de probabilidad binomial.

### Ejemplo

Supongamos que se lanza 51 veces un dado de 6 caras y queremos calcular la probabilidad de que el número 3 salga 20 veces.
En este problema un ensayo consiste en lanzar el dado una vez. Consideramos un éxito si obtenemos un 3 pero si no sale 3 lo consideramos como un fracaso. Defínase $X$ como el número de veces que se obtiene un 3 en 51 lanzamientos.
En este caso tenemos $X\sim \operatorname {Bin} (51,1/6)$ por lo que la probabilidad buscada es $\operatorname {P} [X=20]$
    $\operatorname {P} [X=20]={51 \choose 20}(1/6)^{20}(1-1/6)^{51-20}=0.0000744\,\!$

## Propiedades

Si $X$ es una variable aleatoria discreta tal que $X\sim \operatorname {Bin} (n,p)$ entonces
* $\operatorname {E} [X]=np$
* $\operatorname {Var} [X]=np(1-p)$
La primera de ellas es fácil de demostrar, por definición de Esperanza
    $\operatorname {E} [X]=\sum _{x=0}^{n}x{\binom {n}{x}}p^{x}(1-p)^{n-x}\,$
el primer término de la suma, es decir, para $x=0$ el término vale cero por lo que podemos iniciar la suma en $x=1$
    $\operatorname {E} [X]=\sum _{x=1}^{n}x{\binom {n}{x}}p^{x}(1-p)^{n-x}\,$
Dado que
    ${\begin{aligned}{\binom {n}{x}}&={\frac {n!}{x!(n-x)!}}\\&={\frac {n}{x}}{\frac {(n-1)!}{(x-1)!((n-1)-(x-1))!}}\\&={\frac {n}{x}}{\binom {n-1}{x-1}}\end{aligned}}$
para $x\geq 1$.
Reemplazando lo anterior en la expresión de $\operatorname {E} [X]$ obtenemos
    ${\begin{aligned}\operatorname {E} [X]&=\sum _{x=1}^{n}x{\frac {n}{x}}{\binom {n-1}{x-1}}p^{x}(1-p)^{n-x}\\&=n\sum _{x=1}^{n}{\binom {n-1}{x-1}}p^{x}(1-p)^{n-x}\end{aligned}}$
Haciendo el cambio de índice $k=x-1$ obtenemos
    ${\begin{aligned}\operatorname {E} [X]&=n\sum _{k=0}^{n-1}{\binom {n-1}{k}}p^{k+1}(1-p)^{n-(k+1)}\\&=np\sum _{k=0}^{n-1}{\binom {n-1}{k}}p^{k}(1-p)^{(n-1)-k}\end{aligned}}$
Finalmente por la fórmula de Newton (Teorema del binomio)
    $(a+b)^{n}=\sum _{k=0}^{n}{\binom {n}{k}}a^{k}b^{n-k}$
Obtenemos
    $\operatorname {E} [X]=np(p+1-p)^{n-1}=np$.

## Distribuciones Relacionadas
### Suma de Binomiales

Si $X\sim \operatorname {Bin} (n,p)$ y $Y\sim \operatorname {Bin} (m,p)$ son variables aleatorias independientes con la misma probabilidad $p$ entonces la variable aleatoria $Z=X+Y$ también es una variable aleatoria con distribución binomial con parámetros $n+m$ y $p$, es decir $Z=X+Y\sim \operatorname {Bin} (n+m,p)$
    ${\begin{aligned}\operatorname {P} [Z=z]&=\sum _{k=0}^{z}{\binom {n}{k}}p^{k}(1-p)^{n-k}{\binom {m}{z-k}}p^{z-k}(1-p)^{m-z+k}\\&={\binom {n+m}{z}}p^{z}(1-p)^{n+m-z}\end{aligned}}$

### Distribución Bernoulli

Si $X_{1},X_{2},\dots ,X_{n}$ son $n$ variables aleatorias independientes e idénticamente distribuidas tales que $X_{i}\sim \operatorname {Bernoulli} (p)$ entonces
    $\sum _{i=1}^{n}X_{i}\sim \operatorname {Bin} (n,p)$
Lo anterior es equivalente a decir que la distribución Bernoulli es un caso particular de la distribución Binomial cuando $n=1$, es decir, si $X\sim \operatorname {Bin} (1,p)$ entonces $X\sim \operatorname {Bernoulli} (p)$.

### Distribuciones limitantes
#### Teorema límite de Poisson

Si $n\to \infty$ y $p$ es tal que el producto entre ambos parámetros tiende a $\lambda \,\!$, entonces la distribución de la variable aleatoria binomial tiende a una distribución de Poisson de parámetro $\lambda$.

#### Teorema de De Moivre-Laplace

Si $X$ es una variable aleatoria con media $np$ y varianza $np(1-p)$ entonces
    $Z={\frac {X-np}{\sqrt {np(1-p)}}}\sim N(0,1)$
conforme $n\to \infty$, esta aproximación es buena si $np\geq 5$ y $n(1-p)\geq 5$.

## Propiedades reproductivas

Si $X_{1},X_{2},\dots ,X_{k}$ son variables aleatorias independientes tales que $X_{i}\sim \operatorname {Bin} (n_{i},p)$ con $i=1,2,\dots ,k$ entonces
    $\sum _{i=1}^{k}X_{i}\sim \operatorname {Bin} \left(\sum _{i=1}^{k}n_{i},p\right)$

## Inferencia estadística
### Estimación de parámetros

Cuando se conoce n, el parámetro p puede estimarse utilizando la proporción de aciertos:
    ${\widehat {p}}={\frac {x}{n}}.$
Este estimador se encuentra utilizando estimador de máxima verosimilitud y también el método de los momentos. Este estimador es insesgado y uniforme con mínima varianza, demostrado mediante el Teorema de Lehmann–Scheffé, ya que se basa en un estadístico mínimo suficiente y completo (es decir: x). También es consistente tanto en probabilidad como en MSE.
También existe un estimador de Bayes de forma cerrada para p cuando se utiliza la distribución Beta como conjugada de la probabilidad ''a priori''. Cuando se utiliza un $\operatorname {Beta} (\alpha ,\beta )$ a priori, el estimador medio posterior es:
    ${\widehat {p}}_{b}={\frac {x+\alpha }{n+\alpha +\beta }}.$
El estimador de Bayes es asintóticamente eficiente y a medida que el tamaño de la muestra se aproxima a infinito (n → ∞), se aproxima a la solución de máxima verosimilitud. El estimador de Bayes es sesgado, cuánto depende de los priores, admisible y consistente en probabilidad.
Para el caso especial de utilizar la distribución uniforme estándar como a priori no informativo, $\operatorname {Beta} (\alpha =1,\beta =1)=U(0,1)$, el estimador de la media a posteriori se convierte en:
    ${\widehat {p}}_{b}={\frac {x+1}{n+2}}.$
Este método se denomina regla de sucesión, que fue introducido en el siglo XVIII por Pierre-Simon Laplace.
Cuando se estima p con sucesos muy raros y un n pequeño (por ejemplo: si x=0), entonces utilizar el estimador estándar conduce a ${\widehat {p}}=0,$ lo que a veces es poco realista y poco deseable. En estos casos existen varios estimadores alternativos. Una forma es utilizar el estimador de Bayes, lo que lleva a:
    ${\widehat {p}}_{b}={\frac {1}{n+2}}.$
Otro método consiste en utilizar el límite superior del intervalo de confianza obtenido mediante la regla de tres:
    ${\widehat {p}}_{\text{rule of 3}}={\frac {3}{n}}.$

### Intervalos de confianza

Incluso para valores bastante grandes de n, la distribución real de la media es significativamente no normal. Debido a este problema se han propuesto varios métodos para estimar intervalos de confianza.
En las ecuaciones para intervalos de confianza que se presentan a continuación, las variables tienen el siguiente significado:
* n $_1$ es el número de aciertos de n, el número total de ensayos.
* ${\widehat {p\,}}={\frac {n_{1}}{n}}$ es la proporción de aciertos
* $z=1-{\tfrac {1}{2}}\alpha$ es el cuantil de una distribución normal estándar (es decir, probit) correspondiente a la tasa de error objetivo $\alpha$. Por ejemplo, para un nivel de confianza del 95% el error $\alpha$ = 0,05, por lo que $1-{\tfrac {1}{2}}\alpha$ = 0,975 y $z$ = 1,96.