PRÁCTICA CALIFICADA 2
Código:20196323
Sheccid Yanet Gomez Mamani:
PARTE 2
library(rio)
data = import("ENARES2019_FINAL.csv")PREGUNTA 1
Construya un índice aditivo sobre la magnitud de violencia contra las mujeres casadas o en convivencia (con un mínimo de 13 y un máximo de 29), describa e interprete. Para ello, deberá utilizar las 19 variables de violencia psicológica (Pregunta 402) , las 7 variables de autonomía económica (Pregunta 403), y las 15 variables de violencia física (Pregunta 404). Estas variables han sido recodificadas, siendo 0 = “No”, y 1 = “Sí” en respuesta a cada pregunta. (3 puntos).
data$violenciapsicologica=data$C1P402_1+data$C1P402_2+data$C1P402_3+data$C1P402_4+data$C1P402_5+data$C1P402_6+data$C1P402_7+data$C1P402_8+data$C1P402_9+data$C1P402_10+data$C1P402_11+data$C1P402_12+data$C1P402_13+data$C1P402_14+data$C1P402_15+data$C1P402_16+data$C1P402_17+data$C1P402A_18+data$C1P402_19str(data$violenciapsicologica)## int [1:317] NA NA 7 11 NA NA NA NA NA NA ...data$autonomia=data$C1P403_1+data$C1P403_2+data$C1P403_3+data$C1P403_4+data$C1P403_5+data$C1P403_6+data$C1P403_7str(data$autonomia)## int [1:317] 5 7 5 2 7 7 7 7 6 6 ...data$violenciafísica=data$C1P404_1+data$C1P404_2+data$C1P404_3+data$C1P404_4+data$C1P404_5+data$C1P404_6+data$C1P404_7+data$C1P404_8+data$C1P404_9+data$C1P404_10+data$C1P404_11+data$C1P404_12+data$C1P404_13+data$C1P404_14+data$C1P404_15data$violenciafísica## [1] 8 1 9 8 1 7 3 9 2 4 1 1 1 1 2 1 3 2 9 3 3 1 1 5 1
## [26] 1 2 1 5 8 2 1 1 1 6 6 6 1 9 1 7 1 6 1 2 5 6 3 9 5
## [51] 2 1 1 3 6 2 6 1 6 2 1 1 1 1 1 2 1 1 2 1 1 1 2 5 6
## [76] 4 3 1 2 2 1 7 1 3 2 3 1 6 1 1 1 4 8 4 6 9 2 7 2 4
## [101] 5 5 1 5 3 1 4 1 4 1 1 1 1 1 2 3 3 1 2 2 7 2 2 1 3
## [126] 1 3 1 1 1 1 11 1 2 4 10 1 3 5 1 1 2 11 1 1 1 2 1 3 4
## [151] 1 1 3 6 1 2 2 2 9 4 2 1 1 1 2 3 2 4 2 1 2 1 5 5 1
## [176] 3 1 1 1 4 1 3 1 3 1 2 2 1 3 1 1 1 1 1 1 1 7 2 1 3
## [201] 2 6 5 4 2 1 1 7 1 3 4 1 5 2 2 2 2 5 5 3 9 4 2 3 7
## [226] 5 2 2 2 4 9 3 2 3 1 1 2 2 7 1 1 2 1 7 2 8 1 2 6 1
## [251] 4 5 10 9 2 11 7 2 6 2 2 2 2 2 7 1 3 2 1 2 5 5 7 5 9
## [276] 5 5 4 7 1 1 3 3 7 2 3 3 4 2 2 2 1 10 7 1 1 2 4 1 2
## [301] 3 5 1 4 2 1 2 4 5 2 1 2 1 2 1 1 2data$violenciapsicologica=as.numeric(data$violenciapsicologica)
data$autonomia=as.numeric(data$autonomia)
data$violenciafísica=as.numeric(data$violenciafísica)str(data$violenciapsicologica)## num [1:317] NA NA 7 11 NA NA NA NA NA NA ...str(data$autonomia)## num [1:317] 5 7 5 2 7 7 7 7 6 6 ...str(data$violenciafísica)## num [1:317] 8 1 9 8 1 7 3 9 2 4 ...data$violenciaindice=data$violenciapsicologica+data$autonomia+data$violenciafísica-2
summary(data$violenciaindice)## Min. 1st Qu. Median Mean 3rd Qu. Max. NA's
## 13.00 15.00 18.00 17.92 20.00 24.00 281data = data[complete.cases(data$violenciaindice),]PREGUNTA 2
Determine si existen diferencias significativas del promedio del índice de magnitud de la violencia contra la mujer (índice creado en la pregunta 1), entre la frecuencia en la que el esposo toma o tomaba bebidas alcohólicas (C1P408) (2 puntos)
str(data$C1P408)## int [1:36] 2 2 1 2 1 1 1 1 2 1 ...data$C1P408 = factor(data$C1P408,
levels = c(1:3),
labels = c("siempre/casi siempre","algunas veces","rara vez"))str(data$C1P408)## Factor w/ 3 levels "siempre/casi siempre",..: 2 2 1 2 1 1 1 1 2 1 ...table(data$C1P408)##
## siempre/casi siempre algunas veces rara vez
## 15 12 6sum(is.na(data$C1P408))## [1] 3data = data[complete.cases(data$C1P408),]1er PASO:
Ho: No hay diferencia entre las media del índice de magnitud de la violencia contra la mujer y la frecuencia en la que el esposo toma o tomaba bebidas alcohólicas
H1: Si existen diferencia entre las media del índice de magnitud de la violencia contra la mujer y la frecuencia en la que el esposo toma o tomaba bebidas alcohólicas
library(gplots)##
## Attaching package: 'gplots'## The following object is masked from 'package:stats':
##
## lowessplotmeans(data$violenciaindice ~ data$C1P408, connect=F, barwidth=6, xlab="frecuencia en la que esposa toma", ylab="promedio de violencia recibida por la mujer",
main="Promedio de violencia contra la mujer según frecuencia en la que el esposo toma")
anova <- aov(data$violenciaindice~data$C1P408)
summary(anova)## Df Sum Sq Mean Sq F value Pr(>F)
## data$C1P408 2 35.5 17.75 1.627 0.213
## Residuals 30 327.4 10.91INTERPRETACIÓN: Ahora, confirmemos la observación visual con la prueba ANOVA. Obsevamos que el p-value es mayor a 0.05, por lo que podemos aceptar la hipotesis nula y afirmamos que no hay una diferencia significativa entre la frecuencia en que el esposo toma y el promeido que recibe la mujer. Esto coincide con lo reportado por el gráfico anterior.
PREGUNTA 3
¿Existen diferencias significativas entre la media del índice de magnitud de violencia física creado, entre mujeres que tienen al menos un trabajo de una hora para obtener ingresos (C1P307)? Utilice la prueba pertinente, interprete y realice un gráfico de barras de error. Sugerencia: recuerde configurar adecuadamente la variable (3 puntos)
str(data$C1P307)## int [1:33] 1 1 1 1 1 2 1 1 1 1 ...data$C1P307 = factor(data$C1P307,
levels = c(1:2),
labels = c("si","no"))str(data$C1P307)## Factor w/ 2 levels "si","no": 1 1 1 1 1 2 1 1 1 1 ...table(data$C1P307)##
## si no
## 25 8t.test(violenciaindice ~ C1P307, data = data,conf.level = 0.95)##
## Welch Two Sample t-test
##
## data: violenciaindice by C1P307
## t = -0.37506, df = 10.367, p-value = 0.7152
## alternative hypothesis: true difference in means between group si and group no is not equal to 0
## 95 percent confidence interval:
## -3.940059 2.800059
## sample estimates:
## mean in group si mean in group no
## 17.68 18.25PREGUNTA 4
Determine si existen diferencias significativas del promedio del índice de magnitud de la violencia contra la mujer (índice creado en la pregunta 1), entre la frecuencia en la que el esposo la golpeó o agredió bajo los efectos del alcohol en los últimos 12 meses (C1P410) (2 puntos)