Displaying 1 result from an estimated 1 matches for "ngadmin2".
Did you mean:
ladmin2
2011 Jan 13
0
Colouring areas on a map with updated maptools functions
...hich had the same
number of rows as the shapefile and one column populated with numbers from 1
to 8 which represented an attribute for each of those states. The code then
plotted those numbers as colours on the map so each state had the
appropriate colour. This was the code:
x<-read.shape("NGADMIN2.SHP")
xpolys<-Map2poly(x)
coverage<-read.table("coveragemar08.txt",header=T)
colour1<-coverage[,2]
palette(c("red4","red3","red2","orangered","orange","gold1","yellow","white"))
plot(xpolys, axes...