¿Qué tal van las clases, Bartolo? Me pregunta mi barbero. Bien... Dando probabilidad

1. 5. Distribuciones discretas • ¿Qué tal van las clases, • Bartolo? Me pregunta mi barbero. • Bien... Dando probabilidad • y estadística... Respondo. • ¡Ah! Probabilidad... Yo • suelo jugar a la lotería... • Dice mientras me pasa la cuchilla. • Cuando compro un número, tal y como yo lo veo, hay dos posibilidades: ganar o perder. De modo que tengo un 50% de • probabilidad de ganar y un 50% de perder. • -¡Muy bien, Ricardo! Respondo, mientras pienso que no es bueno contradecir a nadie que tenga una navaja en mi cuello...

2. Distribución de Bernoulli Experimento de Bernoulli: solo son posibles dos resultados: éxito o fracaso. Podemos definir una variable aleatoria discreta X tal que: éxito  1 fracaso  0 Si la probabilidad de éxito es p y la de fracaso 1 - p, podemos construir una función de probabilidad: Un típico experimento de Bernoulli es el lanzamiento de una moneda con probabilidad p para cara y (1-p) para cruz.

3. Función de distribución:

4. Ejercicio: Calcular la esperanza y la varianza de la distribución de Bernoulli.

5. Distribución binomial La distribución binomial aparece cuando estamos interesados en el número de veces que un suceso A ocurre (éxitos) en n intentos independientes de un experimento. P. ej.: # de caras en n lanzamientos de una moneda. Si A tiene probabilidad p (probabilidad de éxito) en un intento, entonces 1-p es la probabilidad de que A no ocurra (probabilidad de fracaso).

6. Experimento aleatorio: n = 3 lanzamientos de una moneda. Probabilidad de éxito en cada lanzamiento (cara) = p. Probabilidad de fracaso en cada lanzamiento (cruz) = 1- p = q.

7. Supongamos que el experimento consta de n intentos y definamos la variable aleatoria: X = Número de veces que ocurre A. En nuestro ejemplo: X = Número de veces que sale cara. Entonces X puede tomar los valores 0, 1, 2, ... n. Si consideramos uno de estos valores, digamos el valor x , i.e. en x de los n intentos ocurre A y en n - x no. Entonces la probabilidad de cada posible ordenación es pxqn-x y existen idénticas ordenaciones.

8. Distribución binomial para n = 5 y distintos valores de p, B(5, p) La función de probabilidad P(X = x) será la distribución binomial:

10. Tablero de Galton o quincunx Sir Francis Galton(1822-1911) Quincunx

11. Ejercicio: ¿Cuál es la probabilidad de que en una familia de 4 hijos exactamente 2 sean niñas?

12. Ejercicio: Si una décima parte de personas tiene cierto grupo sanguíneo, ¿cuál es la probabilidad de que entre 100 personas escogidas al azar exactamente 8 de ellas pertenezcan a este grupo sanguíneo?

13. ¿Y si la pregunta es 8 como máximo?

14. Calcula la probabilidad de obtener al menos dos seises al lanzar un dado cuatro veces. p = 1/6, q = 5/6, n = 4 Al menos dos seises, implica que nos valen k = 2, 3, 4. P(2) +P(3) + P (4)

15. Ejercicio: Supongamos que la probabilidad de encontrar una estrella de masa m* >10 M en un cúmulo estelar joven es del 4%. ¿Cuál es la probabilidad de que en una muestra escogida al azar, entre 10 miembros del cúmulo encontremos 3 estrellas con m* >10 M?

16. Chuck-a-luck: Elige un número entre 1 y 6. Lanzas 3 dados. Si el número que has elegido sale en los 3 dados cobras 3 euros. Si sale en 2 cobras 2 euros. Si sale en un dado cobras 1 euro. Y si no sale en ninguno, pagas 1 euro. ¿Es un juego justo?

17. Características de la distribución binomial • Media • = E(X) = n p • = 5 · 0.1 = 0.5 • = 5 · 0.5 = 0.25 n = 5p = 0.1 P(X) .6 .4 .2 0 X 0 1 2 3 4 5 Desviación estándar n = 5p = 0.5 P(X) .6 .4 .2 X 0 0 1 2 3 4 5

19. Distribución multinomial Cuando hay más de dos acontecimientos posibles (A1, A2, A3 ...) con probabilidades p1, p2 , p3 ... constantes y tales que:

20. Un método de diagnóstico tiene 3 resultados posibles: positivo (P), negativo (N) y dudoso (D). Se sabe que, en la población, el 10% de los sujetos son positivos, el 70% negativos y el resto dudosos. ¿Qué probabilidad hay de, en una muestra de 5 individuos, obtener exactamente 1 positivo, 1 negativo y 3 dudosos ?

21. Distribución geométrica Consideremos el siguiente experimento: Partimos de un experimento de Bernoulli donde la probabilidad de que ocurra un suceso es p (éxito) y la probabilidad de que no ocurra q = 1- p (fracaso). Repetimos nuestro experimento hasta conseguir el primer éxito. Definimos la variable aleatoria X, como el número de fracasos hasta que se obtiene el primer éxito. Entonces:

22. p(x) x Función de distribución:

23. ¿Cómo simular de manera sencilla en el ordenador una variable aleatoria binomial X? Sumando n variables aleatorias independientes cuyos valores pueden ser 1 o 0, con probabilidad p y 1-p respectivamente. ¿Y cómo simular una distribución geométrica de parámetro p? Una manera es generar una secuencia de números aleatorios en [0, 1) con la función rnd, y paramos cuando obtengamos un número que no exceda a p, que es el equivalente al primer éxito. El problema es que si p es pequeño, en promedio se necesitan 1/p pasos de tiempo y se consume mucho tiempo de cómputo.

24. Una forma alternativa con tiempo de cómputo independiente del valor de p sería: Sea q = 1-p y definamos la variable Y como el menor entero que satisface: Entonces tenemos: De modo que Y está distribuida geométricamente con parámetro p.

25. Para generar Y, basta con que despejemos de:

27. Un acontecimiento ocurre, en la población, en el 10% de los casos. ¿Qué tamaño de muestra debo tomar para tener una probabilidad del 95% de obtener al menos un éxito ?

28. Distribución binomial negativa(dePascal o de Pólya) Consideremos el siguiente experimento: Partimos de un experimento de Bernoulli donde la probabilidad de que ocurra un suceso es p (éxito) y la probabilidad de que no ocurra q = 1- p (fracaso). Repetimos nuestro experimento hasta conseguir el r-ésimo éxito. Definimos la variable aleatoria X, como el número de fracasos x hasta que se obtiene el r-ésimo éxito. Entonces: El último tiene que ser un éxito. Se denomina binomial negativa porque los coeficiente provienen de la serie binomial negativa:

29. Distribución binomial negativa(dePascal o de Pólya) La distribución binomial negativa también se puede definir como el número de pruebas x hasta la aparición de r éxitos. Como el número de pruebas x, en este caso, contabiliza tanto los éxitos como los fracasos se tendría según ésta definición que:

30. Disponemos de una moneda trucada con probabilidad de cara igual a p=0.25. La lanzamos hasta que obtenemos 2 caras. La distribución del número de lanzamientos x será: P(x) x

31. Elegir al azar con reemplazo Elegir al azar con reemplazo significa que escogemos al azar un elemento de un conjunto y lo regresamos para elegir de nuevo al azar. Esto garantiza la independencia de las elecciones y nos lleva a una distribución binomial. Si una caja contiene N bolas de las cuales A son rojas, entonces la probabilidad de escoger al azar una bola roja es: p = A/N. Si repetimos el experimento sacando n bolas con reemplazo la probabilidad de que x sean rojas es: (Una distribución binomial)

32. Elegir al azar sin reemplazo Elegir al azar sin reemplazo significa que no devolvemos el elemento elegido al azar al conjunto. De modo que las probabilidades de la siguiente elección dependen de las anteriores. Si repetimos el experimento anterior sacando n bolas sin reemplazo, ¿cuál será ahora la probabilidad de que x sean rojas? Para calcular los casos favorables observa que: N = A + (N – A). De las A bolas rojas tomaremos x y de las N – A bolas no rojas tomaremos n – x.

33. Distribución hipergeométrica

34. Queremos seleccionar al azar dos bolas de una caja que contiene 10 bolas, tres de las cuales son rojas. Encuentra la función de probabilidad de la variable aleatoria : X = Número de bolas rojas en cada elección (con y sin reemplazo). Tenemos N = 10, A = 3, N - A = 7, n = 2 Escogemos con reemplazo: Escogemos sin reemplazo:

35. Hipergeométrica Binomial N = 24 n = 5 X = 8 p = 8/24 =1/3 Observa que si N, A, N-A son grandes comparados con n no hay gran diferencia en qué distribución empleemos. La distribución binomial es una aproximación aceptable a la hipergeométrica si n < 5% de N. n = 5 Error P(x) P(x) x -0.0289 0 0.1028 0.1317 0.0133 1 0.3426 0.3292 0.0397 2 0.3689 0.3292 -0.0065 3 0.1581 0.1646 -0.0148 4 0.0264 0.0412 -0.0028 5 0.0013 0.0041 N = 240 n = 5 X = 80 p = 80/240 =1/3 n = 5 x P(x) P(x) Error -0.0028 0 0.1289 0.1317 0.0014 1 0.3306 0.3292 0.0035 2 0.3327 0.3292 -0.0004 3 0.1642 0.1646 -0.0014 4 0.0398 0.0412 -0.0003 5 0.0038 0.0041

36. Distribución de Poisson Cuando en una distribución binomial el número de intentos (n) es grande y la probabilidad de éxito (p) es pequeña, la distribución binomial converge a la distribución de Poisson: donde np =  Observa que si p es pequeña, el éxito es un “suceso raro”. La distribución de Poisson, junto con la uniforme y la binomial, son las distribuciones más utilizadas.

37. Un proceso poissoniano es aquél compuesto de eventos discretos que son independientes en el espacio y/o en el tiempo. Por ejemplo la llegada de fotones a un detector. Usemos la distribución binomial para modelar el proceso. Podemos dividir el intervalo de tiempo en el que ocurre el proceso en n subintervalos suficientemente pequeños, como para asegurarnos que a lo sumo se produce un evento en cada subintervalo. De modo que en cada subintervalo, o se producen 0 o 1 ocurrencias. A lo sumo llega un fotón en cada subintervalo o ninguno. De modo que podemos entender el proceso como un experimento de Bernoulli. Para determinar p, podemos razonar de la siguiente manera:

38. En promedio se producirán λt ocurrencias en un intervalo de tiempo t. Si este intervalo se divide en n subintervalos, entonces esperaríamos en promedio (usando Bernoulli): np ocurrencias. Así: λt = np, p = λt / n. Sin pérdida de generalidad supongamos que t = 1 y que X es la variable aleatoria = número total de ocurrencias. Sabemos que: Observa que para n grande P(X = 0) es aproximadamente e-λ. Además para n grande (y por tanto p muy pequeño):

39. Tenemos entonces la siguiente ecuación iterada: Que nos proporciona:

40. Bombas sobre Londres en la II Guerra Mundial (Feller) Supón que vivías en uno de los 100 bloques que aparecen en la gráfica inferior. La probabilidad de que una bomba cayera en tu bloque era 1/100. Como cayeron 400 bombas, podemos entender el número de impactos en tu bloque como el número de éxitos en un experimento de Bernoulli con n = 400 y p = 1/100. Podemos usar una Poisson con λ=400 1/100=4: 400 bombas Observado Predicho 10 x 10

41. Características de la distribución de Poisson Media = 0.5 P(X) .6   E ( X )   .4 .2 0 X Desviación estándar 0 1 2 3 4 5    = 6 P(X) .6 .4 Nota: el máximo de la distribución se encuentra en x  .2 0 X 0 2 4 6 8 10

42. La distribución de Poisson se obtiene como aproximación de una distribución binomial con la misma media, para ‘n grande’ (n > 30) y ‘p pequeño’ (p < 0,1). Queda caracterizada por un único parámetro μ (que es a su vez su media y varianza).  =  = n p =  Distribución de Poisson para varios valores de .

43. Si la probabilidad de fabricar un televisor defectuoso es p = 0.01, ¿cuál es la probabilidad de que en un lote de 100 televisores contenga más de 2 televisores defectuosos? La distribución binomial nos daría el resultado exacto: El suceso complementario Ac: No más de 2 televisores defectuosos puede aproximarse con una distribución de Poisson con  = np = 1, sumando p(0) + p(1) + p(2).

44. Una distribución de Poisson con μ = 10. P(7) = 107 e−10 / 7! = 0.09, es decir 9% Parece muy baja. Comparemos con el valor de máxima probabilidad que ocurrirá para x = 10: μ = 10 P(10) = 1010 x e−10 / 10! = 0.125, es decir 12.5% Las probabilidades poissonianas para un número de eventos dado, son siempre pequeñas, incluso en el máximo de la distribución de probabilidad. La señal promedio recibida en un telescopio de una fuente celeste es de 10 fotones por segundo. Calcular la probabilidad de recibir 7 fotones en un segundo dado.

45. Si asumimos que un minuto puede dividirse en muchos intervalos cortos de tiempo independientes y que la probabilidad de que un coche entre en uno de esos intervalos es p – que para un intervalo pequeño será también pequeño – podemos aproximar la distribución a una Poisson con  = np = 2. El suceso complementario “entran 3 coches o menos” tiene probabilidad: y la respuesta es 1 – 0.857 = 0.143 Si en promedio, entran 2 coches por minuto en un garaje, ¿cuál es la probabilidad de que durante un minuto entren 4 o más coches?

47. Ley de Benford

48. Primer dígito significativo 0,0174 299.959 1,74 ·10-2 [1,74] = 1 2,99959 ·105 [2,99959] = 2

49. Las barras negras representan las frecuencias de aparición como primer dígito significativo (d = 1,2,3,...,9) en una lista de N = 201 constantes físicas.

50. En barras blancas aparecen las frecuencias de aparición como primer dígito de los números 1 a 9 en el tamaño en bytes de N = 1.295.777 ficheros.