Curso CEC 4-6 Cap
Curso CEC 4-6 Cap
Curso CEC 4-6 Cap
Distribuciones de probabilidad
oas
Distribución binomial
k x = 0,1,..., n.
• A la variable aleatoria X sigue una ley binomial de parámetros 𝑛 y
𝑝 se la notará como 𝑋~𝐵(𝑛, 𝑝).
𝜇 = 𝐸 𝑋 = 𝑛𝑝
𝜎2 = 𝑉(𝑋) = 𝑛𝑝𝑞
oas
Distribución binomial
Ejemplo 4.1:
• En un proceso de fabricación que produce gran cantidad de
artículos, se sabe que en promedio 2% de ellos son defectuosos.
Los artículos son empacados en cajas de 10, y se quiere saber
cuál es la probabilidad de que no haya ningún artículo
defectuoso en cada caja. Si 𝑋 es el número de artículos
defectuosos por caja, entonces se quiere obtener 𝑃(𝑋 = 0), lo
cual es:
Ejemplo 4.1:
• Si se requiere saber cuál es la probabilidad de que cada caja
tenga exactamente un artículo defectuoso 𝑃(𝑋 = 1), entonces:
oas
Distribución binomial
Ejercicio 4.1:
1. La probabilidad de un artículo que sea defectuoso en un proceso
de manufactura continua es 0,06 ¿Cuáles son las respectivas
probabilidades de que 0,1,2,3,4 y 5 artículos sean defectuosos
de una muestra al azar de 5 artículos extraídos del proceso?
2. Supóngase que la tasa de infección de una enfermedad
contagiosa es del 25%. En una oficina hay 10 personas que se
vacunaron contra la enfermedad y ninguna se contagió.
a) Determine la probabilidad de que ninguna persona se hubiera
contagiado a pesar de que no se hubiera vacunado;
b) De este resultado, ¿deduce usted que la vacuna es efectiva?
3. Se conoce que en promedio la mortalidad de transporte de pollos
en pie desde granjas hasta la planta faenadora es del 0.30%.
Los pollos son transportados en camiones con capacidad de
5,000 unidades, se quiere saber cuál es la probabilidad de que
no haya ningún pollo muerto por cada viaje.
oas
Distribución de poisson
𝜇=𝐸 𝑋 =𝜆
𝜎2 = 𝑉(𝑋) = 𝜆
oas
Distribución de poisson
Ejemplo 4.2:
• La probabilidad de que un producto salga defectuoso es de
0,012. ¿Cuál es la probabilidad de que entre 800 productos ya
fabricados hayan 5 defectuosos?
9,65 −9, 6
P ( X = 5) = e = 0,0460
5!
oas
Distribución de poisson
Ejercicio 4.2:
1. Un ingeniero que labora en el departamento de control de
calidad de una empresa eléctrica inspecciona una muestra al
azar de 200 alternadores de un lote. Si el 2% de los
alternadores del lote están defectuosos. Cuál es la probabilidad
de que, en la muestra,
a) ninguno esté defectuoso,
b) uno salga defectuoso,
c) al menos dos salgan defectuosos,
d) más de tres estén con defectos.
f (x ) =
1 − ( x − )2 / 2 2
e
2 x (− , )
• La función de distribución correspondiente es la integral
x
F (x ) =
1 − (t − )2 / 2 2
2 e
−
dt
𝜎2 = 𝑉(𝑋)
oas
Distribución normal
f(x)
x F(x)
1
x
oas
Distribución normal
Algunas Propiedades:
1. Es simétrica respecto a su media, ;
2. La moda y la medina son ambas iguales a la media, ;
3. Los puntos de inflexión de la curva se dan para x=- y x=+.
4. Distribución de probabilidad en un entorno de la media:
• en un intervalo [-, +] se encuentra comprendida,
aproximadamente el 68,27% de la distribución;
• en un intervalo [-2, +2] se encuentra comprendida,
aproximadamente el 95,45% de la distribución;
• en un intervalo [-3, +3] se encuentra comprendida,
aproximadamente el 99,73% de la distribución.
Cálculo de probabilidades:
• Cuando es una distribución normal con parámetros 𝜇 = 0 y 𝜎2 =
1 , entonces a la distribución se le conoce como distribución
normal estándar 𝑋~𝑁(0,1). Estandarizar una variable es fácil,
puesto que si 𝑋 tiene una distribución normal con 𝜇 = 𝐸(𝑋) y 𝜎2 =
𝑉(𝑋), entonces la variable (estandarizada)
X −
Z=
se distribuye de manera normal estándar. La ventaja de
estandarizar es que cualquier probabilidad de interés, por ejemplo,
𝑃(𝑋 ≤ 𝑥) , se puede escribir en términos de la variable
estandarizada 𝑍 como:
X − x−
P ( X x ) = P = P(Z z )
oas
Distribución normal
Ejemplo 4.3:
• Se sabe que la dimensión de una pieza se distribuye normal con
media 𝜇 = 82 𝑚𝑚, y 𝜎 = 0,5. se desea calcular el porcentaje de
piezas que cumplen con especificaciones 82 ± 1. Lo que se
puede obtener calculando la siguiente probabilidad:
83 − 83 − 82
P( X 83) = P Z = P Z = P(Z 2,0)
0,5
oas
Distribución normal
Ejemplo 4.3:
• Para calcular 𝑃(𝑋 < 81) se procede de la misma manera.
oas
Distribución normal
Test de normalidad
• Los test de normalidad se aplican a conjuntos de datos para
determinar su similitud con una distribución normal. La hipótesis
nula es, en estos casos, si el conjunto de datos es similar a una
distribución normal, por lo que un valor- 𝑝 suficientemente
pequeño indica datos no normales.
• Prueba de Kolmogórov-Smirnov
• Test de Lilliefors
• Test de Anderson–Darling
• Test de Ryan–Joiner
• Test de Shapiro–Wilk
• Normal probability plot (rankit plot)
• Test de Jarque–Bera
• Test omnibús de Spiegelhalter
oas
Test de normalidad
Test de Anderson–Darling:
Donde:
• F es la función de distribución acumulada de la distribución
normal.
• Xi son las observaciones ordenadas.
Hipótesis:
Test de Anderson–Darling:
0.75 2.25
𝐴𝐷∗ = 𝐴𝐷 ∗ 1 + + 2
𝑛 𝑛
Rango Fórmula
2)
0.6 ≤ 𝐴𝐷 ∗ < 13 𝑒 (1.2937−5.709∗𝐴𝐷+0.0186∗𝐴𝐷
2)
0.34 ≤ 𝐴𝐷 ∗ < 0.6 𝑒 (0.9177−4.279∗𝐴𝐷−1.38∗𝐴𝐷
2)
0.2 ≤ 𝐴𝐷∗ < 0.34 1 − 𝑒 (−8.318+42.796∗𝐴𝐷−59.938∗𝐴𝐷
2)
0 ≤ 𝐴𝐷 ∗ < 0.2 1 − 𝑒 (−13.436+101.14∗𝐴𝐷−223.73∗𝐴𝐷
oas
Test de normalidad
Ejercicio 4.3:
• Generar una variable con 100 números aleatorios con una
distribución normal con parámetros 25 y 0.83 𝑁~(25; 0.70 ).
Realizar el test de normalidad de Anderson–Darling.
oas
TEMA V
Índices de capacidad
oas
Proceso con doble especificación:
Ejemplo 5.1:
• Una característica de calidad importante en la fabricación de una
salsa de tomate es su peso en el momento del embotellado, que
para cierto tipo de salsa debe ser de 395 cm3 con una tolerancia
de ± 10 cm3. El peso es el resultado del proceso de embotellado,
por lo que este proceso debe garantizar que el peso esté entre la
especificación inferior EI = 385 y la superior ES = 405, con un
valor ideal o nominal de N = 395. Para monitorear el correcto
funcionamiento del proceso de embotellado, cada media hora se
toman cinco botellas y se pesan. De acuerdo con las mediciones
hechas en el último mes, en el que el proceso ha estado
trabajando de manera estable, se tiene que la media y la
desviación estándar del proceso son = 395.88 y = 4.14,
respectivamente. Con base en esto se requiere tener una
evaluación de la manera que el proceso ha estado cumpliendo
con especificaciones. Haremos esto con base en los índices de
capacidad.
oas
Capacidad de un proceso
ES − EI
C p=
6
donde representa la desviación estándar del proceso, y ES y EI
son las especificaciones superior e inferior para la característica
de calidad.
oas
Proceso 6
• Un proceso 6 tiene Cp = 2.
Tabla 5.2 Índices Cp, Cpi y Cps en términos de la cantidad de piezas malas; bajo normalidad
y proceso centrado en el caso de doble especificación. oas
Cálculo del índice Cp
405 − 385 20
C p= = = 0,80
6 * 4,14 24,87
− EI ES −
C pi = C ps =
3 3
• El índice Cpk está definido por:
C pk = mín C pi , C ps
oas
Cálculo de los índices Cpi, Cps, Cpk
C pk = mín0,88;0,73 = 0,73
oas
Índice K
2 * ( − N )
K= *100
ES − EI
• De esta forma el índice K es igual a la media del proceso menos
el valor nominal para la característica de calidad N, dividida entre
la mitad de la distancia de las especificaciones. Es decir, el índice
K mide qué tan descentrado está el proceso en función de la
mitad de la amplitud de las especificaciones, y al multiplicar por
100 se convierte en una medida porcentual.
oas
Índice K
oas
Cálculo del índice K
2 * (395,88 − 395)
K= *100 = 8,80%
405 − 385
oas
Índice Cpm
ES − EI
C pm = = + ( − N )
2 2
6
• El índice Cpm compara el ancho de las especificaciones con 6;
pero no sólo toma en cuenta la variabilidad del proceso, a través
de ², sino que también toma en cuenta el centrado a través de
(-N)². De esta forma, si el proceso está centrado, es decir, si =
N, entonces el Cp, Cpk y Cpm son iguales.
oas
Cálculo del índice Cpm
405 − 385
C pm = = 0,79
6 * 4,24
oas
Capacidad de largo plazo
R
ˆ =
d2
• La capacidad de largo plazo se calcula con muchos datos
tomados en un período largo para que los factores externos
influyen en el proceso, y S se estima mediante la desviación
estándar de todos los datos:
ˆ = S oas
Índices Pp y Ppk
Índice Pp
• Indicador del desempeño potencial del proceso, que se calcula en
forma similar al índice Cp pero usando la desviación estándar de
largo plazo.
ES − EI
P p=
6 L
Índice Ppk
• Indicador del desempeño real del proceso, que se calcula en
forma similar al índice Cpk pero usando la desviación estándar de
largo plazo.
− EI ES −
P pk = mín ,
3 L 3 L
oas
Índice Z
ES − − EI
Z s= Z i=
• La capacidad de un proceso medida en términos del índice Z es
igual al valor más pequeño de entre Zs y Zi, es decir:
Z = mínZ s , Z i
oas
Índice Z
Z m= Z C − Z L
• El índice Zm representa la habilidad para controlar la tecnología. Si
no se conoce Zm, se asume un valor de 1,5.
3C pk = Z C 3P pk = Z L
oas
Índice Z
Ejemplo 5.2:
• En un proceso de envasado de azúcar en un ingenio azucarero
se tiene como especificación del contenido de los costales 50 kg,
con una tolerancia de 0,6 kg. De esta manera forma, la
especificación inferior es EI = 49,4 kg, y la superior ES = 50,6 kg.
De acuerdo con los datos históricos se tiene que la media del
proceso es = 50,01 y la desviación estándar es 0,2 kg. De aquí
que,
Z = mín2,95;3,05 = 2,95
Tabla 5.3 Calidad de corto y largo plazo en términos de Cp, ZC, ZL y PPM.
oas
Métrica 6 para atributos
oas
Métrica 6 para atributos
− EI ES −
C pi = C ps =
3 3
donde y son la media y la desviación estándar de la
característica de calidad, (EI) y (EI) son las especificaciones
inferior y superior. La interpretación de los índices Cpi y Cps es
parecida al Cp y Cpk.
30 − 10,2
C ps = = 1,29
3 * 5,1
• De acuerdo a la tabla 5.2, tiene un porcentaje fuera de
especificación cercano a 0,005%.
oas
Proceso con sólo una especificación
105 − 100
C pi = = 0,26
3 * 6,5
que es un valor pequeño, por lo que el proceso es incapaz de
cumplir con la especificación, el porcentaje aproximado de
productos que tienen espesor menor a 100 micras está entre
18,406% y 27,425%, por lo que es necesario hacer esfuerzos
serios para mejorar ese proceso.
oas
Estudio real (integral) de capacidad
oas
Estudio real (integral) de capacidad
Datos:
Medidas de dispersión:
Medidas de forma:
Porcentaje fuera de
especificaciones:
Histograma:
Índice Z:
Conclusiones finales:
Datos: oas
TEMA VI
Cartas de control
oas
Variación:
Leyes de variación:
oas
Variación:
C15
A10
A11
A12
B10
C10
C11
C12
C13
C14
A6
C5
A1
A2
A3
A4
A5
A7
A8
A9
B1
B2
B3
B4
B5
B6
B7
B8
B9
C1
C2
C3
C4
C6
C7
C8
C9
Figura 6.1 Gráfico de ventas semestrales.
Entendiendo la variación:
¿Cómo?
Maquinaría Métodos
Variable de salida
Figura 6.3 La variabilidad de un proceso. oas
Variación:
oas
Variación:
Ejercicio 6.1:
• Entendiendo la variación: Caligrafía.
oas
Variación:
Causas comunes:
• La variación por causas comunes (o por azar) es aquella que
permanece día a día, lote a lote; y es aportada de forma natural
por las condiciones de las 6M. Esta variación es inherente a las
actuales características del proceso y es resultado de la
acumulación y combinación de diferentes causas que son
difíciles de identificar y eliminar, ya que son inherentes al sistema
y la contribución individual de cada causa es pequeña; no
obstante, a largo plazo representan la mayor oportunidad de
mejora.
oas
Variación:
Ejercicio 6.2:
• Previniendo resultados (proceso “bajo control”): Suma en el
lanzamiento de dos dados.
oas
Variación:
80
60
Línea Central (LC)
40
20
0
Límite de Control Inferior (LIC)
Línea de acción:
• El mejoramiento, en el caso de causas comunes es más compleja.
oas
Variación:
Causas especiales:
• La variación por causas especiales (o atribuibles) es causada por
situaciones o circunstancias especiales que no están de manera
permanente en el proceso. Las causas especiales por su
naturaleza relativamente discreta a menudo pueden ser
identificadas y eliminadas si se cuenta con los conocimientos y
condiciones para ello.
oas
Variación:
80
60
Línea Central (LC)
40
20
0
Límite de Control Inferior (LIC)
Línea de acción:
• Deben ser atacadas inmediatamente.
Entendiendo la variación:
• “El 94% de los problemas y oportunidades de mejora, se deben
a causas comunes. Apenas 6% se deben a causas especiales”
(W. E. Deming1).
Error 1 / Error 2:
Variación de W Causas
Límite de Control Superior (LSC) especiales
Causas
Línea Central (LC) comunes
oas
Límites de control
LCI = W - 3 W
LC = W
LCS = W + 3 W
oas
Límites de control
oas
Tipos de cartas de control
oas
Tipos de cartas de control
X1 R1 X2 R2 X3 R3
ത
Figura 6.8 Operación de una carta 𝑋-R.
oas
ഥ -𝑹
Carta de control 𝑿
Figura 6.9 La carta 𝑋ത detecta cambios significativos en la media del proceso. Cuando la curva se desplaza
la carta manda una o varias señales de fuera de control.
oas
ഥ -𝑹
Carta de control 𝑿
X 3 X
R d2 3
3 X = 3 = R = A2 R
n d2 n
R 3 R
R d 3
LCI = R − 3d 3 = 1 − 3 R = D3 R
d2 d 2
LC = R
R d 3
LCS = R + 3d 3 = 1 + 3 R = D4 R
d2 d 2
Ejemplo 6.1:
• Una característica de calidad en el proceso de fabricación de
bebidas gaseosas es el contenido de sólidos disueltos expresado
en grados Brix; en particular una bebida gaseosa “negra” debe
tener un contenido de sólidos disueltos de 12.0 expresado en
grados Brix, con una tolerancia de ± 1.5 grados Brix. De esta
forma la especificación o tolerancia inferior es EI = 10.5 grados
Brix y la superior ES = 13.5 grados Brix.
oas
ഥ -𝑹
Carta de control 𝑿
• ധ
Se tiene que 𝑋=12.08, ത
𝑅=1.70 y para n =5, 𝐴2 =0.577; entonces:
LCS = 13.06
𝑿
ഥ
LCI = 11.10
Subgrupo
Figura 6.11 Carta 𝑋ത del ejemplo 6.1. oas
ഥ -𝑹
Carta de control 𝑿
LCS = 3.58
Rango
LCI = 0.00
R 1.70
ˆ = = = 0.73
d 2 2.326
• Los límites reales son: LRI = 9.9 y LRS = 14.3, que por cierto
exceden las especificaciones (EI=10.5, ES=13.5), y por tanto el
proceso es incapaz.
oas
ഥ -𝑹
Carta de control 𝑿
ES − EI 13.5 − 10.5
Cp = = = 0.69
6 6 * 0.73
oas
ഥ -𝑹
Carta de control 𝑿
Ejercicio 6.3:
• En una industria alimenticia se quiere garantizar que la
concentración mínima de grasa de un producto sea de 1.8%. En
la tabla 6.2 se muestran los datos obtenidos para un estudio
inicial, con tamaño de subgrupo de 4.
ത
a. Realice un estudio de estabilidad mediante la carta 𝑋-𝑅.
oas
ഥ -𝑹
Carta de control 𝑿
Datos:
LCS
A
B
C
C
B
A
LCI
LCS
LCI
LCS
LCI
LCS
LCI
oas
Interpretación de las cartas de control
LCS
LCI
LCI
Ejemplo 6.2:
• En una fábrica de bolsas de plástico un aspecto importante de
calidad es la dimensión de las bolsas. En una fase del proceso
una máquina debe cortar automáticamente las bolsas, la medida
ideal es 300 mm, con una tolerancia de más menos 5 mm. Utilice
una carta 𝑋ത para los últimos 40 subgrupos de tamaño n = 5 y
calcule el índice de inestabilidad, St.
303
LCS = 302.3
302
301
𝑿
ഥ
300
299
LCI = 229.3
298
297
1 2 3 4 5 6 7 8 9 10111213141516171819202122232425262728293031323334353637383940
Subgrupo
Figura 6.19 Carta 𝑋ത para el ejemplo 6.2.
4
St = *100 = 10%
40 oas
Índice de inestabilidad, St
• En cuanto a capacidad:
ES − EI 305 − 295
Cp = = = 1.5
6 6 *1.12
oas
ഥ -𝑹
Carta de lecturas individuales 𝑿
oas
ഥ -𝑹
Límites de control para la carta 𝑿
• ത
Los límites de control para una carta de lecturas individuales 𝑋-R,
se obtienen de la siguiente manera:
• Donde:
• 𝑋ത es la media de las mediciones de los subgrupos;
• 𝑅ത es la media de los rangos móviles de orden 2;
• D3 y D4 son constantes que dependen del tamaño de
subgrupo o muestra.
oas
ഥ -𝑹
Carta de lecturas individuales 𝑿
Ejemplo 6.3:
• En la producción de tequila se miden varias características de
calidad a lo largo del proceso, ya sea por regulaciones oficiales o
por criterios propios de cada empresa. Una de esas
características es el grado Brix residual después de la molienda
(cantidad de azúcar que queda en el bagazo del agave), que
mide la eficacia del proceso de molienda. Ésta es una variable
del tipo entre más pequeña mejor, y el valor máximo tolerado que
se tiene en una empresa en particular es ES = 3,5%. Después de
moler cada lote se determina el grado Brix residual, por lo que se
considera un proceso lento, que será más apropiado analizarlo
con una carta de individuales. En la tabla se aprecian los datos
para los últimos 40 lotes molidos, y se ha agregado la columna
para el rango móvil de orden 2, que se obtienen del rango entre
los dos datos consecutivos más recientes. A continuación, se
hará un estudio inicial y se analizará la estabilidad del proceso de
molienda.
oas
ഥ -𝑹
Carta de lecturas individuales 𝑿
Datos:
• ത
Se tiene que 𝑋=1.95 ത
y 𝑅=0.43. De aquí que los límites de control
preliminares para el Brix son:
LCS = 1.95 + 2.66(0.43) = 3.10
LC = 1.95
LCI = 1.95 − 2.66(0.43) = 0.81
LCS = 3.10
Brix
LCI = 0.81
LCS = 1.41
Rango móvil
LCI = 0.00
X = 1.92 R = 0.39
con lo que la desviación estándar del proceso es = 0.35.
oas
ഥ -𝑹
Carta de lecturas individuales 𝑿
LCS = 2.97
Brix
LCI = 0.87
Lote
Figura 6.22 Carta 𝑋ത individual ajustada del ejemplo 6.3. oas
ഥ -𝑹
Carta de lecturas individuales 𝑿
Estabilidad:
St =
1
(100) = 2.5
40
• Que se puede considerar relativamente bajo y por tanto la
estabilidad de proceso de molienda se puede considerar
moderadamente estable.
Capacidad:
ES − 3.5 − 1.92
C ps = = = 1.51
3 3(0.39)
Ejercicio 6.5:
• La pureza de un proceso químico es medida para cada lote, y los
datos obtenidos se registran en una carta de individuales.
Actualmente se tiene que los límites de control para dicha carta
son los siguientes:
LCS = 0.92
LC = 0.86
LCI = 0.80
oas
ഥ -𝑹
Carta de lecturas individuales 𝑿
Conteste:
1. Bajo el supuesto de que el proceso está en control estadístico,
explique de manera sencilla el significado práctico de estos
límites.
2. Obtenga la deviación estándar del proceso.
3. Estime los límites reales, ¿coinciden en este caso con los límites
de control? ¿Por qué?
4. Se sugiere que el límite de control inferior sea igual a 0.84, ya
que se tiene la exigencia por parte de la administración de que
ésa sea la pureza mínima tolerable del proceso. ¿Es correcta
esta sugerencia? Argumente.
5. Si la pureza de los últimos 10 lotes es la siguiente: 0.90, 0.85,
0.83, 0.82, 0.84, 0.84, 0.85, 0.81, 0.83, 0.82. Grafique estos
datos en la carta y señale si ha pasado algo.
oas
Cartas de control para atributos
Atributos
Binomial Poisson
(0,1) (0,1,2,3,…)
p np u c
oas
Carta p – proporción de defectos
p (1 − p )
LCS = p + 3
n
LC = p
p (1 − p )
LCI = p − 3
n
• Donde:
n es el tamaño de subgrupo;
𝑝ҧ es la proporción promedio de artículos defectuosos en el
proceso.
oas
Carta p – proporción de defectos
Ejemplo 6.4:
• En una empresa del ramo alimenticio, mediante ciertas máquinas
se empaquetan salchichas en sobres o paquetes. La forma de
evaluar si el subproceso de empaquetado se hizo de manera
adecuada, es haciendo una inspección visual de los paquetes
para determinar que se satisfagan diferentes atributos de calidad,
por ejemplo: que la cantidad de salchichas dentro del sobre sea
la correcta o que dentro del sobre no quede aire. Se han tenido
problemas con este último atributo, por lo que hay un operador
que a la vez que está pendiente de la máquina empaquetadora,
hace la inspección para separar los paquetes con aire. Los que
son destruidos y las salchichas se regresan a la máquina para
que vuelvan a ser empaquetados.
Datos:
275 23942
p= = 0.011 n= = 599
23942 40
0.011(1 − 0.011)
LCS = 0.011 + 3 = 0.023
599
LC = 0.011
0.011(1 − 0.011)
LCI = 0.011 − 3 = 0.000
599 oas
Carta p – proporción de defectos
LCS = 0.023
𝒑
ഥ
LCI = 0
Subgrupo
oas
Proceso tipo D (inestable e incapaz)