In jcpernias/ec1047: Code and data for the EC1047 course

knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)

library(tibble)
library(ggplot2)
library(dplyr, warn.conflicts = FALSE)
library(ec1047)

The data

We are going to use two vectors:

y1 <- c(2417, 7800, 8489, 10072, 12957)
y2 <- c(4417, 7800, 9989, 10572, 8957)

Compute the coordinates of the Lorenz curves:

lorenz_y1 <- lorenz(y1)
lorenz_y2 <- lorenz(y2)

Merge the databases:

lorenz_all <- bind_rows(lorenz_y1 %>% add_column(Variable = "y1"), 
                        lorenz_y2 %>% add_column(Variable = "y2"))

The plot

Use ggplot2 to plot the Lorenz curves:

lorenz_all %>% 
  ggplot(aes(x = x, y = y, color = Variable)) + 
  geom_line() + 
  geom_line(data = tibble(x = c(0, 1), y = c(0, 1)), color = 'black') + 
  coord_fixed() +
  scale_x_continuous(limits = c(0, 1), breaks = (0:5)/5) + 
  scale_y_continuous(limits = c(0, 1), breaks = (0:5)/5) +
  theme_classic() +
  theme(axis.title.y = element_text(angle = 0),
        panel.grid.major = element_line(color = "gray95"),
        panel.grid.minor = element_line(color = "gray85", size = 0.1))