La semana pasada iniciamos con la primera parte de la explicación del Gráfico G. En resumen, cuando hay pocas ocurrencias de un dato de interés (muchas veces solamente una en un período de tiempo), dibujar el tiempo entre eventos en lugar del evento como tal.
En esta segunda parte vamos a usar los datos de las erupciones del Volcán Turrialba (Costa Rica) desde el 19 de setiembre del 2016 hasta el 10 de julio del 2017. A continuación, el gráfico.
| Fecha | Días entre erupciones |
| 9/19/16 | 56 |
| 9/20/16 | 1 |
| 9/22/16 | 2 |
| 9/26/16 | 4 |
| 10/3/16 | 7 |
| 10/4/16 | 1 |
| 10/5/16 | 1 |
| 10/6/16 | 1 |
| 10/7/16 | 1 |
| 10/8/16 | 1 |
| 10/9/16 | 1 |
| 10/10/16 | 1 |
| 10/11/16 | 1 |
| 10/12/16 | 1 |
| 10/13/16 | 1 |
| 10/14/16 | 1 |
| 10/15/16 | 1 |
| 10/16/16 | 1 |
| 10/17/16 | 1 |
| 10/23/16 | 6 |
| 11/9/16 | 17 |
| 1/27/17 | 79 |
| 2/9/17 | 13 |
| 2/13/17 | 4 |
| 3/27/17 | 42 |
| 4/2/17 | 6 |
| 4/17/17 | 15 |
| 5/8/17 | 21 |
| 5/27/17 | 19 |
| 6/16/17 | 20 |
| 7/10/17 | 24 |
Nuevamente surge la interrogante ¿está disminuyendo la actividad del volcán? Esta vez vamos a agregar límites de control (dos propuestas) utilizando R y Minitab.
Primero vamos a usar el paquete qicharts de R.
El código para generar el gráfico es el siguiente:
volcan <- read.csv(“volcan.csv”)
install.packages(“qcc”)
library(qcc)
qcc(volcan$Dias, type = “g”)
Este código genera el siguiente gráfico:
El paquete simplemente calcula el promedio y la desviación estándar de los días entre erupciones y calcula el centro y límites de la forma tradicional +/- 3 desviaciones estándar. Como en otros gráficos, el límite inferior es generalmente negativo y por lo tanto se hace cero.
Minitab usa límites probabilísticos para evitar falsas alarmas. Como la distribución geométrica es muy sesgada a la derecha el límite superior sea muy bajo con respecto a la distribución. Del mismo modo que en el gráfico generado en R, el límite inferior es cero. Para una explicación completa ver aquí.
El resultado en Minitab, considerando límites probabilísticos es el siguiente:
La lectura, o lecturas sobre el límite superior parecen indicar que la frecuencia de erupciones estaría empezando a disminuir, el Ovsicori (Observatorio Vulcanológico y Sismológico de Costa Rica) se mantiene alerta y no baja la guardia, y todavía no recomienda reanudar actividad turísticas y agrícolas normales en la zona.
El Lean y Six Sigma constantemente estamos agregando herramientas a nuestra enorme caja, y esa es una de las principales razones por las que ambas metodologías siguen creciendo después de tantos años. Con el Gráfico G tenemos a disposición una nueva forma de investigar patrones cuando tenemos muy pocos datos, hagámoslo parte de nuestra caja de herramientas.
