Skip to main content

Modifying Some Plots in R

Grant Data: 500 Households that get Social Grant in a certain region, the data comes from Social and Economic Survey. Download the dummy data: DATA
grant <- read.csv("grant.csv")

#PFOODEXP: Proportion of Food Expenditure to Total Expenditure
#HH_Income: Household Income ($)
#HH_FOOD: Household Food Expenditure ($)
#HH_Loc: Household Location (0: Rural, 1: Urban)
#Educ_H: Household Head Education (Year)
#HH_Size: Household Family Member
#Gender_H: Household Head Gender (0: Female, 1: Male )
#Age_H: Household Head Age

#Remove 1st column
grant[1] <- NULL  
#change variable type to factor
grant$Gender_H <- as.factor(grant$Gender_H) 
grant$HH_Loc <- as.factor(grant$HH_Loc)

head(grant)
##   PFOODEXP  HH_Food HH_Loc Gender_H Age_H Educ_H HH_Size HH_Income
## 1 73.16933 675.1766      1        0    32     22       4  922.7590
## 2 69.77417 596.2286      1        1    27     22      13  854.5119
## 3 63.31992 585.4946      0        1    42     15       8  924.6610
## 4 71.77146 559.9470      0        1    26     22       5  780.1806
## 5 52.36861 554.8954      1        0    32      9       3 1059.5954
## 6 79.74131 551.4373      0        1    21     15       5  691.5328
summary(grant)
##     PFOODEXP        HH_Food      HH_Loc  Gender_H     Age_H      
##  Min.   :29.12   Min.   :165.9   0:380   0:167    Min.   :15.00  
##  1st Qu.:56.43   1st Qu.:189.0   1:120   1:333    1st Qu.:25.00  
##  Median :65.35   Median :218.3                    Median :34.00  
##  Mean   :64.74   Mean   :251.1                    Mean   :34.16  
##  3rd Qu.:75.59   3rd Qu.:286.6                    3rd Qu.:42.00  
##  Max.   :86.98   Max.   :675.2                    Max.   :80.00  
##      Educ_H         HH_Size         HH_Income     
##  Min.   : 3.00   Min.   : 1.000   Min.   : 198.6  
##  1st Qu.: 6.00   1st Qu.: 4.000   1st Qu.: 280.7  
##  Median :15.00   Median : 5.000   Median : 351.0  
##  Mean   :12.57   Mean   : 5.392   Mean   : 406.6  
##  3rd Qu.:16.00   3rd Qu.: 6.000   3rd Qu.: 489.0  
##  Max.   :23.00   Max.   :20.000   Max.   :1059.6
#Scatter Plot with Linear Line

library(ggplot2)
library(ggthemes)

ggplot(data=grant, aes(HH_Income,PFOODEXP, colour = Gender_H, size = HH_Size)) + 
  geom_point(alpha=0.8) + geom_smooth(method = "lm", se=FALSE) +
  ylab("% Food Expenditure")+ xlab("Household Income per Month ($)") +
  guides(color = guide_legend(override.aes = list(size=5, linetype = c(0,0)), title = "HH Gender"), 
         size = guide_legend(override.aes = list(linetype = c(0,0)), title = "H Size")) +
  scale_size_continuous(range = c(1, 8), breaks = c(1, 2, 4, 8))+
  scale_color_manual(labels = c("female","male"), values = c("hotpink","deepskyblue"))+
  labs( col = "Gender") +
  ggtitle("Scatterplot of Percentage of Food Expenditure vs Household Income per Month ($)") +
  scale_x_log10()+
  theme_bw()+
  theme(plot.title = element_text(size=10, face= "bold"))
#Scatter Plot with LOESS

library(ggplot2)
library(ggthemes)

ggplot(data=grant, aes(HH_Income,PFOODEXP, shape = Gender_H, colour = Gender_H, size = HH_Size)) + 
  geom_point(alpha=0.8) + geom_smooth(method = "loess") +
  ylab("% Food Expenditure")+ xlab("Household Income per month ($)") +
  guides(colour = FALSE, 
         size = FALSE,
         shape = guide_legend(override.aes = 
                 list(size=5, linetype = c(0,0), 
                      colour = c("azure4","gold")), title = "HH Gender")) +
  scale_size_continuous(range = c(1, 8), breaks = c(1, 2, 4, 8))+
  scale_shape_manual(labels = c("female","male"), values = c("f","m"))+
  scale_color_manual(labels = c("female","male"), values = c("azure4","gold"))+
  ggtitle("Scatterplot of Percentage of Food Expenditure vs Household Income per Month ($)") +
  scale_x_log10()+
  theme_bw()+
  theme(plot.title = element_text(size=10, face= "bold"))
Scatter Plot with LOESS: One of advantage of LOESS method, it doesn’t require the specification of a function to fit a model to all of the data in the sample. One of disadvantage, it doesn’t generate a regression function that is easily represented by mathematical formula.
#Hexbin Plot

library(hexbin)
library(RColorBrewer)

# Create data
  y<-grant%>% pull("HH_Food")
  x<-grant%>% pull("HH_Income")
  
# Make the plot
  bin<-hexbin(x, y)
  rf=colorRampPalette(rev(brewer.pal(10,'Spectral')))
  
  hexbinplot(y~x, data=bin, main="Income vs Food Expenditure",
             colramp=rf, trans=log, inv=exp,mincnt=1, maxcnt=70,
             ylab="food expenditure ($)",
             xlab="income ($)", cex.label=0.7)
Hexbin Plot: Scatterplots can get very hard to interpret when displaying large datasets, as points inevitably overplot. Hexbinplot helps discerning the data individually. Code source: www.everydayanalytics.ca
#Density Plot

qplot(HH_Income,data = grant, geom="density", fill = HH_Loc ,alpha=I(.5), 
      ylab="Density",
      xlab= "Household Income($)", 
      main = "Distribution of Household (HH) Income per Month by Household Location") + 
      scale_fill_manual(labels = c("rural","urban"), values = c("tomato","mediumspringgreen"))+
      labs( fill = "HH Location") + geom_density(alpha= 0.2,aes(HH_Income), colour = "grey85")+
      theme_minimal() + theme(plot.title = element_text(size=10, face= "bold"))
Income distribution is often right skewed, this shows income inequality. The hypothesized reasons are differences in talents, skills, and opportunities. It is not surprising, household income distribution in urban area is more skewed than in rural area.

Comments

Post a Comment

Popular posts from this blog

How to Create Indonesia Map in R

Creating the Map In this article, I will try to explain how to make Indonesia Map in R. I will assume that you are already familiar with the basic codes in R. First, we need the required libraries : require (maps) #loading maps package require (mapdata) #loading mapdata package library(ggplot2) #ggplot2 package library(readxl) #package for read .xlsx file library(ggthemes) #package for ggplot2 theme library(ggrepel) #extendig the plotting package ggplot2 for maps Then, we prepare the data that contains the information of provinces name, latitude, and longitude of every province in Indonesia, e.g. : You can download the data in here:  Data Now open the file and create the polygon: setwd( "your file's path" ) #set your own directory mydata<- read _xlsx( "dummy.xlsx" ) #assign the data to "mydata" View(mydata) #view the data, notice the column of "latitude","longitude", "woe_label" glo...
1 comment
Read more

What Can We Learn from Greek Debt Dramas?

Greek Debt Dramas Before the Global Financial Crisis (GFC) in 2008, the Greek had positive economic growth and it was considered high among countries in eurozone. Average economic growth reached almost four per cent between 1999 and 2007. Then the crisis hit in 2007 where housing bubble burst and made the subprime mortgage market in the United State collapsed. The crisis in the U.S. created a chain reaction which causing global banking crisis and credit crunch that lasts through 2009. The crisis made Lehman Brothers, big financial company, collapsed and the government in the United States and Europe prepared to bail out their banks. Greece failed to pay their huge debt since borrowing costs rose and financing dried up.  The financial crisis affected the Greek economy by reducing financial liquidity and business activity. Greece had been fortunate enough to face the crisis with the euro instead of its national currency, if they were using their national currency the crisis wo...
Post a Comment
Read more

Empirical Evidence of Engel’s Law Among Social Grant Recievers

Engel's law is an observation in economics stating that as income increases, the proportion of income spent on food decreases, even if absolute expenditure on food increases. The law was named after the statistician  Ernst Engel (1821–1896). One application of this statistic is treating it as a reflection of the living standard of a country. As this proportion — or "Engel coefficient" — increases, the country is by nature poorer; conversely a low Engel coefficient indicates a higher standard of living. Engel's Law image source: Wikipedia Using data collected through National Social and Economic Survey (NSES) by BPS-Statistics Indonesia, I tried to examine the existence of Engel's Law among households that received social grants in West Papua-Indonesia. Some studies found that giving additional money to the low-income households resulted in an increase in overall expenditure on food (on absolute) but the proportion  of income spent on food would decrea...
Post a Comment
Read more