在R中使用ggplot2画散点图，并在每个数据点上添加椭圆形，以显示每个数据点的分布情况。

Question

我的数据集由4列组成，如下所示： 左侧的两列表示地理结构的坐标XY，右侧的两列表示“每个”地理单元的大小（南北直径和东西直径）。

我想绘制一个散点图，其中绘制所有坐标并在每个点上绘制一个包括每个地理单元直径的椭圆。

手动地，仅使用两个点，图像应该像这样： 如何使用ggplot2实现它？

您可以在此处下载数据。

- antecessor

---

如果您提供了一个完整的最小可复制示例，很可能会得到一个好的和及时的回答。这个示例将成为我们解决问题的基础，并且能够帮助我们向您展示如何回答您的问题。请参考这里获取示例编写指南。 - Eric Fail

希望我的手动图像有所帮助！ - antecessor

它确实不同。这很棒！你看了我提供的文章吗？接下来要做的是制作一些可用的数据，即 df <- data.frame(X = c(-4275, -4855, ...), Y = (1145, 1330, ...), ....) 再次强调，所有细节都在我提供的文章中详细解释了。 - Eric Fail

如果有人想要帮助你，他们需要你的数据。从图像中获取数据很困难。请发布有效的R代码以生成数据，以便可以复制/粘贴。 - Gregor Thomas

我添加了一个下载数据的链接。 - antecessor

Answer 1

使用ggforce中的geom_ellipse()：

library(ggplot2)
library(ggforce)

d <- data.frame(
  x = c(10, 20),
  y = c(10, 20),
  ns = c(5, 8),
  ew = c(4, 4)
)

ggplot(d, aes(x0 = x, y0 = y, a = ew/2, b = ns/2, angle = 0)) + 
  geom_ellipse() +
  coord_fixed()

^{2019年6月1日创建，使用reprex包（v0.2.1）}

- Claus Wilke

Answer 2

我不会对Claus Wilke已经发布的内容进行任何新代码的添加。所有的荣誉应该归于Claus。我只是用实际数据测试它，并向OP展示如何发布数据。

加载所需包。

# install.packages(c("tidyverse"), dependencies = TRUE)
library(tidyverse)

读取数据，

tbl <- read.table(
text = "
X Y Diameter_N_S Diameter_E_W
-4275   1145    77  96
-4855   1330    30  25
-4850   1612    45  90
-4990   1410    15  15
-5055   1230    60  50
-5065   1503    43  45
-5135   1305    40  50
-5505   1190    55  70
-5705   1430    90  40
-5645   1535    52  60
", header = TRUE, stringsAsFactors = FALSE) %>% as_tibble()

展示数据。

tbl
#> # A tibble: 10 x 4
#>        X     Y Diameter_N_S Diameter_E_W
#>    <int> <int>        <int>        <int>
#>  1 -4275  1145           77           96
#>  2 -4855  1330           30           25
#>  3 -4850  1612           45           90
#>  4 -4990  1410           15           15
#>  5 -5055  1230           60           50
#>  6 -5065  1503           43           45
#>  7 -5135  1305           40           50
#>  8 -5505  1190           55           70
#>  9 -5705  1430           90           40
#> 10 -5645  1535           52           60

正在加载所需的更多软件包

library(ggforce) # devtools::install_github("thomasp85/ggforce")

执行中

ggplot(tbl, aes(x0 = X, y0 = Y, a = Diameter_E_W, b = Diameter_N_S, angle = 0)) +
        geom_ellipsis() + geom_point(aes(X, Y), size = .5) + coord_fixed() + theme_bw()

- Eric Fail