#Cambio en la pobreza
jpeg('grafica1.jpg',width=16,height=7,units='cm',res=200) par(mfrow=c(1,2),mar=c(4,4,3,1)) graf_pob(binom_exp(pobreza~year,pob),tlim=c(0,800000),plim=c(0,45), col=colores[1:2],'','') grid.text('Elaboración propia con base en los resultados de la ENIGH (2010-2012)', x=1,y=0.05,draw=1,just='right',gp=gpar(col='grey',cex=0.7)) grid.text('Figura 1: Poblacion en situación de Pobreza Multidimensional', x=0.12,y=0.92,draw=1,just='left',gp=gpar(cex=0.8,font=2)) dev.off()
binom_exp(pobreza~year,pob)
#El ingreso y la plb #Plb jpeg('grafica2.jpg',width=16,height=7,units='cm',res=200) par(mfrow=c(1,3),mar=c(4,4,4,1)) graf_pob(binom_exp(plb~year,pob),tlim=c(0,1000000),plim=c(0,50), col=colores[3:4],'','') barplot(svyby(~ictpcD,~year,FUN=svymean,pob)[,2],horiz=0,col=colores[3:4], names.arg=c(2010,2012),ylim=c(0,4000),cex.main=0.8,ylab='Pesos', cex.lab=0.8,cex.axis=0.8,cex.names=0.7,axis.lty=1,space=0) segments(y0=confint(svyby(~ictpcD,~year,FUN=svymean,pob))[,1],y1=confint(svyby(~ictpcD,~year,FUN=svymean,pob))[,2], x0=c(0.5,1.5)) segments(y0=confint(svyby(~ictpcD,~year,FUN=svymean,pob))[,1],x0=c(0.5,1.5)-0.2,x1=c(0.5,1.5)+0.2) segments(y0=confint(svyby(~ictpcD,~year,FUN=svymean,pob))[,2],x0=c(0.5,1.5)-0.2,x1=c(0.5,1.5)+0.2) grid.text('Figura 2: Población vulnerable por ingreso', x=0.12,y=0.92,draw=1,just='left',gp=gpar(cex=0.8,font=2)) grid.text('Elaboración propia con base en los resultados de la ENIGH (2010-2012)', x=1,y=0.05,draw=1,just='right',gp=gpar(col='grey',cex=0.7)) dev.off()
binom_exp(plb~year,pob)
#Carencias
jpeg('grafica3.jpg',width=16,height=7,units='cm',res=200) par(mfrow=c(1,2),mar=c(4,4,2.5,1)) graf_pob(binom_exp(carencia~year,pob),tlim=c(0,1400000),plim=c(0,75), col=colores[5:6],'','') grid.text('Elaboración propia con base en los resultados de la ENIGH (2010-2012)', x=1,y=0.05,draw=1,just='right',gp=gpar(col='grey',cex=0.7)) grid.text('Figura 3: Poblacion vulnerable por carencia', x=0.12,y=0.92,draw=1,just='left',gp=gpar(cex=0.8,font=2)) dev.off()
binom_exp(carencia~year,pob)
#Número de carencias
jpeg('grafica4.jpg',width=16,height=7,units='cm',res=200) par(mfrow=c(1,1),mar=c(3.5,4,2,1)) plot(0:6,svytable(~num_car+year,pob)[,1]/1000,type='l',col=colores[5], ylim=c(0,600),lwd=1.7,main='',ylab='miles de personas', xlab='',cex.main=0.8,cex.lab=0.8,cex.axis=0.8) lines(0:6,svytable(~num_car+year,pob)[,2]/1000,type='l',col=colores[6],lwd=1.7) lines(1:6,svytable(~num_car+year,subset(pob,pobreza==1))[,1]/1000,type='l',col=colores[5],lty=2,lwd=1.7) lines(1:6,svytable(~num_car+year,subset(pob,pobreza==1))[,2]/1000,type='l',col=colores[6],lty=2,lwd=1.7) legend('topright',bty='n',lty=c(1,1,2),col=c(colores[5:6],1),legend=c('2010','2012','Población en PM'), cex=0.7,lwd=1.7) grid.text('Elaboración propia con base en los resultados de la ENIGH (2010-2012)', x=1,y=0.05,draw=1,just='right',gp=gpar(col='grey',cex=0.7)) grid.text('Figura 4: Población según número de carencias', x=0.12,y=0.92,draw=1,just='left',gp=gpar(cex=0.8,font=2)) dev.off()
svyby(~num_car,~year,FUN=svymean,pob) svyby(~num_car,~year+pobreza,FUN=svymean,pob)
summary(svyglm(num_car~year,family=quasipoisson,pob)) summary(svyglm(num_car~year,family=quasipoisson,subset(pob,pobreza==1)))
svytable(~num_car+year,pob) svytable(~num_car+year,subset(pob,pobreza==1))
#Calidad y espacios en la vivienda
jpeg('grafica5.jpg',width=16,height=7,units='cm',res=200) par(mfrow=c(1,2),mar=c(4,4,3,1)) graf_pob(binom_exp(ic_cv~year,subset(pob,pobreza==1)),plim=c(0,25),tlim=c(0,200000), col=colores[7:8],'','') grid.text('Elaboración propia con base en los resultados de la ENIGH (2010-2012)', x=1,y=0.05,draw=1,just='right',gp=gpar(col='grey',cex=0.7)) grid.text('Figura 5: Poblacion en PM vulnerable por calidad y espacios en la vivienda', x=0.12,y=0.92,draw=1,just='left',gp=gpar(cex=0.8,font=2)) dev.off()
binom_exp(ic_cv~year,subset(pob,pobreza==1)) binom_exp(ic_cv~year,pob)
binom_exp(icv_mat_pi~year,pob) #proporcional y total binom_exp(icv_mat_te~year,pob) #no hay diferencia binom_exp(icv_mat_mu~year,pob) #no hay diferencia binom_exp(icv_hac~year,pob) #no hay diferencia
#Alimentación
jpeg('grafica6.jpg',width=16,height=7,units='cm',res=200) par(mfrow=c(1,2),mar=c(4,4,3,1)) graf_pob(binom_exp(ic_ali~year,subset(pob,pobreza==1)),plim=c(0,40),tlim=c(0,300000), col=colores[9:10],'','') grid.text('Elaboración propia con base en los resultados de la ENIGH (2010-2012)', x=1,y=0.05,draw=1,just='right',gp=gpar(col='grey',cex=0.7)) grid.text('Figura 6: Poblacion en PM vulnerable por acceso a la alimentación', x=0.12,y=0.92,draw=1,just='left',gp=gpar(cex=0.8,font=2)) dev.off()
binom_exp(ic_ali~year,subset(pob,pobreza==1))
#Educación
jpeg('grafica7.jpg',width=16,height=7,units='cm',res=200) par(mfrow=c(1,2),mar=c(4,4,3,1)) graf_pob(binom_exp(ic_rezedu~year,subset(pob,pobreza==1)),plim=c(0,40),tlim=c(0,300000), col=colores[11:12],'','') grid.text('Elaboración propia con base en los resultados de la ENIGH (2010-2012)', x=1,y=0.05,draw=1,just='right',gp=gpar(col='grey',cex=0.7)) grid.text('Figura 7: Poblacion en PM vulnerable por rezago educativo', x=0.12,y=0.92,draw=1,just='left',gp=gpar(cex=0.8,font=2)) dev.off()
binom_exp(ic_rezedu~year,subset(pob,pobreza==1)) binom_exp(ic_rezedu~year,pob)
#Servicios básicos de vivienda
jpeg('grafica8.jpg',width=16,height=7,units='cm',res=200) par(mfrow=c(1,2),mar=c(4,4,3,1)) graf_pob(binom_exp(ic_sbv~year,subset(pob,pobreza==1)),plim=c(0,40),tlim=c(0,300000), col=colores[7:8],'','') grid.text('Elaboración propia con base en los resultados de la ENIGH (2010-2012)', x=1,y=0.05,draw=1,just='right',gp=gpar(col='grey',cex=0.7)) grid.text('Figura 8: Poblacion en PM vulnerable por acceso a servicios básicos en la vivienda', x=0.12,y=0.92,draw=1,just='left',gp=gpar(cex=0.8,font=2)) dev.off()
binom_exp(ic_sbv~year,subset(pob,pobreza==1))
jpeg('grafica9.jpg',width=16,height=7,units='cm',res=200) par(mfrow=c(1,2),mar=c(4,4,3,1)) graf_pob(binom_exp(ic_segsoc~year,subset(pob,pobreza==1)),plim=c(0,90),tlim=c(0,700000), col=colores[9:10],'','') grid.text('Elaboración propia con base en los resultados de la ENIGH (2010-2012)', x=1,y=0.05,draw=1,just='right',gp=gpar(col='grey',cex=0.7)) grid.text('Figura 9: Poblacion en PM vulnerable por acceso a seguridad social', x=0.12,y=0.92,draw=1,just='left',gp=gpar(cex=0.8,font=2)) dev.off()
binom_exp(ic_segsoc~year,subset(pob,pobreza==1)) binom_exp(ic_segsoc~year,pob) binom_exp(pam~year,subset(pob,pobreza==1))
svytable(~year+pam+pobreza,pob)
jpeg('grafica10.jpg',width=16,height=7,units='cm',res=200) par(mfrow=c(1,2),mar=c(4,4,3,1)) graf_pob(binom_exp(ic_asalud~year,subset(pob,pobreza==1)),plim=c(0,35),tlim=c(0,300000), col=colores[11:12],'','') grid.text('Elaboración propia con base en los resultados de la ENIGH (2010-2012)', x=1,y=0.05,draw=1,just='right',gp=gpar(col='grey',cex=0.7)) grid.text('Figura 10: Poblacion en PM vulnerable por acceso a servicios de salud', x=0.12,y=0.92,draw=1,just='left',gp=gpar(cex=0.8,font=2)) dev.off()
binom_exp(ic_asalud~year,subset(pob,pobreza==1)) binom_exp(SP~year,pob)