为各种stat_ellipse添加图例

Question

我有一个带有不同数据椭圆(stat_ellipse)的图，因此我改变了alpha级别。如何将图例添加到绘图中，以对应于各种椭圆(具有颜色和alpha级别)？

 library(tidyverse)
 library(RColorBrewer)

colors <- brewer.pal(n = 8, name = "Blues")
ggplot(faithful, aes(waiting, eruptions)) +
  geom_point(alpha=0.1, color='navy') +
   stat_ellipse(level=0.1, color=colors[1], size=1.1) +
   stat_ellipse(level=0.2, color=colors[2], size=1.1) +
   stat_ellipse(level=0.3, color=colors[3], size=1.1) +
   stat_ellipse(level=0.4, color=colors[4], size=1.1) +
   stat_ellipse(level=0.5, color=colors[5], size=1.1) +
   stat_ellipse(level=0.6, color=colors[6], size=1.1) +
   stat_ellipse(level=0.7, color=colors[7], size=1.1) +
   stat_ellipse(level=0.8, color=colors[8], size=1.1)  +
   stat_ellipse(level=0.9, color='navy', size=1.1)

Answer 1

我不知道stat_ellipse是否可以做到这一点。这是一个解决方案。它包括使用car包计算椭圆，然后使用geom_path。

library(car)
levels <- seq(0.1, 0.9, by = 0.2)
ell <- dataEllipse(faithful$waiting, faithful$eruptions, draw = FALSE, 
                   levels = levels, segments = 200)

library(purrr)
paths_list <- imap(ell, function(mat, level){
  cbind(mat, level=as.numeric(level))
})
paths <- as.data.frame(do.call(rbind, paths_list))
paths$level <- as.factor(paths$level)

ggplot(faithful, aes(waiting, eruptions)) +
  geom_point(alpha=0.5, color='navy') +
  geom_path(aes(x=x, y=y, group=level, color=level), data=paths, size = 2)

​

​