ANÁLISIS ESTADÍSTICO DE VOLCANES ACTIVOS A NIVEL GLOBAL
VARIABLE NOMINAL: REGIÓN TECTÓNICA

CARRERA DE GEOLOGÍA

GRUPO N°2

ANÁLISIS ESTADÍSTICO INFERENCIAL

1. CARGA DE LIBRERIAS Y DATOS

2. TABLA DE DISTRIBUCIÓN DE PROBABILIDAD

3. GRÁFICA DE DISTRIBUCIÓN DE PROBABILIDAD

4. CALCULO DE PROBABILIDAD

5. CONCLUSION

CARGA DE LIBRERIAS Y DATOS

## Installing package into '/cloud/lib/x86_64-pc-linux-gnu-library/4.6'
## (as 'lib' is unspecified)
## Installing package into '/cloud/lib/x86_64-pc-linux-gnu-library/4.6'
## (as 'lib' is unspecified)
## Installing package into '/cloud/lib/x86_64-pc-linux-gnu-library/4.6'
## (as 'lib' is unspecified)
## 
## Attaching package: 'dplyr'
## The following objects are masked from 'package:stats':
## 
##     filter, lag
## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union
#CARGA DE DATASET
Volcanes_Globales <- read.csv("global_volcano_eruption_intelligence.csv", header = T, sep = ";", dec = ".")

#SELECCION VARIABLE
reg_tec <- Volcanes_Globales$tectonic_region

#LIMPIEZA DE DATOS
sum(is.na(reg_tec))
## [1] 0

TABLA DE DISTRIBUCIÓN DE PROBABILIDAD

CALCULO DE FRECUENCIAS 

# Tabla de frecuencias
TDFTectonica <- table(reg_tec)

# Convertir a dataframe
TDFTectonica <- as.data.frame(TDFTectonica)  

# Cambiar nombres de columnas
colnames(TDFTectonica) <- c("Grupo", "Freq")

# Calcular frecuencia absoluta (ni)
TDFTectonicaFinal <- TDFTectonica %>%
  group_by(Grupo) %>%
  summarise(
    ni = sum(Freq),
    
#Calcular Frecuencia relativa (hi)
    
 # hi porcentual 
  hi = (ni / sum(TDFTectonica$Freq)) * 100,
    
  # hi decimal
  hi_decimal = ni / sum(TDFTectonica$Freq)
  )

#Añadir fila de totales 
TDFTectonicaFinal <- TDFTectonicaFinal %>%
  add_row(
    Grupo = "Total",
    ni = sum(TDFTectonicaFinal$ni),
    hi = sum(TDFTectonicaFinal$hi),
    hi_decimal = sum(TDFTectonicaFinal$hi_decimal)
  )

# Redondear Frecuencias
TDFTectonicaFinal <- TDFTectonicaFinal %>%
  mutate(
    hi = round(hi, 2),
    hi_decimal = round(hi_decimal, 4)
  )

# Reparar categorias invalidas
TDFTectonicaFinal <- TDFTectonicaFinal %>%
  mutate(
    Grupo = iconv(Grupo, from = "CP1252", to = "UTF-8"),
    Grupo = gsub("–", "-", Grupo)
  )

# Ver tabla
TDFTectonicaFinal
## # A tibble: 11 × 4
##    Grupo                                 ni     hi hi_decimal
##    <chr>                              <int>  <dbl>      <dbl>
##  1 Antarctica                             1   0.11     0.0011
##  2 Canary Islands - Atlantic Hot Spot     5   0.56     0.0056
##  3 East African Rift                     28   3.12     0.0312
##  4 Intraplate / Other                    41   4.57     0.0457
##  5 Japan Arc                            149  16.6      0.166 
##  6 Mediterranean Volcanic Belt           71   7.91     0.0791
##  7 Mid-Atlantic Ridge - Iceland          54   6.01     0.0601
##  8 Philippine Arc                        48   5.35     0.0535
##  9 Ring of Fire - Pacific               279  31.1      0.311 
## 10 Sunda Arc - Indonesia                222  24.7      0.247 
## 11 Total                                898 100        1

TABLA DE FRECUENCIAS Y PROBABILIDAD

Distribución de Probabilidad de las Regiones Tectónicas
Análisis de frecuencias globales según la región tectónica de los volcanes
Región Tectónica Frecuencia Absoluta (ni) Frecuencia Relativa (hi %) Probabilidad (pi)
Antarctica 1 0.11 0.0011
Canary Islands - Atlantic Hot Spot 5 0.56 0.0056
East African Rift 28 3.12 0.0312
Intraplate / Other 41 4.57 0.0457
Japan Arc 149 16.59 0.1659
Mediterranean Volcanic Belt 71 7.91 0.0791
Mid-Atlantic Ridge - Iceland 54 6.01 0.0601
Philippine Arc 48 5.35 0.0535
Ring of Fire - Pacific 279 31.07 0.3107
Sunda Arc - Indonesia 222 24.72 0.2472
Total 898 100.00 1.0000

GRÁFICA DE DISTRIBUCIÓN DE PROBABILIDAD

DIAGRAMA DE DISTRIBUCIÓN DE FRECUENCIA RELATIVA 

# Eliminar fila total
TDFTectonicaPlot <- TDFTectonicaFinal[
  TDFTectonicaFinal$Grupo != "Total", ]

#FRECUENCIA RELATIVA LOCAL
par(mar = c(10, 6, 4, 2))

barplot(
  TDFTectonicaPlot$hi,
  
  main = "Gráfica Nº1: Frecuencia relativa de las Regiones Tectónicas 
  de los Volcanes a Nivel Mundial (local)",
  
  ylab = "Porcentaje (%)",
  
  col = "orange",
  
  names.arg = TDFTectonicaPlot$Grupo,
  
  cex.names = 0.8,
  
  las = 2
)

mtext(
  "Región tectónica",
  side = 1,
  line = 8
)

DIAGRAMA DE DENSIDAD DE PROBABILIDAD 

#GRAFICA DE DENSIDAD DE PROBABILIDAD
par(mar = c(10, 6, 4, 2))

barplot(
  TDFTectonicaPlot$hi,
  
  main = "Gráfica Nº2: Densidad de probabilidad de las Regiones Tectónicas 
  de los Volcanes activos a Nivel Mundial ",
  
  ylab = "Probabilidad %",
  
  col = "red",
  
  names.arg = TDFTectonicaPlot$Grupo,
  
  cex.names = 0.8,
  
  las = 2
)

mtext(
  "Región tectónica",
  side = 1,
  line = 8
)

CALCULO DE PROBABILIDAD

# 1. Calculamos el valor de la probabilidad dinámicamente
porcentaje <- round(
  (TDFTectonicaFinal$ni[TDFTectonicaFinal$Grupo == "Japan Arc"] /
     TDFTectonicaFinal$ni[TDFTectonicaFinal$Grupo == "Total"]) * 100,
  1
)

porcentaje
## [1] 16.6

CONCLUSION