R Notebook for Mapping Choropleth map of full India

This is an R Markdown Notebook. When you execute code within the notebook, the results appear beneath the code.

Try executing this chunk by clicking the Run button within the chunk or by placing your cursor inside it and pressing Cmd+Shift+Enter.

Load required library

library(ggplot2)
library(maptools)

Checking rgeos availability: TRUE

library(rgeos)

rgeos version: 0.3-26, (SVN revision 560)
 GEOS runtime version: 3.6.1-CAPI-1.10.1 r0 
 Linking to sp version: 1.2-5 
 Polygon checking: TRUE 

library(ggmap)

Google Maps API Terms of Service: http://developers.google.com/maps/terms.
Please cite ggmap if you use it: see citation('ggmap') for details.

library(scales)
library(RColorBrewer)
set.seed(8000)

Load the India shapefile from the directory

shp <- readShapeSpatial('~/Dropbox/PhD_Arka/Exploratory Indian Data Analysis/India Shapefile With Kashmir/india_shapefile_git/Admin2.shp')

use rgdal::readOGR or sf::st_readuse rgdal::readOGR or sf::st_read

If you just plot the shapefile, remove the # and run the code, otherwise skip the code

#plot(shp)

If you want to Look into the shapefile attributes and content names remove # and run the code, otherwise skip the code

#names(shp)
#print(shp$ST_NM)

Craete a csv file with attribute state name (here ‘id’) and count of data (here ‘count’) and then read the csv

imr=read.csv('~/Dropbox/PhD_Arka/Exploratory Indian Data Analysis/health stat/IMR_R.csv')

Fortify shape into a dataframe

shp.f <- fortify(shp, region = "ST_NM")

Merge shapefile with csv file

merge.shp.coef<-merge(shp.f,imr, by="id", all.x=TRUE)
final.plot<-merge.shp.coef[order(merge.shp.coef$order), ]

Creating the map

ggplot()+geom_polygon(data = final.plot,aes(x = long, y = lat, group = group, fill = count),color = "black", size = 0.25) + coord_map()
ggsave("India_IMR_2013_BLUE.png",dpi = 300, width = 20, height = 20, units = "cm")

Creating a map with different color

ggplot() +
  geom_polygon(data = final.plot,
               aes(x = long, y = lat, group = group, fill = count), 
               color = "black", size = 0.25) + 
  coord_map()+
  scale_fill_distiller(name="IMR", palette = "YlGn")+
  labs(title="IMR of Different States of India 2013")+
  xlab('Longitude')+
  ylab('Latitude')
ggsave("India_IMR_2013_YlGn.png",dpi = 300, width = 20, height = 20, units = "cm")

Creating a map with gradient fill

ggplot() +
  geom_polygon(data = final.plot,
               aes(x = long, y = lat, group = group, fill = count), 
               color = "black", size = 0.25) + 
  coord_map()+
  scale_fill_gradient(name="IMR", limits=c(0,100), low = 'white', high = 'red')+
  labs(title="IMR of Different States of India 2013")+
  xlab('Longitude')+
  ylab('Latitude')
ggsave("India_IMR_2013_NO_ST_NAME.png",dpi = 300, width = 20, height = 20, units = "cm")

Creating map with states name Aggregating the location for mapping

cnames <- aggregate(cbind(long, lat) ~ id, data=final.plot, FUN=function(x) mean(range(x)))

plotting the location in the map

ggplot() +
  geom_polygon(data = final.plot,
               aes(x = long, y = lat, group = group, fill = count), 
               color = "black", size = 0.25) + 
  coord_map()+
  scale_fill_gradient(name="IMR", limits=c(0,100), low = 'white', high = 'red')+
  labs(title="IMR of Different States of India 2013")+
  xlab('Longitude')+
  ylab('Latitude')+
  geom_text(data=cnames, aes(long, lat, label = id), size=3, fontface="bold")
ggsave("India_IMR_2013.png",dpi = 300, width = 20, height = 20, units = "cm")