Histograma básico

hist(ehpm2017$money)

A veces, para asegurarnos que tenemos los casos, nuestro análisis lo hacemos filtrando a quiénes reportaron estar ocupados, para no cargar con los ingresos cero.

hist(ehpm2017[ehpm2017$actpr2012==10,]$money)

Podemos modificarlo de a poco a poco. Modificar el título

hist(ehpm2017[ehpm2017$actpr2012==10,]$money, main="Histograma de Ingresos por trabajo") 

Le podemos modificar el título del eje de las x y de las y

hist(ehpm2017[ehpm2017$actpr2012==10,]$money, 
     main="Histograma de Ingresos por trabajo", 
     xlab="money", ylab="Frecuencia") 

¡A ponerle colorcitos! Recuerda aquí hay una lista http://www.stat.columbia.edu/~tzheng/files/Rcolor.pdf

hist(ehpm2017[ehpm2017$actpr2012==10,]$money, 
     main="Histograma de Ingresos por trabajo",
     xlab="ingresos por trabajo", 
     ylab="Frecuencia", col="deeppink1") 

con pipes

ehpm2017 %>% 
  filter(ehpm2017$actpr2012==10) %>% 
    with(hist(money))

ehpm2017 %>% 
  filter(ehpm2017$actpr2012==10) %>% 
    with(hist(money))

Cuando usamos pipes, se debe de recordar que no es necesario escribir el nombre del data.frame porque es lo primero que colocamos en uestro “pipe”.

Checa que cualquier aditamiento debe ir en el pipe donde está el comando de hist(). Ten cuidado con los paréntesis.

ehpm2017 %>% 
  filter(ehpm2017$actpr2012==10) %>% 
    with(hist(money, main= "histograma"))

Transformando

hist(log(ehpm2017[ehpm2017$actpr2012==10,]$money)) #ojo con los paréntesis

ehpm2017 %>% 
  filter(ehpm2017$actpr2012==10) %>% 
    with(hist(log(money), main= "histograma"))

¿Qué pasa con el logaritmo?

table(is.infinite(ehpm2017$money))

## 
## FALSE 
## 75133

ehpm2017$log_money<-log(ehpm2017$money)
table(is.infinite(ehpm2017$log_money))

## 
## FALSE  TRUE 
## 26125 49008

Hay varias funciones interesantes para revisar nuestros datos. Busca en la ayuda is.na() is.nan() is.infinite() complete.cases() na.omit()

Iremos utilizándolas a lo largo del curso

Gráficos univariados

Podemos hacer histogramas y gráficos de densidad, de manera fácil. La idea es agregar en nuestro “lienzo” una geometría, un valor para dibujar en él. Esto se agrega con un “+” y con la figura que se añadirá a nuestro gráfico.

ggplot(ehpm2017[ehpm2017$actpr2012==10,], aes(x=money)) + geom_histogram()

## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

## Warning: Removed 6286 rows containing non-finite values (stat_bin).

Si queremos usar los “pipes”, debemos crear el filtro de la base

ehpm2017 %>% 
  filter(ehpm2017$actpr2012==10) %>% 
    ggplot(aes(x=money)) +
       geom_histogram()

## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

Como se puede observar, la lógica de las dos sintaxis de programación son consistentes. Elige tu favorita.

¡Tip!. Puedes guardar el lienzo. Esto puede ayudarte a después ir sumando las geometríasy demás opciones que iremos aprendiendo.

lienzo<-ehpm2017 %>% 
  filter(ehpm2017$actpr2012==10) %>% 
    ggplot(aes(x=money)) 

lienzo

lienzo + geom_histogram()

## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

Podemos ir añadiendo más opciones

lienzo + geom_histogram(color="black", fill="white")

## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

lienzo + geom_histogram(binwidth=50)

lienzo + geom_histogram(bins=20)

Theme Assist Addin

Existe un add-in - es decir un complemento - para que podamos modificar más fácilmente los gráficos que hacemos con ggplot. Lo tenemos que instalar. Así como hemos instalado nuestros otras paqueterías.

#install.packages("ggThemeAssist",  dependencies = TRUE)
library(ggThemeAssist)

Para utilizarlo, tenemos que nombrar el objeto que queremos editar (por eso aprendimos a guardarlo en un objeto anteriormente). Lo escribimos. Lo seleccionamos y buscamos el complemento en la pestaña “Addins”

image

Agregando subgrupos a nuestros gráficos

Con subgrupos podemos agregar otra variable en nuestro gráficos

lienzo2 + geom_jitter(aes(shape=as_label(region)))

lienzo2 + geom_jitter(aes(colour =as_label(region)))

lienzo2 + 
  geom_point()+
  geom_smooth(method=lm)

lienzo2 + 
  geom_point()+
  geom_smooth(method=lm) +
  facet_grid(. ~ as_label(region))

lienzo2 + 
  geom_point()+
  geom_smooth(method=lm) +
  facet_grid(as_label(region) ~ .)

lienzo2 + 
  geom_point()+
  geom_smooth(method=lm) +
  facet_wrap(~as_label(region))

lienzo2 + 
  geom_text(aes(label=as_label(region)))

Combinando varias geometrías

lienzo2 + 
  geom_point(alpha=0.2)+
  geom_smooth(method=lm) +
  facet_wrap(~as_label(region)) +
  theme_light()