0. Librerías

library(gt)
library(dplyr)
## 
## Adjuntando el paquete: 'dplyr'
## The following objects are masked from 'package:stats':
## 
##     filter, lag
## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union

1. Leer datos

datos <- read.csv(
  "waterPollution.csv",
  sep = ",",
  stringsAsFactors = FALSE
)

2. Extracción de la variable

# Variable original
datos$parameterSamplingPeriod

# Extraer el año inicial del período de muestreo
datos$Anio <- as.numeric(substr(datos$parameterSamplingPeriod, 1, 4))

# Agrupar los años en intervalos
datos$Intervalo_Muestreo <- cut(
  datos$Anio,
  breaks = c(1990, 1995, 2000, 2005, 2010, 2015, 2017),
  labels = c(
    "1991-1995",
    "1996-2000",
    "2001-2005",
    "2006-2010",
    "2011-2015",
    "2016-2017"
  ),
  include.lowest = TRUE
)

3. Tabla de distribución de frecuencias

# =========================
# TABLA DE FRECUENCIAS
# =========================

TDF_Periodo <- data.frame(
  table(datos$Intervalo_Muestreo)
)

ni <- TDF_Periodo$Freq

hi <- round((ni / sum(ni)) * 100, 2)

TDF_Periodo <- data.frame(
  Periodo = TDF_Periodo$Var1,
  ni,
  hi
)

# =========================
# FILA TOTAL
# =========================

Summary <- data.frame(
  Periodo = "TOTAL",
  ni = sum(ni),
  hi = 100
)

TDF_Periodo_suma <- rbind(
  TDF_Periodo,
  Summary
)

colnames(TDF_Periodo_suma) <- c(
  "Período de muestreo",
  "ni",
  "hi (%)"
)

# =========================
# TABLA FORMATEADA
# =========================

TDF_Periodo_suma %>%
  gt() %>%
  tab_header(
    title = md("Tabla Nro. 1"),
    subtitle = md("*Distribución de frecuencia del período de muestreo
                  en el estudio de la calidad del agua en 
                  Europa (1991-2017)*")
  ) %>%
  tab_source_note(
    source_note = md("Autor: Grupo 3")
  )
Tabla Nro. 1
Distribución de frecuencia del período de muestreo en el estudio de la calidad del agua en Europa (1991-2017)
Período de muestreo ni hi (%)
1991-1995 5 0.03
1996-2000 7 0.03
2001-2005 20 0.10
2006-2010 18290 91.45
2011-2015 1509 7.54
2016-2017 169 0.84
TOTAL 20000 100.00
Autor: Grupo 3

4. Gráficas

4.1. Diagrama de barras

# =========================
# GRÁFICA DE BARRAS
# Frecuencia absoluta (ni)
# =========================

barplot(
  TDF_Periodo$ni,
  main = "Gráfica N°1: Distribución de frecuencia del período de
  muestreo en el estudio de la calidad del agua en Europa (1991-2017)",
  xlab = "Período de muestreo",
  ylab = "Cantidad",
  col = "skyblue",
  ylim = c(0, max(TDF_Periodo$ni) * 1.1),
  names.arg = TDF_Periodo$Periodo,
  las = 2,
  cex.names = 0.6
)

4.2. Diagrama de barras general

# =========================
# GRÁFICA DE BARRAS (ni)
# =========================

barplot(
  ni,
  main = "Gráfica N°2: Distribución de frecuencia del período de
  muestreo en el estudio de la calidad del agua en Europa (1991-2017)",
  xlab = "Período de muestreo",
  ylab = "Cantidad",
  col = "red",
  ylim = c(0, max(ni) * 1.1),
  las = 2,
  cex.names = 0.6,
  names.arg = TDF_Periodo$Periodo
)

4.3. Diagrama de barras porcentual

# =========================
# GRÁFICA DE BARRAS
# =========================

barplot(
  TDF_Periodo$hi,
  main = "Gráfica N°3: Distribución porcentual del período de muestreo
  en el estudio de la calidad del agua en Europa (1991-2017)",
  xlab = "Período de muestreo",
  ylab = "Porcentaje",
  col = "skyblue",
  ylim = c(0, max(TDF_Periodo$hi) * 1.1),
  names.arg = TDF_Periodo$Periodo,
  las = 2,
  cex.names = 0.6
)

4.4. Diagrama de barras porcentual general

# =========================
# GRÁFICA DE BARRAS (hi)
# =========================

barplot(
  hi,
  main = "Gráfica N°4: Distribución porcentual del período de muestreo
  en el estudio de la calidad del agua en Europa (1991-2017)",
  xlab = "Período de muestreo",
  ylab = "Porcentaje",
  col = "red",
  ylim = c(0, 100),
  las = 2,
  cex.names = 0.6,
  names.arg = TDF_Periodo$Periodo
)

4.5. Diagrama circular

# =========================
# DIAGRAMA CIRCULAR
# =========================

etiquetas <- paste0(hi, " %")

colores <- rainbow(length(hi))

par(mar = c(2, 2, 4, 6))

pie(
  hi,
  labels = etiquetas,
  col = colores,
  main = "Gráfica N°5: Distribución del período de muestreo en
  el estudio de la calidad del agua en Europa (1991-2017)",
  cex = 0.7
)

legend(
  "topright",
  legend = TDF_Periodo$Periodo,
  fill = colores,
  title = "Períodos",
  cex = 0.5,
  xpd = TRUE
)

5. Indicadores Estadísticos

# ===========================================================
# INDICADORES ESTADÍSTICOS
# ===========================================================

library(dplyr)
library(gt)

tabla_indicadores <- data.frame(
  Variable = "Período de muestreo",
  Rango = "-",
  X = "-",
  Me = "-",
  Mo = names(sort(table(datos$Intervalo_Muestreo), decreasing = TRUE))[1],
  V = "-",
  Sd = "-",
  Cv = "-",
  As = "-",
  K = "-",
  Valores_Atipicos = "No aplica",
  stringsAsFactors = FALSE
)

tabla_indicadores_gt <- tabla_indicadores %>%
  gt() %>%
  tab_header(
    title = "Tabla Nro. 2",
    subtitle = "Indicadores estadísticos del período de muestreo en el estudio de la calidad del agua en Europa (1991-2017)"
  ) %>%
  cols_label(
    Variable = "Variable",
    Rango = "Rango",
    X = "X",
    Me = "Me",
    Mo = "Mo",
    V = "V",
    Sd = "Sd",
    Cv = "Cv",
    As = "As",
    K = "K",
    Valores_Atipicos = "Valores Atípicos"
  ) %>%
  tab_source_note(
    source_note = "Autor: Grupo 3"
  ) %>%
  tab_options(
    table.border.top.color = "black",
    table.border.bottom.color = "black",
    column_labels.border.bottom.color = "black",
    row.striping.include_table_body = TRUE,
    table.align = "center"
  )

tabla_indicadores_gt
Tabla Nro. 2
Indicadores estadísticos del período de muestreo en el estudio de la calidad del agua en Europa (1991-2017)
Variable Rango X Me Mo V Sd Cv As K Valores Atípicos
Período de muestreo - - - 2006-2010 - - - - - No aplica
Autor: Grupo 3

6. Conclusión

#El período de muestreo con mayor frecuencia corresponde al intervalo 2006–2010 y coincide con la mediana 2006–2010; por lo tanto, predominan los registros de ese período, dentro del rango de 1991–1995 a 2016–2017.