Бары должны быть нанесены на карту

Question

Я хочу нарисовать диаграмму на карте, как на этой странице: http://docs.oracle.com/cd/E16162_01/web.1112/e16181/dv_map.htm

Я использовал следующие данные и код, но он не работает.

ddf = read.table(header=T, sep=",", text="
country, val1, val2, val3
Nigeria,5,10,8
China,8,12,20
Brazil,15,9,20
Italy,9,25,5
Australia,15,5,25
")

require (rworldmap)

sPDF <- joinCountryData2Map(ddf
    , joinCode = "NAME"
    , nameJoinColumn = "country"
    , verbose = TRUE)

mapBars( sPDF, 
       , nameZs=c(ddf$val1, ddf$val2, ddf$val3)
       , symbolSize=4 )

Как я могу это исправить? Пожалуйста помоги.

---

Изменить: я попытался использовать код в ответе ниже, но я получаю пустое окно окна с ошибкой: Ошибка в номинальной (plt = c(xyx, xyy), new = TRUE): недопустимое значение, указанное для графического параметра "plt"

Answer 1

Вы можете попробовать использовать сюжеты:

# ...
plot(getMap())
df <- merge(x=ddf, y=sPDF@data[sPDF@data$ADMIN, c("ADMIN", "LON", "LAT")], by.x="country", by.y="ADMIN", all.x=TRUE)
require(TeachingDemos)
for (i in 1:nrow(df)) 
  subplot(barplot(height=as.numeric(as.character(unlist(df[i, 2:4], use.names=F))), 
                  axes=F, 
                  col=rainbow(3), ylim=range(df[,2:4])),
          x=df[i, 'LON'], y=df[i, 'LAT'], size=c(.3, .3))
legend("topright", legend=names(df[, 2:4]), fill=rainbow(3))