求助如何用R画地图-经管之家官网！

上课时候老师给了代码如下：
library(raster)
library(rgdal)
CN.map <- readOGR(dsn="bou2_4p.shp",layer="bou2_4p")
projection(CN.map)<- CRS("+proj=longlat +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +no_defs")
#地图代码
chinamaps <- function (mapObj, a, grayscale=FALSE, lo, hi,...){
require(maps)
if (length(a)!=34){
stop ("wrong number of provinces")
}
getColor <- function(mapdata, provname, provcol, othercol){
FUN <- function(x, y){
ifelse(x %in% y, which(y==x), 0)
}
colIndex <- sapply(mapdata@data$NAME, FUN=FUN, provname)
cols <- c(othercol, provcol)[colIndex+1];
return(cols)
}
a.long <-getColor(mapObj, levels(mapObj$NAME), a, "white")
b.long <- ifelse(is.na(a.long), "black", "gray85")
a.long <- ifelse(is.na(a.long), "white", a.long)
plot(mapObj, col=a.long, border=b.long, lwd=0.5)
}
# parameter 'a' follows the order as below:
levels(CN.map$NAME)
GDP <- rnorm(34, mean=5000, sd=1000) #
cols1 <- heat.colors(34, alpha = 1)
cols2 <- terrain.colors(34, alpha = 1)
chinamaps(mapObj=CN.map, a=cols1[order(GDP)])
chinamaps(mapObj=CN.map, a=rev(cols1)[order(GDP)])#rev条换颜色
chinamaps(mapObj=CN.map, a=cols2[order(GDP)])
然后作业要求的数据如下，要求是绘制GDP地图
set.seed(1)
GDP <- rnorm(34, mean=5000, sd=50)
问题是，我安装了各种包老师给的代码走了一遍，没有报错可是也看不到地图，也照着一些帖子的代码做了，可是都不能画出来。
请教一下问题出在哪里？我该用什么代码画这个作业的地图呢（最好能附上代码）？