ggplot2 绘制单个椭圆，但按组绘制颜色

Question

我的数据由多个个人的多个数据点组成，每个人都来自特定的研究地点。我想绘制所有点，为每个人绘制 95% 的椭圆，然后按研究地点为椭圆着色。不幸的是，当我指定按站点着色时，似乎为聚合组绘制了椭圆。

数据如下：

dat1 <- data.frame(X=rnorm(21),Y=rnorm(21),indiv_id=rep(c(1,2,3),7),group_id=rep(1,21))
dat2 <- data.frame(X=rnorm(21,5),Y=rnorm(21,5),indiv_id=rep(c(4,5,6),7),group_id=rep(2,21))
dat3 <- data.frame(X=rnorm(21,10),Y=rnorm(21,10),indiv_id=rep(c(7,8,9),7),group_id=rep(3,21))
ggdat <- rbind(dat1,dat2,dat3)
ggdat$indiv_id <- as.factor(ggdat$indiv_id)
ggdat$group_id <- as.factor(ggdat$group_id)

如果我逐个绘制椭圆，我可以单独看到所有的椭圆：

ggplot(ggdat) +
  geom_point(aes(x=X, y=Y,color=indiv_id),size=1) + # 
  stat_ellipse(aes(x=X, y=Y,color=indiv_id),type = "norm")

但如果我按组绘制，则每组只画一个椭圆：

ggplot(ggdat) +
  geom_point(aes(x=X, y=Y,color=indiv_id),size=1) + # 
  stat_ellipse(aes(x=X, y=Y,color=group_id),type = "norm") + #, linetype = 2
  theme(legend.position='none')

如何绘制所有 9 个椭圆并按组着色？感谢您的帮助！

Answer 1

明确定义组：

ggplot(ggdat) +
  geom_point(aes(x=X, y=Y,color=indiv_id),size=1) + # 
  stat_ellipse(aes(x=X, y=Y,color=group_id, group=indiv_id),type = "norm") +
  theme(legend.position='none')