4

我希望能够在与文档数量相关的图表上创建一个带有圆圈的世界地图。

这是我尝试使用 mapBubbles 来完成的，但是在绘制图表时出现了问题。

我拥有的数据文件如下：

ISO3V10   Country   No of Documents    Lat  Lon
ARG     Argentina   41          -64 -34
AUS     Australia   224         133 -27
CAN     Canada      426         -95 60
IRL     Ireland 68           -8 53
ITA     Italy             583           12.8333 42.8333
NLD     Netherlands 327          5.75   52.5
NZL     New Zealand 26           174    -41
ESP     Spain             325            -4  40
GBR     United Kingdom  2849             -2 54
USA     United States   3162            -97 38

我写的代码如下：

# Thanks to Andrie for the reproducible code and Paolo for the suggested edit

zz <-"ISO3V10   Country No.of.Documents Lat  Lon
ARG Argentina   41  -64 -34
AUS Australia   224 133 -27
CAN Canada  426 -95 60
IRL Ireland 68  -8  53
ITA Italy   583 12.8333 42.8333
NLD Netherlands 327 5.75    52.5
NZL 'New Zealand' 26  174 -41
ESP Spain   325 -4  40
GBR 'United Kingdom'  2849    -2  54
USA 'United States'   3162    -97 38
"

dF2 <- read.table(textConnection(zz), header = TRUE)

# dF2 <- read.delim(file="C:\\Users\\js207963\\Documents\\noofpublications_AllUpdated.txt", header = TRUE, sep = "\t")
 dF2[]
 par(mai=c(0,0,0.2,0),xaxs="i",yaxs="i")
 mapBubbles(dF2=getMap(), nameZSize="No.of.Documents", nameZColour="Country",oceanCol="lightblue", landCol="wheat", addLegend=FALSE

所以问题是你能帮我修复代码以正确绘制数据吗？

干杯，杰斯

4

2 回答 2

6

您可以使用以下方法绘制它ggplot：

在此处输入图像描述

重新创建您的数据：

dat <- read.table(text="
ISO3V10   Country   'No of Documents'    Lat  Lon
ARG     Argentina   41          -64 -34
AUS     Australia   224         133 -27
CAN     Canada      426         -95 60
IRL     Ireland 68           -8 53
ITA     Italy             583           12.8333 42.8333
NLD     Netherlands 327          5.75   52.5
NZL     'New Zealand' 26           174    -41
ESP     Spain             325            -4  40
GBR     'United Kingdom'  2849             -2 54
USA     'United States'   3162            -97 38
", header=TRUE)

加载包和情节：

library(ggplot2)
library(maps)

mdat <- map_data('world')

str(mdat)
ggplot() + 
  geom_polygon(dat=mdat, aes(long, lat, group=group), fill="grey50") +
  geom_point(data=dat, 
             aes(x=Lat, y=Lon, map_id=Country, size=`No.of.Documents`), col="red")
于 2012-05-25T09:25:05.737 回答
1

您需要向函数传递更多信息mapBubbles，因为您的数据框不是 SpatialPolygonsDataFrame。以下应该有效（您现在Lat并且Lon可能被贴错标签）：

 mapBubbles(dF=dF2, nameZSize="No.of.Documents",
  nameZColour="Country",oceanCol="lightblue", landCol="wheat",
  addLegend=FALSE, nameX = "Lat", nameY = "Lon")

上面，nameX并被nameY传递给函数以指示在哪里绘制气泡。此外，将数据框传递dF2给参数dF，而不是调用getMap().

编辑：

我也喜欢 Andrie 的回答ggplot2（除其他外，空间比下图中的空间使用得更好），但是由于您专门使用 发布了其他问题rworldmap，我认为坚持使用该软件包是合适的。

带气泡的世界地图

于 2012-05-25T09:13:44.597 回答