modo-del-pozo.R

setwd("/cloud/project/datos")
datos <- read.csv("Petroleo_Ontaro.csv", header=T, dec=".", sep=";")

# modo de pozo

c_modo <- datos$WELL_MODE
TDFc_modo <- table(c_modo)
TDFc_modo

## c_modo
##                               Abandoned and Junked (Lost) 
##                           31                          109 
##               Abandoned Well                  Active Well 
##                        15038                         3554 
##                  Capped Well                No Well Found 
##                            6                          814 
## Plugged back and whipstocked                    Potential 
##                           26                            7 
##               Suspended Well                      Unknown 
##                          873                         6208

Tabla <- as.data.frame(TDFc_modo)
hi <- Tabla$Freq/sum(Tabla$Freq)
hi

##  [1] 0.0011625291 0.0040876022 0.5639390985 0.1332783320 0.0002250056
##  [6] 0.0305257631 0.0009750244 0.0002625066 0.0327383185 0.2328058201

sum(hi)

## [1] 1

hi_porc <- hi*100
sum(hi_porc)

## [1] 100

Tabla_modo <- data.frame(Tabla, hi_porc)
total <- data.frame(c_modo="Total",
                    Freq=sum(Tabla_modo$Freq),
                    hi_porc="100%")

Tabla_modofinal <- rbind(Tabla_modo, total)

colnames(Tabla_modofinal)[colnames(Tabla_modofinal)=="c_modo"] <- "Modo de Pozo"
names(Tabla_modofinal)

## [1] "Modo de Pozo" "Freq"         "hi_porc"

colnames(Tabla_modofinal)[colnames(Tabla_modofinal)=="Freq"] <- "ni"
names(Tabla_modofinal)

## [1] "Modo de Pozo" "ni"           "hi_porc"

colnames(Tabla_modofinal)[colnames(Tabla_modofinal)=="hi_porc"] <- "hi"
names(Tabla_modofinal)

## [1] "Modo de Pozo" "ni"           "hi"

Tabla_modo <- Tabla_modo[order(-Tabla_modo$Freq), ]

Tabla_top6tipo <- head(Tabla_modo, 6)

barplot(Tabla_top6tipo$Freq,
        main = "Gráfica No3.01: Distribucion de Modo de Pozo",
        xlab = "Modo",
        ylab = "Cantidad",
        col = "aquamarine4",
        cex.names = 0.5,
        cex.axis = 0.9,
        names.arg = Tabla_top6tipo$c_modo,
        las = 2
)

barplot(Tabla_top6tipo$Freq,
        main = "Gráfica No3.02: Distribucion de Modo de Pozo (Global)",
        xlab = "Modo",
        ylab = "Cantidad",
        col = "aquamarine4",
        cex.names = 0.5,
        cex.axis = 0.9,
        names.arg = Tabla_top6tipo$c_modo,
        las = 2,
        ylim = c(0, 26666)
)

barplot(Tabla_top6tipo$hi_porc,
        main = "Gráfica No3.03: Distribucion de Modo de Pozo",
        xlab = "Modo",
        ylab = "Porcentaje",
        col = "aquamarine4",
        cex.names = 0.5,
        cex.axis = 0.9,
        names.arg = Tabla_top6tipo$c_modo,
        las = 2
)

barplot(Tabla_top6tipo$hi_porc,
        main = "Gráfica No3.04: Distribucion de Modo de Pozo (Global)",
        xlab = "Modo",
        ylab = "Porcentaje",
        col = "aquamarine4",
        cex.names = 0.5,
        cex.axis = 0.9,
        names.arg = Tabla_top6tipo$c_modo,
        las = 2,
        ylim = c(0, 100)
)

Colores <- colorRampPalette(rainbow(11))
etiquetas <- paste(round(Tabla_top6tipo$hi_porc),"%", sep = "")
pie(Tabla_top6tipo$hi_porc,
    radius = 0.5,
    col = Colores(10),
    main = "Gráfica No3.05: Distribucion de Modo de Pozo",
    labels = etiquetas)
legend("topright", title = "Leyenda", legend = Tabla_top6tipo$c_modo,
       cex = 0.5, fill = Colores(10), title.cex = 0.8)