Gráficos usando comando plot

Comando plot

El comando plot nos permite generar gráficos en 2D de los datos que nosotros le entreguemos, y que necesitamos visualizar en gráficos.

Para este tutorial vamos a usar la base de datos de la universidad de Princenton, en donde analizaremos las variables:

Para esto utilizaremos el siguiente código para leer la base de datos. 

bdata= read.csv("http://www.princeton.edu/~otorres/mydataviz.csv", header = TRUE, stringsAsFactors = FALSE)
bdata<-na.omit(bdata)
#opcional si desean visualizar la data
View(bdata)

Para realizar el análisis vamos a responder las siguientes preguntas, que se encuentran en blog de Cesar Huamaní.

Pregunta 1: Un gráfico de tasa de desempleo para el partido demócrata y otro para el partido republicano.

Para desarrollar esta pregunta primero debemos separar la data por demócratas y republicanos de la siguiente forma. 

democrata<-bdata[bdata$party=="Democrat",]
#creamos un vector democrata para realizar el analisis y haremos otro para republicanos
republicano<-bdata[bdata$party=="Republican",]

Una vez que se realizo la separacion de datos, pasamos a realizar los graficos de la siguiente manera 

#para el gobierno democrata
plot(democrata$unemp, type = "l", col="blue", main="Gobierno Democrata", adj=0, ylab="Tasa de Desempleo (%)")
abline(h=mean(democrata$unemp), col="red", lty=2, lwd=3)
legend("topleft", legend = c("Tasa de desempleo medio 5.75%"), col = c("red"),lty=2, cex = 0.8)

#para el gobierno republicano
republicano<-bdata[bdata$party=="Republican",]
plot(republicano$unemp, type = "l", col="blue", main = "Gobierno Republicano",adj=0, ylab="Tasa de Desempleo (%)")
abline(h=mean(republicano$unemp), col="red", lty=2, lwd=3)
legend("topright", legend = c("Tasa de desempleo medio 5.76%"), col = c("red"),lty=2, cex = 0.8)

Pregunta 2: Realizar la gráfica del modelo lineal entre el nivel de desempleo y el nivel de aprobación

Con los datos separados se puede calcular y graficar el modelo lineal de la siguiente forma. 

#demorcratas
plot(democrata$unemp, democrata$approve, main="Gobierno Democrata", col="blue", ylab = "Nivel de Aprobación", xlab="Nivel de Desempleo")
abline(lm(democrata$approve~democrata$unemp, data = democrata),col="red", lwd=4)
legend("topright", legend = c("Linea de regresion"), col = c("red"),lty=1, cex = 0.8)

#republicanos
plot(republicano$unemp, republicano$approve,main="Gobierno Republicano", col="blue", ylab = "Nivel de Aprobación", xlab="Nivel de Desempleo")
abline(lm(republicano$approve~republicano$unemp, data = republicano),col="red", lwd=4)
legend("topright", legend = c("Linea de regresion"), col = c("red"),lty=1, cex = 0.8)

Finalmente quedaria de esta forma nuestros gráficos y el análisis sobre los datos.