UNIVERSIDAD CENTRAL DEL ECUADOR

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

FECHA: 24/12/2025

#Estadistica Descriptiva

#24/12/2025

#Cargar Datos
library(gt)
library(dplyr)

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

#Tablas Cualitativas 

#Textura Suelo

Textura_Suelo <- datos$Soil_Texture

#Tabla de Distribucion de Frecuencia

TDF_Textura_Suelo <-data.frame(table(Textura_Suelo))

ni <- TDF_Textura_Suelo$Freq
hi <- round((ni / sum(ni)) * 100, 2)
Pi <- hi 

Textura_Suelo <- TDF_Textura_Suelo$Textura_Suelo

TDF_Textura_Suelo <- data.frame(Textura_Suelo, ni, hi, Pi)

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

TDF_Textura_Suelo_suma <- rbind(TDF_Textura_Suelo, Summary)

colnames(TDF_Textura_Suelo_suma) <- c("Textura Suelo", "ni", "hi(%)", "Pi(%)")


# TABLA 

TDF_Textura_Suelo_suma %>%
  gt() %>%
  tab_header(
    title = md("*Tabla Nro. 1*"),
    subtitle = md("Tabla de distribución y probabilidades de las texturas de suelo")
  ) %>%
  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
Tabla de distribución y probabilidades de las texturas de suelo
Textura Suelo ni hi(%) Pi(%)
Clay 718 23.93 23.93
Loamy 770 25.67 25.67
Sandy 756 25.20 25.20
Silty 756 25.20 25.20
TOTAL 3000 100.00 100.00
Autor: Grupo 3
barplot(hi, main = "Gráfica N°1: Distribución porcentual de 
        las texturas de suelo",
        xlab = "Textura de Suelo",
        ylab = "Porcentaje",
        col = "green",
        ylim = c(0,26),
        las = 1,
        cex.names = 0.6,
        names.arg = TDF_Textura_Suelo$Textura_Suelo)

barplot(hi, main = "Gráfica N°2: Distribución para la probabilidad de 
         las texturas de suelo",
        xlab = "Textura de Suelo",
        ylab = "Probabilidad",
        col = "green",
        ylim = c(0,26),
        las = 1,
        cex.names = 0.6,
        names.arg = TDF_Textura_Suelo$Textura_Suelo)

# ¿ Cuál es la probabilidad de encontrar una textura arcillosa en el suelo?

prob_arcillosa <- TDF_Textura_Suelo$Pi[
  TDF_Textura_Suelo$Textura_Suelo == "Loamy"
]

cat(
  "¿Cuál es la probabilidad de encontrar una textura arcillosa en el suelo?\n",
  "La probabilidad de encontrar una textura arcillosa en el suelo es de",
  prob_arcillosa, "%.\n"
)
## ¿Cuál es la probabilidad de encontrar una textura arcillosa en el suelo?
##  La probabilidad de encontrar una textura arcillosa en el suelo es de 25.67 %.