1. CARGA DE DATOS Y LIBRERÍAS
library(dplyr)
##
## 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
library(readxl)
library(gt)
datos <- read_excel("datos_deslizamientos.xlsx")
3. CONTEO
3.1 Consolidación de categorías
datos <- datos %>%
mutate(
landslide_setting_consol = case_when(
landslide_setting == "above_coast" ~ "Above coast",
landslide_setting == "above_river" ~ "Above river",
landslide_setting == "above_road" ~ "Above road",
landslide_setting == "below_road" ~ "Below road",
landslide_setting == "bluff" ~ "Bluff",
landslide_setting == "burned_area" ~ "Burned area",
landslide_setting == "deforested_slope" ~ "Deforested slope",
landslide_setting == "engineered_slope" ~ "Engineered slope",
landslide_setting == "mine" ~ "Mine",
landslide_setting == "natural_slope" ~ "Natural slope",
landslide_setting == "retaining_wall" ~ "Retaining wall",
landslide_setting == "urban" ~ "Urban",
landslide_setting %in% c("unknown", "Unknown", "", "NA") ~ NA_character_,
TRUE ~ NA_character_
)
)
3.2 Conteo de datos
variable <- datos$landslide_setting_consol
variable <- variable[!is.na(variable)]
N <- length(variable)
3.3 Cálculo de frecuencias
TDFSetting <- datos %>%
filter(!is.na(landslide_setting_consol)) %>%
count(landslide_setting_consol, name = "ni") %>%
mutate(
hi = round(ni / sum(ni), 2),
hi_porcentaje = round(hi * 100, 2)
)
TDFSetting_total <- TDFSetting %>%
add_row(
landslide_setting_consol = "TOTAL",
ni = sum(TDFSetting$ni),
hi = round(sum(TDFSetting$hi), 2),
hi_porcentaje = round(sum(TDFSetting$hi_porcentaje), 2)
)
5. GRÁFICAS
5.1 Distribución local del entorno de deslizamiento a nivel
mundial
par(mar = c(8, 5, 4, 2))
max_ni <- max(tabla_graficos$ni)
pos_x <- barplot(
tabla_graficos$ni,
names.arg = tabla_graficos$landslide_setting_consol,
col = "steelblue",
border = "black",
space = 0.2,
las = 2,
ylim = c(0, max_ni + max_ni*0.10),
yaxt = "n",
main = "Gráfica 1: Distribución local del entorno de\ndeslizamiento a nivel mundial",
xlab = "",
ylab = "Cantidad",
cex.names = 0.65
)
mtext("Entorno de deslizamiento", side = 1, line = 6, cex = 1)
ticks_y <- round(
seq(0, max_ni, length.out = 5),
0
)
axis(
side = 2,
at = ticks_y,
labels = ticks_y,
las = 1
)
text(
x = pos_x,
y = tabla_graficos$ni,
labels = tabla_graficos$ni,
pos = 3,
font = 2,
cex = 0.6
)

5.2 Distribución global del entorno de deslizamiento a nivel
mundial
par(mar = c(6, 5, 4, 2))
N_total <- sum(tabla_graficos$ni)
pos_x <- barplot(
tabla_graficos$ni,
names.arg = tabla_graficos$landslide_setting_consol,
col = "steelblue",
border = "black",
space = 0.2,
las = 2,
ylim = c(0, N_total),
yaxt = "n",
main = "Gráfica 2: Distribución global del entorno de\ndeslizamiento a nivel mundial",
xlab = "",
ylab = "Cantidad",
cex.names = 0.65
)
mtext("Entorno de deslizamiento", side = 1, line = 5, cex = 1)
ticks_y <- round(
seq(0, N_total, length.out = 6),
0
)
axis(
side = 2,
at = ticks_y,
labels = ticks_y,
las = 1
)
abline(
h = N_total,
col = "red",
lty = 2,
lwd = 2
)
text(
x = pos_x,
y = tabla_graficos$ni,
labels = tabla_graficos$ni,
pos = 3,
font = 2,
cex = 0.6
)

5.3 Distribución local en porcentaje del entorno de
deslizamiento a nivel mundial
par(mar = c(6, 5, 4, 2))
pos_x <- barplot(
tabla_graficos$hi_porcentaje,
names.arg = tabla_graficos$landslide_setting_consol,
col = "skyblue",
border = "black",
space = 0.2,
las = 2,
ylim = c(0, max(tabla_graficos$hi_porcentaje)),
yaxt = "n",
main = "Gráfica 3: Distribución local en porcentaje del\nentorno de deslizamiento a nivel mundial",
xlab = "",
ylab = "Porcentaje (%)",
cex.names = 0.65
)
mtext("Entorno de deslizamiento", side = 1, line = 5, cex = 1)
ticks_y <- pretty(
c(0, max(tabla_graficos$hi_porcentaje)),
n = 6
)
axis(
side = 2,
at = ticks_y,
labels = ticks_y,
las = 1
)
text(
x = pos_x,
y = tabla_graficos$hi_porcentaje,
labels = round(tabla_graficos$hi_porcentaje, 2),
pos = 3,
font = 2,
cex = 0.6
)

5.4 Distribución global en porcentaje del entorno de
deslizamiento a nivel mundial
par(mar = c(6, 5, 4, 2))
pos_x <- barplot(
tabla_graficos$hi_porcentaje,
names.arg = tabla_graficos$landslide_setting_consol,
col = "skyblue",
border = "black",
space = 0.2,
las = 2,
ylim = c(0, 100),
yaxt = "n",
main = "Gráfica 4: Distribución global en porcentaje del\nentorno de deslizamiento a nivel mundial",
xlab = "",
ylab = "Frecuencia relativa (%)",
cex.names = 0.65
)
mtext("Entorno de deslizamiento", side = 1, line = 5, cex = 1)
ticks_y <- seq(0, 100, by = 20)
axis(
side = 2,
at = ticks_y,
labels = ticks_y,
las = 1
)
abline(
h = 100,
col = "red",
lty = 2,
lwd = 2
)
text(
x = pos_x,
y = tabla_graficos$hi_porcentaje,
labels = round(tabla_graficos$hi_porcentaje, 2),
pos = 3,
font = 2,
cex = 0.6
)

5.5 Diagrama circular
par(mar = c(5, 4, 5, 8), xpd = TRUE)
colores <- rainbow(nrow(tabla_graficos))
pie(
tabla_graficos$hi_porcentaje,
labels = NA,
col = colores,
radius = 0.80,
main = "Gráfica 5. Distribución porcentual del entorno de\ndeslizamiento a nivel mundial"
)
legend(
x = 0.85,
y = 0,
legend = paste0(
tabla_graficos$landslide_setting_consol,
" (", round(tabla_graficos$hi_porcentaje,2), "%)"
),
fill = colores,
cex = 0.45,
bty = "n"
)

par(xpd = FALSE)