E06 DPoisson R01
E06 DPoisson R01
ISBN 9968-9634-0-3
[email protected]
MENÚ 2
Introducción
2 3
P (0) e P (1) e P ( 2) e
P (3)
(2)(3)
e
2
5
Procesos Poisson.
En general, un proceso aleatorio es un procedimiento físico controlado
completamente o en partes por algún tipo de mecanismo aleatorio. Puede
ser una sucesión de lanzamientos de dados, mediciones de la calidad de
productos que salen de una línea de producción, las vibraciones de las alas
del Concorde, el ruido de una señal de radio, el movimiento de partículas en
una solución o muchos otros fenómenos.
Lo que caracteriza a tales procesos, es su dependencia del tiempo. El hecho
de que ciertos eventos suceden o no, a intervalos regulares o en un intervalo
continuo de tiempo.
Situaciones como: la ocurrencia de imperfecciones en un rollo de casimir, la
medición de la radiación de un contador Geiger, el número de artículos
defectuosos de una línea de producción, las entradas de correos electrónicos
en un computador servidor.
El modelo matemático que aproxima la distribución de estos datos es la
Distribución Poisson.
6
Probabilidad de x éxitos en tiempo T.
Problema 1. De Llegadas.
Una empresa que ofrece el servicio de conexiones de
computadoras particulares a la Red Mundial De Computadoras
(INTERNET) tiene muchos problemas con los Virus Cibernéticos.
El esfuerzo es de tal magnitud que lo sobreasa y decide contratar
los servicios de una empresa dedicada a la “Cacería de Virus
Cibernéticos”
El primer paso de esta empresa fue determinar el número de nuevas
agresiones por día. Después de 10 días de monitoreo decidió iniciar
el estudio.
Los programadores y analistas de esta empresa dedicada a la
“Cacería de Virus”, como sistema, suelen hacerse preguntas que
involucran probabilidades, para enfrentar el problema al tenor de
¿cuáles son las probabilidades?
a) De que entren cuatro nuevos virus por día.
b) Que entren 10 virus en dos días consecutivos.
8
Obteniendo el Promedio
Problema 2. Fallas.
Soldadura
En los procesos industriales mecanizados la cantidad de
producto fabricado es muy regular y podría usarse como
variable x la cantidad de envases en lugar del tiempo y decir
3.000 envases, 5.000 envases o la producción de 15.000 envases,
aunque sea preferible usar el tiempo como unidad.
La tabla de frecuencias, promedios y proporciones es:
Evento Frecuencia Frecuencia 3
2 1 2
3 0 0 x 0,23
p 0,023
n 10
Suma de Días 100
Suma Total Inperfecciones 23 Estadísticos obtenidos del
Promedio 0,23
Tamaño de muestra 10 ejercicio de la empresa.
Proporción 0,023
Las Respuestas del Problema
12
2
a) Si la proporción de fallas es de 0,023 por 3 n p 3 0,023 0,069
hora, en 3 horas se espera un promedio de:
λ 1e λ 0,0691 e 0,069
P( x 1)
Se quiere saber la probabilidad de obtener 1! 1
POISSON(1;0,069;0) 0,0644
una falla en ese tiempo. Para esto se aplica
la probabilidad Poisson:
b) Al menos dos imperfecciones en 5 horas. 5 n p 5 0,023 0,069
Primero se obtiene el promedio para 5 0,1150 e 0,115 0,1151 e 0,115
P ( x 2) 1
horas: 0! 1!
1 - 0,8914 0,1025
Al solicitar la probabilidad de 2 o más 1 - POISSON(1;0,115;1) 0,0061
soldaduras malas es más fácil obtenerla por
probabilidades complementarias:
15 15 0,023 0,345
c) Para esta debe obtenerse el promedio de
0,3450 e 0, 345 0,3451 e 0,345
fallas en 15 horas: P( x 2)
0! 1!
Y probabilidad de: 0,8914 0,1025
POISSON(1;0,345;1) 0,9526
13
La Distribución de Frecuencias
El análisis de los datos históricos de un año indicó que la entrada de casos
a emergencias puede considerarse aleatoria.
De acuerdo a los médicos responsables el ingeniero dividió el día
en tres clases: de 6 a 10 de la mañana, de las 16 a las 21 horas y el resto del
día. La tabulación de emergencias se hace en el siguiente cuadro así como
la estimación de promedios por hora: 4,107; 5,408 y 2,022 emergencias por
hora, para la mañana, tarde y resto del día respectivamente.
Emergencias FRECUENCIAS OBSERVADAS EXPANSIÓN DE TOTALES = f * x
x Mañana Tarde Resto Día Mañana Tarde Resto Día
0 5 0 48 0 0 0
1 24 5 97 24 5 97
2 44 34 97 88 68 194
3 86 41 70 258 123 210
4 56 66 36 224 264 144
5 61 50 9 305 250 45
6 48 46 8 288 276 48
7 25 48 0 175 336 0
8 10 40 0 80 320 0
9 3 18 0 27 162 0
10 3 17 0 30 170 0
Sumas 365 365 365 1499 1974 738
Promedios 4,107 5,408 2,022
17
Los Promedios de Emergencias por Hora
fx
i 0
i i
50 241 ... 3(10) 1.499
Para la mañana: xm 10
4,107
5 24 ... 3 365
f
i 0
i
10
fx i i
00 51 ... 17(10) 1.974
Para la tarde: xt i 0
10
5,408
0 5 ... 17 365
f
i 0
i
10
fx i i
480 971 ... 0(10) 738
xr i 0
2,022
Resto del día: 10
48 97 ... 0 365
f
i 0
i
18
Prueba de Bondad de Ajuste para la
Mañana.
En la prueba de Bondad de Ajuste se comparan mediante la prueba de 2
las frecuencias observadas contra las esperadas obtenidas usando la
distribución Poisson.
Emergencias Frecuencias Prueba
x Observada Esperada Ajuste Chi-Cuadr.
0 5 6,0 -1,0 0,1690
En donde la frecuencia esperada
1 24 24,7 -0,7 0,0184 se obtiene multiplicando la
2 44 50,7 -6,7 0,8765
3 86 69,4 16,6 3,9941 probabilidad Poisson por el
4 56 71,2 -15,2 3,2484
5 61 58,5 2,5 0,1078 número de observaciones
6 48 40,0 8,0 1,5850
7 25 23,5 1,5 0,0974
8 10 12,1 -2,1 0,3511
9 3 5,5 -2,5 1,1378
10 3 3,5 -0,5 0,0767 2
10
foi fei
5 6,0 24 24,7
2 2
...
3 3,5 11,6622 2
11 1
Sumas 365 365 11,6622 i 0 fei 6,0 24,7 3,5
Probabilidad 0,3083
1 1
Día.
Procediendo de la misma forma en la distribución de emergencias
para el resto del día.
La prueba de Chi-cuadrada se
Emergencias Frecuencias Prueba
x Observada Esperada Ajuste Chi-Cuadr. estimó en:
0 48 48,3 -0,3 0,0022
1 97 97,7 -0,7 0,0052
10
fo i fe i
χ 211 1
2 97 98,8 -1,8 0,0322 i 0 fe i
3 70 66,6 3,4 0,1760
4 36 33,7 2,3 0,1636
48 48,3
97 97,7 ... 0 0,1 6,2314
2 2 2
Probabilidad.
Se ha probado que la distribución Poisson aproxima apropiadamente la
distribución de frecuencias de las emergencias en el día. Las estimaciones
incluyendo el promedio se muestran en el siguiente cuadro y gráfico.
PROMEDIOS Y ESTIMACIONES POISSON PARA
Emergencia 4,107 5,408 2,022 3,846
Estimación Poisson
x Mañanas Tardes R. Día Promedio
0 0,0165 0,0045 0,1324 0,0214 0,30
1 0,0676 0,0242 0,2677 0,0822
0,25
2 0,1388 0,0655 0,2706 0,1580
Probabilidades
3 0,1900 0,1181 0,1824 0,2026 0,20
4 0,1951 0,1597 0,0922 0,1948
5 0,1602 0,1727 0,0373 0,1498 0,15
6 0,1097 0,1557 0,0126 0,0960
0,10
7 0,0644 0,1203 0,0036 0,0527
8 0,0330 0,0813 0,0009 0,0254 0,05
9 0,0151 0,0489 0,0002 0,0108
10 0,0062 0,0264 0,0000 0,0042 0,00
11 0,0023 0,0130 0,0000 0,0015 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17
12 0,0008 0,0059 0,0000 0,0005 Emergencias por hora
13 0,0002 0,0024 0,0000 0,0001
14 0,0001 0,0009 0,0000 0,0000 Mañanas Tardes R. Día Promedio
15 0,0000 0,0003 0,0000 0,0000
16 0,0000 0,0001 0,0000 0,0000
17 0,0000 0,0000 0,0000 0,0000
respuestas
Para la mañana la probabilidad de que lleguen más de 2,5 pacientes por
hora es:
P(x 2,5)= 0,2229+ [0,4149 – 0,2229](0,05) = 0,3179
Y
P(x > 2,5) = 1 – 0,3179 = 0,6821
Particulares
En los hospitales hay muchas maneras de cubrir estas situaciones. Nuevos
datos indican que el tiempo promedio que necesita un doctor de 24 minutos
para atender a un paciente se convierte en 12 cuando tiene 6 pacientes (5
esperando turno).
En estos momentos, la enfermera de recepción clasifica a los
pacientes en tres grupos, los que son verdaderas emergencias, los que
presentan una emergencia media y los que pueden esperar por tiempo
largo. Pasa a los pacientes del grupo 1, después a los del dos, así hasta
finalizar la cola. Con esta nueva información calcule la probabilidad de que
lleguen mas pacientes de los 60 / 12 = 5 que pueden atenderse por hora.
Las probabilidades se obtiene fácilmente de las acumulativas o las
individuales siendo de:
23,18% para la mañana; 45,53% para la tarde; 1,74% para el resto del
día y de 19,12% en promedio.
Información suficiente para tomar decisiones en lo tocante a las
emergencias del hospital.
Planes de Control de la 25
Lata.
El departamento de control de la calidad de la empresa ENVASES
HL, S, A. Tiene por costumbre tomar tres muestras de 10 envases
por jornada de 8 horas de trabajo en un plan de muestreo
secuencial con iniciación aleatoria (ver sección de técnicas de
muestreo). Al final del día hace un recuento de las fallas que
acumula en un reporte que todos los días actualiza con los últimos
33 reportes en una hoja de control de la calidad que incluye,
además, un control del proceso recopilando el conjunto de factores
que concurren cuando el proceso de pone Fuera de Control
Estadístico.
Al inicial el día, el supervisor de la producción tiene en su mano las
dos cartas de control de las máquinas de soldado para que efectúe
los ajustes pertinentes y tome las acciones que considere necesarias.
El reporte de fallas en la HE.
Las Distribuciones: Binomial y 27
Poisson.
En el control de fallas están involucradas dos distribuciones que aunque
son de la misma familia su utilidad es diferente:
La Binomial: Que se utiliza para establecer planes de control sobre
las proporciones o porcentajes;
La Poisson: Que se utiliza para establecer planes de control sobre el
número de fallas.
Usualmente se utiliza sólo una de ellas dependiendo de los objetivos del
control.
De manera general las cartas de control son alternativas pictóricas para
mostrar el comportamiento estadístico de los errores mediante intervalos
confiables con probabilidad específica en una sucesión de muestras para
determinar el momento en que el sistema de producción deja su variación
aleatoria debido a la presencia de un factor que provoca la falla. Siendo la
ubicación y control de estos factores el objetivo del sistema.
La Binomial para Proporciones. 28
0,07
allá de la norma al
0,06
0,05
0,04
0,03
0,02
sobrepasar el Límite
0,01
0,00
Promedio
LIC Superior de Control.
1 3 5 7 9 11 13 15 17 19 21 23 25 27 29 31 33
Muestras
La Poisson para Número de 29
Fallas.
Los límites para la Carta de Control 99% Confiable para
observaciones se define como:
LIC λ z 0,7 2,576 0,70 0,30 LSC λ z 0,7 2,576 0,70 2,85
Número de Fallas
2,5
línea de color café oscuro. 2
Es evidente que el control 1,5
1
de fallas es más estricto que el de 0,5
proporciones pues 5 muestras 0
1 3 5 7 9 11 13 15 17 19 21 23 25 27 29 31 33
presentan situaciones fuera de Muestras
control estadístico.
El Límite Inferior de Control únicamente tiene fines informativos para este
caso, pues es obvio que debe ser 0.
Propiedades de la Distribución 30
Poisson
A estas alturas el estudiante está en capacidad de deducir las dos
propiedades de las variables con distribución Poisson:
Primera: El valor de la varianza de la distribución es muy próxima al
valor de la media ;
Segunda: Si una serie de variables independientes X1, X2, X3,,, cada
una con distribución Poisson con medias 1, 2, 3,,, su suma tendrá
una distribución Poisson con una media igual a la suma de las medias
(1 + 2 + 3,,,).
Debe verse la distribución Poisson como una variante de la Binomial, sólo
qué ésta tiene dos alternativas para el promedio de medias y varianzas:
Para números:
np; 2 np(1 p )
Para proporciones y porcentajes:
p; 2 p(1 p )
Por estas dos propiedades se ha visto que la distribución de Poisson es útil para
aproximarse al fenómeno cuando el promedio es pequeño, aún cuando el valor
de n no está bien definido y aún, sí los valores n y p se presume que varían de una
muestra a otra.
El ejemplo que ha servido desde hace mucho tiempo para ejemplificar el uso de la
distribución de Poisson se debe a Ladislaus von Bortkewitch (1.898) que se refiere
al número de hombres en la caballería Prusiana que mueren en un año por coces
de caballos. Si se tienen N = 200 observaciones una por cada 10 batallones en 20
años, Cualquier día un hombre está expuesto a la pequeña probabilidad de que
sea coceado, pero no está claro que valor tiene n y también p que puede ser
variable para cada individuo.
Problema 5: Tradicional del Proceso 32
Poisson
Una empresa de la construcción reporta los siguientes accidentes
laborales de 155 ingenieros entre 30 y 35 años durante su primer año de
trabajo:
Accidentes de Trabajo 0 1 2 3 4 o más
Nº de Ingenieros 80 61 13 1 0
Frecuencias
Nº Accidentes Frecuencia Frecuencia
x Observada por evento
0 80 0
1 61 61
2 13 26
3 1 3 4
4 0 0 fx i i
80 0 61 1 ... 0 4 90
Estadísticos i 0
4
0,5806
80 61 ... 4 150
Nº Observaciones
Nº de Accidentes
155
90
f
i 0
i
95%.
Frecuencias Esperadas:
x e 0,5806 2 e 0,5806
Nº Accidentes Probabilidad Frecuencias Bondad fei n ; Ej : fe2 155 14,6
x! 2 !
x Poisson Observada Esperada Ajuste
0 0,5595 80 86,7 0,5220
1 0,3249 61 50,4 2,2488 Y la prueba de 2:
2 0,0943 13 14,6 0,1795 80 86,7 2 61 50,4
2
1 1
4, 5909 8 1 2 4 ,5909
F4,5909; 8-1 Y0
0
4,5909 2 e d 0,7090
= 2
Nº Accidentes Probabilidad Estimadores El Promedio:
x Poisson Promedio Varianza 7
0 0,559537 0,0000 0,1886 pi xi 0,559537 0 0,324893 1 0,000002 7
1 0,324893 0,3249 0,0571 i 0
2 0,094324 0,1886 0,1900 0,5806
3 0,018256 0,0548 0,1069
4
5
0,002650
0,000308
0,0106
0,0015
0,0310
0,0060
La varianza:
6 0,000030 0,0002 0,0009 n
i 0
Estimadores 0,5806 0,5806
0,324893(1 0,5806) 2 ... 0,000002(7 - 0,5806) 2
0,5806
Normal
Una empresa procesadora E.vaca ha decidido cambiar
de leche de sus
máquinas envasadoras en cartón sanitario inducidos por los vendedores
de maquinas que aseguran una mayor precisión de llenado necesaria para
cubrir las normas internacionales que piden una variación promedio
inferior a 5% y cuando más un 2% de envases por debajo de la norma. Se
estableció un programa de control de la calidad en el que se muestrean 5
unidades por cuatro veces, en una jornada de trabajo de 8 horas
iniciando la toma de la primera muestra de forma aleatoria eligiendo un
número entre el minuto 0 y el 120 para la primera muestra, las siguientes
se toman en secuencia de 2 horas. Los resultados de 100 muestreos, o 20
días se presentan en la HE. Los estadísticos obtenidos directamente:
N
x i
9
Promedio: μ i 1
0,018
N 500
N
x μ
2
i
Varianza: σ 2 i 1 p 1 p 0,018
N
El Promedio es Suficiente 37
95%
Aún corriendo el riesgo de ser reiterativo, es necesario valuar, sí las
frecuencias esperadas se pueden aproximar mediante la
distribución Poisson para n = 5.
x e 0,09 2 e 0, 09
Evento Probabilidad Frecuencias Prueba P ( xi ) ; P ( x 2) 0,037
x Poisson Observada Esperada B. Ajuste x! 2 1
0 0,9139 92 91,4 0,0040
1 0,0823 7 8,2 0,1826 fei P( xi )m; fe2 0,037 100 0,4
2 0,0037 1 0,4 1,0718
3 0,0001 0 0,0 0,0111
26 1
92 91,4 7 8,2
2
2
...
0 0,00.. 1,2697
2
También
Proporción.
se puede probar la 1 1
pˆ P 0,018 - 0,02
2 N ; z 2 500 -0,1597
hipótesis usando directamente z
pq
c
0,02 0,98
las proporciones: N 500
Control
La muestra 2: Alcanzan la zona A
Carta de Control de Envasadora
por la parte superior. No se valora
1.001,5
Zona A
como falla; acción, buscar la
1.001,0
Zona B
causa).
Volumen Envasado
1.000,5
Zona C
La muestra 26: Alcanza la zona A
1.000,0
Zona C
inferior y hay una observación
999,5 considerada falla; acción, buscar
999,0
Zona B la causa.
Zona A
998,5 Las muestras 3 y 4 se presentan en
1
3
5
7
9
13
15
17
19
21
23
25
27
29
11
Prefacio
Desde el primer momento debe imaginar que la charla está preparada para
ser expuesta por el profesor en el aula y que el estudiante cuenta con una
copia del “capítulo” en su computadora personal.
La charla le irá preparando en bases teóricas para enfrentar problemas
que en su momento fueron parte de una investigación y han sido resueltos,
algunos de ellos se han modificado ligeramente con fines didácticos.
El problema que se comenta en la CHARLA se ha resuelto paso a paso en
la Hoja Electrónica que además le proporcionará un conjunto de resultados
obtenidos mediante modelos de simulación que el estudiante puede definir
de una manera parcial y con los que puede “entrenar” el uso de la
herramienta computacional e intentar conclusiones y recomendaciones.
La presentación y los comentarios están ligados de manera que puede
interactuar con ellos.
Además puede obtener una copia impresa del capítulo o de todo el libro en
un formato de pantalla que se edita en una hoja tamaño carta.
Manuel Pontigo Alvarado:
ITCR, Diciembre 2004.