1 Introducción

Se presentan algunas distribuciones probabilísticas discretas de uso común que son de especial importancia porque representan los modelos teóricos de los fenómenos aleatorios más frecuentes.

2 Distribución Binomial

La distribución de Bernoulli es una distribución de probabilidad que consiste en la observación de un experimento aleatorio con dos posibles resultados. A uno de los resultados del experimento se le denomina “éxito”, con probabilidad de ocurrencia \(\pi\), y al otro se le llama “fracaso”, con probabilidad de ocurrencia \(1 - \pi\), donde \(\pi\) es un número real tal que \(0<\pi<1\).

La distribución binomial está conformada por \(n\) repeticiones independientes de un experimento de Bernoulli.

Se dice que una v.a.d. \(X\) que asume valores \(0,1,\ldots,n\) tiene una distribución binomial con parámetros \(n\) y \(\pi\), lo que se escribe \(X\sim Bin(n,\pi)\), si la f.m.p. de \(X\) está dada por: \[ f_X(x;n,\pi) = \left\{ \begin{array}{ll} \binom{n}{x} \pi^x (1-\pi)^{n-x}, & \hbox{si $x=0,1,\ldots,n$;} \\ 0, & \hbox{en otro caso.} \end{array} \right. \] donde \(n\) es un entero positivo y \(\pi\) es un número real tal que \(0 < \pi <1\).

Cuando \(n=1\), la distribución binomial coincide con la distribución Bernoulli de parámetro \(\pi\), lo que se escribe \(X \sim Ber(\pi)\).

2.1 Propiedades

Si \(X\) es una v.a. tal que \(X \sim Bin(n,\pi)\), entonces:

  • \(\textsf{E}[X] = n\pi\).
  • \(\textsf{V}[X] = n\pi(1-\pi)\).

2.2 Ejemplo

Una empresa farmacéutica desarrolló un nuevo medicamento y lo suministró a 10 enfermos elegidos aleatoriamente. La experiencia ha demostrado que 30% de las personas que padecen la enfermedad se recupera al tomar dicho medicamento. ¿Cuál es la probabilidad de que por lo menos nueve de las 10 personas que toman el medicamento se recuperen?

En este caso se tiene que la v.a. de estudio es el “número de personas en la muestra de 10 pacientes que se recupera de la enfermedad” y el éxito consiste en recuperarse de la enfermedad y esto ocurre con una probabilidad de \(0.3\). En consecuencia, la f.m.p. de \(X\) es \[ f_X(x;10,0.3) = \left\{ \begin{array}{ll} \binom{10}{x} (0.3)^x (0.7)^{10-x}, & \hbox{si $x=0,1,\ldots,10$;} \\ 0, & \hbox{en otro caso.} \end{array} \right. \]

Además, se pide calcular la probabilidad de que por lo menos nueve de las 10 personas que toman el medicamento se recuperen, esto es, \(\textsf{P}(X \geq 9)\). Así, \[\begin{align*} \textsf{P}(X \geq 9) &= \textsf{P}(X=9;X=10)\\ &= \textsf{P}(X=9) + \textsf{P}(X=10) \\ &= \binom{10}{9} (0.3)^9 (0.7)^{10-9} + \binom{10}{10} (0.3)^{10} (0.7)^{10-10} \\ %&= 0.000138 + 0.000006 \\ &= 0.000144. \end{align*}\]

De otra parte, se observa que:

  • \(\textsf{E}[X] = 10(0.3)=3\). Este valor indica que se espera la recuperación de 3 enfermos de una muestra aleatoria de 10 pacientes.
  • \(\textsf{V}[X] = 10(0.3)(0.7) = 2.1\)

En la siguiente figura se presenta el gráfico de la f.m.p. y de la f.d.a. de una variable con distribución binomial con parámetros \(n=10\) y \(\pi=0.3\).

# parametros
p <- 0.3
n <- 10
x <- 0:n
# P(X >= 9)
sum(dbinom(x = c(9, 10), size = n, prob = p))
## [1] 0.0001436859
# valor esperado
n*p
## [1] 3
# varianza 
n*p*(1-p)
## [1] 2.1
# f.m.p.
fx <- dbinom(x = x, size = n, prob = p)
# f.d.a.
Fx <- pbinom(q = x, size = n, prob = p)
# graficos
par(mfrow = c(1,2))
# f.m.p
plot(x = x, y = fx, xlab = "x", ylab = "f(x)", pch = 15, col = "blue", main="Función de probabilidad")
segments(x0 = x, y0 = 0, x1 = x, y1 = fx, lwd = 2, col = "blue")
# f.d.a.
plot(x = c(0, x), y = c(0, Fx), type = "s", xlab = "x", ylab = "F(x)", col = "blue", lwd = 2, main="Función de distribución")
points(x, Fx, col = "blue", pch = 15)

3 Distribución Poisson

La distribución Poisson se utiliza para caracterizar probabilísticamente el número de veces que ocurre un evento en relación con una unidad de medida bien definida (como una unidad de tiempo o espacio, por ejemplo), de forma que:

  1. La probabilidad de que el evento ocurra en una unidad de medida dada es igual para todas las unidades.
  2. El número de eventos que ocurren en una unidad de medida es independiente del número de eventos que ocurren en otras unidades.

Se dice que una v.a. \(X\) tiene distribución de Poisson de parámetro \(\lambda\), se escribe \(X\sim Pois(\lambda).\), si la f.m.p. de \(X\) está dada por: \[ f_X(x;\lambda) = \left\{ \begin{array}{ll} \cfrac{e^{-\lambda} \lambda^x}{x!}, & \hbox{si $x=0,1,2,\ldots$;} \\ 0, & \hbox{en otro caso.} \end{array} \right. \] donde \(\lambda\) es un número real positivo.

3.1 Propiedades

Si \(X\) es una v.a. tal que \(X \sim P(\lambda)\), entonces:

  • \(\textsf{E}[X] = \lambda\).
  • \(\textsf{V}[X] = \lambda\).

3.2 Ejemplo

Los pacientes que entran a un centro de salud lo hacen a una tasa esperada de 0.50 clientes por minuto. Hallar la probabilidad de que el número de clientes que entran en un intervalo específico de 10 minutos sea a lo más 3.

Las hipótesis del proceso de Poisson parecen ser razonables en este contexto. Se da por sentado que los pacientes no llegan en grupos (o es posible contar al grupo entero como un solo paciente) y que la entrada de un paciente no aumenta ni disminuye la probabilidad de que llegue otro.

Para obtener \(\lambda\), se observa que a una tasa media de 0.50 por minuto durante un periodo de 10 minutos, se sigue que \(\lambda = (0.50)(10) = 5\) entradas. Sea \(X\) la v.a. dada por el “número de pacientes que entran en un intervalo de 10 minutos”. Así, se tiene que \(X \sim Pois(5)\), por lo que la f.m.p de \(X\) es \[ f_X(x;5) = \left\{ \begin{array}{ll} \frac{e^{-5} 5^x}{x!}, & \hbox{si $x=0,1,2,\ldots$;} \\ 0, & \hbox{en otro caso.} \end{array} \right. \]

Se pide calcular

\[\begin{align*} \textsf{P}(X \leq 3) &= \textsf{P}(X=0;X=1;X=2;X=3) \\ &= \textsf{P}(X=0) + \textsf{P}(X=1) + \textsf{P}(X=2) + \textsf{P}(X=3) \\ &= \frac{e^{-5} 5^0}{0!} + \frac{e^{-5} 5^1}{1!} + \frac{e^{-5} 5^2}{2!} + \frac{e^{-5} 5^3}{3!} \\ &= 0.2650. \end{align*}\] Se observa además que \[ \textsf{P}(X > 3) = 1 - \textsf{P}(X \leq 3) = 1 - 0.2650 = 0.7350. \]

En siguiente figura presenta el gráfico de la f.m.p. y la f.d.a. de una variable con distribución de poisson de parámetro \(\lambda=5\).

# parametros
lambda <- 5
x <- 0:20
# P(X <= 3)
sum(dpois(x = 0:3, lambda = lambda))
## [1] 0.2650259
# otra manera
ppois(q = 3, lambda = lambda)
## [1] 0.2650259
# f.m.p.
fx <- dpois(x = x, lambda = lambda)
# f.d.a.
Fx <- ppois(q = x, lambda = lambda)
# graficos
par(mfrow = c(1,2))
# f.m.p
plot(x = x, y = fx, xlab = "x", ylab = "f(x)", pch = 15, col = "blue", main="Función de probabilidad")
segments(x0 = x, y0 = 0, x1 = x, y1 = fx, lwd = 2, col = "blue")
# f.d.a.
plot(x = c(0, x), y = c(0, Fx), type = "s", xlab = "x", ylab = "F(x)", col = "blue", lwd = 2, main="Función de distribución")
points(x, Fx, col = "blue", pch = 15)