```{r} # importar datos datos <- read.csv(“Focos de Calor 2021.csv”, header = T, sep = “,”, dec = “.”) datos <- read.csv(“Focos de Calor 2021.csv”, header = TRUE, sep = “,”, dec = “.”, fileEncoding = “latin1”)
```{r}
#Estructura de los datos
str(datos)
#Extraer la variable nominal
DPA_DESPRO <- datos$ DPA_DESPRO
#EDA variable nominal
#Tabla de Distribucion de Frecuencia
TDFDPA_DESPRO <- table(DPA_DESPRO)
TDFDPA_DESPRO
TDFfinalDPA_DESPRO <- as.data.frame(TDFDPA_DESPRO)
colnames(TDFfinalDPA_DESPRO) <- c("variable", "ni")
TDFfinalDPA_DESPRO$ni
sum(TDFfinalDPA_DESPRO$ni)
hi <- TDFfinalDPA_DESPRO$ni / sum(TDFfinalDPA_DESPRO$ni)
hi <- hi*100
sum(hi)
TDFfinalDPA_DESPRO$hi <- (TDFfinalDPA_DESPRO$ni / sum(TDFfinalDPA_DESPRO$ni)) * 100
TDFfinalDPA_DESPRO <- subset(TDFfinalDPA_DESPRO, variable != "Total")
ni <- TDFfinalDPA_DESPRO$ni
sum(ni)
```{r} variable <- TDFfinalDPA_DESPRO$variable TDFfinalDPA_DESPRO <- data.frame(variable, ni, hi) TDFfinalDPA_DESPRO
total_ni <- sum(ni) total_hi <- sum(hi) TDFfinalDPA_DESPRO <- TDFfinalDPA_DESPRO[TDFfinalDPA_DESPRO$variable != “Total”, ]
str(TDFfinalDPA_DESPRO) str(data.frame(variable = “Total”, ni = total_ni, hi = total_hi)) str(data.frame(variable = “Total”, ni = total_ni, hi = total_hi)) names(TDFfinalDPA_DESPRO)
nueva_fila <- data.frame( variable = “Total”, ni = total_ni, hi = total_hi, Grupo = NA # o puedes usar “” si prefieres dejarlo en blanco )
TDFfinalDPA_DESPRO <- rbind(TDFfinalDPA_DESPRO, nueva_fila)
install.packages(“gt”) library(gt) install.packages(“stringi”) library(stringi)
if(any(!stri_enc_isutf8(TDFfinalDPA_DESPRO$variable))){ warning(“Quedan cadenas con problemas de codificación”) } else { message(“Codificación UTF-8 correcta”) }
library(gt) TDFfinalDPA_DESPRO\(variable <- iconv(TDFfinalDPA_DESPRO\)variable, from = “latin1”, to = “UTF-8”) tabla_Provincia <- TDFfinalDPA_DESPRO %>% gt() %>% tab_header( title = md(“Tabla Nro. 1.1”), subtitle = md(“Tabla de distribución de frecuencias simples y acumuladas de Provincias del Ecuador”) ) %>% 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” ) %>% tab_style( style = cell_text(weight = “bold”), locations = cells_body( rows = variable == “Total” ) )
tabla_Provincia
## Including Plots
```{r pressure, echo=FALSE}
#Tabla No.1
#Distribuccion de provincias extendido
#GDF 1
TDFfinalDPA_DESPRO$variable <- iconv(TDFfinalDPA_DESPRO$variable, from = "latin1", to = "UTF-8", sub = "")
barplot(
height = TDFfinalDPA_DESPRO$ni,
names.arg = TDFfinalDPA_DESPRO$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)
```{r} # 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)
```{r}
# DIAGRAMA DE BARRAS (Porcentaje)
barplot(TDFfinalDPA_DESPRO$hi,
main = "Grafica No 1.3: Distribución de frecuencias (porcentaje) de
Provincias",
xlab = "",
ylab = "Porcentaje (%)",
col = colores,
names.arg = TDFfinalDPA_DESPRO$variable,
ylim = c(0, 30),
las = 2,
cex.names = 0.7
)
mtext("Provincias", side = 1, line = 4, cex = 1)
```{r} #Crear agrupacion 1
if(!“ni” %in% colnames(TDFfinalDPA_DESPRO)) stop(“No existe la columna ‘ni’ en TDFfinalDPA_DESPRO”) TDFfinalDPA_DESPRO$Grupo <- NA
for(i in 1:nrow(TDFfinalDPA_DESPRO)) { if(TDFfinalDPA_DESPRO\(ni[i] < 2000) { TDFfinalDPA_DESPRO\)Grupo[i] <- “Zona Occidental” } else if(TDFfinalDPA_DESPRO\(ni[i] >= 2000 && TDFfinalDPA_DESPRO\)ni[i] <= 5000) { TDFfinalDPA_DESPRO\(Grupo[i] <- "Zona Central Andina" } else { TDFfinalDPA_DESPRO\)Grupo[i] <- “Zona Amazónica” } } print(TDFfinalDPA_DESPRO) Grupo <- table(TDFfinalDPA_DESPRO\(Grupo) conteo_grupos <- table(TDFfinalDPA_DESPRO\)Grupo)
```{r}
# 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
```{r} # Diagrama circular
colores <- rev(heat.colors(length(hi_grupo))) pie(hi_grupo, 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”)
```