La mayoría de los experimentos no son de carácter deterministico, es decir, no conocemos de antemano el resultado del mismo, razón por la cual, es necesario medir la posibilidad de ocurrencia de un evento en presencia de incertidumbre. A esta medición se le llama probabilidad.
Además de medir la incertidumbre de un evento, a partir de la probabilidad es posible construir modelos que permitan describir las características de una variable en una población, y a partir de éstos, hacer inferencias sobre los parámetros del modelo con base en los datos de una muestra.
Por lo tanto, resulta de vital importancia el conocimiento de los conceptos básicos de probabilidad, pues los modelos de probabilidad, son la base de la construcción de los intervalos de confianza y de los estadísticos de prueba en el juzgamiento de hipótesis.
A continuación una breve descripción de las operaciones entre conjuntos y su notación, dado que son de vital importancia en la teoría de probabilidad.
Sean \(\Omega\) el conjunto universal, \(A\in\Omega\) y \(B\in\Omega\):
Dado que el cálculo de probabilidades implica contar en algunas situaciones prácticas, a continuación se presentan algunos principios del conteo. En particular vamos a contestar la siguiente pregunta: para \(k<n\), ¿de cuántas formas se pueden seleccionar \(k\) elementos de \(n\)?
Sean \(E_1,E_2,\ldots,E_m\) conjuntos (procesos). Si:
entonces, el número de formas en que puede ocurrir \(E_1\cap E_2\cap... E_m\) es \(\prod_{j=1}^{m}n_j\).
En consulta, un médico determina que su paciente debe tomar medicación para dos condiciones, para la primera tiene 3 opciones de medicamentos y para la segunda 5, ¿cuántas posibles prescripciones diferentes puede hacer el médico?
Así, el número de prescripciones es: \(n_1*n_2=3*5=15\).
Nota: al resulado de la multiplicación de los \(n\) primeros enteros se le llama \(n!\), \[1\times 2\times ...\times n=n!\] Por definición \(0!=1\).
La operación resultante se llama permutación sin repetición: \[P^n_k=n\times (n-1)\times ...\times (n-(k-1))=\frac{n!}{(n-k)!}\]
La operación resultante se llama permutación con repetición: \[n^k=n\times n\times....\times n\]
La operación resultante se llama combinación:
\[\binom{n}{k} = \frac{n!}{k!(n-k)!}\]
# permutacion
factorial(4)/factorial(0)
## [1] 24
# lista de posibilidades
library(gtools)
lista <- permutations(n = 4, r = 4, v = c("A", "T", "G", "C"))
head(lista)
## [,1] [,2] [,3] [,4]
## [1,] "A" "C" "G" "T"
## [2,] "A" "C" "T" "G"
## [3,] "A" "G" "C" "T"
## [4,] "A" "G" "T" "C"
## [5,] "A" "T" "C" "G"
## [6,] "A" "T" "G" "C"
dim(lista )
## [1] 24 4
# permutacion con repeticion
5^3
## [1] 125
# lista de posibilidades
lista <- permutations(n = 5, r = 3, v = 1:5, repeats.allowed = T)
head(lista)
## [,1] [,2] [,3]
## [1,] 1 1 1
## [2,] 1 1 2
## [3,] 1 1 3
## [4,] 1 1 4
## [5,] 1 1 5
## [6,] 1 2 1
dim(lista)
## [1] 125 3
# combinacion
choose(n = 20, k = 3)
## [1] 1140
# otra manera
factorial(20)/(factorial(3)*factorial(17))
## [1] 1140
# lista de posibilidades
lista <- combinations(n = 20, r = 3, v = 1:20)
head(lista)
## [,1] [,2] [,3]
## [1,] 1 2 3
## [2,] 1 2 4
## [3,] 1 2 5
## [4,] 1 2 6
## [5,] 1 2 7
## [6,] 1 2 8
dim(lista)
## [1] 1140 3
Ejercicios:
Para definir un espacio de probabilidad es preciso definir algunos conceptos previamente.
Un experimento aleatorio es cualquier experimento que satisface las siguientes condiciones:
El espacio muestral es el conjunto de todos los posibles resultados de un experimento aleatorio.
Cualquier subconjunto del espacio muestral se llama evento aleatorio.
Cuando un evento tiene un solo elemento se denomina evento elemental o simple.
El conjunto \(\Phi\) se llama evento imposible, que nunca sucede, y el conjunto \(\Omega\) se llama evento seguro, que siempre sucede.
Si \(A\) es un evento y el resultado observado del experimento aleatorio es un elemento de \(A\) significa que el evento \(A\) ha sucedido.
Considere un estudio longitudinal con 3 pacientes.
Una medida de probabilidad es una función que le asigna un número entre \(0\) y \(1\) a los eventos de un experimento aleatorio: \[ A\longrightarrow P(A) \] y que satisface:
Observe que la probabilidad de un evento \(A\), denotada con \(P(A)\), es una medida de la incertidumbre relacionada con la posibilidad de la ocurrencia del evento \(A\).
Se supone que un experimento aleatorio se repite \(n\) veces y que un evento \(A\) asociado con estos experimentos ocurre exactamente \(k\) veces. Entonces, la frecuencia relativa de \(A\), \(h_n(A)\), se define como la proporción entre la cantidad de veces que ocurre el evento \(A\) y el número total de repeticiones del experimento aleatorio: \[ h_n(A) = \frac{k}{n} \]
La frecuencia relativa de un evento \(A\) correspondiente a un “gran número” de repeticiones de un experimento aleatorio es igual a la probabilidad del evento \(A\), es decir \[ P(A) = \lim_{n\to\infty} h_n(A) \]
Considerar las siguientes situaciones:
Una máquina produce 100 tubos de ensayo cada 5 minutos. Esta máquina empieza su funcionamiento a las 8:00 a.m. y termina a las 8:00 p.m., hora en la que se toma una muestra de tamaño \(n\) del lote y se revisa el número de tubos de ensayo defectuosos. Si el número de tubos defectuosos es menor que una cantidad predeterminada por el departamento de control de calidad, entonces la producción del día se puede distribuir, de otra manera no es posible. Teniendo en cuenta que no hay cambios en el programa de producción de un día a otro: ¿Cuál es la probabilidad de que la producción de un día no se pueda distribuir? ¿Es correcto estimar la probabilidad mediante la frecuencia relativa?
Un complejo de oficinas se jacta de tener uno de los mejores controles de seguridad. Se quiere corroborar lo anterior estimando la probabilidad de que una persona que no trabaje en la oficina pueda pasar sin tener que anunciarse en recepción. ¿Cómo se calcula esta probabilidad? ¿Es correcto estimar la probabilidad mediante la frecuencia relativa?
En la situación a., si se toman los registros de producción durante cierto número de días, entonces se puede calcular la probabilidad de que la producción de un día no se pueda distribuir mediante: \[ \frac{\text{No. de días en los que la producción no se pudo distribuir}}{\text{No. de días en los que se tomaron los registros}}\,. \] Además, es correcto estimar esta probabilidad mediante el método frecuentista, porque se considera que no hay cambios en el programa de producción de un día a otro.
En la situación b., si se toman los registros de seguridad de un periodo determinado de tiempo, entonces se puede calcular la probabilidad de que una persona que no trabaje en la oficina pueda pasar sin tener que anunciarse en recepción, como sigue: \[ \frac{\text{No. de personas que ingresaron al complejo sin tener que anunciarse y que no trabajan allí}}{\text{No. de personas que entraron al complejo sin trabajar en él}}\,. \] En este caso, es correcto calcular esta probabilidad por medio de la frecuencia relativa siempre y cuando las condiciones de seguridad no varíen drásticamente de un día a otro.
Sea \(\Omega\) un espacio muestral no vacío, \(A\) y \(B\) eventos aleatorios incluidos en \(\Omega\) y \(P(.)\) una medida de probabilidad sobre \(\Omega\). Entonces se satisfacen las siguientes propiedades:
Sean \(A\), \(B\) y \(C\) eventos tales que \(P(A) = 0.50\), \(P(B) = 0.26\), \(P(C) = 0.55\), \(P(A \cap B) = 0.15\), \(P(A \cap C) = 0.25\), \(P(B \cap C) = 0.15\) y \(P(A \cap B \cap C) = 0.05\). Con base en esta información, calcular las siguientes probabilidades: \(P(A \cup B)\), \(P\left(A \cap C^C\right)\), \(P\left(A^C \cup C\right)\).
Aplicando las propiedades se obtiene que \[ P(A \cup B) = P(A) + P(B) - P(A \cap B) = 0.50 + 0.26 - 0.15 = 0.61. \] Además, \(P(A)=P(A\cap C)+P\left(A \cap C^C\right)\) de donde \[\begin{equation*} P\left(A \cap C^C\right) = P(A) - P(A \cap C) = 0.50 - 0.25 = 0.25. \end{equation*}\] Por otra parte, utilizando las leyes de Morgan se concluye que \[\begin{equation*} P\left(A^C \cup C\right) = 1 - P\left(\left(A^C \cup C\right)^C\right) = 1 - P\left(A \cap C^C\right) = 1 - 0.25 = 0.75. \end{equation*}\]
El teorema aditivo está relacionado con la probabilidad de la unión de eventos, así:
\[P\left(\bigcup_{i=1}^{m} A_i\right)=\sum_{i=1}^{m}P(A_i)-\sum_{i>j}P(A_i\cap A_j)+\sum_{i>j>k}P(A_i\cap A_j\cap A_k)-\ldots-(-1)^{m-1}P\left(\bigcap_{i=1}^{m} A_i \right)\]
En particular, para \(m=3\):
\[P(A_1\cup A_2\cup A_3)=P(A_1)+P(A_2)+P(A_3)-P(A_1\cap A_2)-P(A_1\cap A_3)-P(A_2\cap A_3)+P(A_1\cap A_2\cap A_3)\]
La probabilidad condicional es una medida de probabilidad definida por:
\[P(A_i\mid A_j)=\frac{P(A_i\cap A_j)}{P(A_j)}\] siempre que \(P(A_j)\neq 0\).
Dado que la probabilidad condicional también constituye una medida de probabilidad, cumple las propiedades de una medida de probabilidad. Por ejemplo:
El teorema aditivo está relacionado con la probabilidad de la intersección de eventos.
El teorema de la probabilidad total se utiliza cuando se quiere encontrar la probabilidad de un evento que se encuentra “repartido” en las partes de una patición.
Sea \(E_1, E_2,\ldots, E_m\) una partición de \(\Omega\), es decir, \(E_1, E_2,\ldots, E_m\) es una colección de subconjuntos de \(\Omega\) tales que: \[\bigcup_{i=1}^m E_i=\Omega\] \[E_i\cap E_j=\phi, \qquad i\neq j\] Entonces, la probabilidad de ocurrencia de un evento \(A\) se reparte en la partición de la siguiente forma:
\[\begin{align*} P(A)&=P(A\cap E_1)+P(A\cap E_2)+\ldots+P(A\cap E_m)\\ &=P(A\mid E_1)P(E_1)+P(A\mid E_2)P(E_2)+\ldots+P(A\mid E_m)P(E_m)\\ &=\sum_{i=1}^{m}P(A\mid E_i)P(E_i) \end{align*}\]
El terorema de Bayes es ampliamente utilizado en epidemiología, especialmente en la evaluación de las características de pruebas diagnósticas. En el teorema parte de la probabilidad a priori de la ocurrencia de un evento, \(P(E_k)\), para calcular su probabilidad a posteriori, \(P(E_k\mid A)\).
\[\begin{align*} P(E_k\mid A)&=\frac{P(A\mid E_k)P(E_k)}{\sum_{i=1}^{m}P(A\mid E_i)P(E_i)} \end{align*}\]
Una empresa elabora sus productos en cuatro fábricas, a saber, fábrica 1, fábrica 2, fábrica 3 y fábrica 4. El porcentaje de producción total que se elabora en cada fábrica es del 40%, 30%, 20% y 10%, respectivamente, y además el porcentaje de envasado incorrecto en cada fábrica es del 1%, 2%, 7% y 4%. Tomando un producto de la empresa al azar, ¿cuál es la probabilidad de que se encuentre envasado de forma incorrecta?
Con la información del Ejempo 30, calcular la probabilidad de que un producto envasado de forma incorrecta provenga de la primera fábrica.
Según el teorema de Bayes y teniendo en cuenta las probabilidades dadas, se tiene que: \[\begin{align*} P(F_1|A) &=\frac{P(F_1)P(A |F_1)}{\sum_{i=1}^{4} P(F_i)P(A|F_i)}\\ &= \frac{P(F_1)P(A|F_1)}{P(F_1)P(A|F_1) + P(F_2)P(A|F_2) + P(F_3)P(A|F_3)+ P(F_4)P(A|F_4)} \\ &= \frac{(0.4)(0.01)}{ (0.4)(0.01) + (0.3)(0.02) + (0.2)(0.07) + (0.1)(0.04)} \\ &= 0.14285. \end{align*}\] Por lo tanto, la probabilidad de que un producto envasado de forma incorrecta provenga de la primera fábrica es 0.14285.
Se ha observado que los hombres y las mujeres reaccionan de una manera diferente en ciertas circunstancias; el \(70\) por ciento de las mujeres reaccionan positivamente mientras que el porcentaje de los hombres es sólo del \(40\) por ciento. Se sometió a prueba a un grupo de 20 personas, 15 mujeres y 5 hombres, y se les pidió llenar un cuestionario para descubrir sus reacciones. Una respuesta escogida al azar de las 20 resultó negativa. ¿Cuál es la probabilidad de que haya sido contestada por un hombre?.
Un médico aplica una prueba diagnóstica a 10 alumnos de un colegio para detectar una enfermedad cuya incidencia sobre una población de niños es del 10 por ciento. La sensibilidad del test es del 80 por ciento y la especificidad del 75 por ciento.