4.4: La distribución binomial

Como se puede imaginar, las distribuciones de probabilidad varían enormemente, y hay una enorme gama de distribuciones por ahí. Sin embargo, no todos son igualmente importantes. De hecho, la gran mayoría del contenido de este libro se basa en una de cinco distribuciones: la distribución binomial, la distribución normal, la\(t\) distribución, la distribución\(\chi^2\) (“chi-cuadrada”) y la\(F\) distribución. Ante esto, lo que voy a hacer en las próximas secciones es dar una breve introducción a las cinco, prestando especial atención al binomio y a lo normal. Comenzaré con la distribución binomial, ya que es la más simple de las cinco.

 

Presentamos el Binomio

La teoría de la probabilidad se originó en el intento de describir cómo funcionan los juegos de azar, por lo que parece adecuado que nuestra discusión sobre la distribución binomial implique una discusión sobre rodar dados y voltear monedas. Imaginemos un simple “experimento”: en mi manita caliente sostengo 20 dados idénticos de seis caras. En una cara de cada dado hay una imagen de una calavera; las otras cinco caras están todas en blanco. Si procedo a tirar los 20 dados, ¿cuál es la probabilidad de que obtenga exactamente 4 cráneos? Asumiendo que los dados son justos, sabemos que la posibilidad de que alguien muera subiendo cráneos es de 1 en 6; para decirlo de otra manera, la probabilidad de calavera para un solo dado es aproximadamente\(.167\). Esta es información suficiente para responder a nuestra pregunta, así que echemos un vistazo a cómo se hace.

Como es habitual, vamos a querer introducir algunos nombres y alguna notación. Dejaremos\(N\) denotar el número de tiradas de dados en nuestro experimento; que a menudo se conoce como el parámetro de tamaño de nuestra distribución binomial. También usaremos\(\theta\) para referirnos a la probabilidad de que un solo dado surja cráneos, una cantidad que suele llamarse la probabilidad de éxito del binomio. Por último, usaremos\(X\) para referirnos a los resultados de nuestro experimento, es decir, el número de cráneos que obtengo cuando arroje los dados. Dado que el valor real de\(X\) se debe al azar, nos referimos a él como una variable aleatoria. En cualquier caso, ahora que tenemos toda esta terminología y notación, podemos utilizarla para exponer el problema un poco más precisamente. La cantidad que queremos calcular es la probabilidad que\(X = 4\) dado que sabemos que\(\theta = 0.167\) y\(N=20\). La “forma” general de la cosa que me interesa calcular podría escribirse como

\[P(X \ | \ \theta, N) \nonumber \]

y nos interesa el caso especial donde\(X=4\),\(\theta = .167\) y\(N=20\). Sólo hay una pieza más de notación a la que quiero referirme antes de pasar a discutir la solución al problema. Si quiero decir que\(X\) se genera aleatoriamente a partir de una distribución binomial con parámetros\(\theta\) y\(N\), la notación que usaría es la siguiente:

\[X \sim \mbox{Binomial}(\theta, N) \nonumber \]

Sí, sí. Sé lo que estás pensando: notación, notación, notación. En serio, ¿a quién le importa? Muy pocos lectores de este libro están aquí por la notación, así que probablemente debería seguir adelante y hablar sobre cómo usar la distribución binomial. He incluido la fórmula para la distribución binomial en la Tabla 4.3.1, ya que algunos lectores pueden querer jugar con ella ellos mismos, pero como a la mayoría de la gente probablemente no le importa tanto y porque no necesitamos la fórmula en este libro, no voy a hablar de ello con ningún detalle. En cambio, solo quiero mostrarte cómo es la distribución binomial. Para ello, Figura\(\PageIndex{1}\) traza las probabilidades binomiales para todos los valores posibles de\(X\) para nuestro experimento de lanzamiento de dados, desde\(X=0\) (sin cráneos) hasta\(X=20\) (todos los cráneos). Tenga en cuenta que esto es básicamente un gráfico de barras, y no es diferente a la gráfica de “probabilidad de pantalones” que dibujé en la Figura 4.3.1. En el eje horizontal tenemos todos los eventos posibles, y en el eje vertical podemos leer la probabilidad de cada uno de esos eventos. Entonces, la probabilidad de rodar 4 cráneos de 20 veces es de aproximadamente 0.20 (la respuesta real es 0.2022036, como veremos en un momento). En otras palabras, esperarías que eso suceda alrededor del 20% de las veces que repitiste este experimento.

Figura\(\PageIndex{1}\): La distribución binomial con parámetro de tamaño de N = 20 y una probabilidad de éxito subyacente de 1/6. Cada barra vertical representa la probabilidad de un resultado específico (es decir, un posible valor de X). Debido a que esta es una distribución de probabilidad, cada una de las probabilidades debe ser un número entre 0 y 1, y las alturas de las barras deben sumar a 1 también.

 

Trabajando con la distribución binomial en R

R tiene una función llamada dbinom que calcula probabilidades binomiales para nosotros. Los principales argumentos a la función son

• x Este es un número, o vector de números, especificando los resultados cuya probabilidad estás tratando de calcular.
• tamaño Este es un número que le dice a R el tamaño del experimento.
• prob Esta es la probabilidad de éxito para cualquier ensayo en el experimento.

Entonces, para calcular la probabilidad de obtener cráneos, a partir de un experimento de ensayos, en el que la probabilidad de obtener una calavera en cualquier prueba es... bueno, el comando que usaría es simplemente este:

dbinom( x = 4, size = 20, prob = 1/6 )

0.202203581217173

Para darle una idea de cómo cambia la distribución binomial cuando alteramos los valores de\(\theta\) y\(N\), supongamos que en lugar de tirar dados, en realidad estoy volteando monedas. En esta ocasión, mi experimento consiste en voltear una moneda justa repetidamente, y el resultado que me interesa es el número de cabezas que observo. En este escenario, la probabilidad de éxito es ahora\(\theta = 1/2\). Supongamos que iba a voltear la moneda\(N=20\) veces. En este ejemplo, he cambiado la probabilidad de éxito, pero mantuve igual el tamaño del experimento. ¿Qué le hace esto a nuestra distribución binomial?

Figura\(\PageIndex{2}\): Dos distribuciones binomiales, que involucran un escenario en el que estoy volteando una moneda justa, por lo que la probabilidad de éxito subyacente es 1/2. En el panel (a), suponemos que estoy volteando la moneda N = 20 veces. En el panel (b) asumimos que la moneda es volteada N = 100 veces.

Bueno, como\(\PageIndex{2} (a)\) muestra Figura, el principal efecto de esto es desplazar toda la distribución, como cabría esperar. Bien, ¿y si volteamos una moneda\(N=100\) veces? Bueno, en ese caso, obtenemos Figura\(\PageIndex{2} (b)\). La distribución se mantiene aproximadamente en el medio, pero hay un poco más de variabilidad en los posibles resultados.

En este punto, probablemente debería explicar el nombre de la función dbinom. Obviamente, la parte “binom” viene del hecho de que estamos trabajando con la distribución binomial, pero el prefijo “d” probablemente sea un poco misterioso. En esta sección voy a dar una explicación parcial: específicamente, voy a explicar por qué hay un prefijo. En cuanto a por qué es una “d” específicamente, tendrás que esperar hasta la siguiente sección. Lo que está pasando aquí es que R en realidad proporciona cuatro funciones en relación con la distribución binomial. Estas cuatro funciones son dbinom, pbinom, rbinom y qbinom, y cada una calcula una cantidad diferente de interés. No sólo eso, R hace lo mismo por cada distribución de probabilidad que implemente. No importa de qué distribución estés hablando, hay una función d, una función p, una función r a y una función q.

Echemos un vistazo a lo que hacen las cuatro funciones. En primer lugar, las cuatro versiones de la función requieren que especifiques los argumentos size y prob: no importa lo que estés tratando de obtener R para calcular, necesita saber cuáles son los parámetros. No obstante, difieren en cuanto a lo que es el otro argumento, y cuál es la salida. Así que vamos a verlos uno a la vez.

• La forma d que ya hemos visto: se especifica un resultado particular x, y la salida es la probabilidad de obtener exactamente ese resultado. (la “d” es la abreviatura de densidad, pero ignora eso por ahora).
• La forma p calcula la probabilidad acumulada. Se especifica un cuantil particular q, y éste le indica la probabilidad de obtener un resultado menor o igual a q.
• La forma q calcula los cuantiles de la distribución. Se especifica un valor de probabilidad p, y éste le da el percentil correspondiente. Es decir, el valor de la variable para la que existe una probabilidad p de obtener un resultado menor que ese valor.
• La forma r es un generador de números aleatorios: específicamente, genera n resultados aleatorios a partir de la distribución.

Esto es un poco abstracto, así que veamos algunos ejemplos concretos. Nuevamente, ya hemos cubierto dbinom así que centrémonos en las otras tres versiones. Empezaremos con pbinom, y volveremos al ejemplo cráneo-dado. Nuevamente, estoy rodando 20 dados, y cada dado tiene una probabilidad de 1 en 6 de subir cráneos. Supongamos, sin embargo, que quiero saber la probabilidad de rodar 4 o menos cráneos. Si quisiera, podría usar la función dbinom para calcular la probabilidad exacta de rodar 0 cráneos, 1 calavera, 2 cráneos, 3 cráneos y 4 cráneos y luego sumar estos, pero hay una manera más rápida. En cambio, puedo calcular esto usando la función pbinom. Aquí está el comando:

pbinom( q= 4, size = 20, prob = 1/6)

0.768749218992842

Es decir, hay un 76.9% de posibilidades de que ruede 4 o menos cráneos. O, para decirlo de otra manera, R nos está diciendo que un valor de 4 es en realidad el percentil 76.9 de esta distribución binomial.

A continuación, consideremos la función qbinom. Digamos que quiero calcular el percentil 75 de la distribución binomial. Si nos quedamos con nuestro ejemplo de cráneos, usaría el siguiente comando para hacer esto:

qbinom( p = 0.75, size = 20, prob = 1/6 )

4

Hm. Aquí está pasando algo extraño. Vamos a pensarlo bien. Lo que la función qbinom parece estar diciéndonos es que el percentil 75 de la distribución binomial es 4, aunque vimos por la función que 4 es en realidad el percentil 76.9. Y definitivamente es la función pbinom la que es correcta. Te lo prometo. La rareza aquí viene del hecho de que nuestra distribución binomial realmente no tiene un percentil 75. En realidad no. ¿Por qué no? Bueno, hay un 56.7% de probabilidad de rodar 3 o menos cráneos (puedes escribir pbinom (3, 20, 1/6) para confirmar esto si quieres), y un 76.9% de probabilidad de rodar 4 o menos cráneos. Entonces hay un sentido en el que el percentil 75 debería estar “entre” 3 y 4 cráneos. ¡Pero eso no tiene ningún sentido en absoluto! No se puede tirar 20 dados y obtener 3.9 de ellos subir cráneos. Este problema se puede manejar de diferentes maneras: podrías reportar un valor intermedio (o valor interpolado, para usar el nombre técnico) como 3.9, podrías redondear hacia abajo (a 3) o podrías redondear hacia arriba (a 4).

La función qbinom redondea hacia arriba: si pides un percentil que en realidad no existe (como el 75 en este ejemplo), R encuentra el valor más pequeño para el que el rango del percentil es al menos lo que pediste. En este caso, dado que el “verdadero” percentil 75 (lo que sea que eso signifique) se encuentra en algún lugar entre 3 y 4 cráneos, R redondea y te da una respuesta de 4. Esta sutileza es tediosa, lo admito, pero por suerte solo es un problema para distribuciones discretas como el binomio. Las otras distribuciones de las que voy a hablar (normal,\(t\),\(\chi^2\) y\(F\)) son todas continuas, y así R siempre puede devolver un cuantil exacto cada vez que lo pides.

Por último, tenemos el generador de números aleatorios. Para usar la función rbinom, especifica cuántas veces R debe “simular” el experimento usando el argumento n, y generará resultados aleatorios a partir de la distribución binomial. Entonces, por ejemplo, supongamos que iba a repetir mi experimento de laminación de troqueles 100 veces. Podría obtener R para simular los resultados de estos experimentos usando el siguiente comando:

rbinom( n = 100, size = 20, prob = 1/6 )

1. 2
2. 4
3. 3
4. 5
5. 5
6. 1
7. 2
8. 7
9. 2
10. 1
11. 4
12. 3
13. 4
14. 1
15. 1
16. 4
17. 5
18. 2
19. 1
20. 2
21. 7
22. 6
23. 1
24. 1
25. 1
26. 2
27. 4
28. 4
29. 3
30. 6
31. 4
32. 3
33. 3
34. 3
35. 4
36. 2
37. 4
38. 4
39. 2
40. 1
41. 3
42. 2
43. 5
44. 4
45. 2
46. 4
47. 1
48. 1
49. 6
50. 2
51. 3
52. 1
53. 4
54. 1
55. 1
56. 3
57. 4
58. 4
59. 3
60. 6
61. 5
62. 5
63. 5
64. 3
65. 2
66. 4
67. 5
68. 5
69. 4
70. 4
71. 2
72. 1
73. 1
74. 5
75. 3
76. 6
77. 1
78. 6
79. 3
80. 4
81. 3
82. 5
83. 2
84. 7
85. 1
86. 5
87. 3
88. 2
89. 2
90. 0
91. 3
92. 4
93. 3
94. 6
95. 6
96. 4
97. 5
98. 1
99. 4
100. 3

Como puede ver, estos números son más o menos lo que esperaría dada la distribución que se muestra en la Figura\(\PageIndex{1}\). La mayoría de las veces ruedo en algún lugar entre 1 y 5 cráneos. Hay muchas sutilezas asociadas con la generación de números aleatorios usando una computadora, pero para los fines de este libro no necesitamos preocuparnos demasiado por ellas.