E.I_Variable_CodRegional.R

#UNIVERSIDAD CENTRAL DEL ECUADOR
#FACULTAD DE INGENIRÍA EN MINAS, PETRÓLEOS Y AMBIENTAL 
#TEMA:VARIABLE NOMINAL CÓDIGO REGIONAL 
#AUTOR: DOMÉNICA YEPEZ

#ESTADISTICA INFERENCIAL  
#DIRECTORIO DE TRABAJO 
setwd("C:/Users/HP/OneDrive - Universidad Central del Ecuador/SEMESTRE III/Estadistica/Incendios en Chile/Datos")
# CARGAR DATOS 
library(readxl)
datos <- read_excel("iNCENDIOS_FORESTALES_CHILE_DATOS.xlsx")
######################-ESTADÍSTICA INFERENCIAL-#################################
options(scipen = 999)
plot(1, type = "n", axes = FALSE, xlab = "", ylab = "") # Crear un gráfico vacío
text(x = 1, y = 1, 
     labels = "ESTADÍSTICA INFERENCIAL",
     cex = 2,  # Tamaño del texto (ajustable)
     col = "blue",  # Color del texto
     font = 6) #tipo

#Extraer la variable Aleatoria
cod_reg <-datos$codreg
#Frecuencia la Variable nominal 
TDFcod_reg<- table(cod_reg)
tablacod_reg<- as.data.frame(TDFcod_reg)
hi<- tablacod_reg$Freq/sum(tablacod_reg$Freq)
hi_porc<- hi*100
sum(hi)

## [1] 1

TDFcod_reg<- data.frame(tablacod_reg,round(hi_porc,2))
colnames(TDFcod_reg)<- c("Código Regional","ni","hi (%)")
TDFcod_reg

##    Código Regional   ni hi (%)
## 1               10   92   1.76
## 2               11   21   0.40
## 3               12   28   0.53
## 4               13  388   7.41
## 5               14   78   1.49
## 6               15    1   0.02
## 7                3   14   0.27
## 8                4  101   1.93
## 9                5  929  17.75
## 10               6  255   4.87
## 11               7  631  12.06
## 12               8 1946  37.18
## 13               9  750  14.33

Tabla_probabilidad<-data.frame(TDFcod_reg$`Código Regional`,TDFcod_reg$`hi (%)`)
colnames(Tabla_probabilidad)<-c("Código Regional","Probabilidad (%)")
Tabla_probabilidad

##    Código Regional Probabilidad (%)
## 1               10             1.76
## 2               11             0.40
## 3               12             0.53
## 4               13             7.41
## 5               14             1.49
## 6               15             0.02
## 7                3             0.27
## 8                4             1.93
## 9                5            17.75
## 10               6             4.87
## 11               7            12.06
## 12               8            37.18
## 13               9            14.33

#Grafica 
barplot(Tabla_probabilidad$`Probabilidad (%)`, main = "Gráfica N°1: Distribucion de Probabilidad de Cada Código Regional",
        ylab="Probabilidad (%)",xlab = "Código Regional",cex.names =0.8, names.arg = Tabla_probabilidad$`Código Regional`,
        col = "brown")

#calculo de probabilidad 
plot(1, type = "n", axes = FALSE, xlab = "", ylab = "") # Crear un gráfico vacío
text(x = 1, y = 1, 
     labels = "Cálculos de probabilidad\n(Estimación general)\n
     ¿Cuál es la probabilidad de que\nel código regional \n sea 8?\n
     R: 37.18 (%)",
     cex = 1.4,  # Tamaño del texto (ajustable)
     col = "blue",  # Color del texto
     font = 6) #tipo

# Crear vector de colores: todas marrón, excepto región 8 en rojo
colores <- ifelse(Tabla_probabilidad$`Código Regional` == 8, "red", "brown")

# Graficar todas las barras, resaltando la de la región 8
barplot(Tabla_probabilidad$`Probabilidad (%)`,
        main = "Gráfica N°1.1: Distribución de Probabilidad del código regional 8",
        ylab = "Probabilidad (%)",
        xlab = "Código Regional",
        cex.names = 0.8,
        names.arg = Tabla_probabilidad$`Código Regional`,
        col = colores, ylim = c(0,100))

# Agregar texto de probabilidad sobre la barra de la región 8
region_8_index <- which(Tabla_probabilidad$`Código Regional` == 8)
prob_8 <- Tabla_probabilidad$`Probabilidad (%)`[region_8_index]

# Añadir texto sobre la barra de región 8
text(x = region_8_index,
     y = prob_8 + 1,  # un poco arriba de la barra
     labels = paste0(prob_8, "%"),
     col = "red",
     font = 2,
     cex = 1.1)