R Notebook
library(haven)
datos <- read_stata("~/Downloads/MexicoLAPOP")
library(ggplot2)
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(magrittr)
#Descriptivos
grafico1 <- datos %>% filter(e5!= "NA") %>% ggplot(aes(x=factor(e5))) +
geom_bar() +
xlab("") + ylab("Numero de casos") +
labs(title="Grafico 1: Participacion en manifestaciones legales",
subtitle="Hasta que punto aprueba o desaprueba", caption="Fuente: LAPOP
(2017)")
grafico1
library(psych)
##
## Attaching package: 'psych'
## The following objects are masked from 'package:ggplot2':
##
## %+%, alpha
describe(datos$e5)
## vars n mean sd median trimmed mad min max range skew kurtosis
## X1 1 1522 6.33 3.03 7 6.54 4.45 1 10 9 -0.36 -1.06
## se
## X1 0.08
[Link](datos$e5, paired=FALSE, mu = 6, [Link] = 0.95)
##
## One Sample t-test
##
## data: datos$e5
## t = 4.2491, df = 1521, p-value = 2.276e-05
## alternative hypothesis: true mean is not equal to 6
## 95 percent confidence interval:
## 6.177570 6.482089
## sample estimates:
## mean of x
## 6.329829
datos$estratopri <- factor(datos$estratopri)
levels(datos$estratopri) <- c("Norte", "Centro Occidente", "Centro", "Sur")
datos$estratosec <- factor(datos$estratosec)
levels(datos$estratosec) <- c("Grande (mas de 100,000)", "Mediana (Entre
25,000 - 100,000)", "Pequena (< 25,000)")
datos$q1 <- factor(datos$q1)
levels(datos$q1) <- c("Hobre", "Mujer")
datos$etid <- factor(datos$etid)
levels(datos$etid) <- c("Blanco", "Mestizo", "Indigena", "Negro", "Mulato",
"Otro")
datos$q10newf <- cut(datos$q10new, c(0,4,12,16))
levels(datos$q10newf) <- c("Menos de 2900 pesos mensuales", "Entre 2900 y
6200 pesos mensuales", "Mas de 6200 mensuales")
datos$pol1 <- factor(datos$pol1)
levels(datos$pol1) <- c("Mucho", "Algo", "Poco", "Nada")
regresion1 <- lm(e5~ estratopri + estratosec + q1 + l1 + factor(etid) + q2,
data=datos)
regresion2 <- lm(e5~ estratopri + estratosec + q1 + l1 + factor(etid) + q2 +
q10new, data=datos)
regresion3 <- lm(e5~ estratopri + estratosec + q1 + l1 + factor(etid) + q2 +
q10newf + ed + factor(pol1), data=datos)
library(sjPlot)
summary(regresion3)
##
## Call:
## lm(formula = e5 ~ estratopri + estratosec + q1 + l1 + factor(etid) +
## q2 + q10newf + ed + factor(pol1), data = datos)
##
## Residuals:
## <Labelled double>
## Min 1Q Median 3Q Max
## -7.6308 -1.9618 0.1846 2.3058 6.4666
##
## Labels:
## value label
## 1 Strongly disapprove
## 10 Strongly approve
## NA(a) Don't Know
## NA(b) No Response
##
## Coefficients:
## Estimate Std. Error t value
## (Intercept) 5.829448 0.646117 9.022
## estratopriCentro Occidente -0.074058 0.267509 -0.277
## estratopriCentro -0.918779 0.236733 -3.881
## estratopriSur -0.870634 0.274880 -3.167
## estratosecMediana (Entre 25,000 - 100,000) 0.194577 0.241223 0.807
## estratosecPequena (< 25,000) -0.486723 0.301241 -1.616
## q1Mujer -0.290492 0.184059 -1.578
## l1 -0.007288 0.035028 -0.208
## factor(etid)Mestizo 0.413828 0.231738 1.786
## factor(etid)Indigena 0.867380 0.335235 2.587
## factor(etid)Negro 1.049278 0.574909 1.825
## factor(etid)Mulato 1.207575 0.610348 1.979
## factor(etid)Otro 0.285118 0.392986 0.726
## q2 0.008809 0.006332 1.391
## q10newfEntre 2900 y 6200 pesos mensuales 0.264546 0.215288 1.229
## q10newfMas de 6200 mensuales 0.482040 0.265881 1.813
## ed 0.107997 0.025975 4.158
## factor(pol1)Algo -0.418755 0.323650 -1.294
## factor(pol1)Poco -0.589392 0.318338 -1.851
## factor(pol1)Nada -1.543276 0.324844 -4.751
## Pr(>|t|)
## (Intercept) < 2e-16 ***
## estratopriCentro Occidente 0.781955
## estratopriCentro 0.000111 ***
## estratopriSur 0.001584 **
## estratosecMediana (Entre 25,000 - 100,000) 0.420071
## estratosecPequena (< 25,000) 0.106462
## q1Mujer 0.114816
## l1 0.835228
## factor(etid)Mestizo 0.074436 .
## factor(etid)Indigena 0.009808 **
## factor(etid)Negro 0.068275 .
## factor(etid)Mulato 0.048141 *
## factor(etid)Otro 0.468301
## q2 0.164440
## q10newfEntre 2900 y 6200 pesos mensuales 0.219430
## q10newfMas de 6200 mensuales 0.070128 .
## ed 3.48e-05 ***
## factor(pol1)Algo 0.196010
## factor(pol1)Poco 0.064392 .
## factor(pol1)Nada 2.32e-06 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 2.854 on 1020 degrees of freedom
## (523 observations deleted due to missingness)
## Multiple R-squared: 0.1073, Adjusted R-squared: 0.0907
## F-statistic: 6.455 on 19 and 1020 DF, p-value: 4.455e-16
