UNIVERSIDAD CENTRAL DEL ECUADOR

ESTUDIO ESTADÍSTICO DE LA CONTAMINACIÓN DEL SUELO Y SU IMPACTO EN LA SALUD

FECHA: 20/11/2025

#Estadistica Descriptiva

#19/11/2025

library(dplyr)
library(gt)

#Cargar los datos

setwd("/cloud/project")
datos<-read.csv("soil_pollution_diseases.csv",header = TRUE,dec = ".",
                sep = ",")

#Tablas Cualitativas Nominales

#Medida de Mitigación

Medida_Mitigacion<-datos$Mitigation_Measure

# Tabla de distribución de frecuencia

TDF_Medida_Mitigacion<-data.frame(table(Medida_Mitigacion))

ni <- TDF_Medida_Mitigacion$Freq
hi <- round((ni / sum(ni)) * 100, 2)
Medida_Mitigacion <- TDF_Medida_Mitigacion$Medida_Mitigacion
TDF_Medida_Mitigacion <- data.frame(Medida_Mitigacion, ni, hi)
Summary <- data.frame(Medida_Mitigacion = "TOTAL", ni = sum(ni),hi = 100)

TDF_Medida_Mitigacion_suma<-rbind(TDF_Medida_Mitigacion, Summary)
colnames(TDF_Medida_Mitigacion) <- c("Medida de Mitigación", "ni", "hi(%)")

# TABLA 

TDF_Medida_Mitigacion_suma %>%
  gt() %>%
  tab_header(
    title = md("*Tabla Nro. 10*"),
    subtitle = md("**Tabla de distribución
                  de las Medidas de Mitigación**")
  ) %>%
  tab_source_note(
    source_note = md("Autor: 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. 10
Tabla de distribución de las Medidas de Mitigación
Medida_Mitigacion ni hi
Chemical Neutralization 564 18.80
Community Awareness 592 19.73
Crop Rotation 626 20.87
Government Regulation 622 20.73
Soil Remediation 596 19.87
TOTAL 3000 100.00
Autor: 3
#Graficas

# Diagrama de barrras local ni

barplot(ni, main = "Gráfica N°46: Distribución 
        de las Medidas de Mitigación ",
        xlab = "Medidas de Mitigación",
        ylab = "Cantidad",
        col = "red",
        ylim = c(0,700),
        las = 1,
        cex.names = 0.5,
        names.arg = TDF_Medida_Mitigacion$`Medida de Mitigación`)

#Diagrama de barras global ni 

barplot(ni, main = "Gráfica N°47: Distribución
        de las Medidas de Mitigación",
        xlab = "Medidad de Mitigación",
        ylab = "Cantidad",
        col = "skyblue",
        ylim = c(0,3000),
        las=1,
        cex.names = 0.5,
        names.arg = TDF_Medida_Mitigacion$`Medida de Mitigación`)

#Diagrama de barras local hi(%)

barplot(hi, main = "Gráfica N°48: Distribución porcentual
        de las Medidas de Mitigación",
        xlab = "Medidas de Mitigación",
        ylab = "Porcentaje",
        col = "green",
        ylim = c(0,25),
        las = 1,
        cex.names = 0.5,
        names.arg = TDF_Medida_Mitigacion$`Medida de Mitigación`)

#Diagrama de barras global hi(%)

barplot(hi, main = "Gráfica N°49: Distribución porcentual 
        de las Medidas de Mitigación",
        xlab = "Medidas de Mitigación",
        ylab = "Porcentaje",
        col = "blue",
        ylim = c(0,100),
        las = 1,
        cex.names = 0.5,
        names.arg = TDF_Medida_Mitigacion$`Medida de Mitigación`)

#Diagrama circular

etiqueta_pie <- paste(TDF_Medida_Mitigacion$hi, "%")
pie(hi,
    main = "Gráfica N°50: Distribución porcentual
    de las Medidas de Mitigación",
    radius = 1,
    labels = etiqueta_pie,
    col = rev(heat.colors(length(hi))),
    cex = 0.9,
    cex.main = 1)

legend("topright",
       legend = TDF_Medida_Mitigacion$`Medida de Mitigación`,
       fill = colores <- c(rev(heat.colors(7))),
       cex = 0.5,
       title = "Leyenda")

# INDICADORES

# Función de moda
moda <- function(x) {
  frec <- table(x)         
  names(frec)[which.max(frec)]  
}

# Moda de la variable Medidas de Mitigación
moda(datos$Mitigation_Measure)
## [1] "Crop Rotation"