library(rio)
acoso=import("Acoso_R.Rdata")
library(dslabs)
str(acoso)## tibble [1,576 x 32] (S3: tbl_df/tbl/data.frame)
## $ SERIE : chr [1:1576] "2012" "2012" "2012" "2012" ...
## ..- attr(*, "label")= chr "Serie"
## ..- attr(*, "format.spss")= chr "A4"
## ..- attr(*, "display_width")= int 10
## $ NVONRO : num [1:1576] 20120001 20120002 20120003 20120004 20120005 ...
## ..- attr(*, "label")= chr "Número de encuesta"
## ..- attr(*, "format.spss")= chr "F8.0"
## ..- attr(*, "display_width")= int 10
## $ NSEGRUP: Factor w/ 3 levels "AB","C","DE": 1 1 1 1 1 2 3 2 2 2 ...
## $ NSE : dbl+lbl [1:1576] 2, 2, 1, 2, 2, 3, 4, 3, 3, 3, 3, 4, 3, 5, 3, 2, 3, 2,...
## ..@ label : chr "Niveles socioeconómicos (version larga)"
## ..@ format.spss : chr "F8.2"
## ..@ display_width: int 10
## ..@ labels : Named num [1:5] 1 2 3 4 5
## .. ..- attr(*, "names")= chr [1:5] "A" "B" "C" "D" ...
## $ SEXO : Factor w/ 2 levels "Hombre","Mujer": 2 1 2 2 2 1 1 1 1 2 ...
## $ GEDAD : dbl+lbl [1:1576] 3, 3, 3, 1, 1, 1, 3, 2, 3, 3, 2, 3, 1, 1, 2, 2, 2, 3,...
## ..@ label : chr "Grupo de edad del entrevistado"
## ..@ format.spss : chr "F2.0"
## ..@ display_width: int 5
## ..@ labels : Named num [1:3] 1 2 3
## .. ..- attr(*, "names")= chr [1:3] "18 a 29" "30 a 44" "45 o más"
## $ EDAD : num [1:1576] 55 62 70 20 18 25 46 38 52 53 ...
## $ NIVEDU : Factor w/ 10 levels "Ninguna","Inicial",..: 5 7 7 8 7 6 5 5 5 5 ...
## $ P10A : dbl+lbl [1:1576] 2, 2, 3, 3, 3, 3, 3, 3, 4, 3, 3, 4, 3, 3, 3, 3, 2, 3,...
## ..@ label : chr "10a. ¿Qué tan seguro(a) o inseguro(a) se siente usted caminando en la calle: A cualquier hora del día?"
## ..@ format.spss : chr "F2.0"
## ..@ display_width: int 6
## ..@ labels : Named num [1:8] -1 1 2 3 4 8 9 98
## .. ..- attr(*, "names")= chr [1:8] "Sin respuesta" "Muy seguro(a)" "Seguro(a)" "Inseguro(a)" ...
## $ P10B : dbl+lbl [1:1576] 3, 2, 3, 4, 4, 4, 3, 3, 4, 4, 3, 4, 3, 3, 3, 3, 2, 3,...
## ..@ label : chr "10b. ¿Qué tan seguro(a) o inseguro(a) se siente usted caminando en la calle: En la noche?"
## ..@ format.spss : chr "F2.0"
## ..@ display_width: int 6
## ..@ labels : Named num [1:8] -1 1 2 3 4 8 9 98
## .. ..- attr(*, "names")= chr [1:8] "Sin respuesta" "Muy seguro(a)" "Seguro(a)" "Inseguro(a)" ...
## $ P10C : dbl+lbl [1:1576] 3, 3, 3, 3, 4, 4, 3, 3, 4, 4, 3, 4, 3, 3, 3, 3, 3, 3,...
## ..@ label : chr "10c. ¿Qué tan seguro(a) o inseguro(a) se siente usted caminando en la calle: Por zonas oscuras o descampados?"
## ..@ format.spss : chr "F2.0"
## ..@ display_width: int 6
## ..@ labels : Named num [1:8] -1 1 2 3 4 8 9 98
## .. ..- attr(*, "names")= chr [1:8] "Sin respuesta" "Muy seguro(a)" "Seguro(a)" "Inseguro(a)" ...
## $ P11A : dbl+lbl [1:1576] 3, 3, 3, 3, 3, 2, 3, 3, 3, 3, 3, 3, 3, 2, 3, 3, 3, 2,...
## ..@ label : chr "11a. Diría usted que es más riesgoso para un hombre, para una mujer o es igual de riesgoso para ambos: Recorrer"| __truncated__
## ..@ format.spss : chr "F2.0"
## ..@ display_width: int 6
## ..@ labels : Named num [1:7] -1 1 2 3 8 9 98
## .. ..- attr(*, "names")= chr [1:7] "Sin respuesta" "Hombre" "Mujer" "Ambos" ...
## $ P11B : dbl+lbl [1:1576] 3, 3, 3, 3, 3, 2, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 2,...
## ..@ label : chr "11b. Diría usted que es más riesgoso para un hombre, para una mujer o es igual de riesgoso para ambos: Recorrer"| __truncated__
## ..@ format.spss : chr "F2.0"
## ..@ display_width: int 6
## ..@ labels : Named num [1:7] -1 1 2 3 8 9 98
## .. ..- attr(*, "names")= chr [1:7] "Sin respuesta" "Hombre" "Mujer" "Ambos" ...
## $ P11C : dbl+lbl [1:1576] 3, 3, 3, 3, 3, 2, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 2,...
## ..@ label : chr "11c. Diría usted que es más riesgoso para un hombre, para una mujer o es igual de riesgoso para ambos: Caminar "| __truncated__
## ..@ format.spss : chr "F2.0"
## ..@ display_width: int 6
## ..@ labels : Named num [1:7] -1 1 2 3 8 9 98
## .. ..- attr(*, "names")= chr [1:7] "Sin respuesta" "Hombre" "Mujer" "Ambos" ...
## $ P12A : num [1:1576] 1 0 1 1 1 1 0 1 0 1 ...
## $ P12B : num [1:1576] 0 0 0 1 1 0 0 0 0 0 ...
## $ P12C : num [1:1576] 1 0 0 1 1 1 0 0 0 0 ...
## $ P12D : num [1:1576] 1 0 0 1 1 0 0 1 0 0 ...
## $ P12E : num [1:1576] 0 0 0 1 0 1 0 1 0 0 ...
## $ P12F : num [1:1576] 0 0 0 0 0 1 0 1 0 0 ...
## $ P12G : num [1:1576] 0 0 0 0 0 0 0 0 0 0 ...
## $ P12H : num [1:1576] 0 0 0 1 1 1 0 1 0 0 ...
## $ P13A : dbl+lbl [1:1576] 1, 2, 1, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 2, 1, 2, 2, 1,...
## ..@ label : chr "13a. ¿Usted considera que las personas deberían ser castigadas si realizan los siguientes actos: Miradas persis"| __truncated__
## ..@ format.spss : chr "F2.0"
## ..@ display_width: int 6
## ..@ labels : Named num [1:6] -1 1 2 8 9 98
## .. ..- attr(*, "names")= chr [1:6] "Sin respuesta" "Sí" "No" "No sabe" ...
## $ P13B : dbl+lbl [1:1576] 1, 2, 2, 2, 2, 2, 2, 2, 1, 1, 2, 1, 2, 2, 2, 1, 2, 2,...
## ..@ label : chr "13b. ¿Usted considera que las personas deberían ser castigadas si realizan los siguientes actos: Silbidos?"
## ..@ format.spss : chr "F2.0"
## ..@ display_width: int 6
## ..@ labels : Named num [1:6] -1 1 2 8 9 98
## .. ..- attr(*, "names")= chr [1:6] "Sin respuesta" "Sí" "No" "No sabe" ...
## $ P13C : dbl+lbl [1:1576] 1, 2, 2, 2, 2, 2, 2, 2, 1, 1, 2, 1, 2, 1, 2, 2, 2, 2,...
## ..@ label : chr "13c. ¿Usted considera que las personas deberían ser castigadas si realizan los siguientes actos: Ruidos de besos?"
## ..@ format.spss : chr "F2.0"
## ..@ display_width: int 6
## ..@ labels : Named num [1:6] -1 1 2 8 9 98
## .. ..- attr(*, "names")= chr [1:6] "Sin respuesta" "Sí" "No" "No sabe" ...
## $ P13D : dbl+lbl [1:1576] 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2,...
## ..@ label : chr "13d. ¿Usted considera que las personas deberían ser castigadas si realizan los siguientes actos: Gestos vulgares?"
## ..@ format.spss : chr "F2.0"
## ..@ display_width: int 6
## ..@ labels : Named num [1:6] -1 1 2 8 9 98
## .. ..- attr(*, "names")= chr [1:6] "Sin respuesta" "Sí" "No" "No sabe" ...
## $ P13E : dbl+lbl [1:1576] 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1,...
## ..@ label : chr "13e. ¿Usted considera que las personas deberían ser castigadas si realizan los siguientes actos: Comentarios e "| __truncated__
## ..@ format.spss : chr "F2.0"
## ..@ display_width: int 6
## ..@ labels : Named num [1:6] -1 1 2 8 9 98
## .. ..- attr(*, "names")= chr [1:6] "Sin respuesta" "Sí" "No" "No sabe" ...
## $ P13F : dbl+lbl [1:1576] 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1,...
## ..@ label : chr "13f. ¿Usted considera que las personas deberían ser castigadas si realizan los siguientes actos: Tocamientos o “manoseos”?"
## ..@ format.spss : chr "F2.0"
## ..@ display_width: int 6
## ..@ labels : Named num [1:6] -1 1 2 8 9 98
## .. ..- attr(*, "names")= chr [1:6] "Sin respuesta" "Sí" "No" "No sabe" ...
## $ P13G : dbl+lbl [1:1576] 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1,...
## ..@ label : chr "13g. ¿Usted considera que las personas deberían ser castigadas si realizan los siguientes actos: Exhibicionismo"| __truncated__
## ..@ format.spss : chr "F2.0"
## ..@ display_width: int 6
## ..@ labels : Named num [1:6] -1 1 2 8 9 98
## .. ..- attr(*, "names")= chr [1:6] "Sin respuesta" "Sí" "No" "No sabe" ...
## $ P13H : dbl+lbl [1:1576] 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1,...
## ..@ label : chr "13h. ¿Usted considera que las personas deberían ser castigadas si realizan los siguientes actos: Roces incómodo"| __truncated__
## ..@ format.spss : chr "F2.0"
## ..@ display_width: int 6
## ..@ labels : Named num [1:6] -1 1 2 8 9 98
## .. ..- attr(*, "names")= chr [1:6] "Sin respuesta" "Sí" "No" "No sabe" ...
## $ index_O: num [1:1576] 3 0 1 6 5 5 0 5 0 1 ...
## $ index_B: num [1:1576] 1 0 1 1 1 1 0 1 0 1 ...library(MASS)
library(dslabs)
acoso$index_O = cut(acoso$index_O, breaks=c(0,2,4,6,8), labels=c("muy poca","poca","media", "alta"))acoso$inde1= as.factor(acoso$index_O)
modelo1 <- polr(inde1~ SEXO + EDAD, data = acoso, Hess=TRUE)
summary(modelo1)## Call:
## polr(formula = inde1 ~ SEXO + EDAD, data = acoso, Hess = TRUE)
##
## Coefficients:
## Value Std. Error t value
## SEXOMujer 0.82281 0.165641 4.967
## EDAD -0.03412 0.005818 -5.864
##
## Intercepts:
## Value Std. Error t value
## muy poca|poca -0.8795 0.2448 -3.5925
## poca|media 0.6578 0.2438 2.6976
## media|alta 2.3680 0.2781 8.5165
##
## Residual Deviance: 1626.027
## AIC: 1636.027
## (874 observations deleted due to missingness)INTERPRETACIÓN: Si eres mujer se espera que el logaritmo de odds aumente en 0.8 de ser victima porlo menos una vez de acoso Por cada año de edad que aumente se espera que el logaritmo de odds disminuya en -0.03412 de ser victima por lo menos una vez de acoso.
modelo2=glm(index_B~ SEXO + EDAD, data = acoso)
summary(modelo2)##
## Call:
## glm(formula = index_B ~ SEXO + EDAD, data = acoso)
##
## Deviance Residuals:
## Min 1Q Median 3Q Max
## -0.84344 -0.38478 -0.01412 0.33689 1.17623
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.6753097 0.0313321 21.55 <2e-16 ***
## SEXOMujer 0.3484548 0.0218459 15.95 <2e-16 ***
## EDAD -0.0100181 0.0006568 -15.25 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for gaussian family taken to be 0.1877342)
##
## Null deviance: 389.31 on 1575 degrees of freedom
## Residual deviance: 295.31 on 1573 degrees of freedom
## AIC: 1841.3
##
## Number of Fisher Scoring iterations: 2