#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 DESCRIPTIVA
#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")
#Extraer la variable
cod_reg <-datos$codreg
#Frecuencia la Variable nominal
TDFcod_reg<- table(cod_reg)
TDFcod_reg
## cod_reg
## 10 11 12 13 14 15 3 4 5 6 7 8 9
## 92 21 28 388 78 1 14 101 929 255 631 1946 750
tablacod_reg<- as.data.frame(TDFcod_reg)
tablacod_reg
## cod_reg Freq
## 1 10 92
## 2 11 21
## 3 12 28
## 4 13 388
## 5 14 78
## 6 15 1
## 7 3 14
## 8 4 101
## 9 5 929
## 10 6 255
## 11 7 631
## 12 8 1946
## 13 9 750
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))
TDFcod_reg
## cod_reg Freq round.hi_porc..2.
## 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
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
######################## Gráficas ###########################################
# Diagrama de barras
# Local
barplot(TDFcod_reg$ni, main = "Gráfica N°6- Distribucion de Cantidad de Cada Código Regional",
ylab="Cantidad",xlab = "Código Regional",cex.names =0.8, names.arg = TDFcod_reg$`Código Regional`,
col = "cyan")
# Global
barplot(TDFcod_reg$ni, main = "Gráfica N°7.- Distribucion de Cantidad de Cada Código Regional",
ylab="Cantidad",xlab = "Código Regional",cex.names =0.8, names.arg = TDFcod_reg$`Código Regional`,
col = "cyan", ylim = c(0,length(cod_reg)))
# Frecuencia relativa
barplot(TDFcod_reg$`hi (%)`,main = "Gráfica N°8.- Porcentaje de cada Código Regional", ylab="Porcentaje (%)",
xlab = "Código Regional",cex.names =0.8, names.arg = TDFcod_reg$`Código Regional`,
col = "salmon")
# Frecuencia absoluta
barplot(TDFcod_reg$`hi (%)`,main = "Gráfica N°9.- Porcentaje de cada Código Regional", ylab="Porcentaje (%)",
xlab = "Código Regional",cex.names =0.8, names.arg = TDFcod_reg$`Código Regional`,
col = "salmon", ylim = c(0,100))
# Agrupación de las regiones
for (i in 1:11) {
if (TDFcod_reg$`Código Regional`[i] %in% c(15, 1, 2, 3, 4)) {
TDFcod_reg$Grupo[i] <- "Norte"
} else if (TDFcod_reg$`Código Regional`[i] %in% c(5, 6, 7, 13)) {
TDFcod_reg$Grupo[i] <- "Centro"
} else if (TDFcod_reg$`Código Regional`[i] %in% c(8, 9, 10, 11, 12, 14, 16))
TDFcod_reg$Grupo[i] <- "Sur"
}
hi_grupo <- table(TDFcod_reg$Grupo) / sum(table(TDFcod_reg$Grupo))
hi_grupo <- round(hi_grupo * 100,2)
hi_grupo
##
## Centro Norte Sur
## 30.77 23.08 46.15
# Diagrama Circular de grupos
colores <- c("Norte" = "#87CEEB", "Centro" = "#98FB98", "Sur"="#FFA07A" )
pie(hi_grupo, labels = paste0( " (", hi_grupo, "%)"),
col = colores, main = "Grafica N°10.-Distribución de Porcentaje de las Regiones")
legend("bottomright", legend = unique(TDFcod_reg$Grupo), fill = colores, title = "Leyenda")
# Indicador estadístico
# Moda
mo<- which.max(TDFcod_reg$ni)
moda <- TDFcod_reg$`Código Regional`[mo]
moda
## [1] 8
## Levels: 10 11 12 13 14 15 3 4 5 6 7 8 9
# Tabla de indicadores
tabla_indicadores <- data.frame("Variable" =c("Regiones de Chile"),
"Rango" = "[1-13]",
"Intrumento"="ITU y SUSESO Gobierno de Chile",
"Mo" = "8")
library(knitr)
kable(tabla_indicadores, align = 'c', caption = "Conclusiones de la variable
código regional")
Conclusiones de la variable código regional
| Regiones de Chile |
[1-13] |
ITU y SUSESO Gobierno de Chile |
8 |