Memoria de cálculo

En el análisis de datos, la limpieza y modificación de bases de datos son esenciales para garantizar la calidad y precisión de los resultados. Este documento presenta el código en R utilizado para depurar y transformar algunas de las bases de datos utilizadas. El objetivo es proporcionar una visión clara del proceso de limpieza y las modificaciones realizadas, asegurando la integridad y utilidad de los datos para su análisis posterior.

Primero, cargaremos las librerías necesarias.

library(readxl)
library(openxlsx)
library(lubridate)
library(dplyr)
library(stringr)

También cargaremos las bases de datos correspondientes.

 CASOS <- read_excel("CASOS _SERV_AGR_OP - MAY-JUN-2024.xlsx", + sheet = "CASOS")
 AGRESORES <- read_excel("CASOS _SERV_AGR_OP - MAY-JUN-2024.xlsx", + sheet = "AGRESORES")
 ORDENES <- read_excel("CASOS _SERV_AGR_OP - MAY-JUN-2024.xlsx", + sheet = "OP")
 SERVICIOS <- read_excel("CASOS _SERV_AGR_OP - MAY-JUN-2024.xlsx", + sheet = "SERV")

se agregará a la base de casos y de agresores el día, mes y año correspondiente a la fecha hecho, se realizará la clasificación de la edad en un rango para ahora tener una variable categórica y además se hará reclasificación en los datos relacionados con escolaridad, situación conyugal y ocupación.

CASOS$dia_semana <- NA
CASOS$mes <- NA
CASOS$año <- NA



CASOS$rangoedad <- "Desconocida"
CASOS$edocivil <- "Desconocida"
CASOS$escolar <- "Desconocida"
CASOS$ocupacion2 <- "Desconocida"



for (x in 1:nrow(CASOS)) {

  
   if(!is.na(CASOS$Edad[i])){
    if(CASOS$Edad[i] <= 5){
      CASOS$rangoedad[i] <- "5 y menos"
    }
    if(CASOS$Edad[i] > 5 & CASOS$Edad[i] <= 14){
      CASOS$rangoedad[i] <- "6 a 14"
    }
    if(CASOS$Edad[i] > 14 & CASOS$Edad[i] <= 24){
      CASOS$rangoedad[i] <- "15 a 24"
    }
    if(CASOS$Edad[i] > 24 & CASOS$Edad[i] <= 44){
      CASOS$rangoedad[i] <- "25 a 44"
    }
    if(CASOS$Edad[i] > 44 & CASOS$Edad[i] <= 59){
      CASOS$rangoedad[i] <- "45 a 59"
    }
    if(CASOS$Edad[i] >= 60){
      CASOS$rangoedad[i] <- "60 y más"
    }
   }
  
  
  if(!is.na(CASOS$`Estado Civil`[i])){
    if(CASOS$`Estado Civil`[i] == "Casada" | CASOS$`Estado Civil`[i] == "Unión   libre" |
       CASOS$`Estado Civil`[i] == "Concubinato" | CASOS$`Estado Civil`[i] == "Sociedad de convivencia"){
      CASOS$edocivil[i] <- "Casada o unida"
    }
    if(CASOS$`Estado Civil`[i] == "Soltera"){
      CASOS$edocivil[i] <- "Soltera"
    }
    if(CASOS$`Estado Civil`[i] == "Divorciada" |
       CASOS$`Estado Civil`[i] == "Viuda" |
       CASOS$`Estado Civil`[i] == "Separada"){
      CASOS$edocivil[i] <- "Separada, divorciada o viuda"
    }
    if(CASOS$`Estado Civil`[i] == "No identificada" |
       CASOS$`Estado Civil`[i] == "Desconocido"){
      CASOS$edocivil[i] <- "Desconocida"
    }
  }
  
  
  
  if(!is.na(CASOS$Escolaridad[i])){
    if(CASOS$Escolaridad[i] == "Primaria"){
      CASOS$escolar[i] <- "Primaria"
    }
    if(CASOS$Escolaridad[i] == "No identificado" |
       CASOS$Escolaridad[i] == "NULL"){
      CASOS$escolar[i] <- "Desconocida"
    }
    if(CASOS$Escolaridad[i] == "Secundaria"){
      CASOS$escolar[i] <- "Secundaria"
    }
    if(CASOS$Escolaridad[i] == "Carrera técnica comercial" |
       CASOS$Escolaridad[i] == "Preparatoria"){
      CASOS$escolar[i] <- "Preparatoria"
    }
    if(CASOS$Escolaridad[i] == "Licenciatura" |
       CASOS$Escolaridad[i] == "Posgrado" |
       CASOS$Escolaridad[i] == "Maestría" |
       CASOS$Escolaridad[i] == "Doctorado"){
      CASOS$escolar[i] <- "Superior"
    }
    if(CASOS$Escolaridad[i] == "Preescolar" |
       CASOS$Escolaridad[i] == "Estudios que no requieren validéz oficial"|
       CASOS$Escolaridad[i] == "Ninguna"){
      CASOS$escolar[i] <- "Ninguna"
    }

  }
  
  if(!is.na(CASOS$VictSeDesconoce[i]) | !is.na(CASOS$VictOtro[i]) |
     !is.na(CASOS$VictActIlicita[i]) | !is.na(CASOS$VictNinguna[i])){
    if(CASOS$VictSeDesconoce[i] == 1 | CASOS$VictOtro[i]== 1 |
       CASOS$VictActIlicita[i] == 1 | CASOS$VictNinguna[i] == 1){
      CASOS$ocupacion2[i] <- "Otra"
    }
  }
  if(!is.na(CASOS$VictTrabajaHogar[i]) | !is.na(CASOS$VictTrabajaFHogar[i])){
    if(CASOS$VictTrabajaHogar[i] == 1 | CASOS$VictTrabajaFHogar[i] == 1){
      CASOS$ocupacion2[i] <- "Trabaja"
    }
  }
  if(!is.na(CASOS$VictEstudia[i])){
    if(CASOS$VictEstudia[i] == 1){
      CASOS$ocupacion2[i] <- "Estudia"
    }
  }
  if(!is.na(CASOS$VictJubiladaPensionada[i]) | !is.na(CASOS$VictPensionada[i])){
    if(CASOS$VictJubiladaPensionada[i] == 1 | CASOS$VictPensionada[i] == 1){
      CASOS$ocupacion2[i] <- "Jubilada pensionada"
    }
  }
  
  

  fecha_numeros <- as.numeric(CASOS$`Fecha hechos`[x])
  base<- as.Date("1899-12-30")
  fecha<- base+fecha_numeros

  fecha_hora_objeto <- as.POSIXct(fecha, format = "%d-%m-%Y")

  # Obtén el día de la semana
  CASOS$dia_semana[x] <- weekdays(fecha_hora_objeto)

  # Obtén el mes y el año
  CASOS$mes[x] <- format(fecha_hora_objeto, "%B") # Obtiene el nombre del mes completo
  CASOS$año[x] <- format(fecha_hora_objeto, "%Y")

  
    
  CASOS$Dia <- NA
  dias <- c("2.Lun", "3.Mar", "4.Mie", "5.Jue", "6.Vie", "7.Sab", "1.Dom")
  dias_com <- c("lunes", "martes", "miércoles", "jueves", "viernes", "sábado", "domingo")
  for (i in 1:nrow(CASOS)) {
    for (j in 1:length(dias_com)){
      if(!is.na(CASOS$dia_semana[i]) & CASOS$dia_semana[i]==dias_com[j]){
        CASOS$Dia[i] <- dias[j]
        break
      }
    }
  }
  




}

AGRESORES$dia_semana <- NA
AGRESORES$mes <- NA
AGRESORES$año <- NA




AGRESORES$rangoedad <- NA
AGRESORES$estadocivil <- "Desconocida"
AGRESORES$escolar <- "Desconocida"
AGRESORES$ocupacion2 <- NA

for (x in 1:nrow(AGRESORES)) {
  fecha_numeros <- as.numeric(AGRESORES$fecha_hechos[x])
  base<- as.Date("1899-12-30")
  fecha<- base+fecha_numeros

  fecha_hora_objeto <- as.POSIXct(fecha, format = "%d-%m-%Y")

  # Obtén el día de la semana
  AGRESORES$dia_semana[x] <- weekdays(fecha_hora_objeto)

  # Obtén el mes y el año
  AGRESORES$mes[x] <- format(fecha_hora_objeto, "%B") # Obtiene el nombre del mes completo
  AGRESORES$año[x] <- format(fecha_hora_objeto, "%Y")
  
  

  
  if(!is.na(AGRESORES$`Edad del Agresor`[i])){
    if(AGRESORES$`Edad del Agresor`[i] <= 5){
      AGRESORES$rangoedad[i] <- "5 y menos"
    }
    if(AGRESORES$`Edad del Agresor`[i] > 5 & AGRESORES$`Edad del Agresor`[i] <= 14){
      AGRESORES$rangoedad[i] <- "6 a 14"
    }
    if(AGRESORES$`Edad del Agresor`[i] > 14 & AGRESORES$`Edad del Agresor`[i] <= 24){
      AGRESORES$rangoedad[i] <- "15 a 24"
    }
    if(AGRESORES$`Edad del Agresor`[i] > 24 & AGRESORES$`Edad del Agresor`[i] <= 44){
      AGRESORES$rangoedad[i] <- "25 a 44"
    }
    if(AGRESORES$`Edad del Agresor`[i] > 44 & AGRESORES$`Edad del Agresor`[i] <= 59){
      AGRESORES$rangoedad[i] <- "45 a 59"
    }
    if(AGRESORES$`Edad del Agresor`[i] >= 60){
      AGRESORES$rangoedad[i] <- "60 y más"
    }
  }
  
  
  if(!is.na(AGRESORES$EdoCivil[i])){
    if(AGRESORES$EdoCivil[i] == "Casado" | AGRESORES$EdoCivil[i] == "Unión libre" |
       AGRESORES$EdoCivil[i] == "Concubinato" | AGRESORES$EdoCivil[i] == "Sociedad de convivencia"){
      AGRESORES$estadocivil[i] <- "Casado(a) o unido(a)"
    }
    if(AGRESORES$EdoCivil[i] == "Soltero"){
      AGRESORES$estadocivil[i] <- "Soltero(a)"
    }
    if(AGRESORES$EdoCivil[i] == "Divorciado" |
       AGRESORES$EdoCivil[i] == "Viudo" |
       AGRESORES$EdoCivil[i] == "Separado"){
      AGRESORES$estadocivil[i] <- "Separado(a), divorciado(a) o viudo(a)"
    }
    if(AGRESORES$EdoCivil[i] == "No identificado" |
       AGRESORES$EdoCivil[i] == "NULL"|
       AGRESORES$EdoCivil[i] == "Seleccione"){
      AGRESORES$estadocivil[i] <- "Desconocida"
    }
  }
  
  
  
  if(!is.na(AGRESORES$Escolaridad[i])){
    if(AGRESORES$Escolaridad[i] == "Primaria"){
      AGRESORES$escolar[i] <- "Primaria"
    }
    if(AGRESORES$Escolaridad[i] == "No identificado" |
       AGRESORES$Escolaridad[i] == "NULL"){
      AGRESORES$escolar[i] <- "Desconocida"
    }
    if(AGRESORES$Escolaridad[i] == "Secundaria"){
      AGRESORES$escolar[i] <- "Secundaria"
    }
    if(AGRESORES$Escolaridad[i] == "Carrera técnica comercial" |
       AGRESORES$Escolaridad[i] == "Preparatoria"){
      AGRESORES$escolar[i] <- "Preparatoria"
    }
    if(AGRESORES$Escolaridad[i] == "Licenciatura" |
       AGRESORES$Escolaridad[i] == "Posgrado" |
       AGRESORES$Escolaridad[i] == "Maestría" |
       AGRESORES$Escolaridad[i] == "Doctorado"){
      AGRESORES$escolar[i] <- "Superior"
    }
    if(AGRESORES$Escolaridad[i] == "Preescolar" |
       AGRESORES$Escolaridad[i] == "Estudios que no requieren validéz oficial"|
       AGRESORES$Escolaridad[i] == "Ninguna"){
      AGRESORES$escolar[i] <- "Ninguna"
    }

  }
  
  
  if(!is.na(AGRESORES$PrinActAgSinDato[i]) | !is.na(AGRESORES$PrinActAgSeDesc[i]) |
     !is.na(AGRESORES$PrinActAgOtro[i]) | !is.na(AGRESORES$PrinActAgNing[i])|
     !is.na(AGRESORES$PrinActAgActIlicita[i])){
    if(AGRESORES$PrinActAgSinDato[i] == 1 | AGRESORES$PrinActAgSeDesc[i] == 1 |
       AGRESORES$PrinActAgOtro[i] == 1 | AGRESORES$PrinActAgNing[i] == 1 |
       AGRESORES$PrinActAgActIlicita[i] == 1 ){
      AGRESORES$ocupacion2[i] <- "Otra"
    }
  }
  if(!is.na(AGRESORES$PrinActAgTrabHogar[i]) | !is.na(AGRESORES$PrinActAgTrabFHogar[i])){
    if(AGRESORES$PrinActAgTrabHogar[i] == 1 | AGRESORES$PrinActAgTrabFHogar[i] == 1){
      AGRESORES$ocupacion2[i] <- "Trabaja"
    }
  }
  if(!is.na(AGRESORES$PrinActAgEstudia[i])){
    if(AGRESORES$PrinActAgEstudia[i] == 1){
      AGRESORES$ocupacion2[i] <- "Estudia"
    }
  }
  if(!is.na(AGRESORES$PrinActAgJubPens[i]) | !is.na(AGRESORES$PrinActAgPens[i])){
    if(AGRESORES$PrinActAgJubPens[i] == 1 | AGRESORES$PrinActAgPens[i] == 1){
      AGRESORES$ocupacion2[i] <- "Jubilado(a) pensionado(a)"
    }
  }
  
  
}

De igual manera se realizó una nueva clasificación a los datos relacionados con el vínculo con agresor y además se trabajo con modalidad de violencia, a la cual se realiza limpieza y a partir de ella, obtenemos una nueva columna de datos, la cual nos indica cuando un caso llegó a ser feminicidio.

CASOS$Modalidad <- NA
CASOS$feminicidios <- 0
CASOS$"Vínculo(original)" <- CASOS$`Detalle de Vínculo con victima`
CASOS$"Habla lengua indigena" <- 0
familiares <- c("Hijo(a)", "Padrastro o madrastra", "Sobrino(a)", "Abuelo(a)",
                "Tío(a)", "Primo(a)", "Hermano(a)", "Suegro(a)", "Nieto(a)")



for(x in 1:nrow(CASOS)){
  if(!is.na(CASOS$`Modalidad de la Violencia`[x])){
    if(CASOS$`Modalidad de la Violencia`[x] == "1 - Familiar"){
     CASOS$Modalidad[x] <- "Familiar"
    }
    if(CASOS$`Modalidad de la Violencia`[x] == "2 - Laboral"){
     CASOS$Modalidad[x] <- "Laboral"
    }
    if(CASOS$`Modalidad de la Violencia`[x] == "3 - Institucional"){
     CASOS$Modalidad[x] <-  "Institucional"
    }
    if(CASOS$`Modalidad de la Violencia`[x] == "4 - En la comunidad" |
       CASOS$`Modalidad de la Violencia`[x] == "Enlacomunidad"){
     CASOS$Modalidad[x] <- "Comunitaria"
    }
    if(CASOS$`Modalidad de la Violencia`[x] == "5 - Escolar/Docente"){
     CASOS$Modalidad[x] <- "Escolar"
    }
    if(CASOS$`Modalidad de la Violencia`[x] ==  "6 - Obstétrica"){
     CASOS$Modalidad[x] <- "Obstétrica"
    }
    if(CASOS$`Modalidad de la Violencia`[x] == "7 - Política"){
     CASOS$Modalidad[x] <- "Política"
    }
    if(CASOS$`Modalidad de la Violencia`[x] == "9 - Feminicida"){
     CASOS$Modalidad[x] <- "Feminicida"
    }
    if(CASOS$`Modalidad de la Violencia`[x] == "8 - Digital/Cibernética"){
     CASOS$Modalidad[x] <- "Digital"
    }
  }
  
  
  if(CASOS$Modalidad[i] == "Feminicida"){
    CASOS$feminicidios[i] <- 1
  }
  
  if(CASOS$"Habla Indigena"[i] == "Si"){
    CASOS$"Habla lengua indigena"[i] <- 1
  }
  
  
  if(!is.na(CASOS$`Detalle de Vínculo con victima`[i])){
    for (j in familiares) {
      if(CASOS$`Detalle de Vínculo con victima`[i] == j){
        CASOS$`Detalle de Vínculo con victima`[i] <- "Otros familiares"
        break
      }
    }
    if(CASOS$`Detalle de Vínculo con victima`[i] == "Otro"){
      CASOS$`Detalle de Vínculo con victima`[i] <- "Otra relación"
    }
    if(CASOS$`Detalle de Vínculo con victima`[i] == "NULL"){
      CASOS$`Detalle de Vínculo con victima`[i] <- "Desconocido"
    }
    if(CASOS$`Detalle de Vínculo con victima`[i] == "Concubina"){
      CASOS$`Detalle de Vínculo con victima`[i] <- "Cónyuge o pareja"
    }
  }
  
  
  
}

Se agrega una nueva columna correspondiente al estado de hecho, la cual se obtiene a partir de los datos de municipio hecho.

Municipio = c("Abasolo", "Agualeguas", "Los Aldamas", "Allende", "Anáhuac",
              "Apodaca", "Aramberri", "Bustamante", "Cadereyta Jiménez",
              "El Carmen", "Cerralvo", "Ciénega de Flores", "China",
              "Doctor Arroyo", "Doctor Coss", "Doctor González",
              "Los Herrera", "Galeana", "García", "General Bravo",
              "General Escobedo", "General Terán", "General Treviño",
              "General Zaragoza", "General Zuazua", "Guadalupe",
              "Hidalgo", "Higueras", "Hualahuises", "Iturbide",
              "Juárez", "Lampazos de Naranjo", "Linares", "Marín",
              "Melchor Ocampo", "Mier y Noriega", "Mina", "Montemorelos",
              "Monterrey", "Paras", "Pesquería", "Los Ramones",
              "Rayones", "Sabinas Hidalgo", "Salinas Victoria",
              "San Nicolás de los Garza", "San Pedro Garza García",
              "Santa Catarina", "Santiago", "Vallecillo", "Villaldama")

CASOS$estado <- NA

for (i in 1:nrow(CASOS)) {
  if(!is.na(CASOS$`Municipio Hecho`[i])){
    for (j in Municipio) {
      if(CASOS$`Municipio Hecho`[i] == j ){
        CASOS$estado[i] <- "Nuevo León"
      }
    }
  }
}