Variable cualitativa nominal: Día de la semana

#------------------ UNIVERSIDAD CENTRAL DEL ECUADOR
##-----Facultad de Ingeniería en Geología, Minas, Petróleos y Ambiental
###-----------------Carrera de Ingeniería Ambiental

# Proyecto: Estudio de Contaminación por Plásticos a Nivel Mundial
# Etapa: Estadística Descriptiva
# Integrantes: Jhoao Cueva, Deisy Guasti, Marck Pillajo, María José Samaniego, Carlos Yugsi
# Fecha: 11/11/2026

# VARIABLE CUALITATIVA NOMINAL: DÍA DE LA SEMANA

# Carga de datos
ruta <- "C:/Users/HOME/Downloads/98631dc5bb9a4ea5a8f9c0b4ec433290_0.csv"

datos <- read.csv(ruta, stringsAsFactors = FALSE, strip.white = TRUE,
                  na.strings = c("", "NA"), fileEncoding = "UTF-8")
nrow(datos)

## [1] 54388

# Limpieza de datos: se eliminan filas desalineadas (sin día ni año válidos)
dias <- c("Monday","Tuesday","Wednesday","Thursday","Friday","Saturday","Sunday")
anio <- suppressWarnings(as.numeric(trimws(datos$Year)))
datos <- datos[trimws(datos$DOW) %in% dias & !is.na(anio) & anio >= 2015 & anio <= 2023, ]
nrow(datos)

## [1] 54388

# Día de la semana como factor ordenado (de lunes a domingo)
datos$DOW <- factor(trimws(datos$DOW), levels = dias)

# Extraer variable nominal
dia_semana <- datos$DOW
dia_semana <- droplevels(dia_semana[!is.na(dia_semana)])

# EDA variable nominal
TDFdia_semana <- table(dia_semana)
Tablafinal <- as.data.frame(TDFdia_semana, stringsAsFactors = FALSE)
names(Tablafinal) <- c("Categoria", "Freq")

ni <- Tablafinal$Freq
sum(ni)

## [1] 54388

hi <- (Tablafinal$Freq / sum(Tablafinal$Freq)) * 100
sum(hi)

## [1] 100

Tablafinal <- data.frame(Tablafinal, round(hi, 2))
names(Tablafinal) <- c("Día de la semana", "ni", "hi(%)")

# Creamos la fila del total
totales <- data.frame("Día de la semana" = "TOTAL", "ni" = sum(ni),
                      "hi(%)" = round(sum(hi), 2), check.names = FALSE)

Tablafinal_con_total <- rbind(Tablafinal, totales)
Tablafinal_con_total

##   Día de la semana    ni  hi(%)
## 1           Monday  4004   7.36
## 2          Tuesday  3721   6.84
## 3        Wednesday  3981   7.32
## 4         Thursday  6408  11.78
## 5           Friday  5612  10.32
## 6         Saturday 24608  45.25
## 7           Sunday  6054  11.13
## 8            TOTAL 54388 100.00

# Nombres cortos para el eje de las gráficas
nombres <- ifelse(nchar(Tablafinal[[1]]) > 25,
                  paste0(substr(Tablafinal[[1]], 1, 22), "..."), Tablafinal[[1]])

# GDF LOCAL ni
par(mar = c(9, 6, 4, 2))
barplot(Tablafinal$ni,
        main = "Gráfica No.1: Distribución de Cantidad del día de la semana
        del Estudio de Contaminación por Plásticos",
        ylab = "", col = "blue3",
        names.arg = nombres, las = 2, cex.names = 0.8)
mtext("Día de la semana", side = 1, line = 7)
mtext("Cantidad", side = 2, line = 4.5)

# GDF GLOBAL ni
par(mar = c(9, 6, 4, 2))
barplot(Tablafinal$ni,
        main = "Gráfica No.2: Distribución de Cantidad del día de la semana
        del Estudio de Contaminación por Plásticos",
        ylab = "", col = "purple3",
        names.arg = nombres, las = 2, cex.names = 0.8,
        ylim = c(0, sum(ni)))
mtext("Día de la semana", side = 1, line = 7)
mtext("Cantidad", side = 2, line = 4.5)

# GDF LOCAL hi
par(mar = c(9, 6, 4, 2))
barplot(Tablafinal$`hi(%)`,
        main = "Gráfica No.3: Distribución de Porcentaje del día de la semana
        del Estudio de Contaminación por Plásticos",
        ylab = "", col = "blue3",
        names.arg = nombres, las = 2, cex.names = 0.8)
mtext("Día de la semana", side = 1, line = 7)
mtext("Porcentaje (%)", side = 2, line = 4.5)

# GDF GLOBAL hi
par(mar = c(9, 6, 4, 2))
barplot(Tablafinal$`hi(%)`,
        main = "Gráfica No.4: Distribución de Porcentaje del día de la semana
        del Estudio de Contaminación por Plásticos",
        ylab = "", col = "purple3",
        names.arg = nombres, las = 2, cex.names = 0.8,
        ylim = c(0, 100))
mtext("Día de la semana", side = 1, line = 7)
mtext("Porcentaje (%)", side = 2, line = 4.5)

# DIAGRAMA CIRCULAR
par(mar = c(2, 2, 4, 2))
pie(
  Tablafinal$`hi(%)`,
  labels = ifelse(Tablafinal$`hi(%)` >= 3, paste(Tablafinal$`hi(%)`, "%"), ""),
  col = heat.colors(length(Tablafinal$`hi(%)`)),
  radius = 1,
  main = "Gráfica No.5: Distribución Porcentual del día de la semana
        del Estudio de Contaminación por Plásticos"
)
legend("topright",
       legend = paste0(nombres, " (", Tablafinal$`hi(%)`, "%)"),
       fill = heat.colors(length(Tablafinal$`hi(%)`)),
       title = "Leyenda",
       cex = 0.65)