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 ...
Creacion de Tabla de
Frecuencias
TDFcodigo_area <- table(Codigo_area)
Tablacodigo_area <- as.data.frame(TDFcodigo_area)
names(Tablacodigo_area) <- c("Categoria", "ni")
Frecuencia
relativa, porcentaje y fila total
hi <- Tablacodigo_area$ni / sum(Tablacodigo_area$ni)
hi_porc <- hi * 100
Tablacodigo_area$hi_porcentaje <- round(hi_porc, 2)
TDFcodigoFinal <- rbind(Tablacodigo_area, data.frame(
Categoria = "TOTAL",
ni = sum(Tablacodigo_area$ni),
hi_porcentaje = 100))
Tabla Final
## Categoria ni hi_porcentaje
## 1 1 7285 92.34
## 2 2 366 4.64
## 3 3 78 0.99
## 4 4 155 1.96
## 5 5 3 0.04
## 6 7 2 0.03
## 7 TOTAL 7889 100.00
Gráficos de Barras
GDF versión local
de Distribución de areas afectadas
par(mar = c(6, 6, 4, 2))
barplot(
TDFcodigo_area,
main = "Gráfica No.1: Distribución del tipo de Área afectada",
xlab = "Categoría",
ylab = "Cantidad",
col = "lightgoldenrod1",
las = 1,
cex.main = 1.4,
cex.lab = 1.2,
cex.axis = 0.8,
cex.names = 0.7)
GDF versión global
de Distribución de tipos de estado del derrame
par(mar = c(6, 6, 4, 2))
barplot(
TDFcodigo_area,
main = "Gráfica No.2: Distribución del tipo de Área afectada",
xlab = "Categoría",
ylab = "",
col = "goldenrod1",
las = 1,
ylim = c(0, sum(TDFcodigoFinal$ni)),
cex.lab = 1,
cex.names = 0.7)
mtext("Cantidad", side = 2, line = 4.5, cex = 1.4)
GDF: Porcentaje de
Distribución de tipos de estado del derrame
par(mar = c(6, 6, 4, 2))
barplot(
TDFcodigoFinal$hi_porcentaje[TDFcodigoFinal$Categoria != "TOTAL"],
main = "Gráfica No.3: Porcentaje del Código de Área",
xlab = "Categoría",
ylab = "Porcentaje",
col = "goldenrod3",
names.arg = TDFcodigoFinal$Categoria[TDFcodigoFinal$Categoria != "TOTAL"],
cex.main = 1.5,
las = 1,
cex.lab = 1,
cex.axis = 1,
cex.names = 0.7,
ylim = c(0, 100))
Diagrama circular
par(mar = c(0, 0, 4, 4), xpd = TRUE)
colores_codigo <- colorRampPalette(c("lightgoldenrod1", "goldenrod3"))(nrow(Tablacodigo_area))
etiquetas_porcentaje <- paste0(Tablacodigo_area$hi_porcentaje, "%")
pie(
Tablacodigo_area$ni,
labels = etiquetas_porcentaje,
main = "Gráfica No.4: Distribución del Código de Área",
col = colores_codigo,
cex = 0.7)
legend(
x = 1, y = -0.5,
legend = Tablacodigo_area$Categoria,
title = "Leyenda",
fill = colores_codigo,
cex = 0.7,
y.intersp = 0.7,
bty = "o",
inset = 0.2)
INDICADOR - Moda
moda_codigo <- names(TDFcodigo_area)[which.max(TDFcodigo_area)]
print(paste("MODA (Código de Área):", moda_codigo))
## [1] "MODA (Código de Área): 1"