data = read.csv("Enem2015.csv",header=TRUE) attach(data) pdf(file="Enem2015-graficos.pdf",width=15,height=10) provas = c("Linguagens","Redacao","Matematica","Humanas","Natureza") provas1 = c("L","R","M","H","N") tipoescola = c("Municipal","Estadual","Federal","Privada") tipoescola1 = c("M","E","F","P") par(mfrow=c(1,5)) for (i in 15:19){ boxplot(data[dependenciaadministrativa=="Municipal",i],data[dependenciaadministrativa=="Estadual",i], data[dependenciaadministrativa=="Federal",i],data[dependenciaadministrativa=="Privada",i], names=tipoescola1,ylab="Media",main=provas[i-14],ylim=c(340,950)) } par(mfrow=c(1,4)) for (i in 1:4){ boxplot(data[dependenciaadministrativa==tipoescola[i],15],data[dependenciaadministrativa==tipoescola[i],16], data[dependenciaadministrativa==tipoescola[i],17],data[dependenciaadministrativa==tipoescola[i],18], data[dependenciaadministrativa==tipoescola[i],19],names=provas1,ylab="Media",main=tipoescola[i],ylim=c(340,950)) } #dev.off() notamedia = apply(data[,15:19],1,mean) par(mfrow=c(1,1)) boxplot(notamedia[dependenciaadministrativa=="Municipal"],notamedia[dependenciaadministrativa=="Estadual"], notamedia[dependenciaadministrativa=="Federal"],notamedia[dependenciaadministrativa=="Privada"], names=tipoescola,ylab="Media",main="Media das medias",ylim=c(400,770)) #pdf(file="enem2015-grafico2.pdf",width=12,height=10) par(mfrow=c(2,2)) plot(taxadeparticipacao[dependenciaadministrativa=="Municipal"], notamedia[dependenciaadministrativa=="Municipal"],xlab="Participacao",ylab="Media", ylim=c(400,770),main="Municipal",xlim=c(50,100)) plot(taxadeparticipacao[dependenciaadministrativa=="Estadual"], notamedia[dependenciaadministrativa=="Estadual"],xlab="Participacao",ylab="Media", ylim=c(400,770),main="Estadual",xlim=c(50,100)) plot(taxadeparticipacao[dependenciaadministrativa=="Federal"], notamedia[dependenciaadministrativa=="Federal"],xlab="Participacao",ylab="Media", ylim=c(400,770),main="Federal",xlim=c(50,100)) plot(taxadeparticipacao[dependenciaadministrativa=="Privada"], notamedia[dependenciaadministrativa=="Privada"],xlab="Participacao",ylab="Media", ylim=c(400,770),main="Privada",xlim=c(50,100)) #dev.off() #pdf(file="enem2015-grafico3.pdf",width=12,height=10) par(mfrow=c(2,2)) plot(indicadordeformacaodocente[dependenciaadministrativa=="Municipal"], notamedia[dependenciaadministrativa=="Municipal"],xlab="Formacao docente",ylab="Media", ylim=c(400,770),main="Municipal") plot(indicadordeformacaodocente[dependenciaadministrativa=="Estadual"], notamedia[dependenciaadministrativa=="Estadual"],xlab="Formacao docente",ylab="Media", ylim=c(400,770),main="Estadual") plot(indicadordeformacaodocente[dependenciaadministrativa=="Federal"], notamedia[dependenciaadministrativa=="Federal"],xlab="Formacao docente",ylab="Media", ylim=c(400,770),main="Federal") plot(indicadordeformacaodocente[dependenciaadministrativa=="Privada"], notamedia[dependenciaadministrativa=="Privada"],xlab="Formacao docente",ylab="Media", ylim=c(400,770),main="Privada") #dev.off() #pdf(file="enem2015-grafico4.pdf",width=12,height=10) par(mfrow=c(1,4)) boxplot( notamedia[(dependenciaadministrativa=="Municipal")&(indicadordepermanencianaescola=="Menos de 20%")], notamedia[(dependenciaadministrativa=="Municipal")&(indicadordepermanencianaescola=="De 20% a 40%")], notamedia[(dependenciaadministrativa=="Municipal")&(indicadordepermanencianaescola=="De 40% a 60%")], notamedia[(dependenciaadministrativa=="Municipal")&(indicadordepermanencianaescola=="De 60% a 80%")], notamedia[(dependenciaadministrativa=="Municipal")&(indicadordepermanencianaescola=="80% ou mais")], names=c("<20","20-40","40-60","60-80",">80"),ylim=c(400,770),main="Municipal") boxplot( notamedia[(dependenciaadministrativa=="Estadual")&(indicadordepermanencianaescola=="Menos de 20%")], notamedia[(dependenciaadministrativa=="Estadual")&(indicadordepermanencianaescola=="De 20% a 40%")], notamedia[(dependenciaadministrativa=="Estadual")&(indicadordepermanencianaescola=="De 40% a 60%")], notamedia[(dependenciaadministrativa=="Estadual")&(indicadordepermanencianaescola=="De 60% a 80%")], notamedia[(dependenciaadministrativa=="Estadual")&(indicadordepermanencianaescola=="80% ou mais")], names=c("<20","20-40","40-60","60-80",">80"),ylim=c(400,770),main="Estadual") boxplot( notamedia[(dependenciaadministrativa=="Federal")&(indicadordepermanencianaescola=="Menos de 20%")], notamedia[(dependenciaadministrativa=="Federal")&(indicadordepermanencianaescola=="De 20% a 40%")], notamedia[(dependenciaadministrativa=="Federal")&(indicadordepermanencianaescola=="De 40% a 60%")], notamedia[(dependenciaadministrativa=="Federal")&(indicadordepermanencianaescola=="De 60% a 80%")], notamedia[(dependenciaadministrativa=="Federal")&(indicadordepermanencianaescola=="80% ou mais")], names=c("<20","20-40","40-60","60-80",">80"),ylim=c(400,770),main="Federal") boxplot( notamedia[(dependenciaadministrativa=="Privada")&(indicadordepermanencianaescola=="Menos de 20%")], notamedia[(dependenciaadministrativa=="Privada")&(indicadordepermanencianaescola=="De 20% a 40%")], notamedia[(dependenciaadministrativa=="Privada")&(indicadordepermanencianaescola=="De 40% a 60%")], notamedia[(dependenciaadministrativa=="Privada")&(indicadordepermanencianaescola=="De 60% a 80%")], notamedia[(dependenciaadministrativa=="Privada")&(indicadordepermanencianaescola=="80% ou mais")], names=c("<20","20-40","40-60","60-80",">80"),ylim=c(400,770),main="Privada") #dev.off() #pdf(file="enem2015-grafico5.pdf",width=12,height=10) par(mfrow=c(1,4)) boxplot( notamedia[(dependenciaadministrativa=="Municipal")&(indicadornivelsocioeconomico=="Muito baixo")], notamedia[(dependenciaadministrativa=="Municipal")&(indicadornivelsocioeconomico=="Baixo")], notamedia[(dependenciaadministrativa=="Municipal")&(indicadornivelsocioeconomico=="Medio baixo")], notamedia[(dependenciaadministrativa=="Municipal")&(indicadornivelsocioeconomico=="Medio")], notamedia[(dependenciaadministrativa=="Municipal")&(indicadornivelsocioeconomico=="Medio alto")], notamedia[(dependenciaadministrativa=="Municipal")&(indicadornivelsocioeconomico=="Alto")], notamedia[(dependenciaadministrativa=="Municipal")&(indicadornivelsocioeconomico=="Muito alto")],ylim=c(400,770), names=c("MuB","B","MeB","Me","MeA","A","MuA"),main="Municipal") boxplot( notamedia[(dependenciaadministrativa=="Estadual")&(indicadornivelsocioeconomico=="Muito baixo")], notamedia[(dependenciaadministrativa=="Estadual")&(indicadornivelsocioeconomico=="Baixo")], notamedia[(dependenciaadministrativa=="Estadual")&(indicadornivelsocioeconomico=="Medio baixo")], notamedia[(dependenciaadministrativa=="Estadual")&(indicadornivelsocioeconomico=="Medio")], notamedia[(dependenciaadministrativa=="Estadual")&(indicadornivelsocioeconomico=="Medio alto")], notamedia[(dependenciaadministrativa=="Estadual")&(indicadornivelsocioeconomico=="Alto")], notamedia[(dependenciaadministrativa=="Estadual")&(indicadornivelsocioeconomico=="Muito alto")],ylim=c(400,770), names=c("MuB","B","MeB","Me","MeA","A","MuA"),main="Estadual") boxplot( notamedia[(dependenciaadministrativa=="Federal")&(indicadornivelsocioeconomico=="Muito baixo")], notamedia[(dependenciaadministrativa=="Federal")&(indicadornivelsocioeconomico=="Baixo")], notamedia[(dependenciaadministrativa=="Federal")&(indicadornivelsocioeconomico=="Medio baixo")], notamedia[(dependenciaadministrativa=="Federal")&(indicadornivelsocioeconomico=="Medio")], notamedia[(dependenciaadministrativa=="Federal")&(indicadornivelsocioeconomico=="Medio alto")], notamedia[(dependenciaadministrativa=="Federal")&(indicadornivelsocioeconomico=="Alto")], notamedia[(dependenciaadministrativa=="Federal")&(indicadornivelsocioeconomico=="Muito alto")],ylim=c(400,770), names=c("MuB","B","MeB","Me","MeA","A","MuA"),main="Federal") boxplot( notamedia[(dependenciaadministrativa=="Privada")&(indicadornivelsocioeconomico=="Muito baixo")], notamedia[(dependenciaadministrativa=="Privada")&(indicadornivelsocioeconomico=="Baixo")], notamedia[(dependenciaadministrativa=="Privada")&(indicadornivelsocioeconomico=="Medio baixo")], notamedia[(dependenciaadministrativa=="Privada")&(indicadornivelsocioeconomico=="Medio")], notamedia[(dependenciaadministrativa=="Privada")&(indicadornivelsocioeconomico=="Medio alto")], notamedia[(dependenciaadministrativa=="Privada")&(indicadornivelsocioeconomico=="Alto")], notamedia[(dependenciaadministrativa=="Privada")&(indicadornivelsocioeconomico=="Muito alto")],ylim=c(400,770), names=c("MuB","B","MeB","Me","MeA","A","MuA"),main="Privada") #dev.off() #pdf(file="enem2015-grafico6.pdf",width=12,height=10) par(mfrow=c(2,3)) boxplot( notamedia[(dependenciaadministrativa=="Municipal")&(indicadornivelsocioeconomico=="Medio baixo")], notamedia[(dependenciaadministrativa=="Estadual")&(indicadornivelsocioeconomico=="Medio baixo")], notamedia[(dependenciaadministrativa=="Federal")&(indicadornivelsocioeconomico=="Medio baixo")], notamedia[(dependenciaadministrativa=="Privada")&(indicadornivelsocioeconomico=="Medio baixo")], names=c("Mun","Est","Fed","Pri"),main="Medio baixo") boxplot( notamedia[(dependenciaadministrativa=="Municipal")&(indicadornivelsocioeconomico=="Medio")], notamedia[(dependenciaadministrativa=="Estadual")&(indicadornivelsocioeconomico=="Medio")], notamedia[(dependenciaadministrativa=="Federal")&(indicadornivelsocioeconomico=="Medio")], notamedia[(dependenciaadministrativa=="Privada")&(indicadornivelsocioeconomico=="Medio")], names=c("Mun","Est","Fed","Pri"),main="Medio") boxplot( notamedia[(dependenciaadministrativa=="Municipal")&(indicadornivelsocioeconomico=="Medio alto")], notamedia[(dependenciaadministrativa=="Estadual")&(indicadornivelsocioeconomico=="Medio alto")], notamedia[(dependenciaadministrativa=="Federal")&(indicadornivelsocioeconomico=="Medio alto")], notamedia[(dependenciaadministrativa=="Privada")&(indicadornivelsocioeconomico=="Medio alto")], names=c("Mun","Est","Fed","Pri"),main="Medio alto") boxplot( notamedia[(dependenciaadministrativa=="Municipal")&(indicadornivelsocioeconomico=="Alto")], notamedia[(dependenciaadministrativa=="Estadual")&(indicadornivelsocioeconomico=="Alto")], notamedia[(dependenciaadministrativa=="Federal")&(indicadornivelsocioeconomico=="Alto")], notamedia[(dependenciaadministrativa=="Privada")&(indicadornivelsocioeconomico=="Alto")], names=c("Mun","Est","Fed","Pri"),main="Alto") boxplot( notamedia[(dependenciaadministrativa=="Municipal")&(indicadornivelsocioeconomico=="Muito alto")], notamedia[(dependenciaadministrativa=="Estadual")&(indicadornivelsocioeconomico=="Muito alto")], notamedia[(dependenciaadministrativa=="Federal")&(indicadornivelsocioeconomico=="Muito alto")], notamedia[(dependenciaadministrativa=="Privada")&(indicadornivelsocioeconomico=="Muito alto")], names=c("Mun","Est","Fed","Pri"),main="Muito alto") #dev.off() #pdf(file="enem2015-grafico7.pdf",width=12,height=10) par(mfrow=c(1,2)) boxplot( notamedia[(dependenciaadministrativa=="Municipal")&((indicadornivelsocioeconomico=="Medio baixo")|(indicadornivelsocioeconomico=="Medio")|(indicadornivelsocioeconomico=="Medio alto"))], notamedia[(dependenciaadministrativa=="Estadual")&((indicadornivelsocioeconomico=="Medio baixo")|(indicadornivelsocioeconomico=="Medio")|(indicadornivelsocioeconomico=="Medio alto"))], notamedia[(dependenciaadministrativa=="Federal")&((indicadornivelsocioeconomico=="Medio baixo")|(indicadornivelsocioeconomico=="Medio")|(indicadornivelsocioeconomico=="Medio alto"))], notamedia[(dependenciaadministrativa=="Privada")&((indicadornivelsocioeconomico=="Medio baixo")|(indicadornivelsocioeconomico=="Privada")|(indicadornivelsocioeconomico=="Medio alto"))], names=c("Mun","Est","Fed","Pri"),main="Medio baixo/Medio/Medio alto",ylim=c(400,770)) boxplot( notamedia[(dependenciaadministrativa=="Municipal")&((indicadornivelsocioeconomico=="Alto")|(indicadornivelsocioeconomico=="Muito alto"))], notamedia[(dependenciaadministrativa=="Estadual")&((indicadornivelsocioeconomico=="Alto")|(indicadornivelsocioeconomico=="Muito alto"))], notamedia[(dependenciaadministrativa=="Federal")&((indicadornivelsocioeconomico=="Alto")|(indicadornivelsocioeconomico=="Muito alto"))], notamedia[(dependenciaadministrativa=="Privada")&((indicadornivelsocioeconomico=="Alto")|(indicadornivelsocioeconomico=="Muito alto"))], names=c("Mun","Est","Fed","Pri"),main="Alto/Muito alto",ylim=c(400,770)) #dev.off() par(mfrow=c(1,2)) boxplot(notamedia[uf=="DF"],notamedia[uf=="GO"],notamedia[uf=="MS"],notamedia[uf=="MT"], notamedia[uf=="ES"],notamedia[uf=="MG"],notamedia[uf=="RJ"],notamedia[uf=="SP"], notamedia[uf=="PR"],notamedia[uf=="RS"],notamedia[uf=="SC"], names=c("DF","GO","MS","MT","ES","MG","RJ","SP","PR","RS","SC"),ylim=c(400,770)) abline(h=500,lty=2) boxplot(notamedia[uf=="AC"],notamedia[uf=="AM"],notamedia[uf=="AP"],notamedia[uf=="PA"], notamedia[uf=="RO"],notamedia[uf=="RR"],notamedia[uf=="TO"],notamedia[uf=="AL"], notamedia[uf=="BA"],notamedia[uf=="CE"],notamedia[uf=="MA"], notamedia[uf=="PB"],notamedia[uf=="PE"],notamedia[uf=="PI"], notamedia[uf=="RN"],notamedia[uf=="SE"], names=c("AC","AM","AP","PA","RO","RR","TO","AL","BA","CE","MA","PB","PE","PI","RN","SE"),ylim=c(400,770)) abline(h=500,lty=2) dev.off() n = nrow(data) mun=rep(0,n) est=rep(0,n) fed=rep(0,n) pri=rep(0,n) tipo = rep(0,n) tipo[mun==1]=1 tipo[est==1]=2 tipo[fed==1]=3 tipo[pri==1]=4 mun[dependenciaadministrativa=="Municipal"]=1 est[dependenciaadministrativa=="Estadual"]=1 fed[dependenciaadministrativa=="Federal"]=1 pri[dependenciaadministrativa=="Privada"]=1 per2=rep(0,n) per3=rep(0,n) per4=rep(0,n) per5=rep(0,n) per2[indicadordepermanencianaescola=="De 20% a 40%"]=1 per3[indicadordepermanencianaescola=="De 40% a 60%"]=1 per4[indicadordepermanencianaescola=="De 60% a 80%"]=1 per5[indicadordepermanencianaescola=="80% ou mais"]=1 nivel2 = rep(0,n) nivel3 = rep(0,n) nivel4 = rep(0,n) nivel5 = rep(0,n) nivel6 = rep(0,n) nivel7 = rep(0,n) nivel2[indicadornivelsocioeconomico=="Baixo"]=1 nivel3[indicadornivelsocioeconomico=="Medio baixo"]=1 nivel4[indicadornivelsocioeconomico=="Medio"]=1 nivel5[indicadornivelsocioeconomico=="Medio alto"]=1 nivel6[indicadornivelsocioeconomico=="Alto"]=1 nivel7[indicadornivelsocioeconomico=="Muito alto"]=1 porte2 = rep(0,n) porte3 = rep(0,n) porte4 = rep(0,n) porte2[portedaescola=="De 31 a 60 alunos"]=1 porte3[portedaescola=="De 61 a 90 alunos"]=1 porte4[portedaescola=="Maior que 90 alunos"]=1 part=taxadeparticipacao reg = lm(notamedia~uf+dependenciaadministrativa+taxaparticipacao+indicadornivelsocioeconomico+indicadordepermanencianaescola+portedaescola) plot(reg$fit,reg$res) summary(reg) reg1 = lm(notamedia~uf+mun+fed+pri+per2+per3+per4+per5+nivel2+nivel3+nivel4+nivel5+nivel6+nivel7+porte2+porte3+porte4+part) plot(reg1$fit,reg1$res) summary(reg1) reg2 = lm(media~part+doc+pri+doc*pri) plot(reg2$fit,reg2$res) summary(reg2)