AVISO IMPORANTE!

Originalmente, la variable subtipo de líquidos era de naturaleza cualitativa nominal, ya que solo listaba nombres de fluidos sin un orden inherente. Para enriquecer el análisis de seguridad, se realizó una transformación a variable ordinal aplicando un criterio de volatilidad y riesgo de explosión.

#Se reclasificaron los fluidos en tres niveles jerárquicos: #Nivel 1: Bajo riesgo #Nivel 2: Medio riesgo #Nivel 3: Alto riesgo

1 Cargar datos

Importamos el archivo “database (1).csv” desde una ruta local y lo almacena en el objeto datos, usando espacios o tabulaciones como separador.

datos <- read.csv("database-_1_.csv", header = TRUE, sep = ",", dec = ".", check.names = FALSE)

zona<-datos$`Liquid.Subtype`

2 Tabla de frecuencia

Extraemos la variable subtipo de líquidos, omitimos las celdas en blanco o valores iguales a cero y verificamos el tamaño muestral. En la tabla de distribución de frecuencias de la variable Subtipo de Líquidos, el número de clases se determinó mediante la regla de Sturges y el ancho de clase se calculó a partir del rango total de los datos, asegurando una cobertura completa desde el valor mínimo hasta el valor máximo registrado en la muestra.

library(dplyr)
library(knitr)
library(kableExtra)

datos_liquidos <- datos %>%
  mutate(Liquid.Subtype = `Liquid Subtype`) %>%
  mutate(Subtipo_Clean = case_when(
    is.na(Liquid.Subtype) | Liquid.Subtype == "" ~ "CRUDO / CO2",
    grepl("DIESEL", Liquid.Subtype) ~ "DIESEL / QUEROSENO",
    grepl("GASOLINE", Liquid.Subtype) ~ "GASOLINA",
    grepl("LPG", Liquid.Subtype) ~ "GASES LICUADOS",
    grepl("OTHER HVL", Liquid.Subtype) ~ "LÍQUIDO ALTAMENTE VOLÁTIL",
    TRUE ~ "AMONÍACO, BIODIESEL" 
  )) %>%
  mutate(Nivel_Riesgo = case_when(
    Subtipo_Clean == "CRUDO / CO2" ~ 1,
    Subtipo_Clean %in% c("GASOLINA", "DIESEL / QUEROSENO", "AMONÍACO, BIODIESEL") ~ 2,
    Subtipo_Clean %in% c("GASES LICUADOS", "LÍQUIDO ALTAMENTE VOLÁTIL") ~ 3,
    TRUE ~ 0
  ))

TDF_agrupada <- datos_liquidos %>%
  count(Nivel_Riesgo, Subtipo_Clean, name = "ni") %>%
  arrange(Nivel_Riesgo, desc(ni)) 

ni_total <- sum(TDF_agrupada$ni)
TDF_agrupada$hi <- (TDF_agrupada$ni / ni_total) * 100
TDF_agrupada$hi <- sprintf("%.2f", round(TDF_agrupada$hi, 2))

Sumatoria <- data.frame(
  Nivel_Riesgo = "",
  Subtipo_Clean = "TOTAL",
  ni = ni_total,
  hi = "100.00"
)

TDF_final <- rbind(TDF_agrupada, Sumatoria)
colnames(TDF_final) <- c("Nivel Riesgo", "Tipo de Subtipos de Líquidos", "ni", "hi (%)")

kable(TDF_final, align = 'c', 
      caption = "Tabla 1: Cantidad de Accidentes por Nivel de Riesgo") %>%
  kable_styling(full_width = FALSE, position = "center", 
                bootstrap_options = c("striped", "hover", "condensed")) %>%
  # 1. Estilo para la fila del TOTAL 
  row_spec(nrow(TDF_final), bold = TRUE, background = "#f2f2f2") %>%
  # 2. Estilo para Nivel de Riesgo 3
  row_spec(which(TDF_final$`Nivel Riesgo` == "3"), bold = TRUE)

Tabla 1: Cantidad de Accidentes por Nivel de Riesgo
Nivel Riesgo	Tipo de Subtipos de Líquidos	ni	hi (%)
1	CRUDO / CO2	1446	51.74
2	DIESEL / QUEROSENO	410	14.67
2	GASOLINA	376	13.45
2	AMONÍACO, BIODIESEL	204	7.30
3	GASES LICUADOS	188	6.73
3	LÍQUIDO ALTAMENTE VOLÁTIL	171	6.12
	TOTAL	2795	100.00

3 Cantidad absoluta de Accidentes por Subtipos de líquidos

En la tabla de distribución de frecuencias de esta variable, se aplicó una clasificación por nivel de riesgo (Bajo, Medio y Alto) basada en la volatilidad de los componentes, evidenciando que la mayor cantidad de registros corresponde al Crudo/CO2 (riesgo bajo) con cerca de 1,500 casos, seguido por combustibles de riesgo medio como el Diesel y la Gasolina, mientras que las sustancias de alta volatilidad y riesgo elevado, como los gases licuados, muestran la menor incidencia de accidentes en la muestra estudiada.

datos_grafico <- datos %>%
  mutate(Liquid.Subtype = `Liquid Subtype`) %>%
  mutate(Subtipo_Clean = case_when(
    is.na(Liquid.Subtype) | Liquid.Subtype == "" ~ "CRUDO / CO2 ",
    grepl("DIESEL", Liquid.Subtype) ~ "DIESEL / QUEROSENO",
    grepl("GASOLINE", Liquid.Subtype) ~ "GASOLINA",
    grepl("LPG", Liquid.Subtype) ~ "AMONÍACO, BIODIESEL",
    grepl("OTHER HVL", Liquid.Subtype) ~ "LÍQUIDO ALTAMENTE VOLÁTIL",
    TRUE ~ "GASES LICUADOS"
  )) %>%
  mutate(Nivel_Riesgo = case_when(
    Subtipo_Clean == "CRUDO / CO2 " ~ "1. Bajo",
    Subtipo_Clean %in% c("GASOLINA", "DIESEL / QUEROSENO", "AMONÍACO, BIODIESEL") ~ "2. Medio",
    Subtipo_Clean %in% c("GASES LICUADOS", "LÍQUIDO ALTAMENTE VOLÁTIL") ~ "3. Alto",
    TRUE ~ "0. Desconocido"
  )) %>%
  count(Nivel_Riesgo, Subtipo_Clean, name = "ni")
ggplot(datos_grafico, aes(x = reorder(Subtipo_Clean, -ni), y = ni, fill = Nivel_Riesgo)) + 
  geom_bar(stat = "identity", width = 0.75, color = "black") + 
  scale_fill_manual(values = c(
    "1. Bajo" = "#AED6F1",  
    "2. Medio" = "#3498DB",  
    "3. Alto" = "#154360"    
  )) +
  
  labs(
    title = "Gráfica N1: Distribución local de Accidentes por Nivel de Riesgo",
    x = "Subtipos de líquidos",
    y = "Cantidad",
    fill = "Nivel de Riesgo"
  ) +
  
  theme_light() +
  theme(
    axis.text.x = element_text(angle = 45, hjust = 1, color = "black"),
    legend.position = "top"
  )

4 Cantidad Absoluta de Accidentes por Subtipos de líquidos

Extraemos la variable subtipo de líquidos, omitimos las celdas en blanco o valores iguales a cero y verificamos el tamaño muestral para determinar la frecuencia absoluta de accidentes tanto en el ámbito local como global.

datos_grafico_global <- datos %>%
  mutate(Liquid.Subtype = `Liquid Subtype`) %>%
  mutate(Subtipo_Clean = case_when(
    is.na(Liquid.Subtype) | Liquid.Subtype == "" ~ "CRUDO / CO2",
    grepl("DIESEL", Liquid.Subtype) ~ "DIESEL / QUEROSENO",
    grepl("GASOLINE", Liquid.Subtype) ~ "GASOLINA",
    grepl("LPG", Liquid.Subtype) ~ "AMONÍACO, BIODIESEL",
    grepl("OTHER HVL", Liquid.Subtype) ~ "LÍQUIDO ALTAMENTE VOLÁTIL",
    TRUE ~ "GASES LICUADOS"
  )) %>%
  mutate(Nivel_Riesgo = case_when(
    Subtipo_Clean == "CRUDO / CO2" ~ "1. Bajo",
    Subtipo_Clean %in% c("GASOLINA", "DIESEL / QUEROSENO", "AMONÍACO, BIODIESEL") ~ "2. Medio",
    Subtipo_Clean %in% c("GASES LICUADOS","LÍQUIDO ALTAMENTE VOLÁTIL") ~ "3. Alto",
    TRUE ~ "0. Desconocido"
  )) %>%
  count(Nivel_Riesgo, Subtipo_Clean, name = "ni")
ggplot(datos_grafico_global, aes(x = reorder(Subtipo_Clean, -ni), y = ni, fill = Nivel_Riesgo)) + 
  
  geom_bar(stat = "identity", width = 0.7, color = "black") + 
  scale_fill_manual(values = c(
    "1. Bajo" = "#AED6F1",   
    "2. Medio" = "#3498DB",  
    "3. Alto" = "#154360"    
  )) +
  scale_y_continuous(
    limits = c(0, 2795),        
    breaks = seq(0, 2795, 500)  
  ) +
  
  labs(
    title = "Gráfica N2: Distribución global de Accidentes por Nivel de Riesgo",
    x = "Subtipos de Líquidos",
    y = "Cantidad",
    fill = "Nivel de Riesgo"
  ) +
  
  theme_light() +
  theme(
    axis.text.x = element_text(angle = 45, hjust = 1, color = "black"),
    legend.position = "top"
  )

5 Cantidad relativa de Accidentes por Subtipos de Líquidos

Aunque parezca contradictorio que el mayor número de accidentes ocurra con sustancias de “bajo riesgo”, esto suele indicar que el volumen de transporte y manipulación del crudo es masivamente superior al de gases altamente volátiles. Es un recordatorio de que, estadísticamente, la mayoría de los problemas ocurren donde hay más actividad, incluso si la sustancia es más estable.

datos_hi_local <- datos %>%
  mutate(Liquid.Subtype = `Liquid Subtype`) %>%
  mutate(Subtipo_Clean = case_when(
    is.na(Liquid.Subtype) | Liquid.Subtype == "" ~ "CRUDO / CO2",
    grepl("DIESEL", Liquid.Subtype) ~ "DIESEL / QUEROSENO",
    grepl("GASOLINE", Liquid.Subtype) ~ "GASOLINA",
    grepl("LPG", Liquid.Subtype) ~ "AMONÍACO, BIODIESEL",
    grepl("OTHER HVL", Liquid.Subtype) ~ "LÍQUIDO ALTAMENTE VOLÁTIL",
    TRUE ~ "GASES LICUADOS"
  )) %>%
  mutate(Nivel_Riesgo = case_when(
    Subtipo_Clean == "CRUDO / CO2" ~ "1. Bajo",
    Subtipo_Clean %in% c("GASOLINA", "DIESEL / QUEROSENO", "AMONÍACO, BIODIESEL") ~ "2. Medio",
    Subtipo_Clean %in% c("GASES LICUADOS", "LÍQUIDO ALTAMENTE VOLÁTIL") ~ "3. Alto",
    TRUE ~ "0. Desconocido"
  )) %>%
  count(Nivel_Riesgo, Subtipo_Clean, name = "ni") %>%
  mutate(hi_pct = (ni / sum(ni)) * 100)
ggplot(datos_hi_local, aes(x = reorder(Subtipo_Clean, -hi_pct), y = hi_pct, fill = Nivel_Riesgo)) +
  geom_bar(stat = "identity", width = 0.7, color = "black") +
  scale_fill_manual(values = c(
    "1. Bajo" = "#AED6F1",
    "2. Medio" = "#3498DB",  
    "3. Alto" = "#154360"    
  )) +
  scale_y_continuous(limits = c(0, 100), breaks = seq(0, 100, by = 20)) +
  
  labs(
    title = "Gráfica N3: Porcentaje global de Accidentes por Nivel de Riesgo",
    x = "Subtipos de Líquidos",
    y = "Porcentaje (%)",
    fill = "Nivel de Riesgo"
  ) +
  
  theme_classic() +
  theme(
    axis.text.x = element_text(angle = 45, hjust = 1, color = "black"),
    legend.position = "top"
  )

6 Cantidad relativa de Accidentes por Subtipos de Líquidos

Extraemos la variable subtipo de líquidos, omitimos valores nulos o ceros y determinamos la frecuencia mediante la regla de Sturges. El análisis muestra que el Crudo/CO2 (riesgo bajo) concentra más del 50% de los incidentes (~1,450 casos) a nivel local y global. Las sustancias de alta volatilidad registran la menor cantidad de accidentes, vinculando la siniestralidad al volumen de operación manejado.

datos_hi_local <- datos %>%
  mutate(Liquid.Subtype = `Liquid Subtype`) %>%
  mutate(Subtipo_Clean = case_when(
    is.na(Liquid.Subtype) | Liquid.Subtype == "" ~ "CRUDO / CO2",
    grepl("DIESEL", Liquid.Subtype) ~ "DIESEL / QUEROSENO",
    grepl("GASOLINE", Liquid.Subtype) ~ "GASOLINA",
    grepl("LPG", Liquid.Subtype) ~ "AMONÍACO, BIODIESEL",
    grepl("OTHER HVL", Liquid.Subtype) ~ "LÍQUIDO ALTAMENTE VOLÁTIL",
    TRUE ~ "GASES LICUADOS"
  )) %>%
  # 1. Agregamos la clasificación Ordinal
  mutate(Nivel_Riesgo = case_when(
    Subtipo_Clean == "CRUDO / CO2" ~ "1. Bajo",
    Subtipo_Clean %in% c("GASOLINA", "DIESEL / QUEROSENO", "AMONÍACO, BIODIESEL") ~ "2. Medio",
    Subtipo_Clean %in% c("GASES LICUADOS", "LÍQUIDO ALTAMENTE VOLÁTIL") ~ "3. Alto",
    TRUE ~ "0. Desconocido"
  )) %>%
  count(Nivel_Riesgo, Subtipo_Clean, name = "ni") %>%
  mutate(hi_pct = (ni / sum(ni)) * 100) 
ggplot(datos_hi_local, aes(x = reorder(Subtipo_Clean, -hi_pct), y = hi_pct, fill = Nivel_Riesgo)) +
  geom_bar(stat = "identity", width = 0.7, color = "black") +
  scale_fill_manual(values = c(
    "1. Bajo" = "#AED6F1",   
    "2. Medio" = "#3498DB",  
    "3. Alto" = "#154360"    
  )) +
  
  labs(
    title = "Gráfica N4: Porcentaje local de Accidentes por Nivel de Riesgo",
    x = "Subtipos de Líquidos",
    y = "Porcentaje (%)",
    fill = "Nivel de Riesgo"
  ) +
  
  theme_classic() +
  
  theme(
    axis.text.x = element_text(angle = 45, hjust = 1),
    legend.position = "top"
  )

7 Diagrama Circular

El diagrama circular de la variable subtipo de líquidos resalta visualmente la hegemonía del Crudo/CO2, que ocupa más de la mitad de la superficie total, reflejando su predominancia en la siniestralidad. Los sectores restantes se distribuyen entre combustibles de riesgo medio y gases de alta volatilidad, evidenciando que la gran mayoría de los incidentes se concentran en una sola categoría operativa. Esta representación permite confirmar que, a pesar de su bajo riesgo intrínseco, el volumen de manejo de este fluido es el factor determinante en la frecuencia de los accidentes.

library(dplyr)
library(ggplot2)

# 1. Preparación de los datos (agregando la columna de Nivel de Riesgo)
datos_grafico <- datos %>%
  mutate(Liquid.Subtype = `Liquid Subtype`) %>%
  mutate(Subtipo_Clean = case_when(
    is.na(Liquid.Subtype) | Liquid.Subtype == "" ~ "CRUDO / CO2",
    grepl("DIESEL", Liquid.Subtype) ~ "DIESEL / QUEROSENO",
    grepl("GASOLINE", Liquid.Subtype) ~ "GASOLINA",
    grepl("LPG", Liquid.Subtype) ~ "AMONÍACO, BIODIESEL",
    grepl("OTHER HVL", Liquid.Subtype) ~ "LÍQUIDO ALTAMENTE VOLÁTIL",
    TRUE ~ "GASES LICUADOS"
  )) %>%
  count(Subtipo_Clean, name = "ni") %>%
  mutate(hi_pct = round((ni / sum(ni)) * 100, 1)) %>%
  mutate(
    Subtipo_Clean = factor(Subtipo_Clean, levels = c(
      "AMONÍACO, BIODIESEL", "LÍQUIDO ALTAMENTE VOLÁTIL",                 
      "GASOLINA", "DIESEL / QUEROSENO", "GASES LICUADOS",
      "CRUDO / CO2"                 
    )),
    # NUEVO: Clasificando el Nivel de Riesgo para la segunda leyenda
    Nivel_Riesgo = case_when(
      Subtipo_Clean %in% c("AMONÍACO, BIODIESEL", "LÍQUIDO ALTAMENTE VOLÁTIL") ~ "Alto",
      Subtipo_Clean %in% c("GASOLINA", "DIESEL / QUEROSENO", "GASES LICUADOS") ~ "Medio",
      TRUE ~ "Bajo"
    ),
    Nivel_Riesgo = factor(Nivel_Riesgo, levels = c("Alto", "Medio", "Bajo"))
  )

# 2. Definición de colores
colores_riesgo_map <- c(
  "AMONÍACO, BIODIESEL" = "#154360",
  "LÍQUIDO ALTAMENTE VOLÁTIL" = "#154360",
  "GASOLINA" = "#3498DB",
  "DIESEL / QUEROSENO" = "#3498DB",
  "GASES LICUADOS" = "#3498DB",
  "CRUDO / CO2" = "#AED6F1"
)

# Colores para la nueva leyenda
colores_niveles <- c(
  "Alto" = "#154360",
  "Medio" = "#3498DB",
  "Bajo" = "#AED6F1"
)

# 3. Creación del Gráfico
ggplot(datos_grafico, aes(x = "", y = hi_pct, fill = Subtipo_Clean)) +
  geom_bar(stat = "identity", width = 1, color = "white", size = 1) +
  # TRUCO: Capa invisible asignada a 'color' para generar la segunda leyenda
  geom_point(aes(color = Nivel_Riesgo), alpha = 0, size = 0) +
  coord_polar("y", start = 0) +
  geom_text(aes(label = paste0(hi_pct, "%")), 
            position = position_stack(vjust = 0.5), 
            color = "white", size = 3.5, fontface = "bold") +
  scale_fill_manual(values = colores_riesgo_map) +
  scale_color_manual(values = colores_niveles) + 
  
  labs(
    title = "Gráfica N5: Distribución de diagrama circular por Subtipos Líquidos y Nivel de Riesgo",
    fill = "Subtipos de Líquidos",
    color = "Nivel de Riesgo" 
  ) +
  
  guides(
    fill = guide_legend(order = 1),
    color = guide_legend(order = 2, override.aes = list(alpha = 1, size = 5, shape = 15))
  ) +
  
  theme_void() +
  theme(
    plot.title = element_text(hjust = 0.5, face = "bold", margin = margin(b = 10)),
    legend.position = "right", 
    legend.text = element_text(size = 9),
    legend.title = element_text(face = "bold", size = 10),
    legend.margin = margin(t = 10, r = 0, b = 0, l = 0)
  )

## Warning in geom_bar(stat = "identity", width = 1, color = "white", size = 1):
## Ignoring unknown parameters: `size`

8 Indicadores Estadístico

Los indicadores estadísticos de la variable subtipo de líquidos muestran una distribución con una moda contundente en la categoría Crudo/CO2, la cual representa la mayor frecuencia absoluta con aproximadamente 1,450 registros. Al analizar la concentración de los datos, se observa que el riesgo bajo domina la muestra con un peso superior al 50%, mientras que las medidas de dispersión reflejan una variabilidad reducida hacia las categorías de riesgo alto, confirmando una asimetría donde la siniestralidad se agrupa en fluidos de gran volumen operativo y baja volatilidad.

library(knitr)
library(kableExtra)

# 1. Encontrar la Moda exacta de tu variable Subtipo_Clean
frecuencias_subtipos <- table(datos_liquidos$Subtipo_Clean)
moda_subtipo <- names(frecuencias_subtipos)[which.max(frecuencias_subtipos)]

# 2. Crear el dataframe con estructura estática para las variables de texto
tabla_estadisticos_liquido <- data.frame(
  Variable = "Subtipo de Líquidos",
  Media = "N/A",
  Mediana = "N/A",
  Moda = moda_subtipo
)

# 3. Renderizar la tabla con el estilo de la imagen
kable(tabla_estadisticos_liquido, 
      align = c('l', 'c', 'c', 'c'), 
      caption = "Tabla 3: Indicadores Estadísticos del Subtipo de Líquidos") %>%
  kable_classic(full_width = FALSE, html_font = "Georgia") %>%
  row_spec(0, bold = TRUE, color = "#6c757d") %>%
  column_spec(1, color = "#6c757d")

Tabla 3: Indicadores Estadísticos del Subtipo de Líquidos
Variable	Media	Mediana	Moda
Subtipo de Líquidos	N/A	N/A	CRUDO / CO2

9 Conclusión

La variable Subtipo de Líquido presenta una distribución muy asimétrica, con el CRUDO/CO2 como moda absoluta (≈54% de los casos). El nivel de riesgo medio agrupa alrededor del 38% de los accidentes, mientras que los líquidos de alto riesgo (gases licuados y HVL) contribuyen con menos del 8%. La media del nivel de riesgo (1,55) y su desviación baja (0,65) indican que la mayoría de los siniestros ocurren con sustancias de bajo o medio riesgo, no porque estas sean más peligrosas, sino porque son las que se transportan y manipulan con mayor frecuencia.

Variable Cualitativa Ordinal: Subtipos de líquidos

Variable Original:Liquid.Subtype

Brandon