Diagrama de Dispersión y Regresión en R
Manuel Sigüeñas Gonzales - Universidad Nacional Agraria La Molina
26 de julio de 2015
Packages
En esta sección se trabajará con los siguientes paquetes: “ggplot2”, “foreign” y “psych”.
###Librerias###
library(foreign)
library(ggplot2)
library(psych)##
## Attaching package: 'psych'
##
## The following object is masked from 'package:ggplot2':
##
## %+%Data Interés por la Política (I)
###Cargamos la data BaseLatinobarómetro2013-1.txt###
BaseIntPol <- read.table("D:/Manuel/MAESTRÍA EN ESTADÍSTICA APLICA UNALM/SEMESTRE I/Regresión Aplicada/Capítulo 1 - Diagrama de Dispersión y Regresión/Aplicación/BaseLatinobarómetro2013-1.txt", header=TRUE, quote="\"")
###Mostramos las variables###
names(BaseIntPol)## [1] "País" "DesConf" "IntPol"###Mostramos la matriz de datos###
(BaseIntPol)## País DesConf IntPol
## 1 Argentina 37.7 44.5
## 2 Bolivia 46.8 50.6
## 3 Brasil 48.2 33.6
## 4 Chile 48.3 52.0
## 5 Colombia 45.8 49.2
## 6 Costa.Rica 57.1 49.4
## 7 Ecuador 37.7 41.9
## 8 El.Salvador 45.5 38.4
## 9 Guatemala 49.1 48.0
## 10 Honduras 55.3 40.4
## 11 México 43.4 33.0
## 12 Nicaragua 39.2 40.3
## 13 Panamá 40.7 34.7
## 14 Paraguay 31.5 25.0
## 15 Perú 50.9 53.8
## 16 Rep.Dominicana 48.9 46.3
## 17 Uruguay 29.0 35.1
## 18 Venezuela 31.9 23.5Gráfica 01.1.: Diagrama de Dispersión entre Desconfianza hacia los Partidos Políticos y Desinterés por la Política
Diagrama de Dispersión Simple
###Gráfica simple###
plot(BaseIntPol[,2], BaseIntPol[,3], main="Gráfica 01.1a.: Diagrama de Dispersión")
abline(lm(BaseIntPol[,3]~BaseIntPol[,2]), col="yellow")
Diagrama de Dispersión con Etiqueta de Variables
###Gráfica con etiqueta de variables###
plot(IntPol ~ DesConf, data = BaseIntPol, xlab = "Desconfianza hacia los partidos políticos", ylab = "Desinterés por la política", xlim = c(29, 57.1), ylim = c(25, 53.8), main="Gráfica 01.1b.: Diagrama de Dispersión")
MLatin <- lm(IntPol ~ DesConf, data = BaseIntPol)
abline(MLatin, col = "red")
Diagrama de Dispersión con ggplot2
###Gráfica con Packages ggplot2###
###Además de la línea de regresión se puede apreciar los intervalos de confianza###
library(ggplot2)
ggplot(BaseIntPol, aes(x=DesConf, y=IntPol)) + geom_point() + ggtitle("Gráfica 01.1c.: Diagrama de Dispersión") + xlab("Desconfianza hacia los partidos políticos") + ylab("Desinterés por la política") + geom_smooth(method=lm)
Gráfica 01.2.: Diagrama de Dispersión y Residuos entre Desconfianza Hacia Los Partidos Políticos y Desinterés Por La Política
par(mfrow = c(2, 2))
plot(IntPol ~ DesConf, data = BaseIntPol, xlab = "Desconfianza hacia los partidos políticos (a)", ylab = "Desinterés por la política", main="Gráfico 01.2a")
m0 <- lm(IntPol ~ DesConf, data = BaseIntPol)
abline(m0)
plot(BaseIntPol$DesConf, residuals(m0), xlab = "Desconfianza hacia los partidos políticos (b)", ylab = "Residuales", main="Gráfico 01.2b")
abline(h = 0, lty = 2)
plot(log10(IntPol) ~ DesConf, data = BaseIntPol, xlab = "Desconfianza hacia los partidos políticos (c)", ylab = "log(Desinterés por la política)", main="Gráfico 01.2c")
m1 <- lm(log10(IntPol) ~ DesConf, data = BaseIntPol)
abline(m1)
plot(BaseIntPol$DesConf, residuals(m1), xlab = "Desconfianza hacia los partidos políticos (d)", ylab = "Residuales", main="Gráfico 01.2d")
abline(h = 0, lty = 2)
Gráfica 01.3.: Función Media
###Función media###
plot(IntPol ~ DesConf, data = BaseIntPol, xlab = "Desconfianza hacia los partidos políticos", ylab = "Desinterés por la política", xlim = c(29, 57.1), ylim = c(25, 53.8), main="Gráfico 01.3: Función Media" )
abline(a = 0, b = 1, lty = 2)
MLatin <- lm(IntPol ~ DesConf, data = BaseIntPol)
abline(MLatin, col = "green")
Data Ausentismo en la Primera Vuelta de las Elecciones Presidenciales de 2011
###Cargamos la data BaseAusentismoPres2011-Perú.sav###
library(foreign)
Ausentismo.Pres2011 <- as.data.frame(read.spss("D:/Manuel/MAESTRÍA EN ESTADÍSTICA APLICA UNALM/SEMESTRE I/Regresión Aplicada/Capítulo 1 - Diagrama de Dispersión y Regresión/Aplicación/BaseAusentismoPres2011-Perú.sav"))## Warning in read.spss("D:/Manuel/MAESTRÍA EN ESTADÍSTICA APLICA UNALM/
## SEMESTRE I/Regresión Aplicada/Capítulo 1 - Diagrama de Dispersión y
## Regresión/Aplicación/BaseAusentismoPres2011-Perú.sav"): D:/Manuel/MAESTRÍA
## EN ESTADÍSTICA APLICA UNALM/SEMESTRE I/Regresión Aplicada/Capítulo 1 -
## Diagrama de Dispersión y Regresión/Aplicación/BaseAusentismoPres2011-
## Perú.sav: Unrecognized record type 7, subtype 18 encountered in system file###Mostramos las variables###
names(Ausentismo.Pres2011)## [1] "Ubigeo"
## [2] "Departamento"
## [3] "Provincia"
## [4] "Distrito"
## [5] "NumeroDePersonasTotalElectoral"
## [6] "NumeroDePersonasAusentismo"
## [7] "PorcentajeDeAusentismoRespectoAlTotalElectoral"
## [8] "PorcentajeDeAusentismoOrdinal"
## [9] "PorcentajeDeAusentismoDicotomica"
## [10] "PorcentajeDeHogaresQueNoCuentanConMediosDeComunicacion"
## [11] "IDH"
## [12] "Alfabetismo"
## [13] "PartidoGanador"
## [14] "PorcentajeDeHogaresEnViviendasCon2oMasNBI"
## [15] "GINI"Gráfica 01.4.: Análisis de los Diagramas de Dispersión
###Análisis de los diagramas de dispersión###
###Diagramas de dispersión###
par(mfrow = c(2, 2))
plot(PorcentajeDeAusentismoRespectoAlTotalElectoral ~ PorcentajeDeHogaresQueNoCuentanConMediosDeComunicacion, data = Ausentismo.Pres2011, xlab = "% de hogares sin medios de comunicación", ylab = "% Ausentismo", main="Gráfico 01.4a")
M0 <- lm(PorcentajeDeAusentismoRespectoAlTotalElectoral ~ PorcentajeDeHogaresQueNoCuentanConMediosDeComunicacion, data = Ausentismo.Pres2011)
abline(M0, col = "red")
plot(PorcentajeDeAusentismoRespectoAlTotalElectoral ~ IDH, data = Ausentismo.Pres2011, xlab = "Índice de desarrollo humano", ylab = "% Ausentismo", main="Gráfico 01.4b")
M1 <- lm(PorcentajeDeAusentismoRespectoAlTotalElectoral ~ IDH, data = Ausentismo.Pres2011)
abline(M1, col = "blue")
plot(PorcentajeDeAusentismoRespectoAlTotalElectoral ~ Alfabetismo, data = Ausentismo.Pres2011, xlab = "% de Alfabetismo", ylab = "% de Ausentismo", main="Gráfico 01.4c")
M2 <- lm(PorcentajeDeAusentismoRespectoAlTotalElectoral ~ Alfabetismo, data = Ausentismo.Pres2011)
abline(M2, col = "green")
plot(PorcentajeDeAusentismoRespectoAlTotalElectoral ~ PorcentajeDeHogaresEnViviendasCon2oMasNBI, data = Ausentismo.Pres2011, xlab = "% de hogares en viviendas con 2 ó más NBI", ylab = "% de Ausentismo", main="Gráfico 01.4d")
M3 <- lm(PorcentajeDeAusentismoRespectoAlTotalElectoral ~ PorcentajeDeHogaresEnViviendasCon2oMasNBI, data = Ausentismo.Pres2011)
abline(M2, col = "yellow")
Data Interés por la Política (II)
###Cargamos la data BaseLatinobarómetro2013-2.txt###
BaseIntPol2 <- read.table("D:/Manuel/MAESTRÍA EN ESTADÍSTICA APLICA UNALM/SEMESTRE I/Regresión Aplicada/Capítulo 1 - Diagrama de Dispersión y Regresión/Aplicación/BaseLatinobarómetro2013-2.txt", header=TRUE, quote="\"")
###Mostramos las variables###
names(BaseIntPol2)## [1] "País" "ProbNeg" "DesConf" "DesPre" "IntPol"###Mostramos la matriz de datos2###
(BaseIntPol2)## País ProbNeg DesConf DesPre IntPol
## 1 Argentina 24.2 37.7 54.3 44.5
## 2 Bolivia 17.1 46.8 39.3 50.6
## 3 Brasil 21.8 48.2 39.8 33.6
## 4 Chile 23.5 48.3 67.7 52.0
## 5 Colombia 29.2 45.8 44.4 49.2
## 6 Costa.Rica 35.2 57.1 76.1 49.4
## 7 Ecuador 10.8 37.7 22.8 41.9
## 8 El.Salvador 17.2 45.5 30.6 38.4
## 9 Guatemala 35.1 49.1 48.2 48.0
## 10 Honduras 47.0 55.3 65.8 40.4
## 11 México 31.3 43.4 49.5 33.0
## 12 Nicaragua 9.1 39.2 36.0 40.3
## 13 Panamá 14.8 40.7 52.3 34.7
## 14 Paraguay 9.1 31.5 68.6 25.0
## 15 Perú 23.3 50.9 57.6 53.8
## 16 Rep.Dominicana 9.3 48.9 23.5 46.3
## 17 Uruguay 18.1 29.0 35.8 35.1
## 18 Venezuela 22.3 31.9 49.2 23.5Gráfica 01.5.: Matriz de Dispersión de la Data Interés por la Política (II)
###Gráfico matriz de dispersión simple###
pairs(IntPol ~ DesConf + ProbNeg + DesPre, data=BaseIntPol2, main="Gráfico 01.5: Matriz de Dispersión del Estudio Desinterés Por la Política")
Gráfica 01.5.: Matriz de Dispersión, Histograma y Correlación de la Data Interés por la Política (II)
###Primero se debe excluir la variable País###
var.fuera <- names(BaseIntPol2) %in% c("País")
Nuevo.BaseIntPol2 <- BaseIntPol2[!var.fuera]
###Se elabora mediante la una matriz de gráficos de dispersión donde se muestra de una sola mirada las dispersiones entre dos variables, la distribución de las variables en estudio y las correlaciones entre dos variables###
library(psych)
pairs.panels(Nuevo.BaseIntPol2, pch=21,main="Gráfico 01.6: Matriz de Dispersión, Histograma y Correlación")