Códigos de el Análisis de las tendencias del suicidio y la demanda de servicios de salud mental en Puerto Rico (2020–2024)
Natalia Almeyda y Axel Pereles
2025-12-16
library(readr)
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, unionlibrary(knitr)
edad <- read.csv("Suicidios_Edad_2020_2024.csv")
pas <- read.csv("/cloud/project/Linea_PAS_2020_2024 (2).csv", header=FALSE, row.names=1)
genero <- read_csv("Suicidios_Genero_2020_2024.csv", show_col_types = FALSE)
mes <- read_csv("Suicidios_Mes_2020_2024.csv", show_col_types = FALSE)Bases de datos utilizadas
Para el desarrollo de este estudio se utilizaron las siguientes bases de datos oficiales, todas en formato agregado y correspondientes al periodo 2020–2024:
Suicidios por grupo de edad (2020–2024)
head(edad)## Grupo_de_edad edad_2020 edad_2021 edad_2022 edad_2023 edad_2024 edad_total
## 1 20-24 15 16 9 12 8 60
## 2 25-29 13 16 11 19 10 69
## 3 30-34 14 14 20 12 20 80
## 4 35-39 16 18 10 7 10 61
## 5 40-44 15 21 14 19 17 86
## 6 45-49 17 17 18 13 13 78
## promedio_total
## 1 12
## 2 14
## 3 16
## 4 12
## 5 17
## 6 16Esta base de datos presenta la distribución de los suicidios según grupos de edad, lo que permite identificar los rangos etarios con mayor incidencia durante el periodo de estudio.
Suicidios por mes (2020–2024)
head(mes)## # A tibble: 6 × 7
## Mes `2020` `2021` `2022` `2023` `2024` Promedio
## <chr> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
## 1 Enero 28 27 24 16 23 23.6
## 2 Febrero 13 23 19 11 5 14.2
## 3 Marzo 16 12 17 32 18 19
## 4 Abril 11 23 18 20 18 18
## 5 Mayo 13 13 15 20 21 16.4
## 6 Junio 21 12 13 20 20 17.2Contiene información mensual agregada sobre los suicidios ocurridos en Puerto Rico, facilitando el análisis de patrones temporales y posibles variaciones estacionales.
Suicidios por género (2020–2024)
head(genero)## # A tibble: 5 × 6
## Periodo Hombres Mujeres Total `% Hombres` `% Mujeres`
## <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
## 1 2020 176 35 211 83.4 16.6
## 2 2021 188 27 215 87.4 12.6
## 3 2022 169 26 195 86.7 13.3
## 4 2023 187 39 226 82.7 17.3
## 5 2024 176 28 204 86.3 13.7Incluye el número anual de suicidios desglosado por género, permitiendo examinar diferencias demográficas relevantes entre hombres y mujeres.
Llamadas relacionadas con conducta suicida – Línea PAS (2020–2024)
head(pas)## V2 V3 V4 V5 V6
## year 2020 2021 2022 2023 2024
## total_llamadas 922797 238228 225626 367606 444030
## llamadas_suicidas 15704 22008 13297 18493 18431Registra el número anual de llamadas relacionadas con conducta suicida atendidas por la Línea PAS, ofreciendo un indicador complementario sobre la demanda de servicios de apoyo en salud mental.
Todas las bases de datos contienen información agregada y no incluyen identificadores personales, garantizando la confidencialidad de los individuos y el cumplimiento de criterios éticos en la investigación.
Codigos y Resultados
Distribución de suicidios por género (2020–2024)
knitr::kable(
genero,
caption = "Tabla 1. Suicidios por año y género en Puerto Rico (2020–2024)"
)## Warning: 'xfun::attr()' is deprecated.
## Use 'xfun::attr2()' instead.
## See help("Deprecated")
## Warning: 'xfun::attr()' is deprecated.
## Use 'xfun::attr2()' instead.
## See help("Deprecated")| Periodo | Hombres | Mujeres | Total | % Hombres | % Mujeres |
|---|---|---|---|---|---|
| 2020 | 176 | 35 | 211 | 83.4 | 16.6 |
| 2021 | 188 | 27 | 215 | 87.4 | 12.6 |
| 2022 | 169 | 26 | 195 | 86.7 | 13.3 |
| 2023 | 187 | 39 | 226 | 82.7 | 17.3 |
| 2024 | 176 | 28 | 204 | 86.3 | 13.7 |
genero$Razon_HM <- genero$Hombres / genero$Mujeres
mean(genero$Razon_HM)## [1] 5.914424Los resultados evidencian una diferencia marcada y persistente entre hombres y mujeres en la ocurrencia de suicidios en Puerto Rico durante el periodo 2020–2024. En todos los años analizados, los hombres concentran más del 80 % de los casos, mientras que las mujeres representan una proporción considerablemente menor. En promedio, por cada mujer que muere por suicidio, aproximadamente seis hombres presentan el mismo desenlace, lo que confirma un patrón estructural y constante a lo largo del periodo estudiado.
Variabilidad estadística por género
medidas <- function(x){
c(
Minimo = min(x),
Maximo = max(x),
Mediana = median(x),
Varianza = var(x),
Desv_Estandar = sd(x),
Coef_Variacion = sd(x) / mean(x) * 100
)
}
round(sapply(genero[, c("Hombres", "Mujeres")], medidas), 2)## Hombres Mujeres
## Minimo 169.00 26.00
## Maximo 188.00 39.00
## Mediana 176.00 28.00
## Varianza 65.70 32.50
## Desv_Estandar 8.11 5.70
## Coef_Variacion 4.52 18.39Las medidas de dispersión muestran que, aunque los hombres presentan un mayor número absoluto de suicidios, su variabilidad interanual es relativamente baja. En contraste, las mujeres presentan un coeficiente de variación más alto, lo que indica fluctuaciones relativas mayores entre años. Esto sugiere que el comportamiento suicida masculino es más estable en el tiempo, mientras que el femenino muestra mayor sensibilidad a cambios contextuales.
Tendencia interanual de suicidios por género
plot(genero$Periodo, genero$Hombres, type = "o", pch = 16,
col = "black",
xlab = "Año",
ylab = "Casos",
ylim = c(min(genero$Hombres) - 5, max(genero$Hombres) + 5),
main = "Figura 1. Tendencia de suicidios por género en Puerto Rico (2020–2024)")
par(new = TRUE)
plot(genero$Periodo, genero$Mujeres, type = "o", pch = 17,
col = "red",
axes = FALSE,
xlab = "",
ylab = "",
ylim = c(min(genero$Mujeres) - 2, max(genero$Mujeres) + 2))
axis(side = 4, col = "red", col.axis = "red")
mtext("Casos (Mujeres)", side = 4, line = 3, col = "red")
legend("topright",
legend = c("Hombres", "Mujeres"),
col = c("black", "red"),
pch = c(16, 17),
lty = 1,
bty = "n")
La tendencia interanual muestra que los suicidios en hombres se mantienen elevados y relativamente estables durante todo el periodo analizado, mientras que los casos femeninos presentan variaciones menores en magnitud. A pesar de las fluctuaciones observadas, la brecha entre ambos géneros se mantiene constante, reforzando el género como un factor determinante en la distribución del suicidio en Puerto Rico.
Distribución de suicidios por grupo de edad (2020–2024)
par(mfrow = c(2, 3), mar = c(4, 4, 2, 1), oma = c(0, 0, 2, 0))
edad <- read_csv("Suicidios_Edad_2020_2024.csv", show_col_types = FALSE)
# 2020
barplot(edad$edad_2020, names.arg = edad$Grupo_de_edad, las = 2,
col = "lightgray", border = "black",
xlab = "Grupo de edad", ylab = "Suicidios", main = "2020")
# 2021
barplot(edad$edad_2021, names.arg = edad$Grupo_de_edad, las = 2,
col = "lightgray", border = "black",
xlab = "Grupo de edad", ylab = "Suicidios", main = "2021")
# 2022
barplot(edad$edad_2022, names.arg = edad$Grupo_de_edad, las = 2,
col = "lightgray", border = "black",
xlab = "Grupo de edad", ylab = "Suicidios", main = "2022")
# 2023
barplot(edad$edad_2023, names.arg = edad$Grupo_de_edad, las = 2,
col = "lightgray", border = "black",
xlab = "Grupo de edad", ylab = "Suicidios", main = "2023")
# 2024
barplot(edad$edad_2024, names.arg = edad$Grupo_de_edad, las = 2,
col = "lightgray", border = "black",
xlab = "Grupo de edad", ylab = "Suicidios", main = "2024")
barplot(edad$promedio_total, names.arg = edad$Grupo_de_edad, las = 2,
col = "lightgray", border = "black",
xlab = "Grupo de edad", ylab = "Suicidios", main = "Promedio")
mtext(
"Figura 2. Distribución de suicidios por grupo de edad (2020-2024)",
outer = TRUE,
cex = 1.1,
font = 2
)
Las gráficas evidencia que los mayores totales de suicidios durante el período 2020–2024 se concentran en los grupos de edad de 50 a 54, 55 a 59 y 60 a 64 años, los cuales presentan cifras acumuladas significativamente superiores al resto de los grupos etarios. A partir de los 30 años se observa un aumento progresivo en el número de casos, que alcanza su punto máximo en la adultez media y posteriormente desciende de forma gradual en los grupos de mayor edad. En contraste, los grupos de 20 a 24 y 25 a 29 años registran los valores más bajos del período, aunque mantienen una presencia constante de casos. Este patrón sugiere la influencia de factores asociados a la adultez y el envejecimiento, como condiciones de salud física y mental, aislamiento social y cambios en la estabilidad económica o familiar, al tiempo que evidencia que el suicidio es un fenómeno presente en todas las etapas del ciclo de vida y requiere intervenciones diferenciadas según la edad.
Tendencia del total anual de suicidios
mod_total <- lm(Total ~ Periodo, data = genero)
summary(mod_total)##
## Call:
## lm(formula = Total ~ Periodo, data = genero)
##
## Residuals:
## 1 2 3 4 5
## 0.2 4.5 -15.2 16.1 -5.6
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 816.80 8593.70 0.095 0.930
## Periodo -0.30 4.25 -0.071 0.948
##
## Residual standard error: 13.44 on 3 degrees of freedom
## Multiple R-squared: 0.001658, Adjusted R-squared: -0.3311
## F-statistic: 0.004982 on 1 and 3 DF, p-value: 0.9482plot(genero$Periodo, genero$Total, type = "o",
xlab = "Año",
ylab = "Total de suicidios",
main = "Figura 3. Tendencia del total anual de suicidios (2020–2024)")
abline(mod_total, col = "red", lwd = 2)
El modelo de regresión lineal no muestra evidencia de una tendencia creciente o decreciente significativa en el total anual de suicidios durante el periodo 2020–2024. Los resultados indican un comportamiento fluctuante alrededor de un promedio estable, sin una evolución progresiva clara en una sola dirección.
Análisis del comportamiento mensual mediante serie de tiempo
library(tidyr)
mes_largo <- mes %>%
pivot_longer(
cols = matches("2020|2021|2022|2023|2024"),
names_to = "Año",
values_to = "Casos"
) %>%
mutate(
Año = as.integer(gsub("\\D", "", Año)),
Mes_num = match(Mes, c("Enero","Febrero","Marzo","Abril","Mayo","Junio",
"Julio","Agosto","Septiembre","Octubre","Noviembre","Diciembre"))
) %>%
arrange(Año, Mes_num)
ts_mensual <- ts(
mes_largo$Casos,
start = c(2020, 1),
frequency = 12
)
plot(ts_mensual,
main = "Figura 4. Serie de tiempo mensual de suicidios en Puerto Rico (2020–2024)",
xlab = "Año",
ylab = "Casos")
La serie de tiempo mensual revela fluctuaciones recurrentes en la ocurrencia de suicidios a lo largo del año, sin una tendencia sostenida de aumento o disminución. Este patrón sugiere la posible influencia de factores estacionales o contextuales que varían según el mes, complementando el análisis anual del fenómeno.
Llamadas relacionadas con conducta suicida atendidas por la Línea PAS (2020–2024)
pas_normal <- read.csv("Linea_PAS_2020_2024 (2).csv", header = FALSE)
colnames(pas_normal) <- as.character(pas_normal[1, ])
pas <- pas_normal[-1, ]
pas_nuevo <- data.frame(
Año = as.numeric(colnames(pas)[-1]),
Total_llamadas = as.numeric(pas[pas[,1]=="total_llamadas", -1]),
Llamadas_suicidas = as.numeric(pas[pas[,1]=="llamadas_suicidas", -1])
)
pas_nuevo## Año Total_llamadas Llamadas_suicidas
## 1 2020 922797 15704
## 2 2021 238228 22008
## 3 2022 225626 13297
## 4 2023 367606 18493
## 5 2024 444030 18431plot(pas_nuevo$Año, pas_nuevo$Llamadas_suicidas,
type = "o",
pch = 16,
xlab = "Año",
ylab = "Llamadas relacionadas a suicidio",
main = "Figura 5. Llamadas relacionadas a suicidio – Línea PAS (2020–2024)")
Las llamadas relacionadas con conducta suicida atendidas por la Línea PAS se mantienen en niveles elevados durante todo el periodo 2020–2024. Aunque el volumen total de llamadas aumenta en los años más recientes, las llamadas específicamente asociadas a conducta suicida muestran variaciones moderadas. Estos resultados destacan el rol de la Línea PAS como un mecanismo de apoyo y prevención temprana dentro del sistema de salud mental en Puerto Rico.