Sismos del Subcontinente Indio

Variable cualitativa nominal

Orientación

1 Cargar datos

Importamos el archivo “Orientacion (1).csv” desde una ruta local y usando espacios o tabulaciones como separador. Luego, la función str(Earthquakes) muestra la estructura del data frame, indicando tipos de datos y ejemplos de sus columnas.

Extraemos la variable

datos1 <- read.csv("C:/Users/HP/Desktop/Orientacion (1).csv", header = TRUE, sep = ",", dec = ".")
str(datos1)

## 'data.frame':    14776 obs. of  1 variable:
##  $ Orientacion: chr  "ENE" "NE" "S" "NNE" ...

2 Tabla de distribución de frecuencias

Orientacion <- datos1$Orientacion
TDFOrientacion <- as.data.frame(table(Orientacion))
TDFOrientacion

Agrupación

# Agrupación:
library(dplyr)

## Warning: package 'dplyr' was built under R version 4.4.3

## 
## Adjuntando el paquete: 'dplyr'

## The following objects are masked from 'package:stats':
## 
##     filter, lag

## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union

TDFOrientacion$variable <- ifelse(TDFOrientacion$Orientacion %in% c("N", "NNE", "NE", "NNW", "NW"), "Norte", 
                                  ifelse(TDFOrientacion$Orientacion %in% c("E", "ESE", "SE", "ENE"), "Este",
                                         ifelse(TDFOrientacion$Orientacion %in% c("S", "SSW", "SSE"), "Sur",
                                                ifelse(TDFOrientacion$Orientacion %in% c("W", "WNW", "WSW", "SW"), "Oeste", "otros"))))

head(TDFOrientacion)

TDFOrientacionFinal <- TDFOrientacion %>%
  group_by(variable) %>%
  summarise(
    ni = sum(Freq), 
    hi = round((sum(Freq) / sum(TDFOrientacion$Freq)) * 100)
  )

TDFOrientacionFinal <- data.frame(TDFOrientacionFinal)

Agregar totales

# Agregar Totales:
total_ni <- sum(TDFOrientacionFinal$ni)
total_hi <- 100

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

##   variable    ni  hi
## 1     Este  2759  19
## 2    Norte  5330  36
## 3    Oeste  4343  29
## 4      Sur  2344  16
## 5    Total 14776 100

library(gt)

## Warning: package 'gt' was built under R version 4.4.3

gt(TDFOrientacionFin) %>%
  tab_header(
    title = md("**Cuadro No. 2**"),
    subtitle = md("Distribución de Cantidad de Orientación de Propagación de los Sismos en el Subcontinente Indio (2000–2024)")
  ) %>%
  tab_source_note(
    source_note = md("Autores: Alvarez Jenny, Ibaza Nathali, Llangoma Sonya, Perez Cristian")
  ) %>%
  tab_style(
    style = list(
      cell_text(align = "center")
    ),
    locations = list(
      cells_title()
    )
  )

Cuadro No. 2
Distribución de Cantidad de Orientación de Propagación de los Sismos en el Subcontinente Indio (2000–2024)
variable	ni	hi
Este	2759	19
Norte	5330	36
Oeste	4343	29
Sur	2344	16
Total	14776	100
Autores: Alvarez Jenny, Ibaza Nathali, Llangoma Sonya, Perez Cristian

3 Histogramas

3.1 Histograma de frecuencia absoluta local (ni)

barplot(TDFOrientacionFinal$ni,
        main = "Gr\u00E1fica N\u00BA1: Distribuci\u00F3n de Cantidad seg\u00FAn la Orientaci\u00F3n en Sismos",
        xlab = "Orientaci\u00F3n", ylab = "Cantidad",
        col = "skyblue", names.arg = TDFOrientacionFinal$variable,
        las = 2, cex.names = 1)

3.2 Histograma de frecuencia relativa local (hi).

barplot(TDFOrientacionFinal$hi,
        main = "Gr\u00E1fica N\u00BA2: Distribuci\u00F3n en Porcentaje seg\u00FAn la Orientaci\u00F3n",
        xlab = "Orientaci\u00F3n", ylab = "Porcentaje (%)",
        col = "skyblue", names.arg = TDFOrientacionFinal$variable,
        las = 2, cex.names = 1)

4 Modelo de probabilidad

P_Ori <- TDFOrientacionFinal$hi
barplot(P_Ori,
        main = "Gr\u00E1fica N\u00BA3: Distribuci\u00F3n de Probabilidad de Orientaci\u00F3n",
        xlab = "Orientaci\u00F3n", ylab = "Probabilidad",
        col = "skyblue", names.arg = TDFOrientacionFinal$variable,
        cex.names = 1)

print("El modelo es Poisson")

## [1] "El modelo es Poisson"

5 Probabilidad

¿Cuál es la probabilidad de que un sismo tenga orientación hacia el Norte?

x <- round((TDFOrientacionFinal$ni[TDFOrientacionFinal$variable == "Norte"] / 14776) * 100, 1)
print(paste("La probabilidad es de:", x, "%"))

## [1] "La probabilidad es de: 36.1 %"