SATELLITE

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 = ".")

#Estructura de los datos 
str(datos)
## 'data.frame':    22476 obs. of  17 variables:
##  $ MES_REPORT: int  11 11 8 6 5 6 11 9 3 3 ...
##  $ DIA_REPORT: int  20 20 6 10 28 10 20 29 22 22 ...
##  $ 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" ...
##  $ TXT_1     : chr  "PARROQUIA RURAL" "PARROQUIA RURAL" "PARROQUIA RURAL" "PARROQUIA RURAL" ...
##  $ LATITUDE  : chr  "-4,981720000000000" "-4,969160000000000" "-4,958520000000000" "-4,957820000000000" ...
##  $ LONGITUDE : chr  "-79,041280000000000" "-79,049490000000006" "-79,118430000000004" "-79,111859999999993" ...
##  $ BRIGHTNESS: chr  "354,759999999999990" "342,009999999999990" "331,860000000000010" "331,399999999999980" ...
##  $ SCAN      : chr  "0,510000000000000" "0,510000000000000" "0,150000000000000" "0,540000000000000" ...
##  $ TRACK     : chr  "0,490000000000000" "0,490000000000000" "0,380000000000000" "0,420000000000000" ...
##  $ SATELLITE : chr  "1" "1" "1" "1" ...
##  $ CONFIDENCE: chr  "n" "n" "n" "n" ...
##  $ VERSION   : chr  "2.0NRT" "2.0NRT" "2.0NRT" "2.0NRT" ...
##  $ BRIGHT_T31: chr  "299,420000000000020" "298,149999999999980" "299,160000000000030" "296,800000000000010" ...
##  $ FRP       : chr  "12,100000000000000" "6,870000000000000" "3,770000000000000" "5,500000000000000" ...
##  $ DAYNIGHT  : chr  "D" "D" "D" "D" ...
#Extraer la variable nominal
SATELLITE <- datos$SATELLITE


#EDA variable nominal
#Tabla de Distribucion de Frecuencia
TDFSATELLITE <- table(SATELLITE)
TDFSATELLITE
## SATELLITE
##     1     N 
## 11327 11149
TDFfinalSATELLITE <- as.data.frame(TDFSATELLITE)
colnames(TDFfinalSATELLITE) <- c("variable", "ni")


hi <- TDFfinalSATELLITE$ni / sum(TDFfinalSATELLITE$ni)
hi <- hi*100
sum(hi)
## [1] 100
TDFfinalSATELLITE$hi <- (TDFfinalSATELLITE$ni / sum(TDFfinalSATELLITE$ni)) * 100
TDFfinalSATELLITE <- subset(TDFfinalSATELLITE, variable != "Total")
ni <- TDFfinalSATELLITE$ni
sum(ni)
## [1] 22476
variable <- TDFfinalSATELLITE$variable
TDFfinalSATELLITE <- data.frame(variable, ni, hi)
TDFfinalSATELLITE
##   variable    ni       hi
## 1        1 11327 50.39598
## 2        N 11149 49.60402
total_ni <- sum(ni)
total_hi <- sum(hi)
TDFfinalSATELLITE <- TDFfinalSATELLITE[TDFfinalSATELLITE$variable != "Total", ]

TDFfinalSATELLITE<- rbind(
  TDFfinalSATELLITE,
  data.frame(variable = "Total", ni = total_ni, hi = total_hi)
)


library(gt)
library(stringi)

if(any(!stri_enc_isutf8(TDFfinalSATELLITE$variable))){
  warning("Quedan cadenas con problemas de codificación")
} else {
  message("Codificación UTF-8 correcta")
}
## Codificación UTF-8 correcta
library(gt)

tabla_Satelite <- TDFfinalSATELLITE %>%
  gt() %>%
  tab_header(
    title = md("Tabla Nro. 4"),
    subtitle = md("*Tabla de distribución de frecuencias de Tipo Satélite*")
  ) %>%
  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_Satelite
Tabla Nro. 4
Tabla de distribución de frecuencias de Tipo Satélite
variable ni hi
1 11327 50.39598
N 11149 49.60402
Total 22476 100.00000
Autor: Grupo 3 
#Tabla No.1
#Distribuccion de provincias extendido 
TDFfinalSATELLITE <- subset(TDFfinalSATELLITE, variable != "Total")

#GDF 1

TDFfinalSATELLITE$variable <- iconv(TDFfinalSATELLITE$variable, from = "latin1", to = "UTF-8", sub = "")
barplot(
  height = TDFfinalSATELLITE$ni,
  names.arg = TDFfinalSATELLITE$variable,
  main = "Gráfica N°4.1: Distribución de Tipo de Satélite",
  xlab = "",
  ylab = "Cantidad",
  col = heat.colors(length(TDFfinalSATELLITE$ni)),
  las = 2,
  cex.names = 0.8
)
mtext("Sátelites", side = 1, line = 2, cex = 1) 

# DIAGRAMA DE BARRAS (GLOBAL)

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

barplot(TDFfinalSATELLITE$ni,
        main = "Gráfica N°4.2: Distribución de frecuencias (Global) 
        de Satélites",
        xlab = "",
        ylab = "Cantidad",
        col = colores,
        names.arg = TDFfinalSATELLITE$variable,
        ylim = c(0, 22500),
        las = 2,
        cex.names = 0.7
)
mtext("Satélites", side = 1, line = 2, cex = 1) 

# DIAGRAMA DE BARRAS (Porcentaje)

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

barplot(TDFfinalSATELLITE$hi,
        main = "Grafica N°4.3: Distribución de Porcentaje
        de Satélites",
        xlab = "",
        ylab = "Porcentaje (%)",
        col = colores,
        names.arg = TDFfinalSATELLITE$variable,
        ylim = c(0, 50),
        las = 2,
        cex.names = 0.7
)
mtext("Sátelites", side = 1, line = 2, cex = 1) 

# DIAGRAMA DE BARRAS (Porcentaje)

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

barplot(TDFfinalSATELLITE$hi,
        main = "Grafica N°4.4: Distribución de Frecuencias (Porcentaje)
        de Satélites",
        xlab = "",
        ylab = "Porcentaje (%)",
        col = colores,
        names.arg = TDFfinalSATELLITE$variable,
        ylim = c(0, 100),
        las = 2,
        cex.names = 0.7
)
mtext("Sátelites", side = 1, line = 2, cex = 1) 

# Frecuencias relativas (porcentaje)
hi_grupo <- prop.table(hi) * 100

# Diagrama circular 
etiquetas <- paste0( names(hi_grupo),": ",
  round(hi_grupo, 1), "%")

colores <- rev(heat.colors(length(hi_grupo)))

pie(
  TDFfinalSATELLITE$hi,
  labels = etiquetas,
  main = "Gráfica N°2.5: Distribución Porcentual de Cantones",
  col = colores
)

legend(
  "bottomright",
  legend = TDFfinalSATELLITE$variable,
  title = "Satelites",
  fill = colores,
  cex = 1,
  border = "black"
)

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

# Aplicar moda a DPA_DESPRO
moda_SAT <- moda(datos$SATELLITE)
print(moda_SAT)
## [1] "1"