Graficas de Control Arvelo
Graficas de Control Arvelo
Graficas de Control Arvelo
Estudios realizados:
Ingeniero Industrial. UCAB Caracas 1968
Máster en Estadística Matemática CIENES , Universidad de Chile 1972
Cursos de Especialización en Estadística No Paramétrica Universidad de
Michigan 1982
Doctorado en Gestión Tecnológica: Universidad Politécnica de Madrid 2006 al
Presente
PROLOGO
El uso de métodos estadísticos en la gestión de calidad esta universalmente
aceptado y reconocido, por todas las instituciones y organismos que certifican
la calidad, tanto en la producción de bienes como en la prestación de
servicios.
Así por ejemplo, las normas ISO en todas sus versiones, hacen referencia en
varias de sus cláusulas, al uso de procedimientos de muestreo en el
aseguramiento de la calidad, tanto en la recepción de materias primas, como
en el control del proceso.
La recopilación de todos estos procedimientos constituye una disciplina
conocida como “Control Estadístico de Procesos” ó por sus siglas en inglés
SPC (Statistical Process Control), de obligatorio estudio en todas las
especialidades de Ingeniería.
Algunas de sus aplicaciones son las siguientes:
1. Gráficas de control, que se utilizan para detectar oportunamente
fluctuaciones irregulares que pudiera presentar un proceso.
2. Muestreo para la aceptación, que se utiliza para decidir acerca del
ingreso de lotes de materias primas que llegan a un proceso.
3. Diseño de Experimentos, que se utiliza para identificar las variables más
influyentes dentro de un proceso.
4. Análisis de Regresión, que se utiliza para pronosticar los valores de
calidad del producto terminado, en función de las variables controlables
del proceso.
5. Análisis de Confiabilidad, que se utiliza para el establecimiento de
garantías, y para decidir acerca de las políticas de mantenimiento a
seguir.
En Venezuela, su uso no es reciente, y las normas de calidad venezolanas
también sugieren el empleo de métodos estadísticos, tal como se evidencia en
la norma COVENIN – ISO9001: 1995.
Sin embargo, a pesar de lo difundido de los métodos, y de su elevado grado de
aceptación, ellos presentan desde el punto de vista teórico un punto débil, que
consiste en apoyarse en ciertos supuestos de distribución para modelar el
comportamiento del proceso, especialmente el supuesto de normalidad.
Graficas de Control 3
Angel F Arvelo
GRAFICAS DE CONTROL
En caso de que todos los puntos obtenidos caigan dentro de los límites de
control, estos pasan a ser definitivos; si algunos se salen, es necesario
investigar la causa asignable que actuó, eliminar los puntos, y recalcular los
nuevos límites.
Según el Western Electric Handbook (1956), además de un punto fuera de los
límites de control, existen otros síntomas que pueden revelar que el proceso se
encuentra fuera de control, tales como:
Graficas de Control 8
Angel F Arvelo
Al igual que la anterior, estas dos gráficas como casi todas las de Swehart
suponen que las muestras han sido tomadas como “subgrupos racionales”., y
en ellas los límites de control son obtenidos con el criterio tres sigmas bajo el
supuesto de normalidad.
mediante la aplicación del criterio tres sigmas presenta los siguientes límites de
P(1 − P)
L.S.C= P + 3
n
control para p̂ : Central = P
L.I.C= P − 3 P(1 − P)
n
• Grafica NP: Esta es también otra gráfica por atributos que tiene el
mismo objetivo de la anterior, con la diferencia de que utiliza como estadístico
de control al número “X” de no conformes presentes en el subgrupo racional.
Los límites de control para “X” por criterio tres sigmas y utilizando la
L.S.C =n P + 3 n P(1 − P)
aproximación normal a la binomial son: Central = n P
L.I.C =
n P − 3 n P(1 − P)
Dado que la distribución exacta de “X” es conocida cuando el proceso está bajo
control, pues sigue una Binomial con parámetros “n” y “P ”; en esta gráfica es
posible obtener los límites de control con criterio probabilístico, que usualmente
resultan más estrechos que con criterio tres sigmas.
Se utiliza como alternativa del gráfico “C” cuando la unidad de muestreo esta
C
formada por “n” unidades de producción, y se define U = .
n
U
L.S.C= U + 3
n
Los límites de control para U resultan: Central = U
L.I.C
= U+3
U
n
Es de hacer notar que en todos los gráficos de Shwehart el tamaño “n” del
subgrupo racional pudiera ser variable, lo que ocasionaría limites de control
variables, situación ésta de escaso valor práctico.
que el proceso esté bajo control, aproximarse a cero en la medida en que “m”
aumente.
Si la media del proceso llegara a sufrir un desajuste, se presentará un
desequilibrio entre los desvíos positivos y los negativos, y la suma acumulada
se incrementará positivamente si el deslizamiento es hacia un valor mayor que
“µ o ” , o negativamente en caso contrario.
La siguiente gráfica muestra un gráfico CUSUM en donde la media del proceso
a partir de la muestra No 16, ha sufrido un presunto deslizamiento hacia un
valor mayor que el nominal, y la suma acumulada ha comenzado a crecer
indefinidamente.
Graficas de Control 14
Angel F Arvelo
El principal problema que tropezó este gráfico inicial CUSUM fue la ausencia
de unos límites de control para la suma acumulada que permitieran detectar
cuando el proceso se ha desajustado.
La solución de este problema condujo a dos variantes del gráfico conocidas
como el CUSUM algorítmico y la plantilla o mascarilla V.
Para construir estos gráficos, es necesario definir un valor “k” llamado “valor de
referencia”, que representa el mayor deslizamiento en la media del proceso que
podría tolerarse, y a partir del cual la suma acumulada se considera
significativa y en consecuencia el proceso fuera de control.
El valor de “k” también depende de la probabilidad con que se quiera detectar
un deslizamiento, y su metodología de cálculo es distinta en el CUSUM
algorítmico y para la plantilla V.
El CUSUM algorítmico, considera dos sumas acumuladas, la de los desvíos
positivos llamada C i + y la de los negativos C i - , y la gráfica de control toma la
siguiente apariencia:
Graficas de Control 15
Angel F Arvelo
Si se aplica sucesivamente esta fórmula para Z i-1 , Z i-2 etc., puede verse
fácilmente que “Z i ” resulta un promedio exponencialmente ponderado entre
todas las observaciones, en donde la más reciente resulta con una ponderación
mayor, y las más antiguas con otra mucho menor.
j= i−1
Zi = λ ∑ (1 − λ ) j Xi− j + (1-λ )iµ0
j=0
ν
probabílistico del 99% resultan:
χ0.995;
2
ν
L.I.C = σ0
ν
Para procesos en donde se requiere un control por atributos, la Graficas EWMA
también ofrece una alternativa cuando el parámetro a controlar es el número de
defectos, el cual se quiere mantener controlado en un valor “µ o ”.
En este caso, al igual que en la gráfica “C” de Shwehart, se supone que el
número de defectos sigue una Distribución Poisson,
El valor del EWMA se mantiene al igual que antes: Z i = λX i + (1-λ) Z i-1 ; siendo
X i el número de defectos encontrados en la i-ésima unidad inspeccionada,
Z o = µ o , y 0 < λ ≤ 1.
λµo
L.S.C = µ0 + A 1 − (1 − λ )2
2−λ
Los límites de para “Z i ” en este caso son: Central = µ0
L.I.C = µ − A λµo 1 − (1 − λ )2
2−λ
0
σ
L.S.C =µ0 + 3 w
Para todas las siguientes: Central = µ0
σ
L.I.Ci =µ0 − 3
w
En la mayoría de los casos prácticos se selecciona w = 5, y la apariencia de la
gráfica es como sigue:
Graficas de Control 20
Angel F Arvelo
Estudios hechos por simulación, han demostrado que la gráfica de control con
promedios móviles es más efectiva que la de Shwehart para detectar
corrimientos leves en la media del proceso.
Una variante de este tipo de gráficas consiste en calcular el promedio móvil en
forma ponderada, concediéndole un mayor peso a la última observación, y una
ponderación menor a las anteriores.
la elipse resulta con ejes paralelos a los ejes coordenados, pero si no lo son
entonces σ12
2
≠ 0 y la elipse resulta con ejes oblicuos cuya inclinación depende
Gráfico de Control Bidimensional para dos Gráfico de Control Bidimensional para dos
variables independientes variables con correlación positiva
En caso de que las variables a controlar sean tres ó mas, el cálculo del
estadístico de control χ2 debe hacerse mediante la siguiente expresión
T2 =
n(X − X)´ S-1 (X − X)
Graficas de Control 23
Angel F Arvelo
donde X y S son las estimaciones del vector de medias “µ” y de la matriz “Σ”
de varianzas y covarianzas respectivamente.
Al igual que el anterior, el gráfico de control T2 de Hotelling presenta un solo
límite superior de control sin línea central, y sobre él se representan los
sucesivos valores de T2 para las diferentes muestras examinadas.
Cuando se presentan puntos fuera de control en gráficos multivariantes, la
interpretación es mucho mas difícil que en el caso univariante, pues no se
puede identificar la variable que ha ocasionado esta situación.
Algunos autores sugieren complementar con gráficos individuales para cada
una de las variables a fin de facilitar esta identificación. Este procedimiento
además de largo puede fallar, especialmente cuando la señal fuera de control
no se debe a una sola de las variables, si no a la interacción entre varias de
ellas.
Otro procedimiento útil para analizar un punto fuera de control es descomponer
el valor de la T2 de Hotelling en distintas componentes que reflejen la acción de
cada una de las variables. ( véase referencias bibliográficas anteriores)..
Hayter y Tsui (1994) desarrollaron también un procedimiento basado en
intervalos de confianza simultáneos, que permite identificar las causas que
originaron la señal fuera de control.
EJERCICIOS PROPUESTOS
Los siguientes ejercicios están diseñados para el curso de “Control Estadístico
de la Calidad” que el Prof. Arvelo dicta en la Universidad Central de Venezuela,
y complementan a los que se proponen al final de cada capítulo del texto
“Control Estadístico de la Calidad” de Douglas Montgomery.
1º) El contenido de jabón en unas cajas sigue una distribución normal. Si el 80%
de las cajas contienen menos de 429,20 gramos de jabón , y el 90% de las cajas
contienen más de 418,60 gramos:
a) Determine la media y la desviación típica de la distribución.
b) ¿Cual es el porcentaje de cajas cuyo contenido se encuentra en el intervalo
(420 ± 5 ) gramos?
Solución: a) µ =425 grms , σ =5 grms b) 47,72 % .
2º) Una pieza debe caer dentro de la especificación (200 ± 4) mm., para que sea
considerada como buena. Las piezas que caen por debajo del límite inferior de la
especificación, son consideradas como defectuosas y deben ser desechadas;
mientras que, las que caen por encima del límite superior de la especificación,
pueden ser corregidas, y llevadas a los límites de la especificación.
Para producir estas piezas, se dispone de una máquina que las fabrica según
una Distribución Normal, con media 200 mm., y desviación típica de 5 mm .
El costo de producción para una pieza es de Bs. 7 , y el costo de la corrección es
de Bs.2 ; mientras que el precio de venta de la pieza, es de Bs. 15 .
a) ¿Cual es el ganancia esperada en la producción de una pieza? .
b) ¿Cual es la probabilidad de que en un lote de 10 piezas, alguna caiga fuera
de los límites de especificación?
Solución : a) Bs. 4,40 b) 0,9960
4º) El consumo diario de un cierto producto sigue una distribución gamma con
media 6 y varianza 18. ¿Cuál es la probabilidad de que en un día, el consumo
sea más de 9?.
Solución: 0,1991
6°) El contenido de refresco en una botella sigue una Distribución Normal con
media 200 cc y desviación típica 20 cc. Al final del proceso de producción existe
un dispositivo de control que rechaza las botellas con un contenido inferior a 180
cc.
a) ¿Qué porcentaje de la producción resultará rechazada?
b) ¿Cuál es la probabilidad de que una botella que aprobó el control, contenga
más de 210 cc ?
Solución: a) 15,87 % b) 0,3667
7°) La duración de unas pilas sigue una Distribución Normal con media 20 horas.
Se sabe que el 80% de estas pilas sobrepasan las 15 horas de funcionamiento.
Un radio lleva 4 pilas conectadas en serie, es decir, que si falla alguna, entonces
falla el radio.
a) ¿Cuál es la probabilidad de que el radio falle antes de las 18 horas de uso? .
b) ¿ Por cuánto tiempo podría garantizarse la operación del radio, para tener
probabilidad 0,95 , de cumplir con lo garantizado ? .
Solución: a) 0,8393 b) 6,70 horas aproximadamente
9º) El contenido de unas cajas de cereal sigue una Distribución Normal con una
desviación típica de 25 gramos.
Una muestra aleatoria de 36 cajas arrojó una media de 500 gramos.
¿ Con qué nivel de confianza puede decirse que el contenido medio de todas las
cajas de cereal, está entre 495 y 505 gramos .
Solución: 76,99 %
Graficas de Control 27
Angel F Arvelo
19°) La vida de una componente sigue una distribución exponencial con media
10 horas. Se dispone de 30 componentes que serán usadas una tan pronto falle
la anterior.
a) Calcule la confiabilidad para un servicio de 335 horas.
b) ¿Cuántas componentes se necesitan para que la confiabilidad sea de 0,99 por
lo menos? .
Solución : a) 0,25
20°) Una orden de producción establece que las especificaciones para una pieza
deben ser de ( 20.00 ± 1.00 ) mm.
El ingeniero de producción debe decidir entre dos máquinas para ejecutar esta
orden.
La siguiente tabla da los parámetros de producción para cada máquina y sus
costos:
Media ( mm) Desviación típica (m m) Costo ($ / pieza)
Máquina A 19,95 0,85 1,00
Máquina B 20,02 0,50 1,20
Si la pieza es defectuosa, el fabricante pierde el costo de producirla; y si resulta
conforme a las especificaciones la puede vender en $ 1,80
¿Cuál de las dos máquinas debe elegir para producir las piezas?.
Solución: Conviene la máquina “B” porque da una ganancia esperada mayor. $
0,52 para B y $ 0, 37 para A.
22°) Un Ingeniero Industrial afirma que si una cierta operación se efectúa por un
cierto método, se logrará una economía en tiempo de 10 segundos en promedio,
por lo menos.
Se realizan observaciones de tiempos sobre grupos diferentes de obreros, por
los dos procedimientos, el convencional y el propuesto por él, obteniendo:
Group Statistics
Mann-Whitney Test
Ranks
TIPO N Mean Rank Sum of
Ranks
DURACION 1 10 6,20 62,00
2 10 14,80 148,00
Total 20
Test Statistics
DURACION
Graficas de Control 32
Angel F Arvelo
Mann-Whitney 7,000
U
Wilcoxon W 62,000
Z -3,250
Asymp. Sig. ,001
(2-tailed)
Test Statistics
DURACION
Mann-Whitney 8,500
U
Wilcoxon W 29,500
Z -2,184
Asymp. Sig. (2- ,029
tailed)
Exact Sig. ,026
[2*(1-tailed
Sig.)]
Graficas de Control 33
Angel F Arvelo
4°) Los siguientes datos representan los tiempos de operación en horas, para
tres tipos de calculadoras científicas de bolsillo antes de que requieran recarga:
A 4.9 6.1 4.3 4.6 5.3
Calculadora B 5.5 5.4 6.2 5.8 5.5 5.2 4.8
C 6.4 6.8 5.6 6.5 6.3 6.6
8°) Una compañía de taxis está interesada en identificar cual de dos marcas de
cauchos, le conviene instalar en sus unidades.
Se instalan ocho cauchos marca “A” y nueve marca “B” en taxis del mismo
modelo, y se observa el kilometraje obtenido en miles de kilómetros. Los
resultaron fueron:
Marca
A 37.7 34.8 35.0 28.7 35.4 36.9 32.8 33.9
B 38.6 39.6 40.1 36.2 35.7 38.4 36.7 39.2 35.5
Sin suponer normalidad, analice si existe diferencia significativa entre los
kilometrajes ofrecidos por las dos marcas de cauchos.
¿Cuál marca recomendaría Ud.?
9°) Las especificaciones exigen que la resistencia media de unas cadenas debe
ser de 150 Kgs. como mínimo.
Se examina una muestra de 100 cadenas, encontrándose los siguientes
resultados:
Media Muestral 146,80 Kgs
Desviación estándar muestral 2,80 Kgs.
Número de cadenas con resistencia < 150 Kgs 62
Número de cadenas con resistencia = 150 Kgs 4
Suma de los rangos para las cadenas < 150 3350
Kgs
Analice a un nivel de significación del 1%, si estas cadenas cumplen con lo
especificado.
Haga su análisis de la siguiente manera:
a) Suponiendo normalidad en la resistencia de las cadenas.
b) Sin suponer normalidad en la resistencia de las cadenas, y por dos
procedimientos diferentes.
584166357636 ; [email protected]
i = 40 i = 40
El resumen de los resultados obtenidos es: ∑ X i = 2980
i =1
∑S
i =1
i = 120
a) Analice si este gráfico está correctamente construido con criterio tres sigma,
y en caso de no estarlo, halle los límites correctos.
b) El departamento de control de calidad ha sugerido que es más conveniente
controlar el proceso con un gráfico ( X , S ), y subgrupos de tamaño 20. Halle
los límites de control para este nuevo gráfico.
c) Si las especificaciones para una pieza son: 600 ± 20 , estime el porcentaje
de piezas no conformes que actualmente produce el proceso, halle los
coeficientes C P , C PK , y obtenga sus conclusiones acerca de la capacidad del
proceso.
d) De ocurrir un desajuste que deslice la media hacia 606, ¿Cuál es la
probabilidad de que el nuevo gráfico ( X , S ) lo detecte en la primera muestra
subsecuente?.
X = 139.49 ; R = 1.175
a) Calcule los límites de control para la gráfica ( X , R) , y analice si el proceso
ha estado bajo control estadístico durante este lapso.
b) De ocurrir un deslizamiento en la media del proceso hacia 139.00 , ¿ cuál
es la probabilidad de detectarlo, a lo más en la segunda muestra después
de ocurrido?.
c) Suponga que las especificaciones son 140 ± 2. Estime el porcentaje de
piezas fuera de especificación.
d) Estime los índices de capacidad del proceso, y haga las recomendaciones
que considere conveniente.
Solución: a) 4 y 27 b) 7 y 26 c) 0,8228
584166357636 ; [email protected]
1°) Lotes de tamaño 100 fabricados por un proveedor nuevo, van a ser
examinado con un plan de muestreo simple para la aceptación por atributos.
El AQL convenido es de 1.5%
a) ¿Cuál es el plan sugerido por la Norma MIL STD 105 D, utilizando un nivel
general de inspección II. ?.
b) ¿Cuál es el riesgo del productor para este plan?.
c) Suponga que bajo la vigencia de este plan, las muestras de los últimos 10
lotes examinados han arrojado el siguiente número de defectuosas:
0 0 0 10 0 0 1 1 1
Aplique los criterios establecidos en la norma MIL STD 105 D, para decidir bajo
qué tipo de inspección debe ser examinado el próximo lote.
Solución: a) n = 32 c = 1 b) 0,0830 c) Normal
4o) Un primer lote de tamaño 4000 fabricado por un proveedor nuevo, va a ser
sometido a un plan de muestreo para la aceptación.
El AQL convenido es de 0,15 %, mientras que el LTPD es de 1%.
Las especificaciones exigidas para el producto son 80.00 ± 0.50
a) Si se decide un plan simple por atributos, ¿cuál es el plan sugerido por la
Norma MIL STD 105 D, utilizando un nivel general de inspección II. ?.
b) ¿Cuál es el riesgo del productor para este plan simple?. Escriba la expresión
exacta, y luego aproxímela.
c) Si se decide un plan doble por atributos, ¿cuál es el plan sugerido por la
Norma MIL STD 105 D, utilizando este mismo nivel general de inspección II?
d) ¿Cuál es el riesgo del consumidor para este plan doble?.
e) Suponga que Ud. considera que estos riesgos no resultan convenientes, y
Ud. quisiera que fuesen de 0,05 para el productor, y de 0,10 para el
consumidor.
e.1) Diseñe el plan de muestreo simple correspondiente.
e.2) Diseñe el plan de muestreo doble correspondiente, con n 2 = 2 n 1 .
f) Si se decide un plan por variables, ¿cuál es el plan sugerido por la Norma
MIL STD 414 bajo el nivel de inspección IV?
g) Si la muestra sugerida en f), arroja: X = 79,80, S = 0,16. ¿ Qué decisión
debe tomarse respecto al lote ?
Solución: a) n=315 Ac=1 Re=2 b) 0,0819 c) n 1 =n 2 =200 , Ac 1 =0 Re 1 = 2 Ac 2 =1
Re 2 = 2 d) 0,1702 e.1) n=507 c= 2 e.2) n 1 = 400 n 2 = 800 c 1 = 1 c 2 =3
f) n= 50 M= 0,503 g) Rechazar el lote
e) Suponga que las especificaciones para este producto exigen 10.00 ± 0.08, y
Ud. decide un plan por variables en lugar de uno por atributos, con Nivel
General III, ¿cuál es el plan sugerido por la Norma MIL STD 414
f) Si la muestra sugerida en e) , arroja: X = 10,04 , S = 0,02 . ¿ Qué decisión
debe tomarse respecto al lote ?
Solución: a) n= 200 c=1 b) α= 0,09 β= 0,1969 c) n 1 = 240 c 1 = 1 n 2 = 480 c 2 = 3
e) n=30 M=0,879 f) Rechazar el lote
6º) Un lote de tamaño 5000 va a ser examinado para la aceptación con el plan
simple n= 25, c= 1 e inspección rectificadora.
Escriba la ecuación que permite determinar el AOQL para este plan, y explique
el procedimiento para encontrarlo
Solución: (1-p)2 + 26p (1-p) – 600p2 = 0 Al despejar p= AOQL