您好,我只想在此地图上添加一个图例,代表法国各部门年轻女性的生育力:
使用此代码:
library(raster)
fr <- getData(country='France',level=2)
classes <- c('0-1','1-1.5','1.5-2','2-3')
fr$fertility <- sample(seq(1,3,0.1),length(fr),replace = T)
fr$fertgroups[fr$fertility <=1] <-1
fr$fertgroups[fr$fertility > 1 & fr$fertility <=1.5] <-2
fr$fertgroups[fr$fertility > 1.5 & fr$fertility <=2] <-3
fr$fertgroups[fr$fertility > 2 & fr$fertility <=3] <- 4
fr$cols <- c('red','orange','yellow','green')[fr$fertgroups]
plot(fr,col=fr$cols)
我不知道如何添加这个传奇。