1 Carga de datos

setwd("/cloud/project/")
datos<-read.csv("DATOS.csv", header = TRUE, sep=";" , dec=".")
str(datos)
## 'data.frame':    10190 obs. of  17 variables:
##  $ Distrito_edit                        : chr  "1" "1" "1" "1" ...
##  $ Year_edit_Fecha_del_derrame          : int  2013 2013 2013 2013 2013 2013 2013 2013 2013 2013 ...
##  $ Mes_edit_Fecha_del_derrame           : int  6 3 4 4 6 6 3 9 10 6 ...
##  $ Categoria_Instalaciones              : chr  "Instalacion fija" "Pozos" "Pozos" "Pozos" ...
##  $ Operacion_general                    : chr  "Produccion" "Otro" "Produccion" "Produccion" ...
##  $ Categoria_Fuente                     : chr  NA "Tanques/Almacenamiento" "Lineas/Tuberias" "Infraestructura Fija" ...
##  $ Grupo_causas_probable                : chr  NA "Afectaciones externas" "Factores humanos" "Problemas tecnicos" ...
##  $ Liberacion_petroleo_crudo_edicion    : num  0 0 0 0 0 ...
##  $ Edicion_recuperacion_petroleo_crudo  : num  NA 0 0 0 0 0 0 0 0 NA ...
##  $ Volumen_liberado_Cond_Final          : num  0 0 0 10 0 0 0 1 0 0 ...
##  $ Liberacion_agua_de_produccion_edicion: num  6720 3780 5040 420 10920 ...
##  $ Liberacion_volumen_gas               : num  0 0 0 0 0 0 0 0 0 0 ...
##  $ Volumen_condensado_recuperado        : num  NA 0 0 1 0 0 0 0 0 NA ...
##  $ Edicion_Recuperacion_agua_producida  : num  NA 420 4620 0 10920 ...
##  $ Derrame_sobre_agua_limpio            : chr  "NO" "NO" "NO" "NO" ...
##  $ Estado_general                       : chr  "Observaciones tecnicas" NA NA NA ...
##  $ Codigo_area                          : int  1 1 1 1 1 1 1 1 1 3 ...

2 Extraccion de variable

Estado <- datos$Estado_general
Estado <- na.omit(Estado)

3 Creacion de Tabla de Frecuencias

TDFEstado <- table(Estado)
TablaEstado <- as.data.frame(TDFEstado)
names(TablaEstado) <- c("Categoria", "ni")

3.0.1 Frecuencia relativa, porcentaje y fila total

hi <- TablaEstado$ni / sum(TablaEstado$ni)
hi_porc <- hi * 100
TablaEstado$hi_porcentaje <- round(hi_porc, 2)

TDFFinalEstado <- rbind(TablaEstado, data.frame(
  Categoria = "TOTAL",
  ni = sum(TablaEstado$ni),
  hi_porcentaje = 100))

3.0.2 Tabla Final

print(TDFFinalEstado)
##                  Categoria   ni hi_porcentaje
## 1             Caso abierto  271          2.77
## 2     Cerrado / Controlado 7969         81.47
## 3 Evaluación / Seguimiento 1532         15.66
## 4   Observaciones tecnicas   10          0.10
## 5                    TOTAL 9782        100.00

4 Gráficos de Barras

4.0.1 GDF versión local de Distribución de tipos de estado del derrame

Nuevonombre_Tipoestado<-gsub("/","/\n",TablaEstado$Categoria)

par(mar = c(6, 6, 4, 2)) 
barplot(
  TDFEstado, 
  main = "Gráfica No.1: Distribución del Estado General del Derrame",
  xlab = "Categoría",
  ylab = "Cantidad",
  names.arg = Nuevonombre_Tipoestado,
  col = "darkseagreen2",
  las = 1,
  cex.main = 1.4,    
  cex.lab = 1.2,   
  cex.axis = 0.8,
  cex.names = 0.7)

4.0.2 GDF versión global de Distribución de tipos de estado del derrame

par(mar = c(6, 6, 4, 2)) 
barplot(
  TDFEstado,
  main = "Gráfica No.2: Distribución del Estado General del Derrame",
  xlab = "Categoría",
  ylab = "", 
  col = "darkseagreen3",
  names.arg = Nuevonombre_Tipoestado,
  las = 1, 
  ylim = c(0, sum(TDFFinalEstado$ni)),
  cex.names = 0.7)
mtext("Cantidad", side = 2, line = 4.5, cex = 1)

4.0.3 GDF: Porcentaje de Distribución de tipos de estado del derrame

n_filtro <- TDFFinalEstado$Categoria[TDFFinalEstado$Categoria != "TOTAL"]
Nuevonombre_Tipestado <- gsub("/","/\n", n_filtro)
par(mar = c(6, 6, 4, 2)) 
barplot(
  TDFFinalEstado$hi_porcentaje[TDFFinalEstado$Categoria != "TOTAL"],
  main = "Gráfica No.3: Distribución del Estado General del Derrame",
  xlab = "Categoría",
  ylab = "Porcentaje",
  names.arg = Nuevonombre_Tipestado,
  col = "darkseagreen4",
  cex.main = 1.5,
  las = 1,
  cex.lab = 1,
  cex.axis = 1,
  cex.names = 0.7,
  ylim = c(0, 100))

5 Diagrama circular

par(mar = c(0, 0, 4, 4), xpd = TRUE)
colores_estado <- colorRampPalette(c("darkseagreen1", "darkolivegreen3"))(nrow(TablaEstado))
etiquetas_porcentaje <- paste0(TablaEstado$hi_porcentaje, "%")

pie(
  TablaEstado$ni,
  labels = etiquetas_porcentaje,
  main = "Gráfica No.4: Distribución del Estado General del Derrame",
  col = colores_estado,
  cex = 0.8)

legend(
  x = 1, y = -0.5,
  legend = TablaEstado$Categoria,
  title = "Leyenda",
  fill = colores_estado,
  cex = 0.7,
  y.intersp = 0.7,
  bty = "o",
  inset = 0.2)

6 INDICADOR - Moda

moda_estado <- names(TDFEstado)[which.max(TDFEstado)]
print(paste("MODA (Estado General):", moda_estado))
## [1] "MODA (Estado General): Cerrado / Controlado"