0.- Carga de Librerías

#Estadística descriptiva
#Variable cualitativa ordinal: Fecha de Estatuto Operativo (STAT_DATE)
#Autor: Luis Cruz
#Fecha: 2026-07-05

library(dplyr)
library(gt)

1.- Carga de Datos

datos <- read.csv("~/Estudio/TERCER SEMESTRE/Estadistica/Dataset.csv",
                   sep = ";", stringsAsFactors = FALSE)

2.- Seleccionar la variable

Se extrae el año de STAT_DATE, se filtra al rango 1990-2010 y se agrupa en 3 periodos para mantener consistencia metodológica entre fases del proyecto.

FCH_ESTATUTO <- trimws(as.character(datos$STAT_DATE))
FCH_ESTATUTO <- FCH_ESTATUTO[!is.na(FCH_ESTATUTO) & FCH_ESTATUTO != ""]

anios <- substr(FCH_ESTATUTO, 1, 4)
anios <- anios[as.numeric(anios) >= 1990 & as.numeric(anios) <= 2010]
n     <- length(anios)

periodo <- character(length(anios))
for (i in seq_along(anios)) {
  a <- as.numeric(anios[i])
  if (a >= 1990 & a <= 2001) {
    periodo[i] <- "1990-2001"
  } else if (a >= 2002 & a <= 2004) {
    periodo[i] <- "2002-2004"
  } else {
    periodo[i] <- "2005-2006"
  }
}

# ORDEN JERÁRQUICO (IMPORTANTE)
orden_periodo <- c("1990-2001", "2002-2004", "2005-2006")
periodo <- factor(periodo, levels = orden_periodo, ordered = TRUE)

datos_periodo <- data.frame(periodo)

3.- Frecuencia

# FRECUENCIA RESPETANDO ORDEN
TDF_Periodo <- datos_periodo %>%
  group_by(periodo) %>%
  summarise(fi = n()) %>%
  arrange(periodo)

# frecuencia relativa
TDF_Periodo$fri <- round((TDF_Periodo$fi / sum(TDF_Periodo$fi)) * 100, 2)

# TOTAL
Summary <- data.frame(
  periodo = "TOTAL",
  fi      = sum(TDF_Periodo$fi),
  fri     = 100
)

# TABLA FINAL
TDF_Periodo_suma <- rbind(
  data.frame(periodo = as.character(TDF_Periodo$periodo),
             fi = TDF_Periodo$fi, fri = TDF_Periodo$fri),
  Summary
)
colnames(TDF_Periodo_suma) <- c("Periodo", "fi", "fri(%)")

4.- Tabla de distribución de frecuencia

TDF_Periodo_suma %>%
  gt() %>%
  tab_header(
    title    = md("Tabla Nro. 1"),
    subtitle = md("Distribución de frecuencias de la variable Fecha de Estatuto Operativo (STAT_DATE) para el análisis de la actividad minera y emisión de gasesen Estados Unidos.")
  ) %>%
  cols_align(align = "center", columns = everything()) %>%
  tab_style(
    style = cell_text(weight = "bold"),
    locations = cells_body(rows = Periodo == "TOTAL")
  ) %>%
  tab_source_note(source_note = md("Autor: Luis Cruz"))

Periodo	fi	fri(%)
Tabla Nro. 1
Distribución de frecuencias de la variable Fecha de Estatuto Operativo (STAT_DATE) para el análisis de la actividad minera y emisión de gasesen Estados Unidos.
1990-2001	890	35.56
2002-2004	811	32.40
2005-2006	802	32.04
TOTAL	2503	100.00
Autor: Luis Cruz

5.- Gráfico de distribución de frecuencia

Diagramas de Cantidad

par(mar = c(6, 5, 4, 2))
barplot(
  TDF_Periodo$fi,
  names.arg = TDF_Periodo$periodo,
  col       = "#7030A0",
  main      = "Gráfica Nro. 1\nDistribución de FCH_ESTATUTO en instalaciones\nmineras de Estados Unidos (1990-2006)",
  ylab      = "Cantidad (fi)",
  cex.names = 0.9,
  cex.axis  = 0.9
)

par(mar = c(6, 5, 4, 2))
barplot(
  TDF_Periodo$fi,
  names.arg = TDF_Periodo$periodo,
  col       = "#7030A0",
  main      = "Gráfica Nro. 2\nDistribución de FCH_ESTATUTO en instalaciones\nmineras de Estados Unidos (1990-2006)",
  ylab      = "Cantidad (fi)",
  ylim      = c(0, max(TDF_Periodo$fi) + 500),
  cex.names = 0.9,
  cex.axis  = 0.9
)

Diagrama Porcentual

par(mar = c(6, 5, 4, 2))
barplot(
  TDF_Periodo$fri,
  names.arg = TDF_Periodo$periodo,
  col       = "#C5A3D8",
  main      = "Gráfica Nro. 3\nDistribución de frecuencia de FCH_ESTATUTO en instalaciones\nmineras de Estados Unidos (1990-2006)",
  ylab      = "Porcentaje",
  ylim      = c(0, max(TDF_Periodo$fri) + 5),
  cex.names = 0.9,
  cex.axis  = 0.9
)

par(mar = c(6, 5, 4, 2))
barplot(
  TDF_Periodo$fri,
  names.arg = TDF_Periodo$periodo,
  col       = "#C5A3D8",
  main      = "Gráfica Nro. 4\nDistribución de frecuencia de FCH_ESTATUTO en instalaciones\nmineras de Estados Unidos (1990-2006)",
  ylab      = "Porcentaje",
  ylim      = c(0, 100),
  cex.names = 0.9,
  cex.axis  = 0.9
)

Diagrama Circular

colores <- c("#7030A0","#C5A3D8","#4B0082")

par(mar = c(5, 4, 4, 14), xpd = TRUE)
pie(
  TDF_Periodo$fri,
  labels = paste0(TDF_Periodo$fri, " %"),
  col    = colores[1:nrow(TDF_Periodo)],
  radius = 0.8,
  main   = "Gráfica Nro. 5\nDistribución de FCH_ESTATUTO en instalaciones\nmineras de Estados Unidos (1990-2006)"
)
legend("right",
       inset  = c(-0.35, 0),
       legend = TDF_Periodo$periodo,
       fill   = colores[1:nrow(TDF_Periodo)],
       title  = "Leyenda",
       bty    = "o", xpd = TRUE, cex = 0.8)

par(xpd = FALSE)

6.- Indicadores Estadísticos

# 1. CÁLCULO DE LA MODA (Mo)
# Es el periodo con la mayor frecuencia (fi)
Moda_row <- TDF_Periodo[which.max(TDF_Periodo$fi), ]
Mo_calc  <- as.character(Moda_row$periodo[1])

# 2. CÁLCULO DE LA MEDIANA (Me)
# Frecuencias acumuladas
Ni_calc      <- cumsum(TDF_Periodo$fi)
N_total      <- sum(TDF_Periodo$fi)
Posicion_Me  <- N_total / 2
# Categoría que contiene la posición N/2
Mediana_row_index <- which(Ni_calc >= Posicion_Me)[1]
Me_calc <- as.character(TDF_Periodo$periodo[Mediana_row_index])

# RANGO (VARIABLE ORDINAL)
niveles <- levels(datos_periodo$periodo)
Rango_calc <- paste0(niveles[1], " a ", niveles[length(niveles)])
Rango <- Rango_calc

Variable <- "FCH_ESTATUTO (Periodo de Estatuto Operativo)"
Me <- Me_calc
X  <- "-"
Mo <- Mo_calc
sd <- "-"
CV <- "-"
As <- "-"
K  <- "-"

Tabla de Indicadores

Tabla_indicadores <- data.frame(Variable, Rango, X, Me, Mo, sd, CV, As, K)
colnames(Tabla_indicadores) <- c("Variable","Rango","X","Me","Mo","sd","CV","As","K")

Tabla_indicadores %>%
  gt() %>%
  tab_header(
    title    = md("Tabla Nro. 2"),
    subtitle = md("Indicadores estadísticos de la variable Fecha de Estatuto Operativo (STAT_DATE)")
  ) %>%
  cols_align(align = "center", columns = everything()) %>%
  tab_source_note(source_note = md("Autor: Luis Cruz"))

Variable	Rango	X	Me	Mo	sd	CV	As	K
Tabla Nro. 2
Indicadores estadísticos de la variable Fecha de Estatuto Operativo (STAT_DATE)
FCH_ESTATUTO (Periodo de Estatuto Operativo)	1990-2001 a 2005-2006	-	2002-2004	1990-2001	-	-	-	-
Autor: Luis Cruz

7.- Conclusión

En conclusión:

La variable FCH_ESTATUTO registra el periodo en que se determinó el estado operativo de cada instalación minera. El periodo con mayor número de actualizaciones es 1990-2001, lo que indica que ese periodo concentró la mayor actividad de cambios en el estatus operativo de las instalaciones mineras registradas. La mediana se ubica en 2002-2004, es decir, el 50% de las instalaciones tuvieron su último cambio de estado antes o durante ese periodo.

Variable Cualitativa Ordinal: Fecha de Estatuto Operativo (STAT_DATE)

Análisis de la Actividad Minera y Emisión de Gases

Luis Cruz

2026-07-05

0.- Carga de Librerías

1.- Carga de Datos

2.- Seleccionar la variable

3.- Frecuencia

4.- Tabla de distribución de frecuencia

5.- Gráfico de distribución de frecuencia

Diagramas de Cantidad

Diagrama Porcentual

Diagrama Circular

6.- Indicadores Estadísticos

Tabla de Indicadores

7.- Conclusión