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 Histograma de 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 Histograma de 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 Histograma de 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 Histograma de 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

La titularidad pública presenta la mayor frecuencia dentro de los pozos petrolíferos analizados, por lo que constituye la categoría predominante de la distribución. Esto evidencia una fuerte concentración de los registros en pozos de titularidad pública.