Una fabricante de lapiceros tiene un programa de control de calidad que incluye la inspección de lapiceros recibidos para revisar que no tengan defectos. Supongamos que, en cierto día, recibe lapiceros en lotes de cinco y se seleccionan dos lapiceros de un lote para inspeccionarlos. Podemos representar los posibles resultados del proceso de selección por pares. Por ejemplo, el par (3,4) representa la selección de los lapiceros 3 y 4 para inspeccionarlos.
El espacio muestral es
\[\Omega = \binom{5}{2} = 10\] \[\Omega = \{(1,2),(1,3),(1,4),(1,5),(2,3),(2,4),(2,5),(3,4),(3,5),(4,5)\}\]
combinaciones <- combn(5,2)
resultado <- apply(combinaciones, 2, function(x) paste0("(", x[1], ",", x[2], ")"))
print(resultado)
## [1] "(1,2)" "(1,3)" "(1,4)" "(1,5)" "(2,3)" "(2,4)" "(2,5)" "(3,4)" "(3,5)"
## [10] "(4,5)"
\[ X = \text{número de lapiceros defectuosos} \] \[f(0) = P(X=0) = P\big((1,2),(1,5),(2,5)\big) = \frac{3}{10} = 0.3\] \[f(1) = P(X=1) = P\big((1,3),(1,4),(2,3),(2,4),(3,5),(4,5)\big) = \frac{6}{10} = 0.6\] \[f(2) = P(X=2) = P\big((3,4)\big) = \frac{1}{10} = 0.1\]
## X Probabilidad
## 1 0 0.3
## 2 1 0.6
## 3 2 0.1
\[ P(X=0)=0.3, \quad P(X=1)=0.6, \quad P(X=2)=0.1. \]
Entonces, la función de distribución acumulada \[ F(x) = P(X \leq x) \] es
\[ F(x) = \begin{cases} 0, & x < 0, \\[6pt] 0.3, & 0 \leq x < 1, \\[6pt] 0.9, & 1 \leq x < 2, \\[6pt] 1, & x \geq 2. \end{cases} \]
Usando la formula de esperanza:
\[ \mu = E(X) = \sum_{k=0}^{2} x_k \, f(x_k) = x_0 f(x_0) + x_1 f(x_1) + x_2 f(x_2) \]
\[ = (0)(0.3) + (1)(0.6) + (2)(0.1) = 0 + 0.6 + 0.2 = 0.8 \]
\[ \mathrm{V}(X) = \sum_{k=0}^{2} (k^2 \cdot f(x_k)) - \left( \sum_{k=0}^{2} (k \cdot f(x_k)) \right)^2 \]
\[ = (0^2 \cdot f(x_0) + 1^2 \cdot f(x_1) + 2^2 \cdot f(x_2)) - \left( 0 \cdot f(x_0) + 1 \cdot f(x_1) + 2 \cdot f(x_2) \right)^2 \] \[ E[X] = 0(0.3) + 1(0.6) + 2(0.1) = 0.8 \]
\[ E[X^2] = 0^2(0.3) + 1^2(0.6) + 2^2(0.1) = 1.0 \]
\[ \mathrm{V}(X) = 1.0 - (0.8)^2 = 0.36 \]
x <- c(0, 1, 2)
p <- c(0.3, 0.6, 0.1)
esperanza <- sum(x * p)
esperanza_cuadrado <- sum((x^2) * p)
varianza <- esperanza_cuadrado - (esperanza^2)
esperanza
## [1] 0.8
esperanza_cuadrado
## [1] 1
varianza
## [1] 0.36