Determinación del Tiempo de Garantía Vía Utilidad...
Transcript of Determinación del Tiempo de Garantía Vía Utilidad...
Determinación del Tiempo de Garantía Vía Utilidad Máxima.
Fifth ISBIS, Universidad Nacional Mayor de San Marcos.
Lima, Perú,12 de Enero de 2006.
Dr. Víctor Aguirre TorresDepartamento de Estadística, ITAM.
Trabajo conjunto con Dr. Humberto Gutiérrez (Universidad de Guadalajara) y Dr. Andrés Christen (CIMAT)
Tiempo de Garantía, V. Aguirre, ITAM
2
Bosquejo
Introducción.Función de Utilidad.Densidad Predictiva a Posteriori.Aplicación.Conclusiones.Bibliografí[email protected] http://cursos.itam.mx/vaguirre/Presentaciones
Tiempo de Garantía, V. Aguirre, ITAM
3
Introducción
Uno de los objetivos básicos de muchos de los estudios de confiabilidad es fijar el período de la garantía para un productoPráctica común sólo recomendar el uso de cuantiles pequeños de la distribución del tiempo de fallaLibros sobre confiabilidad: Lawless (1982), Nelson (1982), Ansell y Phillips (1994), Meeker y Escobar (1998).
Tiempo de Garantía, V. Aguirre, ITAM
4
IntroducciónGarantía: una decisión importante
La garantía habla de la calidad y confiabilidad del productoSi no está respaldada por datos:– Consumidores insatisfechos– Gastos adicionales para cumplirLa garantía tiene varias implicaciones: legales, costos, mercadotecnia, sociales, administración, y de diseño e ingeniería.
Tiempo de Garantía, V. Aguirre, ITAM
5
Introducción
Literatura especializada: “Product WarrantyHandbook” Eds. Blischke y Murthy (1996): perspectivas legal, social, de mercado, ingeniería y administración.– Algunos modelos matemáticos para los aspectos
de garantía, como costos y modelos de falla.La aproximación estadística a la decisión deltiempo de garantía prácticamente no se aborda
Tiempo de Garantía, V. Aguirre, ITAM
6
Introducción
Menezes y Currim (1992): Determinar período de garantía usando una función de demanda loglineal que depende de: precio, tiempo de garantía, publicidad, puntos de venta, confiabilidad del producto, y características distintivas del productoRequiere del conocimiento de los mismos aspectos de productos competidores.Necesita tanta información previa que resulta difícil de aplicar.
Tiempo de Garantía, V. Aguirre, ITAM
7
Introducción
Kim and Rao (2000), Suzuki (1985), Lawless(1998), Kalbfleisch y Lawless(1996): Usan reclamaciones y unidades que no han fallado para estimar distribución del tiempo de vida y costos de garantía.No se cuestiona cómo se selccionó el tiempo de garantía.
Tiempo de Garantía, V. Aguirre, ITAM
8
IntroducciónSingpurwalla y Wilson (1993): exponen una visión global del problema de garantía tomando en cuenta en particular la maximización de una función de utilidad esperada pero no llegan a una solución específica, aparte el modelo de falla está indexado por dos escalas (uso y tiempo).Singpurwalla y Wilson (1998): modelo de falla indexado por dos escalas (uso y tiempo). Cálculo de reservas monetarias.
Tiempo de Garantía, V. Aguirre, ITAM
9
Introducción
Objetivo: elección del tiempo de garantía vía maximización de la utilidad del productor tomando en cuenta– beneficio económico– costo de garantía– costo de imagen– calidad del producto
Enfoque Bayesiano da contexto natural.
Tiempo de Garantía, V. Aguirre, ITAM
10
Función de Utilidad
L= Manufacturas vendidasti = tiempo de falla del i-ésimo producto.tw= tiempo de garantía.u(t, tw) = utilidad monetaria cuando el producto falla al tiempo t y se ofreció un tiempo de garantía tw .
Tiempo de Garantía, V. Aguirre, ITAM
11
Función de Utilidad
btw: beneficio económico para tiempo de garantía tw .rt, tw: costo de cumplir con garantía. It, tw: costo de insatisfacción del cliente.ut, tw btw rt, tw It, twEn necesario determinar estas funciones.
Tiempo de Garantía, V. Aguirre, ITAM
12
Utilidad Esperada
Utilidad total
Utilidad esperada:
utw i1
Lut i, tw
Eutw ELbtw E i1
Lrt i, tw E
i1
LIt i, tw
Tiempo de Garantía, V. Aguirre, ITAM
13
Función de UtilidadInformación sobre el producto
ps = precio de ventac = costo de producciónv = ps - c = ganancia bruta del productorte = garantía estándarM = tamaño del mercado para productocr = costo de reparar o reponer producto
Tiempo de Garantía, V. Aguirre, ITAM
14
Función Beneficio Económico.
Debería tener ganancias decrecientes a medida que tw tiende a infinito.
ELbtw A2 1 eA1tw
Tiempo de Garantía, V. Aguirre, ITAM
15
Función de UtilidadInformación prospectiva del productor
π(te , ps ) = porción actual del mercado (%)ta = tiempo de garantía atractivo al clienteIa = incremento en porción del mercado con taπ(ta , ps ) = π(te , ps ) (1+ Ia)L~ binomial M,t e,ps o t e,ps1 Ia
Tiempo de Garantía, V. Aguirre, ITAM
16
Función Beneficio
E(Lb( te))= M π(te , ps ) vE(Lb( ta))= M π(te , ps ) (1+ Ia) vSistema de ecuaciones no lineales. Respecto a A1
Existe solución única para si
btebta
1eA1te
1eA1ta 11Ia
teta
11Ia
ELbtw A2 1 eA1tw
gx 1 exte/1 exta
gA1 11I a
Tiempo de Garantía, V. Aguirre, ITAM
17
Función Beneficio
Respecto a A2
A2 bta
A2 Mvt e,ps1 Ia
Ca Mt e,ps1 Ia A2 v
btw A2Ca1 eA1 tw
ELbtw A2 1 eA1tw
Tiempo de Garantía, V. Aguirre, ITAM
18
Costo Esperado de Garantía.
Tipos de garantía: reemplazo, reparación, prorrateo.En prorrateo A3 es la proporción del precio de venta que recibe el cliente.En reemplazo y reparación A3=0 .
rt, tw cr10,twt1 A 3 ttw
E i1
Lrt i, tw ELErt, tw Ca o
tw cr1 A 3 ttwft|Xdt
Tiempo de Garantía, V. Aguirre, ITAM
19
Costo Esperado de Insatisfacción.
Refleja la insatisfacción del cliente por falla dentro de garantía.qps= Costo del máximo nivel deinsatisfacción cuando el producto falla dentro de la garantía (0<q).
It, tw qps10,twt1 ttw
E i1
LIt i, tw ELEIt, tw Ca o
tw qps1 ttwft|Xdt
Tiempo de Garantía, V. Aguirre, ITAM
20
Selección del Tiempo de Garantía.
Por lo tanto, si
Entonces
Implementación numérica
tw
tw
arg max utw
utw Eutw
utw 1N k1
N ut k, tw
Tiempo de Garantía, V. Aguirre, ITAM
21
Densidad Predictiva a Posteriori.
ft|XRepresenta la calidad del producto. Se obtieneconsiderando:– información a priori sobre t, – los modelos usuales en confiabilidad – los datos de experimentos en confiabilidad (usualmente de
naturaleza censurada).
ft|X ft|θπθ|Xdθ
Teorema de Bayes|X LX|
LX|d
Tiempo de Garantía, V. Aguirre, ITAM
22
Densidad Predictiva a Posteriori.
Tiene forma analítica para los casos sencillos, pero no para las situaciones usuales en confiabilidad: muestras censuradas y modelos como el Weibull.Obtención numérica de ft|X por Sampling-Importance-ResamplingPara especificar π(θ) se procedió de acuerdo al método de Gutiérrez, Aguirre y Christen (2005).
Tiempo de Garantía, V. Aguirre, ITAM
23
Aplicación: Juego de Balatas Automotrices. Parámetros Iniciales.
ps = precio de venta = 900 pesosc = costo de producción = 700 pesosv = ps - c = ganancia del productor = 200 pesoste = garantía estándar = 10 ( x 1000 Kms.)ta = tiempo de garantía aún atractivo al cliente = 15 (x1000 Kms.)Ia = incremento en penetración de mercado con ta = 10% .Garantía de reemplazo, A3=0 .cr = costo de reparar o reponer producto = 700 pesosq = costo de imagen = 5%
Tiempo de Garantía, V. Aguirre, ITAM
24
Modelos
Los datos corresponden a 40 observaciones de naturaleza censurada por la derecha de manera aleatoria con r=17 fallas.Como modelos se consideraron: normal (M1), lognormal(M2), valor extremo(M3) y Weibull(M4) .Utilizando el método de Gutiérrez, Aguirre y Christen (2003) se obtuvo que Pr(M4|X)mayor a los otros casos.
Tiempo de Garantía, V. Aguirre, ITAM
25
Distribución a Priori.
Modelo Weibull con parámetro θ = (η,λ)
Para definir la distribución inicial se usó el método de Gutiérrez, Aguirre y Christen(2005).Requiere del usuario intervalos iniciales de la media [18,24] y desviación estándar [2.5,6.5] para el tiempo de vida.π(θ|α,β,m,k) = Normal-gama (2.16,11.6; 2.8,0.195).
1,2|,,m,k
k
21/22
1/2 exp k22 1 m2 2
Tiempo de Garantía, V. Aguirre, ITAM
26
Distribución Posterior Predictiva f(t|X)
t
f(t|X
)
10 20 30
0.0
0.02
0.04
0.06
t
f(t|X)
Tiempo de Garantía, V. Aguirre, ITAM
27
Utilidad Esperada
tw
u*(tw
)
6 8 10 12 14 16 18
130
140
150
160
170
180
tw
u*(tw)
Óptimo=13.00462Pr(T<13.0|X)=0.085
Tiempo de Garantía, V. Aguirre, ITAM
28
Análisis de Sensibildad
Se hizo un análisis de sensibilidad para evaluar los efectos de la especificación de la a priori: los intervalos para la media y para la desviación estándar, y del valor de IaSe utilizó un diseño 23 haciendo variaciones de ±10% en los tres factores. Los cambios en t* no resultaron significativos.SIR con 100,000 ensayos y efectividad de 25%.
Tiempo de Garantía, V. Aguirre, ITAM
29
Decisión
La mejor decisión para el productor es t*=13,000 Km, tomando en cuenta la confiabilidad del producto, los costos, la satisfacción del cliente y apreciaciones sobre el mercado.La intención de subirlo hasta 15,000 no va en el interés del productor, ya que su utilidad decrece un 5% aprox.
Tiempo de Garantía, V. Aguirre, ITAM
30
ConclusionesSe tiene un método para determinar el tiempo de garantía que incluye: confiabilidad del producto, competitividad de la garantía, costos de garantía y pérdida de imagen por fallas en la garantía.El enfoque se basa en la maximización de una función de utilidad esperada respecto a f(t|X), y no en algo como “un percentil bajo de la distribución del tiempo de falla”.Se mostró un ejemplo de aplicación
Tiempo de Garantía, V. Aguirre, ITAM
31
Desarrollos Posteriores.
Inclusión de competidores.Generalizar a garantías con dos escalas (tiempo, uso)[email protected] http://cursos.itam.mx/vaguirre/Presentaciones
Tiempo de Garantía, V. Aguirre, ITAM
32
ReferenciasAnsell, J.I., and Phillips M.J.(1994). Practical Methods forReliability Data Analysis. Clarendon Press, Oxford. Blischke, W.R. and Murthy, D.N. P. (Eds.) (1996). ProductWarranty Handbook. Marcel Dekker: New York. Gutiérrez-Pulido, H., Aguirre-Torres, V. and Christen, J.A.(2005) "A practical method for obtaining prior distributions in reliability". IEEE Transactions on Reliability 54(2), 262- 269 Gutiérrez-Pulido, H., Aguirre-Torres, V. and Christen, J.A.(2003) "Bayesian model evaluation in reliability". Reporte técnico DE-C03.17, Departamento de Estadística, ITAM, México. Gutiérrez-Pulido, H., Aguirre-Torres, V. and Christen, J.A.(2006) "A bayesian approach for the determination of warranty length". Aceptado para publicación Journal of Quality Technology.Kalbfleisch, J.D., Lawless J.F., and Robinson J.A.(1991)."Methods for the Analysis and Prediction of Warranty Claims." Technometrics, Vol. 33, No. 3, 273-285
Tiempo de Garantía, V. Aguirre, ITAM
33
ReferenciasKim, H-G. and Rao B.M.(2000). "Expected Warranty Cost of Two-attribute Free-Replacement Warranties Based on a Bivariate Exponential Distribution." Computers & Industrial Engineering, 38, 425-434. Lawless, J.F. (1982). Statistical Models and Methods for Lifetime Data. John Wiley, New York. Lawless, J.F. (1998). "Statistical Analysis of Product Warranty data". International Statistical Review, 66, 41-60. Meeker, W.Q. and E. Escobar (1998). Statistical Methods for Reliability Data. John Wiley, New York. Menezes, M.A.J. and Currim, I.S. (1992). "An Approach for Determination of Warranty Length". International Journal of Research in Marketing, 9, 177-195. Nelson, W. (1982). Applied Life Data Analysis. John Wiley, New York.Singpurwalla; N.D. and S.P. Wilson (1993). "The warranty problem: its statistical and game theoretic aspects." SIAM Review., Vol. 35, No. 1, 17-42.
Tiempo de Garantía, V. Aguirre, ITAM
34
Referencias
Singpurwalla; N.D. and S.P. Wilson (1998). "Failure ModelsIndexed by Two Scales." Advances in Applied Probability, 30, 1058-1072. Suzuki, K. (1985). "Estimation Method of Lifetime Based on theRecord of Failure During the Warranty Period". Journal of theAmerican Statistical Association, 80, 66-72.