REPORTE DE DIA
Mónica Puchaicela
2025-06-02
Dia de Adquisición
UNIVERSIDAD CENTRAL DEL ECUADOR
PROYECTO: FOCOS DE CALOR EN EL ECUADOR
AUTORES: GUERRERO MARIA GABRIELA,PUCHAICELA MONICA, ZURITA JOHANNA
FECHA: 14/05/2025
datos <- read.csv("maate_focosdecalor_bdd_2021diciembre.csv",
header = T, sep = ",", dec = ".")
#EDA variable Discreta
Fecha <- datos$DIA_REPORT
datos$Fecha <- as.Date(datos$DIA_REPORT)
Variable <- format(datos$Fecha, "%d")
TDFDia <- as.data.frame(table(Variable))
dias <- as.numeric(Variable)
ni_dia <- TDFDia$Freq
sum(ni_dia)## [1] 22476Dia <-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] 17.17988Me<-median(Variable_num)
Me## [1] 17sd<-sd(Variable_num)
sd## [1] 8.770599cv<-sd/x*100
cv## [1] 51.05157moda <- function(x) {
ux <- unique(x)
ux[which.max(tabulate(match(x, ux)))]
}
moda_val <- moda(Variable_num)
print(moda_val)## [1] 17install.packages("e1071")## Installing package into '/cloud/lib/x86_64-pc-linux-gnu-library/4.4'
## (as 'lib' is unspecified)library(e1071)
As<-skewness(Variable_num)
As## [1] -0.1140943K<-kurtosis(Variable_num)
K## [1] -1.089742#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 6482
## 2 Días 11–20 6948
## 3 Días 21–31 9046#Frecuencias simples
ni <- TDFGruposDia$ni_dia
ni## [1] 6482 6948 9046sum(ni)## [1] 22476hi<-ni/sum(ni)*100
sum(hi)## [1] 100#Frecuencias acumuladas
Ni_asc<-cumsum(ni)
Ni_asc## [1] 6482 13430 22476Ni_dsc<-rev(cumsum(rev(ni)))
Ni_dsc## [1] 22476 15994 9046Hi_asc<-round(cumsum(hi),4)
Hi_asc## [1] 28.8397 59.7526 100.0000Hi_dsc<-round(rev(cumsum(rev(hi))),4)
Hi_dsc## [1] 100.0000 71.1603 40.2474Tramos1 <- 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 6482 28.83965 6482 22476 28.8397 100.0000
## 2 Días 11–20 6948 30.91297 13430 15994 59.7526 71.1603
## 3 Días 21–31 9046 40.24737 22476 9046 100.0000 40.2474total_ni <- sum(ni)
total_ni## [1] 22476total_hi <- 100
total_hi## [1] 100TDFGruposDiaFinCompleto<-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 6482 28.83965 6482 22476 28.8397 100
## 2 Días 11–20 6948 30.91297 13430 15994 59.7526 71.1603
## 3 Días 21–31 9046 40.24737 22476 9046 100 40.2474
## 4 Total 22476 100.00000library(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 | 6482 | 28.83965 | 6482 | 22476 | 28.8397 | 100 |
| Días 11–20 | 6948 | 30.91297 | 13430 | 15994 | 59.7526 | 71.1603 |
| Días 21–31 | 9046 | 40.24737 | 22476 | 9046 | 100 | 40.2474 |
| 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")