ggsubplot (ggplot2) 中的饼图答案

【问题标题】：Pie Charts in ggsubplot (ggplot2)ggsubplot (ggplot2) 中的饼图
【发布时间】：2013-12-01 20:48:03
【问题描述】：

另一个有趣的 ggplot2 挑战！我正在尝试创建一个 ggsubplot 调用，该调用在世界地图上提供饼图。困难的事情似乎是 ggplot2 中的饼图是带有极坐标的堆叠条形图，并且 coord_polar 的添加不仅影响子图几何，而且影响整个地图本身。有谁知道如何仅将某个坐标方案应用于子图调用的一部分？这是我目前所拥有的：

library(ggplot2)
library(maps)
library(mapproj)
#install.packages("devtools")
library(devtools)
# install ggsubplot from github repo, not currently on CRAN
install_github(username="garrettgman", repo="ggsubplot")

world = map_data("world")

loc_pie = structure(list(Region = structure(c(3L, 5L, 7L, 8L, 9L, 10L, 
11L, 12L, 13L, 15L, 16L, 2L, 14L, 2L, 4L, 5L, 6L, 7L, 9L, 10L
), .Label = c("", "ANT/SO", "ARC", "EPR/GAL", "GOM/CAR", "IND", 
"MAR", "MED", "N-ATL", "NE-ATL", "NE-PAC", "NW-ATL", "NW-PAC", 
"SE-ATL", "SE-PAC", "SW-ATL", "SW-PAC"), class = "factor"), Group3 = structure(c(1L, 
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 3L, 3L, 3L, 3L, 
3L, 3L, 3L), .Label = c("Annelida", "Choanoflagellata", "Chordata", 
"Cnidaria", "Crustacea", "Echinodermata", "Foraminifera", "Mollusca", 
"Nematoda", "Other", "Platyhelminthes", "Porifera"), class = "factor"), 
ones = c(1, 1, 5, 1, 1, 1, 18, 3, 1, 4, 8, 1, 1, 2, 1, 1, 
6, 1, 2, 5), tot = c(5, 30, 11, 16, 28, 22, 51, 25, 78, 13, 
32, 57, 61, 57, 15, 30, 20, 11, 28, 22), div = c(0.2, 0.0333333333333333, 
0.454545454545455, 0.0625, 0.0357142857142857, 0.0454545454545455, 
0.352941176470588, 0.12, 0.0128205128205128, 0.307692307692308, 
0.25, 0.0175438596491228, 0.0163934426229508, 0.0350877192982456, 
0.0666666666666667, 0.0333333333333333, 0.3, 0.0909090909090909, 
0.0714285714285714, 0.227272727272727), lat = c(71.4493167, 
19.9897167, 23.5874333, 37.6802167, 55.13365, 36.6889333, 
35.9565333, 35.53935, 30.4266, -30.32195, -33.2038, -65.8756333, 
-17.12415, -65.8756333, 0.1135, 19.9897167, -14.5800667, 
23.5874333, 55.13365, 36.6889333), long = c(-1.0550667, -81.3430667, 
-41.2278667, 15.9298833, -30.4984333, -17.4906167, -149.4363333, 
-63.01795, 156.3570833, -110.23255, -31.20155, -25.4557, 
0.0881833, -25.4557, -101.07455, -81.3430667, 77.4312667, 
-41.2278667, -30.4984333, -17.4906167)), .Names = c("Region", 
"Group3", "ones", "tot", "div", "lat", "long"), row.names = c(NA, 
20L), class = "data.frame")

ggplot(data=loc_pie) + geom_polygon(data=world, aes(x=long, y=lat, group =group),colour="grey40", fill="grey40") + geom_subplot(height=12, aes(long, lat, group=Region, subplot = (geom_bar(aes(x = factor(1),   y=div, fill=factor(Group3)), width =1, height = 2,stat="identity"))))

看起来我不能发布图片，因为我是新用户。

![带有堆叠条形图子图的世界地图][1]

关于如何仅将极坐标合并到子图的任何想法？非常感谢。

【问题讨论】：

嗨！拥有minimal reproducible example 会非常有帮助，这样我们就可以自己重现问题。
是的，感谢您的建议。进行了应该允许更基本复制的编辑。
ggsubplot 不在 CRAN 上。请帮助人们通过使您的代码可重现来帮助您，即添加运行平稳的代码，包括用于安装必要软件包的代码。谢谢。
对，对不起。似乎它已经打开和关闭 CRAN，并且可能很快就会恢复。现在，看起来您可以从 cran 获取以前的版本，或者使用我在上面添加的代码通过 github 安装它。谢谢！
这确实解释了如何在基本图形包中执行此操作，但现在可能在 ggplot2 中添加了 ggsubplot。感谢您的链接！

标签： r map plot ggplot2 subplot

【解决方案1】：

尝试多图功能 http://www.cookbook-r.com/Graphs/Multiple_graphs_on_one_page_(ggplot2)/ i 允许您在一个窗口上绘制多个图

或

ggplot, facet, piechart: placing text in the middle of pie chart slices

这就是您要搜索的内容

【讨论】：

【解决方案2】：

这里是ggplot2和ggtree最新版本的解决方案：

library(maps)
library(mapproj)
library(ggplot2)
library(ggtree)
library(dplyr)
library(magrittr)

world = map_data("world")

p <- ggplot(data=world, aes(x=long, y=lat, group =group)) + geom_polygon(colour="grey40", fill="grey40")

for (name in unique(loc_pie[['Region']])) {
  loc_region <- filter(loc_pie, Region %in% name)
  pie_tmp <- ggplot(data = loc_region, aes(x = factor(1), y = div, fill = Group3)) +
    geom_bar(width = 1, stat = "identity") + coord_polar(theta = "y") +
    scale_fill_discrete(drop = FALSE) +
    xlab(NULL) + ylab(NULL) + theme_tree() +
    theme_transparent()
  lat_region <- loc_region[[1,'lat']]
  long_region <- loc_region[[1,'long']]
  p %<>% subview(pie_tmp, long_region, lat_region, width = .07, height = .07)
}

p

宽度和高度是手动选择的，它们对应于子视图相对于完整视图的比例。

为了完整起见，loc_pie 可以通过以下方式获得：

loc_pie = structure(
  list(
    Region = structure(
      c(
        3L,
        5L,
        7L,
        8L,
        9L,
        10L,
        11L,
        12L,
        13L,
        15L,
        16L,
        2L,
        14L,
        2L,
        4L,
        5L,
        6L,
        7L,
        9L,
        10L
      ),
      .Label = c(
        "",
        "ANT/SO",
        "ARC",
        "EPR/GAL",
        "GOM/CAR",
        "IND",
        "MAR",
        "MED",
        "N-ATL",
        "NE-ATL",
        "NE-PAC",
        "NW-ATL",
        "NW-PAC",
        "SE-ATL",
        "SE-PAC",
        "SW-ATL",
        "SW-PAC"
      ),
      class = "factor"
    ),
    Group3 = structure(
      c(
        1L,
        1L,
        1L,
        1L,
        1L,
        1L,
        1L,
        1L,
        1L,
        1L,
        1L,
        2L,
        2L,
        3L,
        3L,
        3L,
        3L,
        3L,
        3L,
        3L
      ),
      .Label = c(
        "Annelida",
        "Choanoflagellata",
        "Chordata",
        "Cnidaria",
        "Crustacea",
        "Echinodermata",
        "Foraminifera",
        "Mollusca",
        "Nematoda",
        "Other",
        "Platyhelminthes",
        "Porifera"
      ),
      class = "factor"
    ),
    ones = c(1, 1, 5, 1, 1, 1, 18, 3, 1, 4, 8, 1, 1, 2, 1, 1,
             6, 1, 2, 5),
    tot = c(5, 30, 11, 16, 28, 22, 51, 25, 78, 13,
            32, 57, 61, 57, 15, 30, 20, 11, 28, 22),
    div = c(
      0.2,
      0.0333333333333333,
      0.454545454545455,
      0.0625,
      0.0357142857142857,
      0.0454545454545455,
      0.352941176470588,
      0.12,
      0.0128205128205128,
      0.307692307692308,
      0.25,
      0.0175438596491228,
      0.0163934426229508,
      0.0350877192982456,
      0.0666666666666667,
      0.0333333333333333,
      0.3,
      0.0909090909090909,
      0.0714285714285714,
      0.227272727272727
    ),
    lat = c(
      71.4493167,
      19.9897167,
      23.5874333,
      37.6802167,
      55.13365,
      36.6889333,
      35.9565333,
      35.53935,
      30.4266,
      -30.32195,
      -33.2038,
      -65.8756333,-17.12415,
      -65.8756333,
      0.1135,
      19.9897167,
      -14.5800667,
      23.5874333,
      55.13365,
      36.6889333
    ),
    long = c(
      -1.0550667,
      -81.3430667,-41.2278667,
      15.9298833,
      -30.4984333,
      -17.4906167,
      -149.4363333,-63.01795,
      156.3570833,
      -110.23255,
      -31.20155,
      -25.4557,
      0.0881833,
      -25.4557,
      -101.07455,
      -81.3430667,
      77.4312667,-41.2278667,
      -30.4984333,
      -17.4906167
    )
  ),
  .Names = c("Region",
             "Group3", "ones", "tot", "div", "lat", "long"),
  row.names = c(NA,
                20L),
  class = "data.frame"
)

【讨论】：

【解决方案3】：

以下类型的作品，但不使用ggsubplot。

library(sp)
library(plotGoogleMaps)

loc_pie$divc = 1-loc_pie$div
coordinates(loc_pie) <- ~long+lat
proj4string(loc_pie) <- CRS("+proj=longlat")
m <- segmentGoogleMaps(loc_pie,zcol=c('div','divc'), 
                       scalelist=FALSE,
                       max.radius=1000000,
                       mapTypeId = "ROADMAP",
                       legend = FALSE, control = FALSE)

我无法弄清楚如何获得该功能以使所有气泡大小相同。它似乎忽略了一些参数，如 scalelist 和 bounds。 Plotly 也可以，只要您愿意将 div 的单个值编码为颜色或大小：

library(plotly)
loc_pie$hover <- paste(loc_pie$div*100,"%")
g <- list(
  scope = 'world',
  projection = list(type = "mercator"),
  showland = TRUE,
  landcolor = toRGB("gray85"),
  subunitwidth = 1,
  countrywidth = 1,
  subunitcolor = toRGB("white"),
  countrycolor = toRGB("white")
)
plot_ly(loc_pie, lon = long, lat = lat, text = hover,
        marker = list(size = 100),
        type = 'scattergeo', color=div,
        locationmode = 'country names',geo=g)

【讨论】：