setwd ("/cloud/proyect")

importar datos

datos <- read.csv("Focos de Calor 2021.csv", header = T, sep = ",", dec = ".")

Estructura de los datos

str(datos)

Extraer la variable nominal

DPADESPRO <- datos$ DPADESPRO

EDA variable nominal

Tabla de Distribucion de Frecuencia

TDFDPADESPRO <- table(DPADESPRO) TDFDPA_DESPRO

TDFfinalDPADESPRO <- as.data.frame(TDFDPADESPRO) colnames(TDFfinalDPA_DESPRO) <- c("variable", "ni")

TablaDPADESPRO$ni sum(TablaDPADESPRO$ni)

hi <- TDFfinalDPADESPRO$ni / sum(TDFfinalDPADESPRO$ni) hi <- hi*100 sum(hi)

TDFfinalDPADESPRO$hi <- (TDFfinalDPADESPRO$ni / sum(TDFfinalDPADESPRO$ni)) * 100 TDFfinalDPADESPRO <- subset(TDFfinalDPADESPRO, variable != "Total") ni <- TDFfinalDPADESPRO$ni sum(ni)

variable <- TDFfinalDPADESPRO$variable TDFfinalDPADESPRO <- data.frame(variable, ni, hi) TDFfinalDPA_DESPRO

totalni <- sum(ni) totalhi <- sum(hi) TDFfinalDPADESPRO <- TDFfinalDPADESPRO[TDFfinalDPA_DESPRO$variable != "Total", ]

TDFfinalDPADESPRO <- rbind( TDFfinalDPADESPRO, data.frame(variable = "Total", ni = totalni, hi = totalhi) )

install.packages("gt") library(gt) install.packages("stringi") library(stringi)

if(any(!striencisutf8(TDFfinalDPA_DESPRO$variable))){ warning("Quedan cadenas con problemas de codificación") } else { message("Codificación UTF-8 correcta") }

library(gt)

tablaProvincia <- TDFfinalDPADESPRO %>% gt() %>% tabheader( title = md("Tabla Nro. 1.1"), subtitle = md("Tabla de distribución de frecuencias simples y acumuladas de Provincias del Ecuador") ) %>% tabsourcenote( sourcenote = md("Autor: Grupo 3") ) %>% taboptions( table.border.top.color = "black", table.border.bottom.color = "black", table.border.top.style = "solid", table.border.bottom.style = "solid", columnlabels.border.top.color = "black", columnlabels.border.bottom.color = "black", columnlabels.border.bottom.width = px(2), row.striping.includetablebody = TRUE, heading.border.bottom.color = "black", heading.border.bottom.width = px(2), tablebody.hlines.color = "gray", tablebody.border.bottom.color = "black" ) %>% tabstyle( style = celltext(weight = "bold"), locations = cells_body( rows = variable == "Total" ) )

tabla_Provincia

Tabla No.1

Distribuccion de provincias extendido

TablaDPA_DESPRO

GDF 1

TDFfinalDPADESPRO$variable <- iconv(TDFfinalDPADESPRO$variable, from = "latin1", to = "UTF-8", sub = "") barplot( height = TDFfinalDPADESPRO$ni, names.arg = TDFfinalDPADESPRO$variable, main = "Gráfica No. 1.1: Distribución de Provincias", xlab = "", ylab = "Cantidad", col = heat.colors(length(TDFfinalDPA_DESPRO$ni)), las = 2, cex.names = 0.7 ) mtext("Provincias", side = 1, line = 4, cex = 1)

DIAGRAMA DE BARRAS (GLOBAL)

colores <- c("yellow", "orange", "red")

barplot(TDFfinalDPA_DESPRO$ni, main = "Gráfica 1.2: Distribución de frecuencias (Global) de Provincias", xlab = "", ylab = "Cantidad",

    col = colores,
    names.arg = TDFfinalDPA_DESPRO$variable,
    ylim = c(0, 22500),
    las = 2,
    cex.names = 0.7

) mtext("Provincias", side = 1, line = 4, cex = 1)

DIAGRAMA DE BARRAS (Porcentaje)

barplot(TDFfinalDPADESPRO$hi, main = "Grafica No 1.3: Distribución de frecuencias (porcentaje) de Provincias", xlab = "", ylab = "Porcentaje (%)", col = colores, names.arg = TDFfinalDPADESPRO$variable, ylim = c(0, 30), las = 2, cex.names = 0.7 ) mtext("Provincias", side = 1, line = 4, cex = 1)

Crear agrupacion 1

if(!"ni" %in% colnames(TDFfinalDPADESPRO)) stop("No existe la columna 'ni' en TDFfinalDPADESPRO") TDFfinalDPA_DESPRO$Grupo <- NA

for(i in 1:nrow(TDFfinalDPADESPRO)) { if(TDFfinalDPADESPRO$ni[i] < 2000) { TDFfinalDPADESPRO$Grupo[i] <- "Zona Occidental" } else if(TDFfinalDPADESPRO$ni[i] >= 2000 && TDFfinalDPADESPRO$ni[i] <= 5000) { TDFfinalDPADESPRO$Grupo[i] <- "Zona Central Andina" } else { TDFfinalDPADESPRO$Grupo[i] <- "Zona Amazónica" } } print(TDFfinalDPADESPRO)

Grupo <- table(TablaDPADESPRO$Grupo) conteogrupos <- table(TDFfinalDPA_DESPRO$Grupo)

Gráfico de barras

barplot( conteo_grupos, main = "Gráfica No 1.4: Agrupación por frecuencia", xlab = "Agrupación de Provincias", ylab = "Cantidad", col = colores )

Frecuencias relativas (porcentaje)

hi_grupo <- prop.table(Grupo) * 100

Diagrama circular

colores <- rev(heat.colors(length(higrupo))) pie(higrupo, main = "Gráfica No. 1.5: Distribución de Frecuencia de provincias", col = colores)

legend("bottomright", legend = names(hi_grupo), title = "Leyenda", fill = colores, cex = 0.8, # tamaño del texto (más grande) border = "black")

INDICADORES

moda <- function(x) { ux <- unique(x) ux[which.max(tabulate(match(x, ux)))] }

Aplicar moda a DPA_DESPRO

modaDPA <- moda(datos$DPADESPRO) print(moda_DPA)