Distribución binomial

Distribución binomial
Función de masa de probabilidad Función de probabilidad
Función de distribución acumulada Función de distribución de probabilidad
Parámetros	$n\geq 0$ número de ensayos (entero) $0\leq p\leq 1$ probabilidad de éxito (real)
Dominio	$x\in \{0,\dots ,n\}\!$
Función de probabilidad (fp)	${n \choose x}p^{x}(1-p)^{n-x}\!$
Función de distribución (cdf)	$I_{1-p}(n-\lfloor x\rfloor ,\lfloor x\rfloor +1)\!$
Media	$np\!$
Mediana	Uno de $\{\lfloor np\rfloor ,\lceil np\rceil \}$​
Moda	$\lfloor (n+1)\,p\rfloor \!$
Varianza	$np(1-p)\!$
Coeficiente de simetría	${\frac {1-2p}{\sqrt {np(1-p)}}}\!$
Curtosis	${\frac {1-6p(1-p)}{np(1-p)}}+3\!$
Entropía	${\frac {1}{2}}\ln \left(2\pi nep(1-p)\right)+O\left({\frac {1}{n}}\right)$
Función generadora de momentos (mgf)	$(1-p+pe^{t})^{n}\!$
Función característica	$(1-p+pe^{it})^{n}\!$

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.

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}^{x}{n \choose k}p^{k}(1-p)^{n-k}$

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/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 ₁ 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.

Véase también

• Distribución Bernoulli
• Distribución Binomial Negativa
• Distribución de Poisson
• Distribución Normal
• Distribución Gamma
• Distribución Beta

Bibliografía

• Hirsch, Werner Z. (1957). «Binomial Distribution—Success or Failure, How Likely Are They?». Introduction to Modern Statistics. New York: MacMillan. pp. 140-153. 
• Neter, John; Wasserman, William; Whitmore, G. A. (1988). Applied Statistics (Third edición). Boston: Allyn & Bacon. pp. 185-192. ISBN 0-205-10328-6. 
• Patrick Bogaert (2005). Probabilités pour scientifiques et ingénieurs. De Boeck Supérieur. p. 402. ISBN 2-8041-4794-0. 

Enlaces externos

• Wikimedia Commons alberga una categoría multimedia sobre Distribución binomial.
• Tablas de la distribución binomial, hasta n=20, en formato PDF.
• Calculadora Distribución binomial
• Cálculo de la probabilidad de una distribución binomial con R (lenguaje de programación)
• Generación estadística de la distribución binomial con números aleatorios usando Python (lenguaje de programación)
• Interactive graphic: Univariate Distribution Relationships
• Binomial distribution formula calculator
• Difference of two binomial variables: X-Y or |X-Y|
• Querying the binomial probability distribution in WolframAlpha

• Datos: Q185547
• Multimedia: Binomial distributions / Q185547