#Set Libraries
library(maps)
library(mapdata)
library(maptools)
library(scales)
library(gdata)
#Get the shapefile for Afghan Districts
area <- readShapePoly("/Users/Tim/Downloads/af_new/af.shp")
#Get data from Excel file
dataset<-read.table("/Users/Tim/Downloads/Workbook2.csv")
dataset2<-as.numeric(t(dataset))
#Plot outline map of Afghanistan
map("world","Afghanistan")
#Associate information with each of the data points in the shapefile
area@data$count <- dataset2
#Set the two colours you are varying between
rbPal <- colorRampPalette(bias=1,c('green','yellow'))
#Split the colours between the ranges of data between the colour choices
area@data$Col <- rbPal(4)[as.numeric(cut(area@data$count, seq(0, 2500, by = 625)))]
area@data$bin <- cut(area@data$count, seq(0, 2500, by = 625),
include.lowest = TRUE, dig.lab = 4)
lev = levels(area@data$bin)
lev2 <- gsub("\\,", " to ", lev)
lev3 <- gsub("\\]$", "", lev2)
lev4 <- gsub("\\(|\\)", " ", lev3)
lev5 <- gsub("^\\[", " ", lev4)
my.levels <- lev5
#Plot the map
plot(area,col=area@data$Col)
legend("topleft", fill = area@data$Col, legend = my.levels, col = area@data$Col)
