FECHA: 14/05/2025
#EDA variable Discreta
datos <- read.csv("Focos de Calor 2021.csv",
header = T, sep = ",", dec = ".")
Fecha <- datos$ACQ_DATE
datos$Fecha <- as.Date(datos$ACQ_DATE, format = "%d/%m/%Y")
Variable <- format(datos$Fecha, "%d")
TDFDia <- as.data.frame(table(Variable))
dias <- as.numeric(Variable)
ni_dia <- TDFDia$Freq
sum(ni_dia)
## [1] 22476
Dia <-TDFDia$Variable
TDFDiaFin<-data.frame(Dia, ni_dia)
#DIAGRAMA DE BARRAS LOCAL ni
barplot(TDFDiaFin$ni_dia, main = "Grafica Nro 10.1: Distribucion de Dias de Focos de calor",
xlab="Dias", ylab="Cantidad", names.arg =TDFDiaFin$Dia, las=2, cex.names = 0.9 ,col="pink")
#INDICADORES
Variable_num <- as.numeric(Variable)
x<-sum(Variable_num )/length(Variable_num)
x
## [1] 16.89433
Me<-median(Variable_num)
Me
## [1] 17
sd<-sd(Variable_num)
sd
## [1] 8.685702
cv<-sd/x*100
cv
## [1] 51.41193
moda <- function(x) {
ux <- unique(x)
ux[which.max(tabulate(match(x, ux)))]
}
moda_val <- moda(Variable_num)
print(moda_val)
## [1] 16
install.packages("e1071")
library(e1071)
As<-skewness(Variable_num)
As
## [1] -0.1013227
K<-kurtosis(Variable_num)
K
## [1] -1.103433
#Agrupado:#Agrupado:Dia
Dia <- as.numeric(as.character(TDFDia$Variable))
grupo_dia <- cut(Dia,
breaks = c(0, 10, 20, 31),
labels = c("Días 1–10", "Días 11–20", "Días 21–31"),
right = TRUE)
ni_dia2 <- TDFDia$Freq
TDFGruposDia <- aggregate(ni_dia ~ grupo_dia, FUN = sum)
print(TDFGruposDia)
## grupo_dia ni_dia
## 1 Días 1–10 6342
## 2 Días 11–20 8085
## 3 Días 21–31 8049
#Frecuencias simples
ni <- TDFGruposDia$ni_dia
ni
## [1] 6342 8085 8049
sum(ni)
## [1] 22476
hi<-ni/sum(ni)*100
sum(hi)
## [1] 100
#Frecuencias acumuladas
Ni_asc<-cumsum(ni)
Ni_asc
## [1] 6342 14427 22476
Ni_dsc<-rev(cumsum(rev(ni)))
Ni_dsc
## [1] 22476 16134 8049
Hi_asc<-round(cumsum(hi),4)
Hi_asc
## [1] 28.2168 64.1885 100.0000
Hi_dsc<-round(rev(cumsum(rev(hi))),4)
Hi_dsc
## [1] 100.0000 71.7832 35.8115
Tramos1 <- TDFGruposDia$grupo_dia
TDFGruposDiaFin <- data.frame(
Tramos = Tramos1,
ni = ni,
hi = hi,
Ni_asc = Ni_asc,
Ni_dsc = Ni_dsc,
Hi_asc = Hi_asc,
Hi_dsc = Hi_dsc
)
print(TDFGruposDiaFin)
## Tramos ni hi Ni_asc Ni_dsc Hi_asc Hi_dsc
## 1 Días 1–10 6342 28.21676 6342 22476 28.2168 100.0000
## 2 Días 11–20 8085 35.97170 14427 16134 64.1885 71.7832
## 3 Días 21–31 8049 35.81153 22476 8049 100.0000 35.8115
total_ni <- sum(ni)
total_ni
## [1] 22476
total_hi <- 100
total_hi
## [1] 100
TDFGruposDiaFinCompleto<-rbind(TDFGruposDiaFin, data.frame(Tramos= "Total",
ni = total_ni,
hi = total_hi, Ni_asc = " ", Ni_dsc = " ",
Hi_asc = " ", Hi_dsc = " "))
print(TDFGruposDiaFinCompleto)
## Tramos ni hi Ni_asc Ni_dsc Hi_asc Hi_dsc
## 1 Días 1–10 6342 28.21676 6342 22476 28.2168 100
## 2 Días 11–20 8085 35.97170 14427 16134 64.1885 71.7832
## 3 Días 21–31 8049 35.81153 22476 8049 100 35.8115
## 4 Total 22476 100.00000
library(gt)
tabla_GruposDia <- TDFGruposDiaFinCompleto %>%
gt() %>%
tab_header(
title = md("*Tabla Nro. 10*"),
subtitle = md("**Tabla de distribución de frecuencias simples y acumuladas de Días de focos de calor**")
) %>%
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 = Tramos == "Total" # ✅ aquí está el cambio clave
)
)
tabla_GruposDia
| Tabla Nro. 10 |
| Tabla de distribución de frecuencias simples y acumuladas de Días de focos de calor |
| Tramos |
ni |
hi |
Ni_asc |
Ni_dsc |
Hi_asc |
Hi_dsc |
| Días 1–10 |
6342 |
28.21676 |
6342 |
22476 |
28.2168 |
100 |
| Días 11–20 |
8085 |
35.97170 |
14427 |
16134 |
64.1885 |
71.7832 |
| Días 21–31 |
8049 |
35.81153 |
22476 |
8049 |
100 |
35.8115 |
| Total |
22476 |
100.00000 |
|
|
|
|
| Autor: Grupo 3 |
#Graficas
#Diagrama de Barras:
barplot(
height = TDFGruposDiaFin$ni,
main = "Grafica Nro 10.2: Distribucion de tramos de dias de los focos de calor",
xlab = "Tramos",
ylab = "Cantidad",
names.arg = TDFGruposDiaFin$Tramos,
las = 1,
cex.names = 0.9,
col = "pink"
)
#DIAGRAMA DE BARRAS GLOBAL ni
barplot(TDFGruposDiaFin$ni, main = "Grafica Nro 10.3: Distribucion de tramos de focos
de calor ",
xlab="Tramos ", ylab="Cantidad", names.arg = TDFGruposDiaFin$Tramos, las=1, col="pink",
ylim=c(0,22476))
#DIAGRAMA DE BARRAS LOCAL hi
barplot(TDFGruposDiaFin$hi, main = "Grafica Nro 10.4: Distribucion de Tramos de Focos de calor",
xlab="Tramos", ylab="Porcentaje(%)", names.arg = TDFGruposDiaFin$Tramos, las=1, col="pink")
#DIAGRAMA DE BARRAS GLOBAL hi
barplot(TDFGruposDiaFin$hi, main = "Grafica Nro 10.5: Distribucion de Tramos de focos de calor",
xlab="Tramos", ylab="Porcentaje(%)", names.arg =TDFGruposDiaFin$Tramos, las=1, col="pink",
ylim=c(0,100))
#Graficar la ojiva ascendente (Hi_asc/dsc)
plot(x=1:length(TDFGruposDiaFin$Tramos), y=TDFGruposDiaFin$Hi_asc, type="p",
col="red", xlab="Tramos", ylab="Cantidad",
main ="Gráfica No 10.6: Ojivas Ascendentes y Descendentes de Frecuencias Absolutas
de Tramos", xaxt="n")
axis(side=1, at=1:length(TDFGruposDiaFin$Hi_asc), labels=TDFGruposDiaFin$Tramos, las=1)
points(x=1:length(TDFGruposDiaFin$Tramos), y=TDFGruposDiaFin$Hi_dsc, type="p",col="blue")
#DIAGRAMA DE CAJA Y BIGOTE
boxplot(dias,
horizontal = TRUE,
col = "brown",
main = "Gráfica No 10.7: Diagrama de Caja de Días de Focos de Calor",
xlab = "Día del mes")
