Titularidad de los Pozos Petroleros

1 Carga de Librerias

En esta sección se incorporan las librerías necesarias para facilitar la manipulación y el análisis de los datos.

setwd("C:/Users/Usuario/Desktop/Nueva carpeta")
library(readxl)
library(dplyr)
library(gt)
datos <- read_excel("tabela_de_pocos_janeiro_2018.xlsx")

2 Cargar datos

En esta sección se carga la base de datos que contiene la información de los pozos petrolíferos para su posterior análisis estadístico.

3 Extraer la variable

En esta sección se extrae la variable correspondiente a la titularidad de los pozos petrolíferos para su posterior análisis.

Titularidad <- datos$TITULARIDADE

TDFTitularidad <- as.data.frame(table(Titularidad))

TDFTitularidad
##    Titularidad  Freq
## 1 Confidencial   705
## 2      Público 28870

4 Conteo

Se calculan las frecuencias absolutas y relativas de los pozos de acuerdo con su titularidad.

TDFTitularidad$Freq <- as.numeric(TDFTitularidad$Freq)

TDFTitularidad1 <- TDFTitularidad %>%
  group_by(Titularidad) %>%
  summarise(
    ni = sum(Freq),
    hi = round(sum(Freq) / sum(TDFTitularidad$Freq) * 100, 5)
  )

TDFTitularidad1 <- data.frame(TDFTitularidad1)

TDFTitularidad1
##    Titularidad    ni       hi
## 1 Confidencial   705  2.38377
## 2      Público 28870 97.61623
TDFTitularidad1$fi <- TDFTitularidad1$ni / sum(TDFTitularidad1$ni)

TDFTitularidad1$hi <- TDFTitularidad1$fi * 100

TDFTitularidad1 <- TDFTitularidad1[, c(
  "Titularidad", "ni", "hi", "fi"
)]

total_ni <- sum(TDFTitularidad1$ni)
total_hi <- sum(TDFTitularidad1$hi)
total_fi <- sum(TDFTitularidad1$fi)

TDFTitularidadcompleta <- rbind(
  TDFTitularidad1,
  data.frame(
    Titularidad = "Total",
    ni = total_ni,
    hi = total_hi,
    fi = total_fi
  )
)

print(TDFTitularidadcompleta)
##    Titularidad    ni        hi        fi
## 1 Confidencial   705   2.38377 0.0238377
## 2      Público 28870  97.61623 0.9761623
## 3        Total 29575 100.00000 1.0000000

5 Tabla de frecuencia

Se presenta la distribución de frecuencias de los pozos petrolíferos según su titularidad.

gt(TDFTitularidadcompleta) %>%
  tab_header(
    title = md("**Tabla N°1: DISTRIBUCIÓN DE FRECUENCIAS DE POZOS PETROLEROS DE BRASIL**"),
    subtitle = "Titularidad de los pozos petrolíferos"
  ) %>%
  tab_spanner(
    label = md("**Frecuencia Relativa**"),
    columns = c(hi, fi)
  ) %>%
  cols_label(
    Titularidad = md("**Titularidad**"),
    ni = md("**ni**"),
    hi = md("Porcentual (%)"),
    fi = md("Fracción")
  ) %>%
  fmt_number(columns = hi, decimals = 2) %>%
  fmt_number(columns = fi, decimals = 4) %>%
  cols_align(
    align = "center",
    columns = everything()
  ) %>%
  tab_style(
    style = list(
      cell_fill(color = "#2E4053"),
      cell_text(color = "white", weight = "bold")
    ),
    locations = cells_title()
  ) %>%
  tab_style(
    style = list(
      cell_fill(color = "#F2F3F4"),
      cell_text(weight = "bold", color = "#2E4053")
    ),
    locations = cells_column_labels()
  ) %>%
  tab_style(
    style = list(
      cell_fill(color = "#2E4053"),
      cell_text(color = "white", weight = "bold")
    ),
    locations = cells_column_spanners()
  ) %>%
  tab_style(
    style = list(
      cell_fill(color = "#D5D8DC"),
      cell_text(weight = "bold", color = "#2E4053")
    ),
    locations = cells_body(rows = nrow(TDFTitularidadcompleta))
  ) %>%
  tab_options(
    table.border.top.color = "#2E4053",
    table.border.bottom.color = "#2E4053",
    column_labels.border.bottom.color = "#2E4053",
    data_row.padding = px(6),
    table.font.size = px(13)
  )
Tabla N°1: DISTRIBUCIÓN DE FRECUENCIAS DE POZOS PETROLEROS DE BRASIL
Titularidad de los pozos petrolíferos
Titularidad ni
Frecuencia Relativa
Porcentual (%) Fracción
Confidencial 705 2.38 0.0238
Público 28870 97.62 0.9762
Total 29575 100.00 1.0000
 
TDFTitularidadcompleta_ <- TDFTitularidadcompleta[
  TDFTitularidadcompleta$Titularidad != "Total",
]

6 Gráficas

6.1 Diagrama de barras frecuencia absoluta local

Se representa la cantidad de pozos correspondiente a cada categoría de titularidad.

barplot(
  TDFTitularidadcompleta_$ni,
  main = "Gráfica N°1: Distribución en cantidad según titularidad",
  xlab = "Titularidad",
  ylab = "Cantidad",
  col = "#2E4053",
  names.arg = TDFTitularidadcompleta_$Titularidad,
  las = 1,
  cex.names = 1,
  cex.axis = 0.8,
  cex.main = 1
)

6.2 Diagrama de barras frecuencia absoluta global

Se representa la cantidad de pozos según su titularidad utilizando una escala global.

barplot(
  TDFTitularidadcompleta_$ni,
  main = "Gráfica N°2: Distribución global según titularidad",
  xlab = "Titularidad",
  ylab = "Cantidad",
  col = "#2E4053",
  names.arg = TDFTitularidadcompleta_$Titularidad,
  las = 1,
  cex.names = 1,
  cex.axis = 0.8,
  cex.main = 1,
  ylim = c(0, sum(TDFTitularidadcompleta_$ni))
)

6.3 Diagrama de barras frecuencia relativa local

Se representa la distribución porcentual de los pozos según su titularidad.

barplot(
  TDFTitularidadcompleta_$hi,
  main = "Gráfica N°3: Distribución porcentual según titularidad",
  xlab = "Titularidad",
  ylab = "Porcentaje",
  col = "#2E4053",
  names.arg = TDFTitularidadcompleta_$Titularidad,
  las = 1,
  cex.names = 1,
  cex.axis = 0.8,
  cex.main = 1
)

6.4 Diagrama de barras frecuencia relativa global

Se representa la distribución porcentual utilizando una escala global de 0 a 100 %.

barplot(
  TDFTitularidadcompleta_$hi,
  main = "Gráfica N°4: Distribución porcentual global según titularidad",
  xlab = "Titularidad",
  ylab = "Porcentaje",
  col = "#2E4053",
  names.arg = TDFTitularidadcompleta_$Titularidad,
  las = 1,
  cex.names = 1,
  cex.axis = 0.8,
  cex.main = 1,
  ylim = c(0, 100)
)

6.5 Diagrama circular

Se representa la participación porcentual de cada categoría de titularidad.

pie(
  TDFTitularidadcompleta_$hi,
  main = "Gráfica N°5: Distribución porcentual según titularidad",
  radius = 0.9,
  labels = paste0(round(TDFTitularidadcompleta_$hi, 2), "%"),
  col = c("#264653", "#F4A261"),
  cex = 1,
  cex.main = 1,
  init.angle = 90
)

legend(
  x = -1.7,
  y = 1,
  legend = TDFTitularidadcompleta_$Titularidad,
  fill = c("#264653", "#F4A261"),
  cex = 1,
  title = "Titularidad"
)

7 Indicadores Estádisticos

moda_titularidad <- as.character(
  TDFTitularidad1$Titularidad[
    which.max(TDFTitularidad1$ni)
  ]
)

Conclusiones <- data.frame(
  Variable = "Titularidad",
  `Rango [Min; Max]` = "N/A",
  `Media (X̄)` = "N/A",
  `Mediana (Me)` = "N/A",
  `Moda (Mo)` = moda_titularidad,
  `Varianza (S²)` = "N/A",
  `Desv. Est. (S)` = "N/A",
  `C.V. (%)` = "N/A",
  `Asimetría (As)` = "N/A",
  `Curtosis (K)` = "N/A",
  `Valores Atípicos` = "N/A",
  check.names = FALSE
)

gt(Conclusiones) %>%
  tab_header(
    title = md("**CONCLUSIONES Y ESTADÍSTICOS**"),
    subtitle = "Resumen de Indicadores de la Titularidad de los Pozos Petrolíferos"
  ) %>%
  tab_source_note(
    source_note = "Autor: Caleb Yanez"
  ) %>%
  cols_align(
    align = "center",
    columns = everything()
  ) %>%
  tab_style(
    style = list(
      cell_fill(color = "#2E4053"),
      cell_text(color = "white", weight = "bold")
    ),
    locations = cells_title()
  ) %>%
  tab_style(
    style = list(
      cell_fill(color = "#F2F3F4"),
      cell_text(weight = "bold", color = "#2E4053")
    ),
    locations = cells_column_labels()
  ) %>%
  tab_options(
    table.border.top.color = "#2E4053",
    table.border.bottom.color = "#2E4053",
    column_labels.border.bottom.color = "#2E4053",
    data_row.padding = px(6)
  )
CONCLUSIONES Y ESTADÍSTICOS
Resumen de Indicadores de la Titularidad de los Pozos Petrolíferos
Variable Rango [Min; Max] Media (X̄) Mediana (Me) Moda (Mo) Varianza (S²) Desv. Est. (S) C.V. (%) Asimetría (As) Curtosis (K) Valores Atípicos
Titularidad N/A N/A N/A Público N/A N/A N/A N/A N/A N/A
Autor: Caleb Yanez

8 Conclusión

El valor más frecuente de Titularidad es Público.