Mes_adquisicion
Monica Puchaicela
2025-06-03
MES_REPORT
UNIVERSIDAD CENTRAL DEL ECUADOR
PROYECTO: FOCOS DE CALOR EN EL ECUADOR
AUTORES: GUERRERO MARIA GABRIELA,PUCHAICELA MONICA, ZURITA JOHANNA
FECHA: 14/05/2025
setwd("/cloud/project")
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")
#Estructura de los datos
str(datos)## 'data.frame': 22476 obs. of 12 variables:
## $ ACQ_DATE : chr "20/11/2021" "20/11/2021" "06/08/2021" "10/06/2021" ...
## $ DPA_DESPRO: chr "ZAMORA CHINCHIPE" "ZAMORA CHINCHIPE" "ZAMORA CHINCHIPE" "ZAMORA CHINCHIPE" ...
## $ DPA_DESCAN: chr "CHINCHIPE" "CHINCHIPE" "CHINCHIPE" "CHINCHIPE" ...
## $ DPA_DESPAR: chr "CHITO" "CHITO" "PUCAPAMBA" "PUCAPAMBA" ...
## $ LATITUDE : num -4.98e+15 -4.97e+15 -4.96e+15 -4.96e+15 -4.96e+15 ...
## $ LONGITUDE : num -7.90e+16 -7.90e+16 -7.91e+16 -7.91e+16 -7.92e+16 ...
## $ BRIGHTNESS: num 3.55e+17 3.42e+17 3.32e+17 3.31e+17 3.28e+17 ...
## $ SCAN : num 5.1e+14 5.1e+14 1.5e+14 5.4e+14 5.0e+14 ...
## $ TRACK : num 4.9e+14 4.9e+14 3.8e+14 4.2e+14 4.9e+14 ...
## $ INSTRUMENT: chr "VIIRS" "VIIRS" "VIIRS" "VIIRS" ...
## $ VERSION : chr "2.0NRT" "2.0NRT" "2.0NRT" "2.0NRT" ...
## $ FRP : num 1.21e+16 6.87e+15 3.77e+15 5.50e+15 2.90e+15 ...#EDA variable Discreta
Fecha <- datos$ACQ_DATE
datos$Fecha <- as.Date(datos$ACQ_DATE, format = "%d/%m/%Y")
Variable2 <- format(datos$Fecha, "%m")
table(Variable2) ## Variable2
## 01 02 03 04 05 06 07 08 09 10 11 12
## 771 418 149 273 298 497 1120 1335 2176 3105 8266 4068head(Variable2)## [1] "11" "11" "08" "06" "05" "06"TDFMes <- as.data.frame(table(Variable2))
TDFMes## Variable2 Freq
## 1 01 771
## 2 02 418
## 3 03 149
## 4 04 273
## 5 05 298
## 6 06 497
## 7 07 1120
## 8 08 1335
## 9 09 2176
## 10 10 3105
## 11 11 8266
## 12 12 4068Meses <- as.numeric(Variable2)
ni_mes <- as.numeric(table(Variable2))
sum(ni_mes)## [1] 22476Mes <-TDFMes$Variable
TDFMesFin<-data.frame(Mes, ni_mes)
#DIAGRAMA DE BARRAS LOCAL ni
barplot(TDFMesFin$ni_mes, main = "Grafica Nro 9.1: Distribucion de Meses de Focos de calor",
xlab="Mes", ylab="Cantidad", names.arg =TDFMesFin$Mes, las=1, cex.names = 0.8
,col="purple") 
#INDICADORES
Variable_num9 <- as.numeric(Variable2)
x9<-sum(Variable_num9 )/length(Variable_num9)
x9## [1] 9.63312Me9<-median(Variable_num9)
Me9## [1] 11sd9<-sd(Variable_num9)
sd9## [1] 2.675167cv9<-sd9/x9*100
cv9## [1] 27.77051moda <- function(x) {
ux <- unique(x)
ux[which.max(tabulate(match(x, ux)))]
}
moda_val9 <- moda(Variable_num9)
print(moda_val9)## [1] 11install.packages("e1071")## Installing package into '/cloud/lib/x86_64-pc-linux-gnu-library/4.4'
## (as 'lib' is unspecified)library(e1071)
As9<-skewness(Variable_num9)
As9## [1] -1.811187K9<-kurtosis(Variable_num9)
K9## [1] 2.828461#Agrupado:#Agrupado:Dia
Mes <- as.numeric(as.character(TDFMes$Variable))
grupo_mes <- cut(Mes,
breaks = c(0, 4, 8, 12),
labels = c("MES 1–4", "MES 5–8", "MES 9–12"),
right = TRUE)
ni_mes2 <- TDFMes$Freq
TDFGruposMes <- aggregate(ni_mes ~ grupo_mes, FUN = sum)
print(TDFGruposMes)## grupo_mes ni_mes
## 1 MES 1–4 1611
## 2 MES 5–8 3250
## 3 MES 9–12 17615#Frecuencias simples
ni9 <- TDFGruposMes$ni_mes
ni9## [1] 1611 3250 17615sum(ni9)## [1] 22476hi9<-ni9/sum(ni9)*100
sum(hi9)## [1] 100#Frecuencias acumuladas
Ni9_asc<-cumsum(ni9)
Ni9_asc## [1] 1611 4861 22476Ni9_dsc<-rev(cumsum(rev(ni9)))
Ni9_dsc## [1] 22476 20865 17615Hi9_asc<-round(cumsum(hi9),4)
Hi9_asc## [1] 7.1676 21.6275 100.0000Hi9_dsc<-round(rev(cumsum(rev(hi9))),4)
Hi9_dsc## [1] 100.0000 92.8324 78.3725Tramos9 <- TDFGruposMes$grupo_mes
TDFGruposMesFin <- data.frame(
Tramos = Tramos9,
ni = ni9,
hi = hi9,
Ni_asc = Ni9_asc,
Ni_dsc = Ni9_dsc,
Hi_asc = Hi9_asc,
Hi_dsc = Hi9_dsc
)
print(TDFGruposMesFin)## Tramos ni hi Ni_asc Ni_dsc Hi_asc Hi_dsc
## 1 MES 1–4 1611 7.167645 1611 22476 7.1676 100.0000
## 2 MES 5–8 3250 14.459868 4861 20865 21.6275 92.8324
## 3 MES 9–12 17615 78.372486 22476 17615 100.0000 78.3725total_ni9 <- sum(ni9)
total_ni9## [1] 22476total_hi9 <- 100
total_hi9## [1] 100TDFGruposMesFinCompleto<-rbind(TDFGruposMesFin, data.frame(Tramos= "Total",
ni = total_ni9,
hi = total_hi9, Ni_asc = " ", Ni_dsc = " ",
Hi_asc = " ", Hi_dsc = " "))
print(TDFGruposMesFinCompleto) ## Tramos ni hi Ni_asc Ni_dsc Hi_asc Hi_dsc
## 1 MES 1–4 1611 7.167645 1611 22476 7.1676 100
## 2 MES 5–8 3250 14.459868 4861 20865 21.6275 92.8324
## 3 MES 9–12 17615 78.372486 22476 17615 100 78.3725
## 4 Total 22476 100.000000library(gt)
tabla_GruposMes <- TDFGruposMesFinCompleto %>%
gt() %>%
tab_header(
title = md("*Tabla Nro. 9*"),
subtitle = md("**Tabla de distribución de frecuencias simples y acumuladas de Meses 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_GruposMes| Tabla Nro. 9 | ||||||
| Tabla de distribución de frecuencias simples y acumuladas de Meses de focos de calor | ||||||
| Tramos | ni | hi | Ni_asc | Ni_dsc | Hi_asc | Hi_dsc |
|---|---|---|---|---|---|---|
| MES 1–4 | 1611 | 7.167645 | 1611 | 22476 | 7.1676 | 100 |
| MES 5–8 | 3250 | 14.459868 | 4861 | 20865 | 21.6275 | 92.8324 |
| MES 9–12 | 17615 | 78.372486 | 22476 | 17615 | 100 | 78.3725 |
| Total | 22476 | 100.000000 | ||||
| Autor: Grupo 3 | ||||||
#Graficas
#Diagrama de Barras:
barplot(
height = TDFGruposMesFin$ni,
main = "Grafica Nro 9.2: Distribucion de tramos de Meses de los focos de calor",
xlab = "Tramos",
ylab = "Cantidad",
names.arg = TDFGruposMesFin$Tramos,
las = 1,
cex.names = 0.9,
col = "purple"
)
#DIAGRAMA DE BARRAS GLOBAL ni
barplot(TDFGruposMesFin$ni, main = "Grafica Nro 9.3: Distribucion de tramos de focos
de calor ",
xlab="Tramos ", ylab="Cantidad", names.arg = TDFGruposMesFin$Tramos, las=1, col="purple",
ylim=c(0,22476))
#DIAGRAMA DE BARRAS LOCAL hi
barplot(TDFGruposMesFin$hi, main = "Grafica Nro 9.4: Distribucion de Tramos de Focos de calor",
xlab="Tramos", ylab="Porcentaje(%)", names.arg = TDFGruposMesFin$Tramos, las=1, col="purple") 
#DIAGRAMA DE BARRAS GLOBAL hi
barplot(TDFGruposMesFin$hi, main = "Grafica Nro 9.5: Distribucion de Tramos de focos de calor",
xlab="Tramos", ylab="Porcentaje(%)", names.arg =TDFGruposMesFin$Tramos, las=1, col="purple",
ylim=c(0,100))
#Graficar la ojiva ascendente (Hi_asc/dsc)
plot(x=1:length(TDFGruposMesFin$Tramos), y=TDFGruposMesFin$Hi_asc, type="p",
col="red", xlab="Tramos", ylab="Cantidad",
main ="Gráfica No 9.6: Ojivas Ascendentes y Descendentes de Frecuencias Absolutas
de Tramos", xaxt="n")
axis(side=1, at=1:length(TDFGruposMesFin$Hi_asc), labels=TDFGruposMesFin$Tramos, las=1)
points(x=1:length(TDFGruposMesFin$Tramos), y=TDFGruposMesFin$Hi_dsc, type="p",col="blue")
#DIAGRAMA DE CAJA Y BIGOTE
boxplot(Meses,
horizontal = TRUE,
col = "brown",
main = "Gráfica No 9.7: Diagrama de Caja de Meses de Focos de Calor",
xlab = "Mes del año")