Seis Sigma BB Analisis
Seis Sigma BB Analisis
Seis Sigma BB Analisis
1
Fase de Análisis
Propósitos:
Establecer hipótesis sobre las posibles Causas Raíz
Refinar, rechazar, o confirmar la Causa Raíz
Seleccionar las Causas Raíz más importantes:
Las pocas Xs vitales
Salidas:
Causas raíz validadas
Factores de variabilidad identificados
2
QFD
FASE DE ANÁLISIS
Diagrama de
relaciones
Diagrama
Causa Efecto Diagrama de
Ishikawa
Diagrama
de Árbol
Definición
Y=X1 + X2+. .Xn
CTQs = Ys Medición Y,
Operatividad X1, X2, Xn
X's
Causas
Análisis del Modo y Efecto de potenciales
Falla (AMEF)
3
VI. Análisis
A. Medición y modelaje de relación entre
variables
B: Pruebas de hipótesis
4
A. Medición y modelaje de
relación entre variables
5
A. Medición y modelaje de relación
entre variables
1. Coeficiente de correlación
2. Regresión
3. Herramientas Multivariadas
4. Estudios Multivari
6
VI.A.1 Coeficiente de correlación
7
Definiciones
Correlación
Establece si existe una relación entre las variables y
responde a la pregunta, ”¿Qué tan evidente es esta
relación?"
Regresión
Describe con más detalle la relación entre las variables.
8
Correlación
Propósito:
Propósito: Estudiar
Estudiarla
laposible
posiblerelación
relación
entre
entredos
dosvariables.
variables.
Accidentes laborales
• • •
•
•
•
• Correlación
• •
• •
•
• • positiva,
•
•
•
•
•
•
posible
• • •
• • •
•
• •
9
Coeficiente de correlación (r )
Mide la fuerza de la relación lineal entre las variables
X y Y en una muestra.
10
Correlación de la información (R ) de las X y las Y
Correlación Positiva Correlación Negativa
25
Evidente 25
Evidente
20 20
15 15
10
Y
Y
10
5
5
0
0 5 10 15 20 25 Sin Correlación 0
0 5 10 15 20 25
X 25
R=1 20
X
R=-1
15
Correlación Y 10
Correlación
5
25
Positiva 0 Negativa
0 5 10 15 20 25 25
20
15
X
R=0 20
15
Y
10
Y
10
5
5
0
0 5 10 15 20 25 0
0 5 10 15 20 25
X
R=>1 X
R=>-1
11
Coeficiente de correlación
El coeficiente de correlación r asume el mismo signo
de la pendiente de la recta 1 siendo cero cuando
1 =0
12
Coeficiente de correlación
Reglas empíricas
13
Correlaciones (Pearson)
Tabla de Correlación mínima
(n-1)
R2(Adj) = 1 – (1 – r2) (n-p)
18
Correlación vs causación
Tener cuidado de no tener variables colineales, por
ejemplo peso de un coche y peso de las personas
que transporta, o que no la relación no tenga
sentido, como si lavo mi coche, llueve.
19
VI.A.2 Regresión
20
Análisis de Regresión
El
Elanálisis
análisisde
deregresión
regresión es esun
unmétodo
método
estandarizado
estandarizadopara
paralocalizar
localizarla
lacorrelación
correlaciónentre
entredos
dos
grupos
grupos de
dedatos,
datos, y,
y, quizá
quizámásmásimportante,
importante,crear
crear un
un
modelo
modelode depredicción.
predicción.
Puede
Puedeser
ser usado
usadopara
paraanalizar
analizar las
lasrelaciones
relacionesentre:
entre:
••Una
Unasola
sola“X”
“X”predictora
predictorayyuna
unasola
sola “Y”
“Y”
••Múltiples
Múltiplespredictores
predictores“X”
“X”yy una
unasola
sola“Y”
“Y”
••Varios
Variospredictores
predictores“X”
“X”entre
entresí
sí
21
Supuestos de la regresión lineal
Los principales supuestos que se hacen en el análisis de regresión
lineal son los siguientes:
La relación entre las variables Y y X es lineal, o al menos bien
aproximada por una línea recta.
y 0 1 X
El término de error tiene media cero.
y 0 1 X
23
Regresión Lineal Simple
La línea de regresión se calcula por el método de mínimos cuadrados.
Un residuo es la diferencia entre un punto de referencia en particular (xi,
yi) y el modelo de predicción ( y = a + bx ). El modelo se define de tal
manera que la suma de los cuadrados de los residuales es un mínimo. La
suma residual de los cuadrados es llamada con frecuencia la suma de los
cuadrados de los errores (SSE) acerca de la línea de regresión
yi y = b0 + b1x
ei
• • •
• •
• • a y b son
• • •
• • • • Estimados de
•
•
•
•
• 0 y 1
•
• • •
• • •
•
• •
xi
SSE = ei2 = yi - yi2
Gráfica de la Línea de Ajuste
Recta de regresión
Y=-.600.858+5738.89X
R2 = .895
600
Retención
500 Regresión
95% Intervalo
de confianza
95% Intervalo
400 de predicción
26
Interpretación de los Resultados
28
Errores residuales
^
Los residuos ei Yi Y i , i 1,2,3..., n pueden ser graficados para:
Checar normalidad.
Checar el efecto del tiempo si su orden es conocido en los
datos.
Checar la constancia de la varianza y la posible necesidad de
transformar los datos en Y.
Checar la curvatura de más alto orden que ajusta en las X’s.
R es idu a l
20
0 10
0 X=0.000
-10 -10
-20
-20 -30
-40 -3.0SL=-43.26
-50
-2 -1 0 1 2 0 5 10
Marcador Normal Número de Observación
Histograma - Histograma de Residuales Residuales vs. Ajustes
¿curva de 3 20
Fre cu e n cia
10
campana?
R es idu a l
2
0
¿Aleatorio
1 -10
Ignórese 0
-20
para grupos
-25 -20 -15 -10 -5 0 5 10 15 450 500
Ajuste
550 alrededor de
pequeños de cero, sin
Buscar
Buscarlas
lasinconsistencias
inconsistencias tendencias?
información
mayores
mayores
(<30)
30
Ejemplo
Considere el problema de predecir las ventas mensuales en
función del costo de publicidad. Calcular el coeficiente de
correlación, el de determinación y la recta.
1 1.2 101
2 0.8 92
3 1.0 110
4 1.3 120
5 0.7 90
6 0.8 82
7 1.0 93
8 0.6 75
9 0.9 91
10 1.1 105 31
Cálculo manual
Calcular columnas para Suma X, Suma Y, Xi2, XiYi y Yi2
Xi Yi
MES Publicidad Ventas Xi2 XiYi Yi2
1 1.2 101 1.44 121.2 10201
2 0.8 92 0.64 73.6 8464
3 1.0 110 1.00 110.0 12100
4 1.3 120 1.69 156 14400
5 0.7 90 0.49 63.0 8100
6 0.8 82 0.64 65.6 6724
7 1.0 93 1.00 93.0 8649
8 0.6 75 0.36 45.0 5625
9 0.9 91 0.81 81.9 8281
10 1.1 105 1.21 115.5 11025
y 2
b0 y b1 x
SS y y
2
x y
SS xy xy
n 33
Análisis de varianza en la regresión
Tabla de Análisis de Varianza .
Fuente df SS MS = SS/df Fc
Regresión 1 SSR b1 S XY MS REG MSreg/s2
n i 1 n
i 1
ei Yi Y i Yi Y i Yi Y (Y i Y ) (Y i Y ) 2 (Y i Y ) 2 (Yi Y i ) 2
34
1
t ( n 2,1 ).S
b1 2
__
i )2
( X X
1
MSE S t ( n 2,1 ). S
se(b1 ) b1 2
S XX S XX __
( X i X )2 35
1
t ( n 2,1 ).S
b1 2
__
i )2
( X X
( n 2) MSE ( n 2) MSE
2
2 / 2,n 2 12 / 2,n 2
Intervalos de confianza para la Y estimada promedio
__
1 (X0 X )
2
^
Y0 t a / 2 , n 2 MSE
n S XX
Intervalo de predicción para un valor particular de Y estimado
__
__
1 ( X X )2 1 ( X X )2
Yˆ0 t / 2 ,n 2 MSE 1 0
Y0 Yˆ0 t / 2,n 2 MSE 1 0
n S XX n S XX
36
1
t ( n 2,1 ).S
b1 2
__
i )2
( X X
b1
t0
MSE
S XX
37
1
t ( n 2,1 ).S
b1 2
__
i )2
( X X
S XY
r
S XX SYY
Coeficiente de determinación: r2
R2 mide la proporción de la variación total respecto a la media que
es explicada por la regresión. Se expresa en porcentaje.
^ __
R2
( SS .de.la.regresión . por.b0 )
(Y Y ) 2
1
SSE
__
( SSTotal .corregido. para.la.media ) SYY
(Y i Y )2
38
1
t ( n 2,1 ).S
b1 2
__
i )2
( X X
t 0 t / 2,n 2
Si se rechaza la hipótesis Ho, indicando que existe una
correlación significativa
39
Riesgos de la regresión
Los modelos de regresión son válidos como ecuaciones de
interpolación sobre el rango de las variables utilizadas en el
modelo. No pueden ser válidas para extrapolación fuera de
este rango.
* *
* * * Sin A y B
* * * *
*B
X 40
Riesgos de la regresión
Los outliers u observaciones aberrantes pueden distorsionar
seriamente el ajuste de mínimos cuadrados.
Y
*A *
* * *
**
* * *
** *
**
* * *
**
* *
S = 6.84
43
Ejemplo (cont..)
Inferencias respecto a la pendiente de la línea b1:
45
Ejemplo (cont..)
46
Análisis de Regresión
1. Teclear los datos para Xi y Yi
Column 1 Column 2
Column 1 1 0.875442
Column 2 0.875442 1
NOTAS:
ANOVA
df SS MS F Significance F
Regression 1 1226.927 1226.927 26.24633 0.000904
Residual 8 373.973 46.74662
Total 9 1600.9
Confidence 95%
Coefficients Standard Error t Stat P-value Lower Upper
140
120
100
80
Y
60
40
20
0
0 20 40 60 80 100
Sample Percentile
50
Gráfica de Residuos vs. X de Excel
20
Residuals
10
0
-10 0 0.5 1 1.5
X Variable 1
51
Ejercicio
Calcular la recta de predicción con sus bandas de confianza, la
correlación y la determinación para la respuesta de un Taxi,
los datos se muestran a continuación:
Distancia Tiempo
0.8 200
2.2 400
1.0 160
0.6 120
1.0 360
1.4 280
2.2 560
0.6 320
52
Relaciones no Lineales
¿Qué pasa si existe una relación causal, no lineal?
¿Cómo describiría
El siguiente es un conjunto de datos
experimentales codificados, sobre
esta relación?
resistencia a la compresión de una
aleación especial:
Resistencia a
Concentración la Compresión
x y
10.0 25.2 27.3 28.7
15.0 29.8 31.1 27.8
20.0 31.2 32.6 29.7
25.0 31.7 30.1 32.3
30.0 29.4 30.8 32.8
53
Resultados del Análisis de Regresión - Modelo Cuadrático
Analysis of Variance
Source DF SS MS F P
Regression 1 8.9296 8.9296 160.26 0.000
Residual Error 23 1.2816 0.0557
Total 24 10.2112
56
Regresión cuadrática
Los residuos
No son normales
Se deben transformar
Las variables
57
Otros Patrones No Lineales
A veces es posible transformar una o ambas variables, para mostrar
mejor la relación entre ambas. La meta es identificar la relación
matemática entre las variables, para que con la variable transformada
se obtenga una línea más recta. Algunas transformaciones comunes
incluyen:
x’ = 1/x
x’ = Raíz cuadrada de (x)
59
Transformación de variables del
ejemplo de regresión cuadrática
Transformando la variable X’ = 1/X se tiene, utilizando Minitab
The regression equation is
Y = 2.98 - 6.93 1/X
60
Transformación de variables del
ejemplo de regresión cuadrática
Transformando la variable X’ = 1/X se tiene, utilizando Minitab
61
Transformación de variables del
ejemplo de regresión cuadrática
Los residuos ahora ya se muestran normales
62
Transformación para homoestacidad
de la varianza
Algunas transformaciones para estabilizar la varianza
Relación de 2 a E(Y) Transformación
2 constante ..............................Y ' Y
2 E (Y ).................................Y ' Y Datos de Poisson
2 E (Y ) 3 ...........................Y ' Y 1 / 2
63
Transformación para homoestacidad
de la varianza
Ejemplo: Se hizo un estudio entre la demanda (Y) y la energía
eléctrica utilizada (X) durante un cierto periodo de tiempo
Obs X Y Fit SE Fit Residual St Resid
1 679 0.790 1.649 0.351 -0.859 -0.61
2 292 0.440 0.308 0.490 0.132 0.10
3 1012 0.560 2.802 0.293 -2.242 -1.57
4 493 0.790 1.004 0.412 -0.214 -0.15
5 582 2.700 1.312 0.381 1.388 0.98
6 1156 3.640 3.301 0.297 0.339 0.24
7 997 4.730 2.750 0.294 1.980 1.38
8 2189 9.500 6.880 0.651 2.620 2.00R
9 1097 5.340 3.097 0.293 2.243 1.57
10 2078 6.850 6.495 0.600 0.355 0.27
11 1818 5.840 5.595 0.488 0.245 0.18
12 1700 5.210 5.186 0.441 0.024 0.02
13 747 3.250 1.884 0.333 1.366 0.96
14 2030 4.430 6.329 0.579 -1.899 -1.42
15 1643 3.160 4.988 0.420 -1.828 -1.31
16 414 0.500 0.730 0.441 -0.230 -0.17
17 354 0.170 0.523 0.465 -0.353 -0.25
18 1276 1.880 3.717 0.313 -1.837 -1.29
19 745 0.770 1.877 0.333 -1.107 -0.78
20 435 1.390 0.803 0.433 0.587 0.42
21 540 0.560 1.167 0.395 -0.607 -0.43
22 874 1.560 2.324 0.307 -0.764 -0.53
23 1543 5.280 4.642 0.384 0.638 0.45
24
25
1029
710
0.640
4.000
2.861
1.756
0.293
0.343
-2.221
2.244
-1.55
1.58 64
Transformación para homoestacidad
de la varianza
Ejemplo: Se hizo un estudio entre la demanda (Y) y la energía
eléctrica utilizada (X) durante un cierto periodo de tiempo
The regression equation is
Y = - 0.704 + 0.00346 X
Predictor Coef SE Coef T P
Constant -0.7038 0.6170 -1.14 0.266
X 0.0034645 0.0005139 6.74 0.000
S = 1.462 R-Sq = 66.4% R-Sq(adj) = 64.9%
Analysis of Variance
Source DF SS MS F P
Regression 1 97.094 97.094 45.45 0.000
Residual Error 23 49.136 2.136
Total 24 146.231
65
Transformación para homoestacidad
de la varianza
Se observa que la varianza se incrementa conforme aumenta X
66
Transformación para homoestacidad
de la varianza
Se observa que la varianza se incrementa conforme aumenta X
67
Transformación para homoestacidad
de la varianza
Transformando a X por su raíz cuadrada se tiene:
Obs X SQR-Y Fit SE Fit Residual St Resid
1 679 0.8888 1.1694 0.1092 -0.2805 -0.64
2 292 0.6633 0.7717 0.1524 -0.1084 -0.25
3 1012 0.7483 1.5115 0.0912 -0.7632 -1.71
4 493 0.8888 0.9783 0.1280 -0.0894 -0.21
5 582 1.6432 1.0697 0.1184 0.5735 1.31
6 1156 1.9079 1.6595 0.0922 0.2484 0.56
7 997 2.1749 1.4961 0.0914 0.6788 1.52
8 2189 3.0822 2.7208 0.2024 0.3614 0.89
9 1097 2.3108 1.5989 0.0911 0.7120 1.60
10 2078 2.6173 2.6068 0.1867 0.0105 0.03
11 1818 2.4166 2.3397 0.1518 0.0770 0.18
12 1700 2.2825 2.2184 0.1371 0.0641 0.15
13 747 1.8028 1.2392 0.1035 0.5635 1.27
14 2030 2.1048 2.5575 0.1800 -0.4527 -1.09
15 1643 1.7776 2.1598 0.1304 -0.3822 -0.88
16 414 0.7071 0.8971 0.1372 -0.1900 -0.44
17 354 0.4123 0.8354 0.1445 -0.4231 -0.98
18 1276 1.3711 1.7828 0.0974 -0.4116 -0.93
19 745 0.8775 1.2372 0.1037 -0.3597 -0.81
20 435 1.1790 0.9187 0.1347 0.2603 0.60
21 540 0.7483 1.0265 0.1228 -0.2782 -0.64
22 874 1.2490 1.3697 0.0955 -0.1207 -0.27
23 1543 2.2978 2.0571 0.1195 0.2407 0.55
68
24 1029 0.8000 1.5290 0.0910 -0.7290 -1.64
25 710 2.0000 1.2012 0.1065 0.7988 1.81
Transformación para homoestacidad
de la varianza
Transformando a X por su raíz cuadrada se tiene:
Regression Analysis: SQR-Y versus X
69
Transformación para homoestacidad
de la varianza
Transformando a X por su raíz cuadrada se tiene:
70
Regresión lineal múltiple
71
Regresión múltiple
Cuando se usa más de una variable independiente para predecir
los valores de una variable dependiente, el proceso se llama
análisis de regresión múltiple, incluye el uso de ecuaciones
lineales.
Yu 0 1 X u1 2 X u 2 ....... k X uk u
72
Regresión múltiple
Estimación de los parámetros del modelo
Se trata de minimizar los errores cuadráticos en:
N
R ( 0 , 1 ,..., k ) (Yu 0 1 X u1 2 X u 2 ..... uk ) 2
u 1
73
Regresión múltiple
Estimación de los parámetros del modelo:
b = (X’X)-1 X’Y
El vector de valores ajustados Yˆ Xb se puede expresar como:
Yˆ Xb X ( X ' X ) 1 X ' Y Hy
SSE (Y Xb )' (Y Xb) Y ' Y b' X ' Y Y ' Xb b' X ' Xb Y ' Y 2b' X ' Y b' X ' Xb
SSE
SSE Y ' Y b' X ' Y s 2 MSE
Np 74
Tamaño de muestra
Tomar 5 observaciones para cada una de las
variables independientes, si esta razón es menor de5
a 1, se tiene el riesgo de “sobreajustar” el modelo
75
Ejemplo de regresión múltiple
Un embotellador está analizando las rutas de servicio de
máquinas dispensadoras, está interesado en predecir la
cantidad de tiempo requerida por el chofer para surtir las
máquinas en el local (Y).
76
X2-Dist X1-CAS Y-TENT Fit SE Fit Residual St Resid
Obs
16.68 7.0 16.680 21.708 1.040 -5.028 -1.63
1
11.50 3.0 11.500 10.354 0.867 1.146 0.36
2
12.03 3.0 12.030 12.080 1.024 -0.050 -0.02
3
1 1 1 1 1 1 1 1 1 1 1 1
6 9 10 6 7 3 17 10 26 9 8 4
462 448 776 200 132 36 770 140 810 450 635 150 ]
Matrix M6 = X'Y
[ 25 219 10232
219 3055 133899
10232 133899 6725688 ]
Matrix M7 = X'Y
[ 560
7375
337072 ]
78
Ejemplo de regresión múltiple
Solución matricial
Matrix M8 = INV(X'X)
2.34123
1.61591
0.01438
79
Ejemplo de regresión múltiple
Solución matricial
Cálculo de la estimación de la varianza:
Data Display
[ 16.68 11.50 12.03 14.88 13.75 18.11 8.00 17.83 79.24 21.50
40.33
21.00 13.50 19.75 24.00 29.00 15.35 19.00 9.50 35.10 17.90 52.32
18.75 19.83 10.75 ]
SS E 233.732
S2 10.624
Np 25 3
80
Ejemplo de regresión múltiple
Solución matricial
Intervalo de confianza para Beta 1
81
Ejemplo de regresión múltiple
Solución matricial
El embotellador desea construir un intervalo de confianza sobre
el tiempo medio de entrega para un local requiriendo:
1 2.34123
X 0 8 Yˆ0 X ' 0 b 1,8,275 1.61591 19.22minutos
275 0.01438
82
Ejemplo de regresión múltiple
Solución matricial
El intervalo de confianza sobre el tiempo medio de entrega para
un local requiriendo es para 95% de nivel de confianza:
83
Ejemplo de regresión múltiple
Solución matricial
El análisis de varianza es:
Analysis of Variance
(559.6) 2
SST = 18,310.629 - = 5784.5426
25
(559.6) 2
SSR = 18,076.930 - = 5,550.8166
25
SSE = SST – SSR = 233.7260
85
Multicolinealidad
La multicolinealidad implica una dependencia cercana entre
regresores (columnas de la matriz X ), de tal forma que si hay
una dependencia lineal exacta hará que la matriz X’X sea
singular.
86
Multicolinealidad
Una prueba fácil de probar si hay multicolinealidad entre dos
variables es que su coeficiente de correlación sea mayor a 0.7
1
VIF j
1 R 2j
Si es mayor a 10 implica que se tienen serios problemas de
multicolinealidad.
87
Análisis de los residuos
Los residuos graficados vs la Y estimada, pueden mostrar
diferentes patrones indicando adecuación o no adecuación del
modelo:
Yˆ 5 X 1 1000 X 2
89
Escalamiento de residuos
Residuos estandarizados
Se obtienen dividiendo cada residuo entre la desviación
estándar de los residuos
ei
di ,
MSE
Después de la estandarización, los residuos tienen una
media de 0 y desviación estándar de 1
94
Matriz sombrero
Los valores en la diagonal de la matriz sombrero miden dos
aspectos:
Para cada observación miden la distancia de la observación
95
Distancia de Mahalanobis
D2 es una medida comparable a los valores sombrero (hat
values) que considera sólo la distancia de una observación del
valor medio de las variables independientes.
96
Influencia en coeficientes
individuales
El impacto de eliminar una observación simple en cada
uno de los coeficientes de la regresión múltiple se muestra
con la DFBETA y su versión estandarizada SDFBETA.
99
Ejemplo de Regresión Múltiple
Cat. (US News) GMAT Salario Inicial ($) % Aceptación
Stanford 1 711 82000 7.4
Harvard 2 670 80000 12.8
Penn (Wharton) 3 662 79000 14.7
MIT (Sloan) 4 650 78000 15.1
Chicago 5 680 65000 25.0
Northwestern 6 660 70000 16.0
Columbia 7 660 83000 14.8
Dartmouth 8 670 70000 12.6
Duke 9 646 67500 20.5
Berkeley 10 653 70000 13.3
Virginia 11 660 66000 18.9
Michigan 12 645 65000 28.0
NYU 13 646 70583 20.9
Carnegie Mellon 14 640 67200 30.8
Yale 15 675 65000 23.5
U.N.C. 16 630 60000 19.8
UCLA 17 651 65000 17.5
Texas-Austin 18 630 60000 27.3
Indiana 19 630 61500 44.7
Cornell 20 637 64000 25.4
Rochester 21 630 58500 36.0
Ohio State 22 611 61000 23.2
Emory 23 626 60000 33.0
Purdue 24 603 63700 20.7
Maryland 25 640 53000 18.9
100
Interpretación de Resultados de Excel- Regresión Multiple
SUMMARY OUTPUT
Regression Statistics
Multiple R 0.8749313 R Square 0.76550478
Adjusted R Square 0.732005463 Standard Error 4050.855918 Observations 25
ANOVA
df SS MS F Significance F
Regression 3 1.12E+09 374977790.1 22.851355 8.17E-07
Residual 21 3.45E+08 16409433.67
Total 24 1.47E+09
Regression Statistics
Multiple R 0.855974
R Square 0.732691
Adjusted R Square 0.721069
Standard Error 4132.688
Observations 25
ANOVA
df SS MS F Significance F
Con
Consólo
sóloX1,
X1,el
elModelo
Modelosesesimplifica
simplificaenormemente
enormemente
poca
poca importancia
importanciapráctica
prácticase
sepierde
pierdeen
enRR2 (ajustada)
2
(ajustada)
Reducción
Vuelva a del Modelo
correr la regresión usando la categoría
US News, como el único agente de predicción (“predictor”)
Análisis de Variancia
Fuente DF SS MS F p
Regresión 1 1076712008 1076712008 63.04 0.000
Error 23 392819470 17079107
Total 24 1469531477
El
El Modelo
Modelo se
se simplifica
simplifica enormemente..…poca
enormemente..…poca
importancia
importancia práctica se pierde en R (ajustada)
práctica se pierde en R 22
(ajustada)
Corrida en Minitab
Se introducen los datos en varias columnas C1 a C5
incluyendo la respuesta Y (heatflux) y las variables
predictoras X’s (North, South, East)
HeatFlux Insolation East South North
271.8 783.35 33.53 40.55 16.66
264.0 748.45 36.50 36.19 16.46
238.8 684.45 34.66 37.31 17.66
230.7 827.80 33.13 32.52 17.50
251.6 860.45 35.75 33.71 16.40
257.9 875.15 34.46 34.14 16.28
104
Corrida en Minitab
Utilzar el archivo de ejemplo Exh_regr.mtw
Opción: Stat > Regression > Regression
Para regresión lineal indicar la columna de respuesta
Y (Score2) y X (Score1)
3.5
Score2
2.5
Regression
1.5 95% CI
95% PI
2 3 4 5 6 7 8 9
Score1
107
Resultados de la regresión Múltiple
The regression equation is
HeatFlux = 389 - 24.1 North + 5.32 South + 2.12 East
Predictor Coef SE Coef T P
Constant 389.17 66.09 5.89 0.000
North -24.132 1.869 -12.92 0.000
South 5.3185 0.9629 5.52 0.000
East 2.125 1.214 1.75 0.092
S = 8.598 R-Sq = 87.4% R-Sq(adj) = 85.9%
Analysis of Variance
Source DF SS MS F P
Regression 3 12833.9 4278.0 57.87 0.000
Residual Error 25 1848.1 73.9
Total 28 14681.9
Source DF Seq SS
North 1 10578.7
South 1 2028.9
108
East 1 226.3
Resumen de la Regresión
• La regresión sólo puede utilizarse con información de variables
continuas.
110
Herramientas multivariadas
1. Introducción
3. Análisis factorial
4. Análisis discriminante
5. MANOVA
111
Introducción
En el análisis multivariado se incluyen dos o más
variables dependientes Y1, Y2, etc. Consideradas
simultáneamente para las variables independientes
X1, X2, …., Xn
113
Análisis de componentes principales
Pasos de análisis en Minitab
Se usa una matriz de correlación para determinar la
relación entre componentes
Las matrices definen cantidades como eigenvalores y
eigenvectores
Se suman los eigenvalores y se calculan las
proporciones de cada componente
Se identifican los PC1, PC2, … que explican la mayor
parte de la varianza
Se puede hacer un diagrama de Pareto como apoyo
114
Ejemplo: Alimentos en Europa
X1 X2 X3 X4 X5 X6 X7 X8 X9
País RMEAT WMEAT EGGS MILK FISH CERL STARCH NUTS FR-VEG
1 10.1 1.4 0.5 8.9 0.2 42.3 0.6 5.5 1.7
2 8.9 14 4.3 19.9 2.1 28 3.6 1.3 4.3
3 13.5 9.3 4.1 17.5 4.5 26.6 5.7 2.1 4
4 7.8 6 1.6 8.3 1.2 56.7 1.1 3.7 4.2
5 9.7 11.4 2.8 12.5 2 34.3 5 1.1 4
6 10.6 10.8 3.7 25 9.9 21.9 4.8 0.7 2.4
7 8.4 11.6 3.7 11.1 5.4 24.6 6.5 0.8 3.6
8 9.5 4.9 2.7 33.7 5.8 26.3 5.1 1 1.4
9 18 9.9 3.3 19.5 5.7 28.1 4.8 2.4 6.5
10 10.2 3 2.8 17.6 5.9 41.7 2.2 7.8 6.5
11 5.3 12.4 2.9 9.7 0.3 40.1 4 5.4 4.2
12 13.9 10 4.7 25.8 2.2 24 6.2 1.6 2.9
13 9 5.1 2.9 13.7 3.4 36.8 2.1 4.3 6.7
14 9.5 13.6 3.6 23.4 2.5 22.4 4.2 1.8 3.7
15 9.4 4.7 2.7 23.3 9.7 23 4.6 1.6 2.7
16 6.9 10.2 2.7 19.3 3 36.1 5.9 2 6.6
17 6.2 3.7 1.1 4.9 14.2 27 5.9 4.7 7.9
18 6.2 6.3 1.5 11.1 1 49.6 3.1 5.3 2.8
19 7.1 3.4 3.1 8.6 7 29.2 5.7 5.9 7.2
20 9.9 7.8 3.5 24.7 7.5 19.5 3.7 1.4 2
21 13.1 10.1 3.1 23.8 2.3 25.6 2.8 2.4 4.9
22 17.4 5.7 4.7 20.6 4.3 24.3 4.7 3.4 3.3
23 9.3 4.6 2.1 16.6 3 43.6 6.4 3.4 2.9
24 11.4 12.5 4.1 18.8 3.4 18.6 5.2 1.5 3.8
25 4.4 5 1.2 9.5 0.6 55.9 3 5.7 3.2
115
Corrida en Minitab
2 Stat > Multivariate > Principal components
3 En Variables, X1, X2, X3, X4, X6, X7, X8, X9
116
Ejemplo: Alimentos en Europa
Scree Plot of RMEAT, ..., FR-VEG
3
Loading Plot of RMEAT, ..., FR-VEG
Eigenvalue
WMEA T C ERL
MILK
2 0.2
0.1 RMEA T
EGGS
0.0
Second Component
1
-0.1 NUTS
-0.2
0
-0.3 STA RC H
1 2 3 4 5 6 7 8 9
Component Number -0.4
-0.5 FR-VEG
-0.6 FISH
117
Ejemplo: Alimentos en Europa
Se tiene la gráfica siguiente de países:
Europa occidental Europa oriental Balcanes
Scatterplot of Z2 vs Z1
2 1
4
18
14 2 11
25
1 12
8
21
5
24 20
22 23
0 6 3
7 13
16
9 15
10
-1
Z2
-2
19
-3
-4 17
-5
-3 -2 -1 0 1 2 3 4
Z1
Península ibérica
118
Análisis factorial
Es una técnica de reducción de variables para
identificar factores que expliquen la variación,
aunque se reiere un juicio subjetivo.
119
Análisis factorial
Se especifican un cierto número de factores comunes
120
Corrida con Minitab
2 Stat > Multivariate > Factor Analysis.
3 En Variables, X1, X2, X3, X4, X6, X7, X8, X9
4 En Number of factors to extract, 4.
En Method of Extraction, seleccionar Principal components
6 En Type of Rotation, seleccionar Varimax.
7 Click Graphs y seleccionar Loading plot for first 2 factors
y Scree Plot.
Click Results y seleccionar Sort loadings.
Seleccionar Storage e indicar columnas para ponderaciones,
coeficientes, Z’s, eigenvalores, etc.
Click OK en cada uno de los cuadros de d
121
Loading Plot of RMEAT, ..., FR-VEG
CERL
0.50
NUTS
Ejemplo
0.25
FR-VEG STA RC H
Second Factor
0.00 FISH
WMEA T
-0.25
-0.50
MILK
EGGS
-0.75
RMEA T
-1.00
-0.50 -0.25 0.00 0.25 0.50 0.75 1.00
First Factor
Scatterplot of Z2 vs Z1
2
Yugoslav ia
Portugal
Rumania Hungría
Bélgica
Grecia
-1 Suiza Irlanda
Francia
-2 Reino Unido
-2 -1 0 1 2
Z1
123
Análisis discriminante
Si se tiene una muestra con grupos conocidos, el
análisis discriminante clasifica las observaciones o
atributos en dos o más grupos
124
Ejemplo de actividades en países
No Grupo Ciudad Agr Min Man Ps Con Ser Fin Sps Tc
1 1 Bélgica 3.3 0.9 27.6 0.9 8.2 19.1 6.2 26.6 7.2
2 1 Dinamarca 9.2 0.1 21.8 0.6 8.3 14.6 6.5 32.2 7.1
3 1 Francia 10.8 0.8 27.5 0.9 8.9 16.8 6.0 22.6 5.7
4 1 Alemania Occ. 6.7 1.3 35.8 0.9 7.3 14.4 5.0 22.3 6.1
5 1 Irlanda 23.2 1.0 20.7 1.3 7.5 16.8 2.8 20.8 6.1
6 1 Italia 15.9 0.6 27.6 0.5 10.0 18.1 1.6 20.1 5.7
7 1 Luxenburgo 7.7 3.1 30.8 0.8 9.2 18.5 4.6 19.2 6.2
8 1 Holanda 6.3 0.1 22.5 1.0 9.9 18.0 6.8 28.5 6.8
9 1 Inglaterra 2.7 1.4 30.2 1.4 6.9 16.9 5.7 28.3 6.4
10 1 Austria 12.7 1.1 30.2 1.4 9.0 16.8 4.9 16.8 7.0
11 1 Finlandia 13.0 0.4 25.9 1.3 7.4 14.7 5.5 24.3 7.6
12 2 Grecia 41.4 0.6 17.6 0.6 8.1 11.5 2.4 11.0 6.7
13 1 Noruega 9.0 0.5 22.4 0.8 8.6 16.9 4.7 27.6 9.4
14 2 Portugal 27.8 0.3 24.5 0.6 8.4 13.3 2.7 16.7 5.7
15 2 España 22.9 0.8 28.5 0.7 11.5 9.7 8.5 11.8 5.5
16 1 Suecia 6.1 0.4 25.9 0.8 7.2 14.4 6.0 32.4 6.8
17 1 Suiza 7.7 0.2 37.8 0.8 9.5 17.5 5.3 15.4 5.7
18 2 Turquía 66.8 0.7 7.9 0.1 2.8 5.2 1.1 11.9 3.2
19 3 Bulgaria 23.6 1.9 32.3 0.6 7.9 8.0 0.7 18.2 6.7
20 3 Checa 16.5 2.9 35.5 1.2 8.7 9.2 0.9 17.9 7.0
21 3 Alemania Ori. 4.2 2.9 41.2 1.3 7.6 11.2 1.2 22.1 8.4
22 3 Hungría 21.7 3.1 29.6 1.9 8.2 9.4 0.9 17.2 8.0
23 3 Polonia 31.1 2.5 25.7 0.9 8.4 7.5 0.9 16.1 6.9
24 3 Rumania 34.7 2.1 30.1 0.6 8.7 5.9 1.3 11.7 5.0
25 3 Rusia 23.7 1.4 25.8 0.6 9.2 6.1 0.5 23.6 9.3
26 3 Yugoslavia 48.7 1.5 16.8 1.1 4.9 6.4 11.3 5.3 4.0
125
Corrida con Minitab
2 Stat > Multivariate > Discriminant Analysis.
126
Corrida con Minitab
Canonical Discriminant Functions
3
1 3
1
0
-1 GRUPO
F u n c t io n 2
2 Group Centroids
-2
3
-3
2
-4 1
-6 -4 -2 0 2 4 6
Function 1
127
Análisis de conglomerados
128
Análisis de conglomerados
Se usa para determinar agrupaciones o
clasificaciones de un conjunto de datos
129
Ejemplo
Suponer que un estudio de Variables V1 V2
mercado trata de determinar
segmentos de mercado en A 3 2
base a los patrones de B 4 5
lealtad de marcas (V1) y
tiendas (V2), medidas del 0 C 4 7
al 10 en 7 personas (A-G).
D 2 7
E 6 6
F 7 7
G 6 4
130
Corrida en Minitab
Stat > Multivariate Análisis > Cluster Observations
Distance Measured Euclidean Seleccionar Show
Dendogram OK
Dendrogram with Single Linkage and Euclidean Distance
3.16
2.11
Distance
1.05
0.00
1 2 3 4 5 6 7
Observations
131
Análisis de correlación canónico
Prueba la hipótesis de que los efectos pueden tener
causas múltiples y de que las causas pueden tener
efectos múltiples (Hotelling 1935)
132
Análisis de correlación canónico
Se usan los más altos valores de correlación para los
conjuntos. Los pares de combinaciones lineales se
denominan variates canónicas con correlaciones
canónicas (Rc con valor mayor a 0.3)
133
MANOVA
(Análisis de varianza múltiple)
Es un modelo para analizar la relación entre una o
más variables independientes y dos o más variables
dependientes
134
MANOVA
(Análisis de varianza múltiple)
135
Diferencias de ANOVA y MANOVA
136
Ejemplo:
Extrusión de película plástica
Se realiza un estudio para determinar las condiciones
óptimas para extruir película plástica.
137
Ejemplo:
Extrusión de película plástica
Tear Gloss Opacity Extrusión Additive
6.5 9.5 4.4 1 1
6.2 9.9 6.4 1 1
5.8 9.6 3 1 1
6.5 9.6 4.1 1 1
6.5 9.2 0.8 1 1
6.9 9.1 5.7 1 2
7.2 10 2 1 2
6.9 9.9 3.9 1 2
6.1 9.5 1.9 1 2
6.3 9.4 5.7 1 2
6.7 9.1 2.8 2 1
6.6 9.3 4.1 2 1
7.2 8.3 3.8 2 1
7.1 8.4 1.6 2 1
6.8 8.5 3.4 2 1
7.1 9.2 8.4 2 2
7 8.8 5.2 2 2
7.2 9.7 6.9 2 2
7.5 10.1 2.7 2 2
7.6 9.2 1.9 2 2
138
Ejemplo:
Extrusión de película plástica
1 Abrir el archivo EXH_MVAR.MTW.
2 Seleccionar Stat > ANOVA > Balanced
MANOVA.
3 En Responses, poner Tear Gloss Opacity.
4 En Model, poner Extrusion | Additive.
5 Click Results. En Display of Results, seleccionar
Matrices (hypothesis, error, partial
correlations) y Eigen analysis.
6 Click OK en cada cuadro de diálogo.
139
Ejemplo
Criterion Statistic F Num Denom P
Wilks' 0.38186 7.554 3 14 0.003
SSCP Matrix for Extrusion
Tear Gloss Opacity
Tear 1.740 -1.505 0.8555
Gloss -1.505 1.301 -0.7395
Opacity 0.855 -0.739 0.4205
SSCP Matrix for Error
Tear Gloss Opacity
Tear 1.764 0.0200 -3.070
Gloss 0.020 2.6280 -0.552
Opacity -3.070 -0.5520 64.924
Partial Correlations for the Error SSCP Matrix
Eigenvector 1 2 3
Tear 0.6541 0.4315 0.0604
Gloss -0.3385 0.5163 0.0012
Opacity 0.0359 0.0302 -0.1209
140
Ejemplo:
Extrusión de película plástica
Las matrices SSCP evalúan la contribución a la
variabilidad de manera similar a la suma de
cuadrados en la ANOVA univariada.
141
VI.A.4 Estudios Multivari
142
Estudios Multivari
La carta multivari permite analizar la variación dentro
de la pieza, de pieza a pieza o de tiempo en tiempo
143
Estudios Multivari
La variación dentro de las muestras (cinco puntos en
cada línea). La variación de muestra a muestra como
posición vertical de las líneas.
E
S
P
E
S
O
R
145
Estudios Multivari
Procedimiento de muestreo:
Seleccionar el proceso y la característica a
investigar
146
Estudios Multivari
Procedimiento de muestreo:
Realizar la carta Multivari
Unir los valores observados con una línea
Variaciones por la localización dentro de un proceso que
produce múltiples unidades al mismo tiempo. Por
ejemplo las diferentes cavidades de un molde
Variaciones de máquina a máquina; operador a
operador; ó planta a planta
149
Cartas Multivari
Ejemplo: Se toman 3 a 5 unidades consecutivas, repitiendo el
proceso tres o más veces a cierto intervalo de tiempo, hasta que al
menos el 80% de la variación en el proceso se ha capturado.
1 2 3 4 5 27 28 29 30 31 55 56 57 58 59
VARIACIÓN POSICIONAL DENTRO DE LA UNIDAD
150
Cartas Multivari
Ejemplo: (cont...)
1 2 3 4 5 27 28 29 30 31 55 56 57 58 59
VARIACIÓN CÍCLICA DE UNIDAD A UNIDAD
151
Cartas Multivari
Ejemplo: (cont...)
1 2 3 4 5 27 28 29 30 31 55 56 57 58 59
152
Cartas Multivari
Ejemplo: Un proceso produce flecha cilíndricas, con un
diámetro especificado de 0.0250” 0.001”.
153
Cartas Multivari
Se tomaron cuatro lecturas en cada flecha, dos a cada
lado. Estas muestran una disminución gradual desde
el lado izquierdo al lado derecho de las flechas,
además de excentricidad en cada lado de la flecha.
154
Cartas Multivari
8 AM 9 AM 10 AM 11 AM 12 AM
.0.2510”
0.2500”
0.2490”
Izquierda
Máximo
Derecha
Mínimo
155
Cartas Multivari
Un análisis rápido revela que la mayor variación es temporal con un cambio mayor entre las 10 AM y las 11 AM.
A las 10 AM se para el equipo para el almuerzo y se arranca a las 11 AM, con lecturas similares a las de las 8 AM. Conforme pasa el tiempo
las lecturas tienden a decrecer más y más, hasta que se invierten a las 10 A.M. en forma drástica.
La variación en temperatura era causada por que la cantidad de refrigerante no era la adecuada, lo cual se notaba más cuando se paraba el
equipo y se volvía a arrancar. Se adicionó, reduciendo la variación en 50% aproximadamente..
156
Cartas Multivari
También se encontró que el acabado cónico era causado por que la herramienta de corte estaba mal
alineada. Se ajustó, contribuyendo a otra reducción del 10% de la variabilidad.
157
Cartas Multivari
Tipo de % var. Causas de Acción % de variación
Variación Total Variación Correctiva Reducida
Temporal 50 Bajo nivel de Adicionar Casi 50
Tiempo a tiempo Refrigerante refrigerante
Dentro de 10 Ajuste no Ajuste de la Casi 10
la flecha no paralelo herramienta de
corte
Dentro de 30 Rodamiento Nuevo Casi 30
la flecha gastado rodamiento
Flecha a 5 -??? - -
flecha
158
Cartas Multivari
Resultados:
La variación total en la siguiente corrida de producción se redujo de 0.0025” a 0.0004”
Como beneficios se redujo a cero el desperdicio y no hubo necesidad de adquirir una nueva máquina.
Se observa que antes de cambiar equipo o máquinas, es conveniente realizar un estudio de variabilidad para identificar las fuentes de
variación y tratar de eliminarlas.
159
Cartas Multivari
Ejemplo: Búsqueda de fuentes de variación con el diagrama sistemático.
Diámetro de Flecha
(0.150" +/- .002)
Variación
Variación de
de
sist. medición
proceso
Operador a
operador
Programa Máquina Accesorios
160
Cartas Multivari
Ejemplo (cont..):
• Al realizar la prueba de homogeneidad de varianza F, se
encontró que había una diferencia significante entre los
operadores.
163
Resultados
Multi-Vari Chart for Strength by SinterTime - MetalType
SinterTime
0.5
23.5
1.0
2.0
22.5
21.5
Strength
20.5
19.5
18.5
17.5
15 18 21
MetalType
164
VI.A.5 Análisis de datos
por atributos
165
Análisis de datos por atributos
Si los CTQ’s son variables continuas, se usa la
regresión, dependiendo de la naturaleza de la
característica crítica para el cliente (CTS’s) como éste
la expresa:
CTS HERRAMIENTA
Nominal (Verde, Rojo, azul) Regresión Logística Nominal
Atributo (Pasa/No pasa) Regresión Logística Binaria
Ordinal (1, 2, 3, 4, 5) Regresión Logística Ordinal
166
Análisis de datos por atributos
El análisis de datos por atributos se organiza en
valores, categorías o grupos dicotómicos
167
Análisis de datos por atributos
Regresión logística
Relaciona variables independientes categóricas a una
variable dependiente (Y). Minitab incluye los modelos
binario, ordinal y nominal
Regresión logit
Es subconjunto del modelo log-lineal. Tiene solo una
variable dependiente, usa determinaciones de
probabilidad o tasa de probabilidad
168
Análisis de datos por atributos
Regresión probit
Es similar a la prueba de vida acelerada, la unidad se
somete a esfuerzo con la respuesta pasa/falla, bueno o
malo. Es una respuesta binaria en un tiempo de falla
futuro
169
Regresión logística o binaria
• En caso de información cualitativa es necesario
traducir las preferencias del cliente expresadas como
atributos a un intervalo de valores aceptables de
variables (Especificaciones).
170
Regresión logística o binaria
Es similar a la regresión múltiple excepto que la
respuesta es binaria (si/no, bueno/malo, etc.) Sus
coeficientes se determinan por el método de máxima
verosimilitud
Yi = 0, 1
171
Regresión logística o binaria
La probabilidad de que el resultado esté en cierta
categoría es:
173
Regresión logística
Relación con ajuste pobre
174
Regresión logística - Procedimiento
Definir el atributo a “traducir” (“y”)
Definir la variable apropiada para el atributo (“x”)
Definir el modelo matemático a probar
Determinar los defectos que está dispuesto a
aceptar
Recolecte información de “x” vs “y”. Asigne 1 si falla
y 0 si es aceptable.
Analice la información mediante Regresión Logística
Binaria
175
Regresión logística- Procedimiento
176
Regresión logística - Procedimiento
P-Value de Deviance
177
Regresión logística - Procedimiento
• Construya el modelo de regresión para la
probabilidad de falla estará dado por :
P(Falla) =
e
b +b x +....
0 1 1
Donde :
1+e
b +b x +....
0
b , b , ... = Coeficientes del modelo
1 1
0 1
178
Ejemplo de riesgo de paro cardiaco
Logistic Regression Table
Odds 95% CI
Predictor Coef SE Coef Z P Ratio Lower Upper
Constant -1.98717 1.67930 -1.18 0.237
Fuma
Si -1.19297 0.552980 -2.16 0.031 0.30 0.10 0.90
Peso 0.0250226 0.0122551 2.04 0.041 1.03 1.00 1.05
180
Regresión logística ordinal -
Procedimiento
Defina la variable de respuesta a “traducir” (“y” ó
CTS)
Defina el CTQ (“x”) ó variable a relacionar con el CTS
Defina el modelo matemático a probar
Determine los defectos que está dispuesto a aceptar
en la categoría de interés
Recolecte información de “x” vs “y”
Analice la información mediante Regresión Logística
Ordinal
181
Regresión logística ordinal -
Procedimiento
Stat > Regression > Ordinal Logistic Regression
Seleccione la respuesta (“y”)
Seleccione los términos que estima tiene el modelo
[3]
Constantes y
Coeficientes
del modelo
182
Regresión logística ordinal -
Procedimiento
Observe el P-Value de “Deviance” en la Sesión, debe
de ser grande (P >0.10)
183
Regresión logística ordinal -
Procedimiento
Constantes y
Coeficientes
del modelo
Expectativas
de Diseño
Matriz de
Importan.
CTQs
Diseño
Producto
(CTS’s)
(DPs)
Tipo
(General) Especificaciones Clientes
LIE LSE Otra
Usuarios
Finales
Producto
(Específico)
185
Análisis Logit
Usa razones para determinar que tanta posibilidad
tiene una observación de pernecer a un grupo que a
otro.
Una posibilidad de 0.8 de estar en el grupo A se
puede expresar como una tasa de posibilidades de
4:1 ( que es p/(1-p)), cuyo logaritmo es el logit.
186
Análisis Logit - ejemplo
50 estudiantes tomaron un examen, donde solo 27 pasaron.
¿Cuáles son las posibilidades de pasar?
Posibilidades = P/(1-P) = 0.54/0.46 = 1.17 o 1.71:1
187
Análisis Probit
Es similar a las pruebas de vida acelerada y análisis
de sobrevivencia. Un artículo sujeto a esfuerzo puede
fallar o sobrevivir. El modelo probit tiene un valor
esperado de 0 y una varianza de 1.
189
VI.B Pruebas de hipótesis
1. Conceptos fundamentales
2. Estimación puntual y por intervalo
3. Pruebas para medias, varianzas y proporciones
4. Pruebas comparativas para varianzas, medias y prop.
5. Bondad de ajustes
6. Análisis de varianza (ANOVA)
7. Tablas de contingencia
8. Pruebas no paramétricas
190
VI.B.1 Conceptos
fundamentales
191
Análisis Estadístico
En CADA prueba estadística, se comparan algunos valores
observados a algunos esperados u otro valor observado
comparando estimaciones de parámetros (media, desviación
estándar, varianza)
Hipótesis alterna Ha
Es la hipótesis que se acepta como verdadera cuando se
193
Conceptos fundamentales
Ejemplos:
Se está investigando si una semilla modificada
proporciona una mayor rendimiento por hectárea, la
hipótesis nula de dos colas asumirá que los
rendimientos no cambian Ho: Ya = Yb
194
Conceptos fundamentales
Estadístico de prueba
Para probar la hipótesis nula se calcula un estadístico
de prueba con la información de la muestra el cual se
compara a un valor crítico apropiado. De esta forma se
toma una decisión sobre rechazar o no rechazar la Ho
inversamente
Incrementando el tamaño de muestra se pueden reducir ambos
riesgos.
Decisión realizada Ho en realidad es Ho en realidad es
Verdadera falsa
No hay evidencia para p = 1- p=
rechazar Ho Decisión correcta Error tipo II
Rechazar Ho p= p=1-
Error tipo I Decisión correcta
196
Conceptos fundamentales
Pruebas de dos colas
Si la Ho: , , , = cte. que un valor poblacional,
entonces el riesgo alfa se reparte en ambos extremos
de la distribución. Por ejemplo si Ho = 10 se tiene:
197
Conceptos fundamentales
Pruebas de una cola
Si la Ho: , , , >= Cte. que un valor poblacional,
entonces el riesgo alfa se coloca en la cola izquierda de
la distribución. Por ejemplo si Ho: >= 10 y Ha: < 10
se tiene una prueba de cola izquierda:
198
Conceptos fundamentales
Pruebas de una cola
Si la Ho: , , , <= Cte. que un valor poblacional,
entonces el riesgo alfa se coloca en la cola derecha de
la distribución. Por ejemplo si Ho: <= 10 y Ha: >
10 se tiene una prueba de cola derecha:
199
Conceptos fundamentales
Tamaño de muestra requerido
Normalmente se determina el error alfa y beta deseado
y después se calcula el tamaño de muestra necesario
para obtener el intervalo de confianza.
Z 2 / 2 2
n
E2
Z 2 / 2 ( p )(1 p )
n
( p ) 2 Z 2 / 2 2
n
(X )2
Z 2 / 2 ( )( 1 )
n
( p )2 201
Conceptos fundamentales
Ejemplo:
¿Cuál es el tamaño de muestra mínimo que al 95% de
nivel de confianza (Z=1.96) confirma la significancia de
una corrida en la media mayor a 4 toneladas/hora (E),
si la desviación estándar (sigma) es de 20 toneladas?
n = (1.96^2)(20^2)/(4)^2 = 96
203
Efecto del tamaño de muestra
204
Efecto del tamaño de muestra
205
Efecto del tamaño de muestra
206
Potencia de la prueba
La potencia de una prueba estadística es su habilidad
para detectar una diferencia crítica
Potencia 1
Si Beta = 0.1 la potencia es del 90%
208
Potencia de la prueba
Para estimar la potencia, Minitab requiere de dos de
los siguientes parámetros:
Tamaños de muestra
Diferencias - un corrimiento significativo de la media
que se desea detectar
Valores de potencia - La probabilidad deseada de
rechazar Ho cuando es falsa
209
Considerando la potencia de prueba
210
VI.2 Significancia estadística
vs práctica
211
Estimación de riesgos
212
Pruebas de Minitab
Permite hacer las siguientes pruebas:
Prueba z de una muestra
Prueba t de una muestra
Prueba t de dos muestras
Prueba de 1 proporción
Prueba de 2 proporciones
ANOVA
Diseños factoriales de dos niveles
Diseños de Packett Burman
213
Calculo manual
214
Calculo manual
215
VI.3 Tamaño de muestra
216
Calculo manual de tamaño de
muestra
217
Calculo manual de tamaño de
muestra – Pruebas de una cola
218
Calculo manual de tamaño de
muestra – Pruebas de una cola
219
Ejemplo con prueba de una media t
Ejemplo: Se tiene una población normal con media de 365 y
límites de especificación de 360 y 370. Si la media se desplaza
2.5 gramos por arriba de la media, el número de defectos sería
inaceptable, la desviación estándar histórica es de 2.403:
CORRIDA DE 2.5 GRS. EN PROMEDIO
Ha: Corrida
0.18 Ho: 367.5 Variable
Meta LIE 370
O riginal
0.16 LIE 360 C orrida
365
0.14
0.12
0.10
Y-Data
0.08
0.06
0.04
0.02
0.00
220
Ejemplo con prueba de una media t
Stat > Power and Sample Size > 1 - Sample t
Completar el diálogo como sigue:
221
Ejemplo con prueba de una media t
Los resultados se muestran a continuación:
Power and Sample Size
1-Sample t Test
Testing mean = null (versus not = null)
Calculating power for mean = null + difference
Alpha = 0.05 Assumed standard deviation = 2.403
222
Ejemplo con prueba de una media t
Stat > Power and Sample Size > 1 - Sample t
Sample Target
Difference Size Power Actual Power
2.5 10 0.80 0.832695
2.5 11 0.85 0.873928
2.5 12 0.90 0.905836
2.5 15 0.95 0.962487
225
Ejemplo con prueba de 1 proporción
Test for One Proportion
Testing proportion = 0.02 (versus > 0.02)
Alpha = 0.05
Alternative Sample Target
Proportion Size Power Actual Power
0.04 391 0.8 0.800388
0.04 580 0.9 0.900226
Si se desea saber la Potencia si se utiliza un tamaño de muestra de 500 se tiene:
Stat > Power and Sample Size > 2 - Sample t
Sample sizes = 500 Alternative values of p = 0.04
Options: Greater Than
Significance Level = 0.05
231
Estimación puntual y por
intervalo
Las medias o desviaciones estándar calculadas de una muestra
se denominan ESTADÍSTICOS, podrían ser consideradas como
un punto estimado de la media y desviación estándar real de
población o de los PARAMETROS.
232
Intervalo de confianza
Error de estimación
234
Estimación puntual y por
intervalo
95% de Nivel de Confianza significa que sólo tenemos un 5% de
oportunidad de obtener un punto fuera de ese intervalo.
p (1 p )
p Z
2
n
236
Para n grande el IC es pequeño
237
Para n grande el IC es pequeño
238
Ejemplo
Dadas las siguientes resistencias a la tensión: 28.7,
27.9, 29.2 y 26.5 psi
239
Ejemplos para la media con
Distribución normal Z
240
Ejemplos para la media y varianza
con Distribución t
241
Ejemplos para proporciones con
Distribución Z
242
Instrucciones con Minitab
Intervalo de confianza para la media
En Options:
OK
243
Instrucciones con Minitab
Intervalo de confianza para proporción
En Options:
Indicar el Confidence Interval -- 90, 95 o 99%
244
VI.B.5 Pruebas de hipótesis para
medias, varianzas y proporciones
245
Elementos de una
Prueba de Hipótesis
246
Elementos de una
Prueba de Hipótesis
247
Pasos en la Prueba de Hipótesis
1. Definir el Problema - Problema Práctico
248
Pasos en la Prueba de Hipótesis
5. Establecer las Hipótesis
- Hipótesis Nula (Ho) - Siempre tiene el signo =, ,
Ha : , 2 , , , parametro de la hipotesis
-Z 0 Z
Pruebas de Hipótesis de cola derecha:
Ho: a b
Ha: a > b Región de
Rechazo
0 Z
Pruebas de Hipótesis cola izquierda:
Ho: a b
Ha: a < b Región de
Rechazo
-Z 0 Z
250
Pasos en la Prueba de Hipótesis
6. Seleccionar el nivel de Alfa (normalmente 0.05 o 5%) o el
nivel de confianza NC = 1 - alfa
d
t ; Pares.de.medias; d i . para.cada. par
Sd / n
( n 1) S 2
X 2
; DF ( n 1); prueba.una.v ar ianza
2
(O E ) 2
X 2
; DF ( r 1)(c 1); bondad .ajuste
E
253
Pasos en la Prueba de Hipótesis
11. Obtener el estadístico correspondiente de tablas o Excel.
-Z 0 Z
Pruebas de Hipótesis de cola derecha:
Ho: a b
Ha: a > b Región de
Rechazo
0 Z
Pruebas de Hipótesis cola izquierda:
Ho: a b
Ha: a < b Región de
Rechazo
-Z 0 Z
255
Prueba de hipótesis para la varianza
Las varianzas de la población se ditribuyen de
acuerdo a la distribución Chi Cuadrada. Por tanto las
inferencias acerca de la varianza poblacional se
basarán en este estadístico
257
Prueba de hipótesis para la varianza
2.17
5.99
261
Ejemplo de Prueba de hipótesis para la media
Para una muestra grande (n>30)probar la hipótesis de una media u
1.) Ho:
2.) Ha:
3.) Calcular el estadístico de prueba
4.) Establecer la región de rechazo
Las regiones de rechazo para prueba de 2 colas: -Z Z
Zcalc= s
n
Región de Región de
Rechazo Rechazo
0
-Z
s
IC. para.estimar. X Z Error estándar
Z alfa/2
101.5
1.95996398
2 n
Intervalo de confianza 2251 198.936344
El intervalo de confianza incluye a la media de la hipótesis
por tanto no se rechaza la Ho. 2052.064 <= <= ### ) 264
Prueba de hipótesis de una población
para muestras pequeñas con t
Se piensa que las ventas promedio de $5,775 se han incrementado gracias a la campaña publicitaria
Probar esta afirmación a un nivel de significancia alfa de 1%
Se inicia con el planteamiento de la hipótesis Alterna
Datos
Semanas n 15 media mu = 5775
Ventas prom X 6012
Desv. Estándar s 977 (Alfa = 0.01 (1-Alfa = 0.99
(Alfa/2 = 0.005 (1-Alfa/2 = 0.995
Paso 1. Establecimiento de hipótesis
s
IC . para.estimar . X t Error estándar 252.260315
2 n Z alfa/2 2.62449406
Como el intervalo de confianza Intervalo de confianza 6012 662.0557002
contiene a la media Hipótesis no se rechaza Ho 5349.94 <=
<= 6674.06 )
266
Prueba de hipótesis
para una proporción con Z
El gerente de mercado considera que el 50% de sus clientes gasta menos de $10 en cada visita a la tienda.
¿Estás de acuerdo con esta afirmación a un nivel de significancia del 5%?
Se inicia con el planteamiento de la hipótesis nula
Datos
Clientes n 50 Proporción media = 0.5
30 gastaron p 0.6
menos de$10 (Alfa = 0.05 (1-Alfa = 0.95
(Alfa/2 = 0.025 (1-Alfa/2 = 0.975
Paso 1. Establecimiento de hipótesis
Ho : c 0.5
Ha : c 0.5 Se trata de una prueba de dos colas
Paso 2. Cálculo del estadístico de prueba Zc
En Options:
Indicar el Confidence Interval -- 90, 95 o 99%
269
Instrucciones con Minitab para la
prueba de hipótesis de una proporción
En Options:
Indicar el Confidence Interval -- 90, 95 o 99%
Indicar la Test Proportion Proporción de la hipótesis
Indicar el signo de la hipótesis alterna: Less Than, Not equal,
Greater than
270
Pruebas de hipótesis para
comparación de varianzas, medias, y
proporciones
271
Prueba de Hipótesis
Supongamos que tenemos muestras de dos reactores que
producen el mismo artículo. Se desea ver si hay diferencia
significativa en el rendimiento de “Reactor a Reactor”.
Reactor A Reactor B
89.7 84.7 Estadísticas Descriptivas
81.4 86.1
84.5 83.2 Variable Reactor N Media Desv.Std
84.8 91.9
87.3 86.3 Rendimiento A 10 84.24 2.90
79.7 79.3 B 10 85.54 3.65
85.1 82.6
81.7 89.1
83.7 83.7
84.5 88.5
272
Prueba de Hipótesis
Pregunta Práctica: Existe diferencia entre los reactores?
Ho: a b
Ha: a b
274
¿Qué representa esto?
Reactor A Reactor B
B B B B B BB BB B
A AA AAAA A A
80.0 82.5 85.0 87.5 90.0 92.5
¿Representan los reactores dos procesos diferentes?
276
Prueba F de dos varianzas
Método 1 Método 2
No. De CDs n1 15 n2 17 Alfa/2 0.025
Desv. Estan. s1 5.4 X2 4.8
2 2
Varianza s1 29.16 s2 23.04
Ho : 12 2
2
Ha : 12 2
2
Por tanto se trata de una prueba de dos colas
2
Tomamos a s1 como el mayor para comparar Fc contra Fexcel (1- Alfa/2)
f(F)
Fe(0.025) = 2.81701784
Ze (0.90)= 1.28
Zc = 1.03099301 Valor p para Zc es igual a
P(-Zc) = 0.149402368
p > Alfa
Como Zc es menor que Zexcel, no cae en el área de rechazo,
y por tanto no hay suficiente evidencia para rechazar Ho
Se concluye que no hay diferencia entre vacas de vacaciones y normales
Paso adicional. El Intervalo de confianza del 90% sobre la diferencia de medias poblacionales,
con sigmas desconocidas es:
s12 s 22
X 1 X 2 = Error estándar 6.11061372
n1 n 2 Z (alfa/2) = 1.64485363
( X 1 X 2 ) Z / 2 s X 1 X 2 = Intervalo de confianza
6.3 + - 10.05106514
La diferencia es del orden de cero,es decir ( -3.75107 < = u < = 16.3511 )
281
Prueba de dos medias
muestras pequeñas
Sigmas descono-
cidas e iguales
Sigmas desconocidas
y desiguales 282
Prueba de hipótesis de dos pob.
Comparando dos medias con t
Investigar si hay diferencia en los promedios de las ventas diarias de dos tiendas
Las varianzas de las dos poblaciones son iguales pero desconocidas 1 2
2 2
Ho : T 1 T 2
Ha : T 1 T 2 Por tanto se trata de una prueba de dos colas
2
s12 (n1 1) s 22 (n 2 1) 19564.25 Sp = 782.57
s
2
p
n1 n 2 2 25
Paso 3. Determinar la te de Excel o de tablas para una alfa de 0.01 que corresponde a alfa/2 = 0.005
Se tienen n1 + n2 - 2 grados de libertad o sean 25
te (0.01) = 2.78743581 DIST.T.INV (0.01, 25)
283
Asi es para dos colas
Paso 4. Comparando los valores tc calculado contra texcel (0.01) se tiene
Zc = 3.39304676
Ze(0.9) = 1.281551566
Valor p para Zc es igual a
P(-Zc) = 0.00034946
p < Alfa
Como Zc es mayo que Zexcel, si cae en el área de rechazo,
y por tanto hay suficiente evidencia para rechazar Ho y aceptar Ha
Se concluye que la diferencia en conv. entre los bonos es significativa
Paso adicional. El Intervalo de confianza del 98% sobre la diferencia de medias poblacionales,
con sigmas desconocidas es:
p1 (1 p1 ) p 2 (1 p 2 )
s p1 p 2 = Error estándar 0.044171193
n1 n2 Zexcel (para alfa/2) 1.644853627
289
Resumen
290
Instrucciones con Minitab para la
comparación de dos varianzas
En Options:
Indicar el Confidence Interval -- 90, 95 o 99%
OK
291
Instrucciones con Minitab para la
comparación de dos medias
En Options:
Indicar el Confidence Interval -- 90, 95 o 99%
Indicar la diferencia a probar Test Difference (normalmente 0)
Indicar el signo de la hipótesis alterna: Less Than, Not equal,
Greater than
En graphs seleccionar las graficas Boxplot e Individual value plot
OK
292
Instrucciones con Minitab para la
comparación de dos medias pareadas
En Options:
Indicar el Confidence Interval -- 90, 95 o 99%
Indicar la diferencia a probar Test Mean (normalmente 0)
Indicar el signo de la hipótesis alterna: Less Than, Not equal,
Greater than
En graphs seleccionar las graficas Boxplot e Individual value plot
OK
293
Instrucciones con Minitab para la
prueba de hipótesis de dos proporciones
En Options:
Indicar el Confidence Interval -- 90, 95 o 99%
Indicar la Test Difference Normalmente 0
Indicar el signo de la hipótesis alterna: Less Than, Not equal,
Greater than
294
VI.B.7 Pruebas de bondad de
ajuste
295
Bondad de ajuste
PRUEBAS DE BONDAD DE AJUSTE
Medidas sobre que tan cerca se ajustan los datos muestrales observados a una forma
de distribución particular planteada como hipótesis
Si el ajuste es razonablemente cercano, puede concluirse que sí exite la forma de distribución
planteada como hipótesis
Por ejemplo:
(Oi Ei) 2
K
2
i 1 Ei
El estadístico Chi cuadrado calculado de 1.17 es menor al de excel de 7.815 por tanto se acepta
la hipótesis nula
PRUEBA.CHI.INV 297
Prueba de Bondad de ajuste
para la distribución de Poisson
1. Plantear la hipótesis nula y alterna
Ho: La población tiene una distribución de prob. De Poisson
Ha: Caso contrario
2. Tomar una muestra aleatoria, anotar la frecuencia observada fi
y calcular la media de ocurrencias
3. Calcular la frecuencia esperada de ocurrencias ei. Multiplicar el
tamaño de muestra con la prob. de Poisson para cada valor de
la variable aleatoria. Si hay menos de 5 combinar las categorías
n
( f i ei ) 2
2
4. Calcular el estadístico de prueba i 1 ei
2 2
5. Rechazar Ho si o si p < alfa. Con gl=k-p-1 y alfa nivel de
significancia
298
Ejemplo:
Distribución de Poisson =5
Ho: No. de clientes que llega en intervalos de 5 min. tiene una distribución
de Poisson Ha: No se sigue una distribución de Poisson
( f i ei ) 2
n
2
i 1 ei
Ho se rechaza si o si p es mayor que alfa.
2
2
301
Prueba de Bondad de ajuste
para la distribución Normal
1. Plantear la hipótesis nula y alterna
Ho: La población tiene una distribución de prob. Normal
Ha: Caso contrario
302
Prueba de Bondad de ajuste
para la distribución Normal
4. Calcular el número de ocurrencias esperado ei, para cada
intervalo de valores. Multiplicar el tamaño de muestra por la
probabilidad de que una variable aleatoria esté en el intervalo.
n
( f i ei ) 2
5. Calcular el estadístico de prueba 2
i 1 ei
de significancia
303
Prueba de Bondad de ajuste
para la distribución Normal
Ejemplo: datos de calificaciones: Media = 68.42; S = 10.41
Calificaciones
71 66 61 65 54 93
60 86 70 70 73 73
55 63 56 62 76 54
82 79 76 68 53 58
85 80 56 61 61 64
65 62 90 69 76 79
77 54 64 74 65 65
61 56 63 80 56 71
79 84 304
Prueba de Bondad de ajuste
para la distribución Normal
Ho: la población tiene una distribución normal con media 68.42 y
S=10.41 Ha: Caso contrario
305
Prueba de Bondad de ajuste
para la distribución Normal
Intervalo Frecuencia Frecuencia Se registran las
observada (fi) esperada (ei) frecuencias de
Menos de 5 5 los datos
55.10 tomados de las
55.10 a 59.68 5 5 calificaciones
59.68 a 63.01 9 5
63.01 a 65.82 6 5
65.82 a 68.42 2 5
68.42 a 71.02 5 5
71.02 a 73.83 2 5
73.83 a 77.16 5 5
77.16 a 81.74 5 5
81.74 o más 6 5
50 50
306
Prueba de Bondad de ajuste
para la distribución Normal
Se determina el estadístico Chi Cuadrado = 7.2
( f i ei ) 2
n
2
i 1 ei
Como 2
2
no se puede rechazar la hipótesis nula de
307
Prueba de Bondad de ajuste
para la distribución Multinomial
1. Enunciar la hipótesis nula y alternativa
Ho: La población sigue una distribución de probabilidad
multinomial con probabilidades especificadas para cada una de
las K categorías Ha: Caso contrario
308
Prueba de Bondad de ajuste
para la distribución Multinomial
4. Se determina el estadístico Chi Cuadrado de prueba
( f i ei ) 2
n
2
i 1 ei
5. Regla de rechazo:
2 2
Si no se puede rechazar la hipótesis nula
309
Prueba de Bondad de ajuste
para la distribución Multinomial
Ejemplo: El año pasado la participación de mercado para la
empresa A fue del 30%, 50% para la empresa B y 20% para la
empresa C. La empresa C hace una prueba con un nuevo
producto para estimar su impacto en las preferencias del
mercado.
310
Prueba de Bondad de ajuste
para la distribución Multinomial
Datos para calcular el estadístico de prueba Chi Cuadrado
Empresa A 0.3 48 60
Empresa C 0.2 54 40
311
Prueba de Bondad de ajuste
para la distribución Multinomial
Chi Cuadrado calculado = 7.34
312
Prueba de Bondad
de ajuste en Minitab
La columna C1 – Observadas contiene las frecuencias observadas
y la C2 – esperadas las frecuencias esperadas
314
Prueba de Bondad
de ajuste en Excel
Ejemplo: investigación de mercado
315
VI.B.6 ANOVA para un factor
principal y una o más variables de
bloqueo
316
Introducción
Cuando es necesario comparar 2 o más medias poblacionales al
mismo tiempo, para lo cual se usa ANOVA.
distribuyen normalmente
El término de error tiene un efecto distribuido normalmente
e independiente
317
Contenido
ANOVA de un factor o dirección
318
ANOVA de un factor
o dirección
319
Introducción
Con el ANOVA las variaciones en la respuesta se dividen
en componentes que reflejan los efectos de una o más
variables independientes
320
ANOVA – Prueba de hipótesis para
probar la igualdad de medias de
varias poblaciones para un factor
Se trata de probar si el efecto de un factor o
Tratamiento en la respuesta de un proceso o sistema es
Significativo, al realizar experimentos variando
Los niveles de ese factor (Temp. 1, Temp. 2, Temp.3, etc.)
Ho : 1 2 3 ......... a
Ha : A lg unas. ' s.son.diferentes
321
ANOVA - Condiciones
Todas las poblaciones son normales
Xij
Gran media
Xij
a b 2
SCT ( Xij X )
i 1 j 1
324
ANOVA – Suma de cuadrados de
renglones (a)-tratamientos
Media Trat. 1 Media Trat. a
a renglones
Gran media
Media trat. 2
SCTr b( X i X ) 2
i 1
325
ANOVA – Suma de cuadrados
del error
X2j X3j
X1j
Media X1.
Media X3.
Media X2.
a b
SCE (X ij X i) 2
i 1 j 1 326
ANOVA – Suma de cuadrados
del error
X2j X3j
X1j
Media X1.
Media X3.
Media X2.
gl.SCT n 1
gl.SCTr a 1
gl.SCE (n 1) (a 1) n a
328
ANOVA – Cuadrados medios:
Total, Tratamiento y Error
MCT SCT /( n 1)
MCTr SCTr /( a 1)
MCE SCE /( n a )
329
ANOVA – Cálculo del estadístico
Fc y Fexcel
MCTr
Fc
MCE
Fexcel FINVALFA , gl .SCTr , gl .SCE
330
Tabla final de ANOVA
TABLA DE ANOVA
Regla: Rechazar Ho si la Fc de la muestra es mayor que la F de Excel para una cierta alfa
o si el valor p correspondiente a la Fc es menor al valor de alfa especificado
331
ANOVA – Toma de decisión
Distribución F Fexcel
Alfa
Fc
332
ANOVA – Toma de decisión
O si el valor de p correspondiente a Fc es
menor de Alfa se rechaza Ho
333
ANOVA – Identificar las medias
diferentes por Prueba de Tukey T
CME
T q , a , n a
b
Para diseños balanceado
(mismo número de columnas
en los tratamientos) el valor
de q se determina por medio de
la tabla en el libro de texto
334
ANOVA – Identificar las medias
diferentes por Prueba de Tukey T
Se calcula la diferencia Di entre cada par de Medias Xi’s:
D1 = X1 – X2 D2 = X1 – X3 D3 = X2 – X3 etc.
335
ANOVA – Identificar las medias
diferentes por Prueba de Diferencia
Mínima Significativa DMS
2(CME ) F ,1,n a
DMS
b
Para diseños balanceados (los tratamientos
tienen igual no. De columnas), se calcula un
factor DMS contra el que se comparan las
diferencias Xi – Xi’. Significativas si lo exceden
336
Prueba DMS para Diseños no
balanceados
1 1
DMS j ,k (CME ) F ,a 1,n a
b j bk
Para diseños no balanceados (los
tratamientos tienen diferente no. De
columnas), se calcula un factor DMS
Para cada una de las diferencias Xi – Xi’
337
Ejemplo:
Considerar un experimento de un factor (máquina)
con tres niveles (máquinas A, B, C). Los datos se
muestran a continuación y debe verificarse si existe
diferencia significativa a un alfa = 0.05
338
Ejemplo:
A B C
5 2 1
7 0 0
6 1 -2
7 -2 -3
6 2 0
340
Ejemplo:
One-way ANOVA: A, B, C
Source DF SS MS F P
341
Corrida en Minitab
Se introducen las respuestas en una columna C1
Se introducen los subíndices de los renglones en una
columna C2 Durability Carpet
18.95 1
12.62 1
11.94 1
14.42 1
10.06 2
7.19 2
7.03 2
14.66 2 342
Corrida en Minitab
Opción: stat>ANOVA – One Way (usar archivo
Exh_aov)
En Response indicar la col. De Respuesta (Durability)
Resultados
Analysis of Variance for Durabili
Source DF SS MS F P
Carpet 3 111.6 37.2 2.60 0.101
Error 12 172.0 14.3
Total 15 283.6
Individual 95% CIs For Mean
Based on Pooled StDev
Level N Mean StDev ---------+---------+---------+-------
1 4 14.483 3.157 (-------*-------)
2 4 9.735 3.566 (-------*--------)
3 4 12.808 1.506 (--------*-------)
4 4 17.005 5.691 (-------*-------)
---------+---------+---------+-------
Pooled StDev = 3.786 10.0 15.0 20.0
Tukey's pairwise comparisons
Family error rate = 0.0500
Individual error rate = 0.0117
Critical value = 4.20
344
ANOVA de dos vías un factor
principal y una variable de
bloqueo
345
ANOVA de 2 vías
Este es un procedimiento extensión de los patrones
del ANOVA de una vía con tres fuentes de variación:
Tratamiento del factor A (columnas), Tratamiento del
factor B (renglones) y Error experimental.
346
ANOVA – Prueba de hipótesis para
probar la igualdad de medias de
varias poblaciones con dos vías
Se trata de probar si el efecto de un factor o
Tratamiento en la respuesta de un proceso o sistema es
Significativo, al realizar experimentos variando
Los niveles de ese factor (Temp.1, Temp.2, etc.)
POR RENGLON
Y
Considerando los niveles de otro factor que se piensa
Que tiene influencia en la prueba – FACTOR DE BLOQUEO
POR COLUMNA
347
ANOVA – 2 vías
Para el tratamiento – en renglones
Ho : 1 2 3 ......... a
Ha : A lg unas. ' s.son.diferentes
Para el factor de bloqueo – en columnas
350
ANOVA de 2 vías
b
SCBl a ( X j X ) 2
j 1
gl.SCBl b 1
CMBl SCBl /(b 1)
351
ANOVA de 2 vías
352
ANOVA –Estadístico Fc y Fexcel
MCTr
Fc
MCE
Fexcel FINVALFA , gl .SCTr , gl .SCE
353
ANOVA – Estadístico Fb
MCBl
Fc
MCE
Fexcel FINVALFA , gl .SCBl , gl .SCE
354
Tabla final ANOVA 2 vías
FUENTE DE VARIACIÓN SUMA DE GRADOS DE CUADRADO VALOR F
CUADRADOS LIBERTAD MEDIO
Regla: No rechazar si la F de la muestra es menor que la F de Excel para una cierta alfa
355
ANOVA – 2 vías: Toma de decisión
Distribución F Fexcel
Alfa
Fc
Tr o Bl 356
ANOVA – 2 vías: Toma de decisión
357
Cálculo de los residuales
ˆ ij yi . y. j y..
y Y estimada
eij yij y
ˆ ij Error o residuo
359
Corrida en Minitab
Se introducen las respuestas en una columna C3 y los
subíndices de renglones en columna C4 y de columnas en C5
Resultados
Suppleme 2 1919 959 9.25 0.015
Lake 1 21 21 0.21 0.666
Interaction 2 561 281 2.71 0.145
Error 6 622 104
Total 11 3123
Individual 95% CI
Suppleme Mean --+---------+---------+---------+---------
1 43.5 (-------*-------)
2 68.3 (--------*-------)
3 39.8 (--------*-------)
--+---------+---------+---------+---------
30.0 45.0 60.0 75.0
Individual 95% CI
Lake Mean ------+---------+---------+---------+-----
Dennison 51.8 (----------------*----------------)
Rose 49.2 (----------------*----------------)
------+---------+---------+---------+-----
42.0 48.0 54.0 60.0
362
ANOVA de un factor y dos o
tres variables de bloqueo
CUADRADO LATINO Y
GRECOLATINO
363
ANOVA – 3 y 4 factores
El diseño de Cuadrado latino utiliza dos factores de
bloqueo adicionales al de Tratamiento
364
Cuadrado Latino
Años exp. Turno
Empleado Mañana Tarde Noche
1 B=15 A=18 C=11
j 1
gl.SCTr a 1 b 1
CMTr SCTr /(b 1)
366
ANOVA – Cuadrado Latino: Cálculo
del error
367
ANOVA – Cálculo del estadístico
Fc y Fexcel
MCTr
Fc
MCE
Fexcel FINVALFA , gl .SCTr , gl .SCE
368
ANOVA – Cuadrado Latino Reng / Col
MC Re ng
Fcreng
MCE
MCCols
Fcols
MCE
Fexcel FINVALFA , gl .SCBl , gl .SCE
369
Tabla final ANOVA 2 Factores
FUENTE DE VARIACIÓN SUMA DE GRADOS DE CUADRADO VALOR F
CUADRADOS LIBERTAD MEDIO
370
Cuadrado latino en Minitab
Se introducen las respuestas en una columna C1
Lotes MP 1 2 3 4 5
373
Cuadrado Greco latino en Minitab
Se introducen las respuestas en una columna C1
Se introducen los subíndices de los renglones en una
columna C2
375
ANOVA – Cuadrado Grecolatino
b
SCG a ( X m X ) 2
m 1
gl.SCG b 1
CMG SCG /(b 1)
376
ANOVA de 2 factores – Suma de
cuadrados, gl. y Cuadrado medio
para el error
377
ANOVA – Cálculo del estadístico
Fc y Fexcel
MCG
Fc
MCE
Fexcel FINVALFA , gl .SCTr , gl .SCE
378
ANOVA – Cuadrado Grecolatino
MCTr
Fc
MCE
Fexcel FINVALFA , gl .SCBl , gl .SCE
379
Tabla final ANOVA 2 Factores
FUENTE DE VARIACIÓN SUMA DE GRADOS DE CUADRADO VALOR F
CUADRADOS LIBERTAD MEDIO
380
ANOVA para diseño factorial AxB
En un experimento factorial involucrando el factor A con (a)
niveles y un factor B con (b) niveles, la suma de cuadrados se
puede dividir en:
SST = SS(A) + SS(B) + SS(AB) + SSE
381
VI.B.8 Tablas de contingencia
382
¿Para qué se utiliza?
2
383
Ejemplo 1: Chi Cuadrada( 2 )
384
Ejemplo 1: Chi Cuadrada( 2 )
Aguilas 63 50 3.38
Soles 37 50 3.38
2 = 3.38 + 3.38
2 = 6.76
385
Ejemplo 1: Chi cuadrada
Función de Distribución Acumulada Chi2 con 1 grado de libertad
(d.f)
2c P(2c > x)
6.7600 p = 1 - 0.9907 = 0.0093
387
Tabla de Valores Críticos Seleccionados de Chi2
df .250 .100 .050 .025 .010 .005 .001
1 1.323 2.706 3.841 5.024 6.635 7.879 10.828
2 2.773 4.605 5.991 7.378 9.210 10.597 13.816
3 4.108 6.251 7.815 9.348 11.345 12.838 16.266
4 5.385 7.779 9.488 11.143 13.277 14.860 18.467
5 6.626 9.236 11.070 12.832 15.086 16.750 20.515
388
100 109.141 118.498 124.342 129.561 135.807 140.169 149.449
Tabla de contingencia
Una tabla de clasificación de dos vías (filas y columnas) que
contiene frecuencias originales, se puede analizar para
determinar si las dos variables (clasificaciones) son
independientes o tienen una asociación significativa.
389
Tabla de contingencia
Para esta prueba se usa la prueba Chi Cuadrada donde:
390
Tabla de contingencia
Ejemplo: Cada una de las 15 celdas hace una contribución al
estadístico Chi Cuadrado (una celda)
391
Ejemplo 2: Chi2 Para comparación de dos
grupos; ¿son las mismas proporciones?)
Ho: No existen diferencias en los índices de defectos de las dos máquinas.
Ha: Existen diferencias en los índices de defectos de las dos máquinas.
392
Ejemplo 2: Chi2 Para comparación de dos
grupos; ¿son las mismas proporciones?)
Cálculo de los valores esperados
2 232 3 235
233.47 1.53
Total 764 5 769
Nota: Chi cuadrada no podrá aplicarse en los casos donde los conteos seas menores a 5 en 20%
de celdas.
Si cualquiera de los conteos esperados en las celdas es menor a uno, no deberá usarse Chi 2.
Si algunas celdas tienen un conteo menor a los esperados, ya sea combinando u omitiendo
renglones y/o columnas, las categorías pueden ser de utilidad.
394
Tabla de Chi2
Tabla de valores críticos seleccionados para Chi 2
395
Problema: Fugas
Beneficios Potenciales: $10,000 de ahorro en retrabajos, y en la
reducción de tiempo de ciclo.
Máquina a máquina
Ho: No existe diferencia en los índices de defecto de las diferentes
máquinas
Ha: Existe diferencia en los índices de defecto de las diferentes
máquinas
Tamaño de la muestra:
5000 + total de oportunidades (172 piezas) 396
Prueba de chi2 (máquina a máquina)
Los conteos esperados están colocados debajo de los conteos observados
Con fugas Sin fugas Total
1 30 610 640
32.11 607.89
3 3 253 256
12.84 243.16
4 18 334 352
17.66 334.34
2 1 127 128
6.61 121.39
4 54 202 256
13.23 242.77
5 5 699 704
36.38 667.62
6 12 116 128
6.61 121.39
Total 278 5102 5380
Chi2 = 0.057 + 0.003 + 4.765 + 0.260 + 0.351 + 0.019 +125.666 + 6.847 + 27.065 + 1.475
+ 4.386 + 0.239 = 171.132
DF= 5; valor P = 0.000 398
¿Qué sucede si los grupos múltiples de variación son estadísticamente significativos?
(en este caso, operador a operador y máquina a máquina)
6 12 116 128
6.61 121.39
(Estos mismos operadores fueron quienes
tuvieron los números más grandes de chi2)
400
Operador a operador: = 0.000
Rechace
Ho y acepte Ha
(Existe una diferencia significativa entre los operadores)
¿Cuál es el próximo paso? Hable con todos los operadores para averiguar qué diferencias
pueden existen en sus técnicas.
El operador 5 encontró un modo de mejor de hacer el ensamble, con lo cual consiguió mejorar
el trabajo de soldadura, aunque esto mostraba un grado de dificultad ergonómica. Se añadió
un colocador para ensamblar la parte en forma segura. (Esto también redujo el tiempo que
requerían los operadores para “acostumbrarse” a trabajar en esta forma)
401
Ejercicios
Radios detectados 27 25 22
Radios no detectados 3 5 8
403
Coeficiente de Contingencia
Coeficiente de contingencia es el grado de relación o
dependencia de las clasificaciones en la tabla de contingencias
es:
X2
C2
X2 N
404
Coeficiente de Contingencia
Para los datos del ejemplo anterior se tiene:
X2 66.22 2
C2 2 0.38
X N
2
66.22 393
2
k 2 82
Max C 0.866
k 8
405
Correlación de atributos
Para tablas de orden k * k, el coeficiente de correlación, r, es :
2
X
r
N (k 1)
Donde 0<= r <= 1
406
VI.B.9 Pruebas de Hipótesis
no paramétricas
407
Pruebas no paramétricas
Las pruebas paramétricas asumen una distribución para la
población, tal como la Normal
408
Prueba de Hipótesis
Variable Atributo
No Normal Tablas de
Contingencia de
Varianza Medianas
Correlación
Correlación
Homogeneidad
Prueba de signos
de la Variación
de Levene Wilcoxon
Normal
Mann-
Whitney Variancia Medias
Kurskal-
Pruebas de t
Wallis
Prueba-F Muestra-1
Residuos
Prueba de Mood Muestra-2
Homogeneidad distribuidos
Friedman de la Variación ANOVA
de Bartlett Una vía normalmente
Dos vías
Correlación
Regresión 409
Resumen de pruebas de Hipótesis
Datos Normales Datos No Normales
Pruebas de Variancias
Pruebas de Varianzas
X2 : Compara la variancia de una Homogeneidad de la varianza de
muestra con una variancia de un Levene : Compara dos o más
universo conocido. varianzas de muestras de la misma
población.
Prueba F : Compara dos varianzas
de muestras.
Homogeneidad de la variancia de
Bartlett: Compara dos o más
varianzas muestras de la misma
población.
410
Resumen de pruebas de Hipótesis
Datos Normales Datos No Normales
Pruebas de los Promedios Pruebas de la Mediana
Mediana: Valor del punto medio de los datos, cuando se ordenan en forma
ascendente (en caso de datos pares, obtener promedio).
Moda : Valor que se repite con más frecuencia sobre el conjunto de datos.
Ejemplo:
Se cuestionó a veinte personas sobre cuánto tiempo les tomaba estar
listas para ir a trabajar, en las mañanas. Sus respuestas (en minutos) se
muestran más adelante. ¿Cuáles son el promedio y la mediana para esta
muestra?
Promedio
Mediana
415
Prueba de Rachas
Considere los siguientes datos (que se muestran aquí en orden cronológico):
325, 210, 400, 72, 150, 145, 110, 507, 56, 120, 99, 144, 110, 110,
320, 290, 101, 0, 80, 500, 201, 50, 140, 80, 220, 180, 240, 309, 80
500
Promedio
Primera 400
"corrida"
300
200
100
Segunda ”racha"
Racha: Un punto o una serie consecutiva de puntos que caen
en un lado del promedio.
Número total de Rachas: 12
Número total de puntos > al promedio: 11
Número total de puntos < al promedio: 18
Prueba de Rachas
Ho: Los datos son aleatorios
Ha:Los datos NO so aleatorios
Prueba de Rachas
Promedio K = 184.4483 Promedio
Z = (G - MediaG) / DesvStG
Runs Test: C1
Runs test for C1
Runs above and below K = 184.448
The observed number of runs = 12
The expected number of runs = 14.6552
11 observations above K, 18 below
P-value = 0.285
P > 0.05
No rechazar
Ho 419
Prueba de Signos de la Mediana
115 144
Después de todo, tal vez
esto SEA lo correcto. 423
Corrida en Minitab
Stat > Nonparametrics > 1-Sample sign Variable C1
Confidence interval 95% Test Median 115 Alternative Not equal
Sign Test for Median: Signos
Sign test of median = 115.0 versus not
= 115.0
424
Prueba de Signos de la Mediana
Para observaciones pareadas
T 0 T
n(n 1)(2n 1)
6
En este caso n = pares eliminando las que son iguales con dif. = 0 para el
trabajador 8.
Como Zc = 2.24 > Z0.025 se rechaza Ho, los métodos son diferentes
429
Prueba en Minitab para prueba de
mediana con Wilconox
File> Open worksheet > Exh_Stat
Stat > Nonparametrics > 1-Sample Wilconox Achievement
Variables C1 Test Median 77 77
Altenative Not equal 88
85
Wilcoxon Signed Rank Test: Achievement
Test of median = 77.00 versus median not = 77.00 74
for Wilcoxon Estimated 75
for Wilcoxon Estimated
N Test Statistic P Median 62
Achievement 9 8 19.5 0.889 77.50 80
70
Ho: Mediana = 77 Ha: Mediana <> 77
Como P de 0.889 >> alfa de 0.05 no se rechaza Ho 83
430
Prueba de Mann-Whitney
Se llevó a cabo un estudio que analiza la frecuencia del pulso en dos
grupos de personas de edades diferentes, después de diez minutos de
ejercicios aeróbicos.
n1 = 10 n2 = 9
Ta = 130.5 Tb = 55.5
432
Prueba de Mann-Whitney
Ho: Las distribuciones de frecuencias relativas de las poblaciones A y B son iguales
Ha: Las distribuciones de frecuencias relativas poblacionales no son idénticas
Ho: 1 = 2 Ha: 1 2 1, 2 = Medianas de las poblaciones
Ordenando los datos y asignándoles su posición relativa se tiene:
Ua = n1*n2 + (n1) * (n1 + 1) /2 - Ta
Ub = n1*n2 + (n2) * (n2 + 1) /2 - Tb
Ua + Ub = n1 * n2
Ua = 14.5 Ub = 79.5
Utilizando el estadístico Z y la distribución normal se tiene:
45 12.24
Z = [ (U - (n1* n2 / 2 ) / Raiz (n1 * n2 * (n1 + n2 + 1) / 12)
Con Ua y Ub se tiene:
Za = (14.5 - 45) / 12.24 = - 2.49 P(Z) = 0.0064 similar a la anterior
Zb = (79.5 -45) / 12.24 = 2.81 P(total) = 2 * 0.0064 = 0.0128 menor = 0.05
El valor crítico de Z para alfa 0.025 por ser prueba de dos colas, es 1.96.
Como Za > Zcrítico se rechaza la Hipótesis Ho de que las medianas son iguales.
135 5 -31 7 9 -5 -1 3 7 11
150 20 -16 22 24 10 14 18 22 26
140 10 -26 12 14 0 4 8 12 16
144 14 -22 16 18 4 8 12 16 20
154 24 -12 26 28 14 18 22 26 30
160 30 -6 32 34 20 24 28 32 36
edad
144 14 -22 16 18 4 8 12 16 20
136 6 -30 8 10 -4 0 4 8 12
148 18 -18 20 22 8 12 16 20 24
Una vez ajustados los "enlaces" (eventos de un mismo valor en ambos grupos de
información), Minitab usa este punto estimado para calcular el valor p.
Corrida en Minitab
Stat > Nonparametrics > Mann Whitney
First Sample C1 Second Sample C2 Conf. Level 95%
Alternative Not equal
n1 = 8 n2 = 10 n3 = 7 N = n1 + n2 + n3 = 25
Ta = 118 Tb = 111.5 Tc = 95.5
437
Prueba de Kruskal Wallis
Como H < 2 crítico, no se rechaza la Hipótesis Ho: Afirmando que no hay diferencia
entre las poblaciones
438
Corrida en Minitab
Stat > Nonparametrics > Kruskal Wallis
Response C1 Factor C2 OK
440
Prueba de Medianas de Mood
Se determina el estadístico Chi Cuadrada con:
(O E ) 2
2
E
Probar Ho: Todas las medianas son iguales
Ha: Al menos una mediana es diferente
441
Corrida con Minitab
Se les da a 179 participantes una conferencia con
dibujos para ilustrar el tema. Después se les da la
prueba OTIS que mide la habilidad intelectual. Los
participantes se clasificaron por nivel educativo 0-No
prof., 1-Prof., 2-Prepa
444
Diseños factoriales aleatorias
bloqueados de Friedman
Resultados de salida:
Se muestra el estadístico de prueba con distribución
445
Diseños factoriales aleatorias
bloqueados de Friedman
Ejemplo:
Se evalúa el efecto del tratamiento de una droga en
446
Diseños factoriales aleatorias
bloqueados de Friedman
Datos: EnzymeActivity Therapy Litter
0.15 1 1
0.26 1 2
0.23 1 3
0.99 1 4
0.55 2 1
0.26 2 2
-0.22 2 3
0.99 2 4
0.55 3 1
0.66 3 2
0.77 3 3
0.99 3 4 447
Diseños factoriales aleatorias
bloqueados de Friedman
Resultados:
Friedman Test: EnzymeActivity versus Therapy
blocked by Litter
449
Diseños factoriales aleatorias
bloqueados de Friedman
Resultados:
450
Diseños factoriales aleatorias
bloqueados de Friedman
Resultados:
451
Diseños factoriales aleatorias
bloqueados de Friedman
Resultados:
452
Prueba de igualdad de
varianzas de Levene
Se usa para probar la hipótesis nula de que las varianzas de k
múltiples poblacionales son iguales
453
Prueba de igualdad de
varianzas de Levene
Para dos muestras el procedimiento es como sigue:
Determinar la media
454
Rot Temp Oxygen
13 10 2
Prueba de igualdad 11 10 2
de Varianzas-Minitab 3
10
10
10
2
6
4 10 6
Se estudian tamaños de papa 7 10 6
inyectando con bacterias y
15 10 10
sujetas a diferentes
2 10 10
temperaturas. Antes del
ANOVA se verifica la 7 10 10
igualdad de varianzas 26 16 2
19 16 2
24 16 2
Stat > ANOVA > Test for
equal variances 15 16 6
22 16 6
Response Rot
18 16 6
Factors Temp Oxigen
20 16 10
Confidence level 95%
24 16 10
8 16 10
455
Resultados
456
Resultados
Test for Equal Variances: Rot versus Temp, Oxygen
95% Bonferroni confidence intervals for standard deviations
Temp Oxygen N Lower StDev Upper
10 2 3 2.26029 5.29150 81.890
10 6 3 1.28146 3.00000 46.427
10 10 3 2.80104 6.55744 101.481
16 2 3 1.54013 3.60555 55.799
16 6 3 1.50012 3.51188 54.349
16 10 3 3.55677 8.32666 128.862
Bartlett's Test (normal distribution)
459
El coeficiente de correlación de
rangos de Spearman (rs)
El coeficiente de correlación es una medida de la asociación que
requiere que ambas variables sean medidas en al menos una
escala ordinal de manera que las muestras u observaciones a
ser analizadas pueden ser clasificadas en rangos en dos series
ordenadas
6 d 2
6(5.5)
rs 1 1 0.03 0.97
990
460
Coeficiente de correlación de
rangos para monotonía de
preferencias
Una persona interesada en adquirir un TV asigna
rangos a modelos de cada uno de 8 fabricantes
Rango
Fab. Preferencia Precio
Di
(rango) Di cuadrada
6 36
1 7 449.50 (1)
-1 1
2 4 525.00 (5)
-1 1
3 2 479.95 (3)
4 6 499.95 (4) 2 4
5 1 580.00 (8) -7 49
6 3 549.95 (7) -4 16
7 8 469.95 (2) 6 36
8 5 532.50 (6) -1 1 461
Coeficiente de correlación de
rangos para monotonía de
preferencias
Ho: No existe asociación entre los rangos
Ha: Existe asociación entre los rangos o es positiva o negativa
Rt = 0.686
462
Tabla de constantes
n Alfa=0.05 Alfa = 0.025
5 0.900 -
6 0.829 0.886
7 0.714 0.786
8 0.643 0.738
9 0.600 0.683
10 0.564 0.648
11 0.523 0.623
12 0.497 0.591
13 0.475 0.566
14 0.457 0.545
15 0.441 0.525
16 0.425 0.507
17 0.412 0.490
18 0.388 0.476
19 0.377 0.462
20 0.368 0.450
21 0.359 0.438
22 0.351 0.428
23 0.343 0.418
24 0.336 0.409
25 0.329 0.400
26 0.329 0.392
27 0.323 0.385
28 0.317 0.377
29 0.311 0.370
30 0.305 0.364 463
Corrida con Minitab
Para la corrida en Minitab primero se Fabric Prefe- Preci
deben determinar los rangos en ante rencia Precio o
forma manual para las variables X
y Y. 1 7 1 449
Stat > Basic statistics > Correlation 2 4 5 525
Variables Preferencia Precio 3 2 3 479
Correlations: Preferencia, Precio 4 6 4 499
Pearson correlation of 5 1 8 580
Preferencia and Precio = -0.714 6 3 7 549
P-Value = 0.047 7 8 2 469
8 5 6 532
464
Ejemplo con Minitab
Se estudia la relación entre colágeno y
Proline en pacientes con cirrosis Paciente Colágeno Proline
Stat > Basic statistics > Correlation
1 7.1 2.8
Variables Colágeno Proline
2 7.1 2.9
3 7.2 2.8
Correlations: Colageno, Proline
4 8.3 2.6
Pearson correlation of Colageno 5 9.4 3.5
and Proline = 0.935
6 10.5 4.6
P-Value = 0.002 7 11.4 5
465
Resumen de pruebas
no paramétricas
Prueba de signos de 1 muestra: Prueba la igualdad de la
mediana a un valor y determina el intervalo de confianza
466
Resumen de pruebas
no paramétricas
467
468
469
VI.C Análisis del Modo y
Efecto de Falla (FMEA)
470
¿ Qué es el FMEA?
El Análisis de del Modo y Efectos de Falla es un grupo sistematizado de
actividades para:
472
Propósitos del FMEA
Soporta el desarrollo de planes de verificación del
desarrollo de diseño robusto
473
Tipos del FMEA
AMEF de concepto (CFMEA)
A nivel de sistema, subsistema y componente
474
Tipos de FMEAs
FMEA de Diseño (AMEFD), su propósito es analizar como
afectan al sistema los modos de falla y minimizar los
efectos de falla en el sistema. Se usan antes de la
liberación de productos o servicios, para corregir las
deficiencias de diseño.
475
PFMEA o AMEF de Proceso
Fecha límite:
Concepto Prototipo Pre-producción /Producción
FMEAD
FMEAP
FMEAD FMEAP
Característica de Diseño Paso de Proceso
Falla Forma en que el Forma en que el proceso falla
producto o servicio falla al producir el requerimiento
que se pretende
Controles Técnicas de Diseño de Controles de Proceso
Verificación/Validación
476
Flujo del FMEA y su rol
en evitar el Modo de Falla
DFMEA
Es un análisis detallado de los modos de falla
477
Beneficios de los tipos de FMEA
FMEA de Diseño
Soporta el proceso de diseño al reducir el riesgo de fallas
(incluyendo las salidas no intencionadas) por:
480
Beneficios de los tipos de FMEA
FMEA de Proceso
Los beneficios de un FMEA de proceso incluyen:
485
Salidas del FMEA de Proceso
Una lista de procesos o acciones de proceso para
reducir la Severidad, eliminar las causas de los modos
de falla del producto o reducir su tasa de ocurrencia, y
mejorar la tasa de Detección de defectos si no se
puede mejorar la capacidad del proceso
486
Modos de fallas vs
Mecanismos de falla
487
Definiciones
Modo
Modode
deFalla
Falla
--La
Laforma
formaenenque
queun unproducto
productoooproceso
procesopuede
puedefallar
fallarpara
paracumplir
cumplir
con
conlas
lasespecificaciones
especificacionesoorequerimientos.
requerimientos.
--Normalmente
Normalmentese
seasocia
asociacon unDefecto,
conun Defecto,falla
fallaooerror.
error.
Diseño
Diseño Proceso
Proceso
Alcance
Alcanceinsuficiente
insuficiente Omisiones
Omisiones
Recursos
Recursosinadecuados
inadecuados Monto
Montoequivocado
equivocado
Servicio
Serviciono
noadecuado
adecuado Tiempo
Tiempodederespuesta
respuestaexcesivo
excesivo
488
Definiciones
Efecto
Efecto
--ElElimpacto enelelCliente
impactoen Clientecuando
cuandoelelModo
Modode
deFalla
Fallano
nose
sepreviene
previene
ninicorrige.
corrige.
--ElElcliente
clienteooelelsiguiente
siguienteproceso
procesopuede
puedeser
serafectado.
afectado.
Ejemplos:
Ejemplos: Diseño
Diseño Proceso
Proceso
Serv.
Serv.incompleto
incompleto Servicio
Serviciodeficiente
deficiente
Operación
Operaciónerrática
errática Claridad
Claridadinsuficiente
insuficiente
Causa
Causa
--Una
Unadeficiencia
deficienciaque
quegenera
generaelelModo
Modode
deFalla.
Falla.
--Las
Lascausas
causasson
sonfuentes deVariabilidad
fuentesde Variabilidadasociada
asociadaconconvariables
variables
de Entrada Claves
de Entrada Claves
Ejemplos:
Ejemplos: Diseño
Diseño Proceso
Proceso
Material
Materialincorrecto
incorrecto Error
Erroren
enservicio
servicio
Demasiado
Demasiadoesfuerzo
esfuerzo No
No cumple
cumplerequerimientos
requerimientos
489
Preparación
Preparación del
del AMEF
AMEF
Se recomienda que sea un equipo
multidisciplinario
490
¿Cuando iniciar un FMEA?
Al diseñar los sistemas, productos y procesos nuevos.
Al cambiar los diseños o procesos existentes o que serán
usados en aplicaciones o ambientes nuevos.
492
AMEF de Diseño
493
AMEF de Diseño
El proceso inicia con un listado de lo que se espera
del diseño (intención) y que no hará el diseño
495
Alcance del DMEA
El alcance se establece en el Diagrama de límites
(Boundary Diagram) por medio de consenso con el
equipo de:
496
Alcance del DMEA
Para determinar la amplitud del alcance, se deben hacer las
decisiones siguientes:
497
Entradas al DFMEA
Herramientas de robustez
Su propósito es reducir la probabilidad de campañas
de calidad, mejorar la imagen, reducir reclamaciones
de calidad e incrementar la satisfacción del cliente
498
Modelo DFMEA – Paso 1
Funciones
499
Modelo DFMEA – Paso 1
Funciones
500
Modelo DFMEA – Paso 1
Funciones
Ejemplos:
Almacenar fluido, X litros sin fugas
de operación
501
Modelo DFMEA – Paso 1
Modos de falla potenciales
502
Modelo DFMEA – Paso 1
Tipos de Modos de falla potenciales
No funciona
Funciona parcialmente / sobre función / degradación
con el tiempo
Función intermitente
A veces causado por los factores ambientales
Función no intencionada
Los limpiadores operan sin haber actuado el switch
El coche va hacia atrás aún con la palanca en Drive
503
Modelo DFMEA – Paso 1
Preguntas para Modos Potenciales de falla
504
Modelo DFMEA – Paso 1
Preguntas para Modos Potenciales de falla
Asumiendo la función:
Almacenar fluido, X litros, 0 fugas, durante 10 años
505
Modelo DFMEA – Paso 1
Efectos Potenciales de falla
506
Modelo DFMEA – Paso 1
Efectos Potenciales de falla
Sistemas
Clientes
Reglamentaciones
507
Modelo DFMEA – Paso 1
Ejemplos de Efectos Potenciales de falla
Ruidos
Operación inestable
Operación intermitente
Fugas
509
Rangos de Severidad (AMEFD)
Efecto Rango Criterio .
No 1 Sin efecto
Muy poco 2 Cliente no molesto. Poco efecto en el desempeño del componente o
servicio.
Poco 3 Cliente algo molesto. Poco efecto en el desempeño del comp. o
servicio.
Menor 4 El cliente se siente un poco fastidiado. Efecto menor en el
desempeño del componente o servicio.
Moderado 5 El cliente se siente algo insatisfecho. Efecto moderado en el
desempeño del componente o servicio.
511
Modelo DFMEA – Paso 1
Acciones recomendadas
Mitigar el efecto
512
Modelo DFMEA – Paso 2
Identificar:
Las Causas asociadas (primer nivel y raíz)
513
Model DFMEA – Paso 2
Causa potencial o mecanismo de falla
514
Model DFMEA – Paso 2
Causa potencial o mecanismo de falla
515
Model DFMEA – Paso 2
Causa potencial o mecanismo de falla
Calor excesivo
516
Model DFMEA – Paso 2
Causa potencial o mecanismo de falla
517
Modelo DFMEA – Paso 2
Causa potencial o mecanismo de falla
Precauciones:
El DFMA no confía en los controles del proceso para subsanar
debilidades del diseño, pero toma en cuenta sus limitaciones
518
Modelo DFMEA – Paso 2
Ocurrencia
519
Modelo DFMEA – Paso 2
Estimación de la Ocurrencia
520
Rangos de Ocurrencia (AMEFD)
Ocurrencia Criterios Rango Probabilidad de Falla
Remota Falla improbable. No existen fallas 1 <1 en 1,500,000 Zlt > 5
asociadas con este producto o con
un producto / Servicio casi idéntico
Muy Poca Sólo fallas aisladas asociadas con 2 1 en 150,000 Zlt > 4.5
este producto / Servicio
casi idéntico
3 1 en 30,000
Poca Fallas aisladas asociadas con
Zlt > 4
productos / Servicios similares
Moderada Este producto / Servicio ha
tenido fallas ocasionales 4 1 en 4,500 Zlt > 3.5
5 1 en 800 Zlt > 3
Alta Este producto / Servicio ha
6 1 en 150 Zlt > 2.5
fallado a menudo
7 1 en 50 Zlt > 2
Muy alta La falla es casi inevitable
8 1 en 15 Zlt > 1.5
9 1 en 6 Zlt > 1
10 >1 en 3 Zlt < 1
522
Modelo DFMEA
Paso 3
Si las causas no se pueden eliminar en paso 1 o 2, Identificar
Controles actuales de prevención usados para establecer la
ocurrencia
Controles actuales de detección (vg. Pruebas) usadas para
establecer la Detección
Determinar la efectividad de los controles de Detección en
escala de 1 a 10
El RPN inicial (Risk Priority Number).
Acciones Recomendadas (Prevenciónn and Detección).
Cuando ya se hayan implementado las acciones recomenddas,
se revisa el formato DFMEA en relación a la Severidad,
Ocurrencia, Detección y RPN
523
Modelo DFMEA – Paso 3
Controles de diseño actuales
Listar las actividades terminadas para prevención,
vaidación/verificación del diseño (DV), u otras actividades que
aseguran la adecuación del diseño para el modo de falla y/o
causa / mecanismo bajo consideración
524
Modelo DFMEA – Paso 3
Controles de diseño actuales
Hay dos tipos de controles de diseño: Prevención y
detección
De prevención:
Previenen la ocurrencia de la causa/mecanismo o Modo
de falla/efecto reduciendo la tasa de Ocurrencia
De detección:
Detectan la causa/mecanismo o Modo de falla/efecto
ya sea por métodos analíticos o físicos antes que el
artículo se libere para Poducción
Si solo se usa una columna indicarlos con P o D
525
Modelo DFMEA – Paso 3
Controles de diseño actuales
Identificación de controles de diseño
Si una causa potencial no fue analizada, el producto
con deficiencia de diseño pasará a Producción. Una
forma de detectarlo es con su Modo de falla
resultante. Se debe tomar acción correctiva
finitos)
Estudios de tolerancias (vg. Tolerancias deométricas
dimensionales)
Estudios de compatibilidad de materiales (vg. Expansión
térmica, corrosión)
Revisión de diseño subjetiva
Ruido)
529
Modelo DFMEA – Paso 3
Detección
Métodos de desarrollo de pruebas (cont…):
Pruebas en muestras de pre-producción o prototipo
Maquetas usando partes similares
532
DFMEA – Acciones recomendadas
Considerar acciones como las siguientes:
Revisión del diseño de la Geometría y/o tolerancias
Revisión de especificación de materiales
Diseños de experimentos (con múltiples causas interactuando)
u otras técnicas de solución de problemas
Revisión de planes de prueba
Sistemas redundantes – dispositivos de aviso – estados de falla
(ON y OFF)
Dar seguimiento:
Desarrollar una lista de características especiales parasu
consideración en el DFMEA
Dar seguimiento a todas las acciones recomendadas y
534
DFMEA – Nivel de riesgo RPN
Después de haber implementado las acciones
preventivas/correctivas, registrar la nueva Severidad,
Ocurrencia y Detección
535
DFMEA – Lista de verificación de
robustez
Es una salida del proceso integrado de robustez:
Resume los atributos de robustez clave y controles de
diseño
536
FMEA de Proceso - PFMEA
537
PFMEA
Equipo
Se inicia por el Ing. responsable de la actividad, en
conjunto con un equipo de personas expertas además
de incluir personas de apoyo
Alcance
Define que es incluido y que es excluido
538
Entradas al PFMEA
Diagrama de flujo del proceso
El equipo debe desarrollar el flujo del proceso,
preguntando ¿Qué se supone que hace el proceso?;
¿Cuál es su propósito?; ¿Cuál es su función?
539
ANALISIS DEL MODO Y EFECTO DE FALLA
AMEF de Diseño / Proceso
Componente ______________________ Responsable del Diseño ____________AMEF Número _________________
Ensamble ________________ Preparó _______________ Pagina _______de _______
Equipo de Trabajo ___________ FECHA (orig.) de FMEA ______(rev.) ______
Resultados de Acción
O Controles de D
Función S Causa(s)
Efecto (s) c Diseño o e R Responsable S O D R
del Producto/ Modos de Falla e Potencial(es) Acción Acción
Potencial (es) c Proceso t P y fecha límite e c e P
Paso del Potenciales v o Mecanismos Sugerida Adoptada
de falla u Actuales e N de Terminación v c t N
proceso . de falla
r c
540
Modelo del PFMEA – Paso 1
Identificar todos los requerimientos funcionales dentro del
alcance
541
Modelo de PFMEA – Paso 1
Ejemplos
Operación No. 20: Hacer perforación de tamaño X de
cierta profundidad
Operación No. 22: Realizar el subensamble X al
ensamble Y
542
ANALISIS DEL MODO Y EFECTO DE FALLA
AMEF de Diseño / Proceso
Componente ______________________ Responsable del Diseño ____________AMEF Número _________________
Ensamble ________________ Preparó _______________ Pagina _______de _______
Equipo de Trabajo ___________ FECHA (orig.) de AMEF ______(rev.) ______
Resultados de Acción
O D
Función S Causa(s) Controles del
Efecto (s) c e R Responsable S O D R
de Modos de Falla e Potencial(es) Diseño / Acción Acción
Potencial (es) c t P y fecha límite e c e P
Componente/Paso Potenciales v de los Mecanismos Proceso Sugerida Adoptada
de falla u e N de Terminación v c t N
de proceso . de falla Actual
r c
Factura correcta
Relacione las
funciones del
diseño del
componente
tiempo
Funcionamiento intermitente
Función no intencionada
544
ANALISIS DEL MODO Y EFECTO DE FALLA
AMEF de Diseño / Proceso
Componente ______________________ Responsable del Diseño ____________AMEF Número _________________
Ensamble ________________ Preparó _______________ Pagina _______de _______
Equipo de Trabajo ___________ FECHA (orig.) de AMEF ______(rev.) ______
Resultados de Acción
Función O D
Causa(s) Controles de
del Efecto (s) D c e R Responsable S O D R
Modos de Falla Potencial(es) Diseño / Acción Acción
componente/ Potencial (es) i c t P y fecha límite e c e P
Potenciales de los Mecanismos Proceso Sugerida Adoptada
Paso del de falla v u e N de Terminación v c t N
de falla Actuales
proceso r c
Factura
correcta
Datos incorrectos Identificar modos
de falla Tipo 1
inherentes al
diseño
Modelo de PFMEA – Paso 1
546
Modelo de PFMEA - Paso 1
547
Modelo de PFMEA –Paso 1
548
ANALISIS DEL MODO Y EFECTO DE FALLA
AMEF de Diseño
Componente ______________________ Responsable del Diseño ____________AMEF Número _________________
Ensamble ________________ Preparó _______________ Pagina _______de _______
Equipo de Trabajo ___________ FECHA (orig.) de AMEF ______(rev.) ______
Resultados de Acción
O Controles de D
Función Causa(s)
Efecto (s) D c Diseño / e R Responsable S O D R
del componente Modos de Falla Potencial(es) Acción Acción
Potencial (es) i c Proceso t P y fecha límite e c e P
/ Paso del Potenciales oMecanismos Sugerida Adoptada
de falla v u Actuales e N de Terminación v c t N
proceso de falla
r c
Esta calificación resulta cuando un modo de falla potencial resulta en un defecto con un cliente final y/o una planta de manufactura
/ ensamble. El cliente final debe ser siempre considerado primero. Si ocurren ambos, use la mayor de las dos severidades
Peligroso Calificación de severidad muy alta cuando un modo potencial de Puede exponer al peligro al operador (máquina o ensamble)
con aviso falla afecta la operación segura del producto y/o involucra un no sin aviso 9
cumplimiento con alguna regulación gubernamental, con aviso
Muy alto El producto / item es inoperable ( pérdida de la función El 100% del producto puede tener que ser desechado op
primaria) reparado con un tiempo o costo infinitamente mayor 8
Alto El producto / item es operable pero con un reducido nivel de El producto tiene que ser seleccionado y un parte desechada
desempeño. Cliente muy insatisfecho o reparada en un tiempo y costo muy alto 7
Modera Producto / item operable, pero un item de confort/conveniencia Una parte del producto puede tener que ser desechado sin
do es inoperable. Cliente insatisfecho selección o reparado con un tiempo y costo alto 6
Bajo Producto / item operable, pero un item de confort/conveniencia El 100% del producto puede tener que ser retrabajado o
son operables a niveles de desempeño bajos reparado fuera de línea pero no necesariamente va al àrea 5
de retrabajo .
Muy bajo No se cumple con el ajuste, acabado o presenta ruidos y El producto puede tener que ser seleccionado, sin desecho,
rechinidos. Defecto notado por el 75% de los clientes y una parte retrabajada 4
Menor No se cumple con el ajuste, acabado o presenta ruidos y El producto puede tener que ser retrabajada, sin desecho,
rechinidos. Defecto notado por el 50% de los clientes en línea, pero fuera de la estación 3
Muy No se cumple con el ajuste, acabado o presenta ruidos, y El producto puede tener que ser retrabajado, sin desecho en
menor rechinidos. Defecto notado por clientes muy críticos (menos del la línea, en la estación 2
25%)
Ninguno Sin efecto perceptible Ligero inconveniente para la operación u operador, o sin
efecto 1
CRITERIO DE EVALUACIÓN DE SEVERIDAD SUGERIDO PARA PFMEA
Esta calificación resulta cuando un modo de falla potencial resulta en un defecto con un cliente final y/o una planta de manufactura
/ ensamble. El cliente final debe ser siempre considerado primero. Si ocurren ambos, use la mayor de las dos severidades
Peligroso Calificación de severidad muy alta cuando un modo potencial de Puede exponer al peligro al operador (máquina o ensamble)
con aviso falla afecta la operación segura del producto y/o involucra un no sin aviso 9
cumplimiento con alguna regulación gubernamental, con aviso
Muy alto El producto / item es inoperable ( pérdida de la función El 100% del producto puede tener que ser desechado op
primaria) reparado con un tiempo o costo infinitamente mayor 8
Alto El producto / item es operable pero con un reducido nivel de El producto tiene que ser seleccionado y un parte desechada
desempeño. Cliente muy insatisfecho o reparada en un tiempo y costo muy alto 7
Modera Producto / item operable, pero un item de confort/conveniencia Una parte del producto puede tener que ser desechado sin
do es inoperable. Cliente insatisfecho selección o reparado con un tiempo y costo alto 6
Bajo Producto / item operable, pero un item de confort/conveniencia El 100% del producto puede tener que ser retrabajado o
son operables a niveles de desempeño bajos reparado fuera de línea pero no necesariamente va al àrea 5
de retrabajo .
Muy bajo No se cumple con el ajuste, acabado o presenta ruidos y El producto puede tener que ser seleccionado, sin desecho,
rechinidos. Defecto notado por el 75% de los clientes y una parte retrabajada 4
Menor No se cumple con el ajuste, acabado o presenta ruidos y El producto puede tener que ser retrabajada, sin desecho,
rechinidos. Defecto notado por el 50% de los clientes en línea, pero fuera de la estación 3
Muy No se cumple con el ajuste, acabado o presenta ruidos, y El producto puede tener que ser retrabajado, sin desecho en
menor rechinidos. Defecto notado por clientes muy críticos (menos del la línea, en la estación 2
25%)
Ninguno Sin efecto perceptible Ligero inconveniente para la operación u operador, o sin
efecto 1
Modelo de PFMEA – Paso 2
Paso 2 identificar:
Las causas asociadas (primer nivel y raíz)
Su tasa de ocurrencia
552
Modelo de PFMEA – Paso 2
553
Efecto(s) Potencial(es) de falla
Evaluar 3 (tres) niveles de Efectos del Modo de Falla
• Efectos Locales
– Impactos Inmediatos
• Efectos Finales
– Efecto en el Usuario Final del producto o Servicio
554
Modelo de PFMEA – Paso 2
555
Modelo de PFMEA – Paso 2
556
Modelo de PFMEA – Paso 2
Ocurrencia:
Es la probabilidad de que una causa/mecanismo ocurra
Se puede reducir o controlar solo a través de un
cambio de diseño
Si la ocurrencia de la causa no puede ser estimada,
entonces estimar la tasa de falla posible
557
CRITERIO DE EVALUACIÓN DE OCURRENCIA SUGERIDO PARA AMEFP
Probabilidad Indices Posibles de ppk Calif.
falla
Muy alta: Fallas 100 por mil piezas < 0.55 10
persistentes
50 por mil piezas > 0.55 9
559
ANALISIS DEL MODO Y EFECTO DE FALLA
AMEF de Diseño / Proceso
Componente ______________________ Responsable del Diseño ____________AMEF Número _________________
Ensamble ________________ Preparó _______________ Pagina _______de _______
Equipo de Trabajo ___________ FECHA (orig.) de AMEF ______(rev.) ______
Resultados de Acción
O D
Función S Causa(s) Controles de
Efecto (s) c e R Responsable S O D R
del componente Modos de Falla e Potencial(es) Diseño / Acción Acción
Potencial (es) c t P y fecha límite e c e P
/ Paso del Potenciales v o Mecanismos Proceso Sugerida Adoptada
de falla u e N de Terminación v c t N
proceso . de falla Actuales
r c
La abertura del
engrane propor La abertura no LOCAL:
ciona una aber- es suficiente Daño a sensor
tura de aire entre de velocidad y
diente y diente engrane
MAXIMO PROXIMO
Usar tabla para
Falla en eje 7 determinar severidad o
gravedad
CON CLIENTE
Equipo
parado
ANALISIS DEL MODO Y EFECTO DE FALLA
AMEF de Diseño / Proceso
Componente ______________________ Responsable del Diseño ____________AMEF Número _________________
Ensamble ________________ Preparó _______________ Pagina _______de _______
Equipo de Trabajo ___________ FECHA (orig.) de AMEF ______(rev.) ______
Resultados de Acción
Función O D
S Causa(s) Controles de
del Efecto (s) c e R Responsable S O D R
Modos de Falla e Potencial(es) Diseño/ Acción Acción
Componente / Potencial (es) c t P y fecha límite e c e P
Potenciales v o Mecanismos Proceso Sugerida Adoptada
Paso del de falla u e N de Terminación v c t N
. de falla Actuales
proceso r c
CON CLIENTE
Molestia
Insatisfacción
Modelo de PFMEA – Paso 3
En el paso 3 identificar:
Controles actuales de prevención del proceso (con
acciones de diseño o proceso) usados para establecer
la ocurrencia
Controles actuales de detección (vg. Inspección)
usados para establecer la tasa de detección
Efectividad de los controles de detección del proceso
en una escala de 1 a 10
El factor de riesgo RPN inicial
Acciones recomendadas (Prevención y Detección)
562
Identificar Causa(s) Potencial(es) de la Falla
– Asignación de recursos
– Equipos planeados
• Mecanismos de Falla
– Rendimiento, tiempo de entrega, información completa
563
ANALISIS DEL MODO Y EFECTO DE FALLA
AMEF de Diseño
Componente ______________________ Responsable del Diseño ____________AMEF Número _________________
Ensamble ________________ Preparó _______________ Pagina _______de _______
Equipo de Trabajo ___________ FECHA (orig.) de AMEF ______(rev.) ______
Resultados de Acción
O D
S Causa(s) Controles de
Función Efecto (s) c e R Responsable S O D R
Modos de Falla e Potencial(es) Diseño/Proces Acción Acción
de Potencial (es) c t P y fecha límite e c e P
Potenciales v de los Mecanismos o Actuales Sugerida Adoptada
Artículo de falla u e N de Terminación v c t N
. de falla
r c
Consideraciones
Incrementar la probabilidad de detección es costosa y
no efectiva
A veces se requiere un cambio en el diseño para
apoyar la detección
El incremento del control de calidad o frecuencia de
inspección sólo debe utilizarse como medida temporal
Se debe hacer énfasis en la prevención de los defectos
565
Identificar Controles de Diseño o de
Proceso Actuales
• Verificación/ Validación de actividades de Diseño o
control de proceso usadas para evitar la causa,
detectar falla anticipadamente, y/o reducir impacto:
Resultados de Acción
Función O D
S Causa(s) Controles de
del Efecto (s) c e R Responsable S O D R
Modos de Falla e Potencial(es) Diseño / Acción Acción
Componente / Potencial (es) c t P y fecha límite e c e P
Potenciales v o Mecanismos Proceso Sugerida Adoptada
Paso del de falla u e N de Terminación v c t N
. de falla Actuales
proceso r c
CON CLIENTE
Molestia
Insatisfacción
Modelo de PFMEA – Paso 3
Moderada Los controles tienen una buena X X Detección de error en operaciones subsiguientes, o medición
Muy Alta Controles seguros para detectar X No se pueden hacer partes discrepantes porque el item ha
pasado a prueba de errores dado el diseño del 1
proceso/producto
Resultados de Acción
Función O D
S Causa(s) Controles de
del Efecto (s) c e R Responsable S O D R
Modos de Falla e Potencial(es) Diseño / Acción Acción
Componente / Potencial (es) c t P y fecha límite e c e P
Potenciales v o Mecanismos Proceso Sugerida Adoptada
Paso del de falla u e N de Terminación v c t N
. de falla Actuales
proceso r c
¿Cuál es la probabilidad
MAXIMO PROXIMO
de detectar la causa de
Contabilidad 7 3 5
falla?
erronea
CON CLIENTE
Molestia
Insatisfacción
Modelo de PFMEA – Paso 3
Acciones recomendadas
Se deben dirigir primero a las de valores altos de
Severidad (9 o 10) o RPNs, después continuar con las
demás
Las acciones se deben orientar a prevenir los defectos
a través de la eliminación o reducción de las causas o
modos de falla
571
Calcular RPN (Número de Prioridad de
Riesgo)
572
Planear Acciones
Acciones tomadas
Identificar al responsable de las acciones
recomendadas y la fecha estimada de terminación
Después de terminar una acción, dar una descripción
breve de la acción real y fecha de efectividad
574
Modelo de PFMEA – Paso 3
RPN resultante
Después de implementadas las acciones
recomendadas, estimar de nuevo los rangos de
Severidad, Ocurrencia y Detección y calcular el nuevo
RPN. Si no se tomaron acciones dejarlo en blanco.
575
ANALISIS DEL MODO Y EFECTO DE FALLA
AMEF de Diseño / Proceso
Componente ______________________ Responsable del Diseño ____________ AMEF Número _________________
Ensamble ________________ Preparó _______________ Pagina _______de _______
Equipo de Trabajo ___________ FECHA (orig.) de AMEF ______(rev.) ______
Resultados de Acción
O D
S Causa(s)
Función Efecto (s) c e R Responsable S O D R
Modos de Falla e Potencial(es) Controles de Acción Acción
de Potencial (es) c t P y fecha límite e c e P
Potenciales v de los Mecanismos Diseño Actual Sugerida Adoptada
Artículo de falla u e N de Terminación v c t N
. de falla
r c
Riesgo = Severidad x
MAXIMO PROXIMO Ocurrencia x Detección
Contabilidad 7 3 5 105
erronea
CON CLIENTE
Molestia
Insatisfacción
Causas probables a
atacar primero
ANALISIS DEL MODO Y EFECTO DE FALLA
AMEF de Diseño / Proceso
Componente ______________________ Responsable del Diseño ____________ AMEF Número _________________
Ensamble ________________ Preparó _______________ Pagina _______de _______
Equipo de Trabajo ___________ FECHA (orig.) de AMEF ______(rev.) ______
Resultados de Acción
O D
Función S Causa(s) Controles de
Efecto (s) c e R Responsable S O D R
del componente Modos de Falla e Potencial(es) Diseño / Acción Acción
Potencial (es) c t P y fecha límite e c e P
/ Paso del Potenciales v o Mecanismos Prcoeso Sugerida Adoptada
de falla u e N de Terminación v c t N
proceso . de falla Actuales
r c
MAXIMO PROXIMO
Contabilidad 7 3 5 105
erronea
CON CLIENTE
Molestia
Insatisfacción Usar RPN para identificar
acciones futuras. Una vez que
se lleva a cabo la acción,
recalcular el RPN.
Herramientas para el FMEA
578
Herramientas
Diagramas de límites
Diagramas de flujo de proceso
Matriz de características
Tormenta de ideas
Árboles de funciones
Lista de efectos: FMEA de diseño
Lista de efectos: FMEA de proceso
Diagrama de Ishikawa
Tecnica de preguntas
579
Herramientas
Análisis de árbol de fallas (FTA)
Análisis del modo de falla (FMA)
Diseño de experimentos (DOE)
Proceso de solución de problemas de 8Ds
Planes de Control
Planeación dinámica de control (DCP)
Despliegue de la función de calidad (QFD)
Análisis de valor/ Ingeniería del valor (VA/VE)
REDPEPR
FMEA Express
FMEA del software
580
Diagrama de límites
Diagramas de límites de funciones
Salida del análisis de funciones para la fase de
concepto CFMEA, ilustran funciones en vez de partes
581
Nombres de verbos útiles
582
Tormenta de ideas
583
Herramientas para el FMEA
Árbol de funciones
Ayuda a que los requerimientos del cliente no
expresados explícitamente sobre el producto o proceso
se cumplan
584
Técnica de preguntas
Hacer una oración con el modo de falla, causa y efecto y ver si
la oración tiene sentido. Un modo de falla es debido a una
causa, el modo de falla podría resultar en efectos, por ejemplo:
MODO DE FALLA: No ajustan los faros delanteros
P: ¿Qué podría ocasionar esta falla?
585
Técnica de preguntas
Paso 3 Paso 1
¿Qué lo causa? Paso 2
Modo de falla ¿Qué efecto tiene?
586
Análisis de árbol de fallas (FTA)
Es una técnica analítica deductiva que usa un árbol para
mostrar las relaciones causa efecto entre un evento indeseable
(falla) y las diversas causas que contribuyen. Se usan símbolos
lógicos para interconectar las ramas
587
Análisis del Modo de Falla (FMA)
Es un enfoque sistemático disciplinado para
cuantificar el modo de falla, tasa de falla, y causa
raíz de fallas o tasas de reparación conocidas (el
FMEA para las desconocidas)
3. Bitácora de preguntas
592
Planeación dinámica de control (DCP)
4. Información de soporte (ES, DFMEAs, DVP&R, PFMEA, etc.)
7. PFMEA
8. Plan de control
593
Despliegue de la función de calidad
(QFD)
El QFD es un método estructurado en el cual los
requerimientos del cliente son traducidos en requerimientos
técnicos para cada una de las etapas del desarrollo del
producto y producción
594
Análisis del valor / Ingeniería del
valor (VA/VE)
Son metodologías usadas comúnmente para despliegue del valor. La
Ingeniería del valor se realiza antes de comprometer el herramental. El
análisis del valor (VA) se realiza después del herramentado. Ambas
técnicas usan la fórmula:
595
REDPEPR (Robust Engineering
Design Product Enhacement Process)
Es una herramienta que proporciona a los equipos de
Diseño:
Un proceso paso a paso para aplicar el RED
Las herramientas necesarias para completar el diagrama P, listas
de verificación de confiabilidad y robustez (RRCL) y la matriz de
demostración de confiabilidad y robustez (RRDM)
Preguntas y tips para guiar al equipo en el proceso
Capacidad para generar reportes en Excel
Un proceso para mejorar la comunicación con el equipo de
ingeneiría
El Web site donde se encuentra el software es www.redpepr.ford.com
596
Aplicaciones del FMEA
Express
Ambiental
De máquinas
597
FMEA Express
Es un proceso que aplica técnicas de FMEA
simultaneamente tanto a los aspectos de diseño
como a los de manufactura de un proyecto:
598
FMEA Express
Desarrollo del FMEA: El equipo de trabajo multidisciplianrio
completa el FMEA utilizando formatos y definiciones estándar
599
E-FMEA ambiental
600
E-FMEA ambiental
601
Matriz de requerimientos ambientales
con criterios múltiples
Para cada alternativa de diseño resumir la siguiente
información
Uso de substancias prohibidas o de uso restringido
Tipo y cantidad de residuos (refleja el nivel de materiales
utilizados)
Consumo de energía por componente
Consumo de agua por componente
Otros objetivos ambientales
602
E-FMEA
Ejemplos de acciones recomendadas (hacer una revisión
previa de efectos secundarios en la vida del
producto):
Sistemas de conexión alternos
Reciclar
603
E-FMEA
Salidas del FMEA ambiental:
Recomendaciones de materiales
604
MFMEA – FMEA de maquinaria
605
FMEA de maquinaria
Su propósito es que a través de un equipo se asegure que los
modos de falla y sus causas/mecanismos asociados se hayan
atendido
606
FMEA de maquinaria
607
FMEA de maquinaria
Ejemplos de descripción de funciones
Proceso de partes – 120 tareas / hora
608
FMEA de maquinaria
Efectos potenciales como consecuencias de falla de subsistemas
en relación a seguridad y “Las 7 grandes pérdidas”
Falla – pérdidas resultado de una pérdida funcional o
609
FMEA de maquinaria
Tiempo de espera y paros menores – pérdidas
resultado de interrupciones menores al flujo del proceso
(como atoramiento de microswitch) requiriendo
intervención del operador. El tiempo de espera sólo se
puede resolver revisando el sistema / línea completa
612
FMEA de maquinaria
Causas potenciales, se asume que la maquinaria se
fabricó, instaló, usó, y se dispuso de acuerdo a sus
especificaciones, preguntarse para identificar causas
potenciales lo siguiente:
¿Cuáles son las circunstancias que pueden orientar al
componente, subsistema y sistema a no cumplir sus
requerimientos funcionales / de desempeño?
614
FMEA de maquinaria
Criterios de Detección
615
Herramientas de la
Fase de Análisis
616
Identificación de causas
potenciales
Tormenta de ideas
Diagrama de Ishikawa
Diagrama de Relaciones
Diagrama de Árbol
Verificación de causas raíz
617
Tormenta de ideas
Técnica para generar ideas creativas cuando la mejor
solución no es obvia.
618
Tormenta de ideas
Permite obtener ideas de los participantes
619
Diagrama de Ishikawa
Anotar el problema en el cuadro de la derecha
Mediciones
Materia Prima
Maquinaria
Personal y
Métodos
o
Las diferentes etapas del proceso de manufactura o
servicio
620
Diagrama de Ishikawa
Medio
ambiente Métodos Personal
Frecuencia Falta de
Rotación de
Clima personal
de visitas supervi
húmedo Falta de
ción
motivación
Posición de Ausentismo
Distancia de exhibidores
la agencia al
Elaboración ¿Qué
changarro de pedidos produce
bajas ventas
Clientes con Calidad del de
ventas bajas Seguimiento producto Tortillinas
Malos
semanal Tía Rosa?
Conocimiento
itinerarios
de los Tipo de
Descompostura mínimos por exhibidor
del camión ruta
repartidor
13/0 2/1
Bandas de
Dancer transmisión
2/4 Taco generador 1/1 Causas a validar
del motor Empaques de arrastre
No strong
commitment
to the group
In = 2 Out = 0
Capacity
may not
meet needs
Planning In = 5 Out = 1
approach not
standardized Outcome
In = 0 Out = 5
External Lack of
Driver factors impact time and
implemen- resources
tation
In = 5 Out = 0
In = 0 Out = 2
Diagrama de árbol o sistemático
Meta Medio
Meta Medio
Meta Medio
Segundo Tercer Cuarto
Primer nivel nivel nivel
nivel
Medios
Medios
Medios
Medios
o planes
Meta u
objetivo
Medios
o planes
625
Diagrama de Arbol- Aplicación Sistema SMED
¿Cómo? ¿Cuándo?
Filmar la preparación 5- 12 - Mar-04
Preparación
para el SMED
Analizar el video 10 y 17 –Mar-04
627
Calendario de las actividades
4 Poleas de 4.1 Puede generar 4.1.1 Verificar alineación, entre poleas de Abril’04 1804 J. R.
transmisión de vibración excesiva ejes principales y polea de transmisión del Embob. U. P.
ejes durante el ciclo de motor.
embobinadores cambio. 4.1.2 Tomar dimensiones de poleas(dientes
. de transmisión).
4.1.3 Tomar dimensiones de bandas (dientes
de transmisión)
4.1.4 Verificar valor de tensión de bandas. 628
Resumen de la validación de las causas
# de Causa
Resultados Causa
Causas Raíz
630
Métodos adicionales de análisis
1. Análisis de brecha
631
VI.D.1 Análisis de brecha
632
El análisis de brecha (Gap Analysis) es una
herramienta de evaluación para comparar el
desempeño actual de la organización, a un
desempeño potencial deseado.
633
Análisis de brecha
Se pueden redirigir los esfuerzos a objetivos como:
Permanecer en el negocio
Mantener o incrementar la participación del mercado
Mejorar el clima laboral
Igualar o exceder a Benchmarks
Igualar o exceder a la competencia
Reducir tiempos de ciclo
Lograr certificaciones
Mejorar la productividad
Mejorar los niveles de calidad
634
Análisis de brecha
Se requieren tres categorías de información
¿Dónde estamos?
¿Dónde queremos ir?
¿Cómo vamos a medir los resultados?
635
Planeación de escenarios
Al elaborar planes estratégicos, los directivos pueden
confiarse o ser orgullosos de aceptar cambios. Por lo
que se sugiere considerar escenarios del mejor y del
peor caso, para evitar errores en la toma de
decisiones
636
Planeación de escenarios
El proceso de planeación es como sigue:
Seleccionar al personal que pueda dar muchas
perspectivas
Desarrollar una lista de cambios percibidos, sociales,
técnicos y económicos
Agrupar estas percepciones en patrones relacionados
Desarrollar una lista de las mejores percepciones
(prioridades)
637
Planeación de escenarios
El proceso de planeación es como sigue:
Desarrollar un escenario grueso del futuro basado en
estas prioridades
Determinar como afectan los escenarios a la
organización
Determinar los cursos de acción potenciales a tomar
Monitorear, evaluar, y revisar los escenarios
638
Planeación de escenarios
Por lo común se perciben de 6 – 10 amenazas u
oportunidades en 2 o 3 escenarios desarrollados.
Evitar las siguientes trampas:
639
Planeación de escenarios
Evitar las siguientes trampas…..:
No incluir a un equipo directivo en el proceso
640
Planeación Hoshin
Es una herramienta de ejecución, usada para
organizar y desplegar planes estratégicos
641
Planeación Hoshin
Tiene seis objetivos:
Alinear las metas organizacionales
Análisis FODA
Análisis PEST
643
Evaluación organizacional
Análisis funcional con datos de colección:
Entrevistas cara a cara
Selección de muestra apropiada
Entradas de grupo de enfoque
Observaciones de visitas a la planta
Datos colectados de fuentes de la industria
644
Evaluación organizacional
Se deben analizar los resultados y presentarlos a la
dirección, quien debe promover e implementar
planes de acción claros
645
Métricas organizacionales
Se establecen metas de desempeño organizacional y
sus métricas en las áreas de:
Utilidades
Tiempos de ciclo
Recursos
Respuestas del mercado
646
Métricas organizacionales
Para los anteriores, las métricas pueden ser:
Utilidades a corto y largo plazo
Tiempos de ciclo
Tiempos de ciclo actuales
Benchmarks internos
Benchmarks externos
Reducción en tiempos de ciclo
647
Métricas organizacionales
Recursos
No. De proyectos de mejora, ROI de proyectos, estudios de
Análisis de devoluciones
Retención de clientes
650
Análisis de causa raíz
651
Análisis de causa raíz
655
5Ws y 1H
El método de las 5Ws y 1H se resume al preguntar
¿quién?, ¿qué?, ¿cuándo?, ¿dónde?, ¿por qué? Y
¿cómo?.
656
Diagrama de causa efecto
Rompe el problema en partes más pequeñas
Muestra muchas causas potenciales gráficamente
Muestra como interactúan las causas
Sigue las reglas de la tormenta de ideas
Las sesiones tienen tres partes:
Tormenta de ideas
Dar prioridades (identificar las tres causas principales)
Desarrollo de un plan de acción
657
Diagrama de Pareto
Sirve para identificar problemas u oportunidades
prioritarias o mayores
658
Método de las 8 disciplinas - Ford
El método de Ford para el análisis de causa raíz es:
D1. Establecer el equipo
D2. Describir el problema
D3. Desarrollar una acción de contención
D4. Identificar la causa raíz
D5. Desarrollar alternativas de solución
D6. Implementar una acción correctiva permanente
D7. Prevenir la recurrencia
D8. Reconocer al equipo y las contribuciones individuales
659
Análisis de árbol de falla - FTA
FTA es un método sistemático deductivo, para definir
un evento singular específico e indeseable, y
determinar todas las posibles razones (fallas) que
pueden hacer que ocurra el evento
660
Análisis de árbol de falla - FTA
Otras áreas de su aplicación son:
Análisis funcional de sistemas complejos
Evaluación de requerimientos de seguridad,
confiabilidad,
defectos de diseño,
riesgos de peligro,
acciones correctivas,
simplificación de mantenimiento y detección de falla,
eliminación lógica de causas de falla
661
Análisis de árbol de falla - FTA
Se prefiere el FTA en vez del FMEA cuando:
La seguridad el personal es importante
Se pueden identificar un número pequeño de eventos
superiores
Hay alto potencial de falla
El problema es cuantificar la evaluación del riesgo
La funcionalidad del producto es altamente compleja
El producto no es reaprables
662
Análisis de árbol de falla - FTA
Se prefiere el FMEA en vez del FTA cuando:
Los eventos superiores no se pueden definir
explícitamente
Son factibles múltiples perfiles potencialmente exitosos
La identificación de todos los modos de falla es
importante
La funcionalidad del producto tiene poca intervención
externa
663
Análisis de árbol de falla - FTA
Símbolos de compuertas lógicas para determinar la
confiabilidad del sistema. Hay símbolos de eventos y
símbolos de compuertas
Símbolos de eventos
664
Análisis de árbol de falla - FTA
Símbolos de compuertas lógicas
665
Análisis de árbol de falla - FTA
Ejemplo: se asume que falla el sistema superior
666
Análisis de árbol de falla - FTA
La probabilidad de falla del sistema es 5.02%. Se
indica que el teclado es prioritario (0.20), después la
CPU (0.015) y el monitor (0.015)
667
VI.D.3 Análisis del Muda
668
Análisis de Muda
Las actividades que no agregan valor se clasifican
como Muda, de acuerdo a Imai son:
Sobreproducción
Inventarios
Reparaciones / rechazos
Movimientos
Transportes
Re – Procesos
Esperas
669
Sobreproducción
Se produce más en cierto momento, por:
Producir más de lo necesario por el siguiente proceso
Producir antes de lo requerido por el siguiente proceso
Producir más rápido de lo requerido por el siguiente
proceso
Sus consecuencias son:
Espacio extra en las instalaciones del cliente
Materias primas adicionales en uso
Utilización de energéticos y transportes adicionales
Costos de programación adicionales
670
Inventario en exceso
Las partes, materias primas, inventario en proceso,
refacciones y productos terminados forman el
inventario, el inventario es Muda ya que requiere:
Espacio en piso, Transporte, Montacargas
Sistemas de transportadores
Interés sobre el costo de los materiales
671
Inventario en exceso
Las partes, materias primas, inventario en proceso,
refacciones y productos terminados forman el
inventario, el inventario es Muda ya que requiere:
Espacio en piso, Transporte, Montacargas
Sistemas de transportadores
Interés sobre el costo de los materiales
672
Reparaciones / defectos
Las reparaciones o el retrabajo de partes defectivas
significa un segundo intento de producirlas bien. Se
rompe el Takt Time
673
Movimientos
Los movimientos adicionales del personal son Muda.
Caminar mucho, cargar pesado, agacharse, estirarse
mucho, repetir movimientos, etc.
674
Movimientos
Algunas reglas de la ergonomía incluyen:
Enfatizar la seguridad todas las veces
Adecuar el empelado a la tarea
Cambiar el lugar de trabajo para que se adecue al
empleado
Mantener posiciones neutrales del cuerpo
Rediseñar las herramientas para reducir esfuerzo y
daños
Variar las tareas con rotación de puestos
Hacer que la máquina sirva al ser humano
675
Reprocesos
Consiste de pasos adicionales en el proceso de
manufactura, por ejemplo:
Remoción de rebabas
Maquinado de partes mal moldeadas
Agregar procesos de manejo adicionales
Realizar procesos de inspección
Repetir cambios al producto innecesarios
Mantener copias adicionales de información
676
Transportes
Todo transporte es Muda excepto la entrega al
cliente. Incluye:
Uso de montacargas
Uso de transportadores
Uso de movedores de pallets y camiones
677
Esperas
Ocurre cuando un operador está listo para realizar su
operación, pero permanece ocioso, por falla de
máquina, falta de partes, paros de línea, etc. El Muda
de espera puede ser por:
Operadores ociosos
Fallas de maquinaria
Tiempos de ajuste y preparación largos
Tareas no programadas a tiempo
Flujo de materiales en lotes
Juntas largas e innecesarias
678
Mudas adicionales
Otros mudas adicionales a los 7 desperdicios son:
Recursos mal utilizados
Recursos poco utilizados
Actividades de conteo
Búsqueda de herramientas o partes
Sistemas múltiples
Manos múltiples
Aprobaciones innecesarias
Fallas de máquinas
Envío de producto defectivo al cliente o mal servicio
679
Salidas de la Fase de Análisis
Causas raíz validadas
680
Preguntas ejemplo
1. En un sentido amplio, cuantas de las siguientes causas de
variación en estudios multi vari pueden incluir elementos de
proceso relacionados con el tiempo:
I. Posicional
II. Cilíndrica
III. Temporal
a. I y II c. II y III
b. I y III d. I, II y III
683
Preguntas ejemplo
7. Una muestra aleatoria de tamaño n se toma de una gran población con
desviación estándar de 1.0”. El tamaño de muestra se determina de
manera que haya un 0.05 de probabilidad de riesgo de exceder 0.1” de
error de tolerancia al usar la media de la muestra para estimar la Mu.
¿Cuál de los siguientes valores es el más cercano al tamaño de muestra
requerido?
8. La diferencia entre poner alfa de 0.05 y alfa igual a 0.01 en una prueba
de hipótesis es:
a. Con alfa de 0.05 se tiene mayor tendencia a cometer un error tipo I
b. Con alfa de 0.05 se tiene más posibilidad de riesgo de cometer un error
tipo II
c. Con alfa de 0.05 es una prueba más “conservadora” de la hipótesis nula
Ho
d. Con alfa de 0.05 se tiene menos posibilidad de cometer un error tipo I
684
Preguntas ejemplo
10. En una muestra aleatoria de 900 vehículos, 80% tienen frenos
ABS. ¿Cuál es el intervalo del 95% para el porcentaje de
vehículos con frenos ABS?
a. 0.778 – 0.821 c. 0.639 – 0.964
b. 0.771 – 0.829 d. 0.774 – 0.826
14. Un análisis de varianza de dos vías tiene r niveles para una variable y c
niveles para la otra, con dos obseraciones por celda. Los grados de
libertad para la interacción son:
a. 2 (r ) (c ) b. (r-1) (c-1) c. rc – 1 d. 2 (r -1) (c – 1)
Preguntas ejemplo
15. Los supuestos básicos del análisis de varianza oncluyen:
I. Las observaciones vienen de poblaciones normales
II. Las observaciones vienen de poblaciones con vaianzas iguales
III. Las observaciones vienen de poblaciones con medias iguales
a. I y II c. II y III
b. I y III d. I, II y III