1. 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

2. Leer datos

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

3. Selección de la variable

# =========================
# SELECCIÓN DE LA VARIABLE
# =========================

RUM <- datos$resultUom

4. Tabla de distribución de frecuencias

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

TDF_RUM <- data.frame(table(RUM))

# Frecuencia absoluta
ni <- TDF_RUM$Freq

# Frecuencia relativa porcentual
hi <- round((ni / sum(ni)) * 100, 2)

# Construcción de la tabla
TDF_RUM <- data.frame(
  RUM = TDF_RUM$RUM,
  ni,
  hi
)

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

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

TDF_RUM_suma <- rbind(TDF_RUM, Summary)

colnames(TDF_RUM_suma) <- c(
  "Unidad de medida",
  "ni",
  "hi (%)"
)

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

TDF_RUM_suma %>%
  gt() %>%
  tab_header(
    title = md("*Tabla Nro. 1*"),
    subtitle = md("**Distribución de frecuencias de las unidades de medida de los resultados en el estudio de la calidad del agua en Europa (1991-2017)**")
  ) %>%
  tab_source_note(
    source_note = md("Autor: Grupo 3")
  ) %>%
  tab_options(
    table.border.top.color = "black",
    table.border.bottom.color = "black",
    table.border.top.style = "solid",
    table.border.bottom.style = "solid",
    column_labels.border.top.color = "black",
    column_labels.border.bottom.color = "black",
    column_labels.border.bottom.width = px(2),
    row.striping.include_table_body = TRUE,
    heading.border.bottom.color = "black",
    heading.border.bottom.width = px(2),
    table_body.hlines.color = "gray",
    table_body.border.bottom.color = "black"
  )
Tabla Nro. 1
Distribución de frecuencias de las unidades de medida de los resultados en el estudio de la calidad del agua en Europa (1991-2017)
Unidad de medida ni hi (%)
% 1145 5.73
[pH] 1157 5.78
{massRatio} 1450 7.25
Cel 898 4.49
mg/L 113 0.56
mg{C}/L 558 2.79
mg{CaCO3}/L 9 0.04
mg{N}/L 3713 18.57
mg{NH3}/L 260 1.30
mg{NH4}/L 1110 5.55
mg{NO2}/L 1266 6.33
mg{NO3}/L 1289 6.44
mg{O2}/L 1609 8.04
mg{P}/L 2477 12.38
mg{PO4}/L 6 0.03
mg{Si}/L 378 1.89
mmol/L 154 0.77
ug/L 1670 8.35
uS/cm 738 3.69
TOTAL 20000 100.00
Autor: Grupo 3

5. Gráficas

5.1 Diagrama de barras de cantidad

# =========================================
# Diagrama de barras generado por RStudio
# =========================================

barplot(
  ni,
  main = "Gráfica N°1: Distribución de las unidades de medida de los resultados en el estudio de la calidad del agua en Europa (1991-2017)",
  xlab = "Unidad de medida",
  ylab = "Cantidad",
  col = "royalblue1",
  ylim = c(0, max(ni) * 1.1),
  las = 2,
  cex.names = 0.8,
  names.arg = TDF_RUM$RUM
)

5.2 Diagrama de barras general

# ===========================================================
# Diagrama de barras con relación a la totalidad de los datos
# ===========================================================

barplot(
  ni,
  main = "Gráfica N°2: Distribución de las unidades de medida de los resultados en el estudio de la calidad del agua en Europa (1991-2017)",
  xlab = "Unidad de medida",
  ylab = "Cantidad",
  col = "royalblue3",
  ylim = c(0, sum(ni)),
  las = 2,
  cex.names = 0.8,
  names.arg = TDF_RUM$RUM
)

5.3 Diagrama de barras porcentual

# ======================================
# Diagrama de barras porcentual
# ======================================

barplot(
  hi,
  main = "Gráfica N°3: Distribución porcentual de las unidades de medida de los resultados en el estudio de la calidad del agua en Europa (1991-2017)",
  xlab = "Unidad de medida",
  ylab = "Porcentaje (%)",
  col = "royalblue1",
  ylim = c(0, max(hi) * 1.05),
  las = 2,
  cex.names = 0.8,
  names.arg = TDF_RUM$RUM
)

5.4 Diagrama de barras porcentual con relación a la totalidad

# ===========================================================
# Diagrama de barras porcentual con relación a la totalidad
# ===========================================================

barplot(
  hi,
  main = "Gráfica N°4: Distribución porcentual de las unidades de medida de los resultados en el estudio de la calidad del agua en Europa (1991-2017)",
  xlab = "Unidad de medida",
  ylab = "Porcentaje (%)",
  col = "royalblue3",
  ylim = c(0, 100),
  las = 2,
  cex.names = 0.8,
  names.arg = TDF_RUM$RUM
)

5.5 Diagrama circular

# ======================================
# Diagrama circular con hi
# ======================================

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 de las unidades de medida de los resultados en el estudio de la calidad del agua en Europa (1991-2017)",
  cex = 0.7
)

legend(
  "topright",
  legend = TDF_RUM$RUM,
  fill = colores,
  title = "Unidad de medida",
  cex = 0.7,
  xpd = TRUE
)

6. Indicadores estadísticos

# ======================================
# Indicadores Estadísticos
# ======================================

# 1. CÁLCULO DE LA MODA (Mo)

TDF_solo_RUM <- TDF_RUM_suma[TDF_RUM_suma$`Unidad de medida` != "TOTAL", ]
Moda_row_RUM <- TDF_solo_RUM[which.max(TDF_solo_RUM$ni), ]
Mo_RUM <- as.character(Moda_row_RUM$`Unidad de medida`[1])

# 2. CREACIÓN DE LA TABLA DE INDICADORES

tabla_indicadores_RUM <- data.frame(
  "Variable" = "Unidad de medida",
  "Rango" = "-",
  "X" = "-",
  "Me" = "-",
  "Mo" = Mo_RUM,
  "V" = "-",
  "Sd" = "-",
  "Cv" = "-",
  "As" = "-",
  "K" = "-",
  "Valores Atípicos" = "-"
)

# 3. FORMATO DE LA TABLA CON GT

tabla_indicadores_RUM %>%
  gt() %>%
  tab_header(
    title = md("*Tabla Nro. 2*"),
    subtitle = md("**Indicadores estadísticos de las unidades de medida de los resultados en el estudio de la calidad del agua en Europa (1991-2017)**")
  ) %>%
  tab_source_note(
    source_note = md("Autor: Grupo 3")
  ) %>%
  tab_options(
    table.border.top.color = "black",
    table.border.bottom.color = "black",
    table.border.top.style = "solid",
    table.border.bottom.style = "solid",
    column_labels.border.top.color = "black",
    column_labels.border.bottom.color = "black",
    column_labels.border.bottom.width = px(2),
    row.striping.include_table_body = TRUE,
    heading.border.bottom.color = "black",
    heading.border.bottom.width = px(2),
    table_body.hlines.color = "gray",
    table_body.border.bottom.color = "black"
  )
Tabla Nro. 2
Indicadores estadísticos de las unidades de medida de los resultados 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
Unidad de medida - - - mg{N}/L - - - - - -
Autor: Grupo 3

7. Conclusión

#La unidad de medida con mayor frecuencia fue mg{N}/L, lo que indica que esta fue la unidad más utilizada para expresar los resultados en el estudio de la calidad del agua en Europa durante el período 1991–2017.