按：由于maptools更新比较大，因此拙作“基于R语言的中国柑橘主要病虫害空间分布图”（参见：http://sciencenet.cn/m/user_content.aspx?id=217688）中的R代码已经失效，这里面提供新的代码：

#新代码：

require(maptools)

readShapeSpatial("pests") -> pestmap
summary(pestmap)

br1 <- c(0.002,0.031,0.068,0.101)
rm <- cut(pestmap$潜叶蛾,breaks=br1,labels=c(1:3))

cols <- c("gray","blue","red")
fgs <- cols[rm]

plot(pestmap,axes=T,col=fgs,main="柑橘潜叶蛾为害率", xlab="经度", ylab="纬度")

legtext <- c("0.003-0.031","0.031-0.068","0.068-0.101")
legend("topright",bg="white", fill=cols,legend=legtext)
#thigmophobe.labels(dat$jd,dat$wd,dat$lb)

dat = read.csv(con <- textConnection("lb,jd,wd
   重 庆 市,106.3,29.6
   浙 江 省,119.2,29
   湖 南 省,111.3,28.1
   江 西 省,115.9,28.2
   湖 北 省,112,31
   四 川 省,103,30.65
   贵 州 省,106.5,27.1
   福 建 省,117.8,26.5
   广 东 省,113.2,23.5
   广西壮族自治区,108.7,24.3
   云 南 省,101.2,24.8
   陕 西 省,108.3,34.1"))
text(dat$jd,dat$wd,dat$lb,col="black")
close(con)

######################################################

#旧代码：

require(maptools)

read.shape("pests") -> pestmap
attributes(pestmap$att.data)

br1 <- c(0.002,0.031,0.068,0.101)
rm <- cut(pestmap$att.data$潜叶蛾,breaks=br1,labels=c(1:3))

cols <- c("gray","blue","red")
fgs <- cols[rm]

plot.Map(pestmap,fg=fgs,ol="black",main="柑橘潜叶蛾为害率", xlab="经度", ylab="纬度")

legtext <- c("0.003-0.031","0.031-0.068","0.068-0.101")
legend("topright",bg="white", fill=cols,legend=legtext)
#thigmophobe.labels(dat$jd,dat$wd,dat$lb)

dat = read.csv(con <- textConnection("lb,jd,wd
   重 庆 市,106.3,29.6
   浙 江 省,119.2,29
   湖 南 省,111.3,28.1
   江 西 省,115.9,28.2
   湖 北 省,112,31
   四 川 省,103,30.65
   贵 州 省,106.5,27.1
   福 建 省,117.8,26.5
   广 东 省,113.2,23.5
   广西壮族自治区,108.7,24.3
   云 南 省,101.2,24.8
   陕 西 省,108.3,34.1"))
text(dat$jd,dat$wd,dat$lb,col="white")
close(con)