### Egyptian cobra world map
# use install.packages("maptools") if you have not already installed it

library(maptools)
data(wrld_simpl)
plot(wrld_simpl)
countries<-wrld_simpl$NAME 
Naja_countries<-c("Tanzania","Kenya","Somalia","Ethiopia","Uganda","South Sudan","Sudan","Cameroon","Nigeria","Niger","Burkina Faso","Mali","Senegal","Mauritania","Morocco","Algeria","Tunisia","Libya","Egypt")
N_countries<-length(Naja_countries)
country_numbers<-NULL
for(i in 1:N_countries) country_numbers<-c(country_numbers,grep(Naja_countries[i],countries))
col.map[country_numbers]<-"red"
plot(wrld_simpl,col=col.map)