-2

df

Date     Score  Team
1/1/2011   3    A
1/2/2011   5    A
1/3/2011   15   A
1/4/2011   39   B
1/5/2011   23   B
1/6/2011  100   B
1/7/2011   4    C
1/8/2011  25    C
1/9/2011   30   C

library(ggplot2)
ggplot(df, aes(Date, Score, group=1)) + geom_point() + geom_smooth(method="loess", se=T, size=1) + facet_wrap(~Team)

我希望能够在 method="loess" 达到最大值的地方画一条水平线。这里有没有人对我将如何做到这一点有任何意见？

4

1 回答 1

8

用 构建绘图ggplot_build()，然后提取数据：

p <- ggplot(df, aes(Date, Score, group=1)) + geom_point() + 
  geom_smooth(method="loess", se=T, size=1)

pp <- ggplot_build(p)
p <- p + geom_hline(yintercept = max(pp$data[[2]]$y), col="red")
p

在此处输入图像描述

结果ggplot_build是一个列表。第二个元素data是您正在寻找的内容：

str(pp$data)
List of 2
 $ :'data.frame':   9 obs. of  4 variables:
  ..$ x    : int [1:9] 1 2 3 4 5 6 7 8 9
  ..$ y    : num [1:9] 3 5 15 39 23 100 4 25 30
  ..$ group: int [1:9] 1 1 1 1 1 1 1 1 1
  ..$ PANEL: int [1:9] 1 1 1 1 1 1 1 1 1
 $ :'data.frame':   9 obs. of  7 variables:
  ..$ x    : int [1:9] 1 2 3 4 5 6 7 8 9
  ..$ y    : num [1:9] 1.29 8.67 19.64 24.54 53.86 ...
  ..$ ymin : num [1:9] -103.1 -60.2 -57.3 -52.4 -23 ...
  ..$ ymax : num [1:9] 105.7 77.6 96.5 101.4 130.8 ...
  ..$ se   : num [1:9] 34.8 23 25.6 25.6 25.6 ...
  ..$ group: int [1:9] 1 1 1 1 1 1 1 1 1
  ..$ PANEL: int [1:9] 1 1 1 1 1 1 1 1 1

请注意，此对象是一个列表，其中包含与每个 geom 对应的数据框。在我的解决方案中，我只是提取最大值y并通过它绘制 a hline。

这也适用于多个方面，但您必须做更多的工作来提取数据：

p <- ggplot(df, aes(Date, Score, group=1)) + geom_point() + 
  geom_smooth(method="loess", se=T, size=1) + facet_wrap(~Team)

pp <- ggplot_build(p)

library(plyr)
hdat <- ddply(pp$data[[2]], .(PANEL), summarize, hline=max(y))
hdat$Team <- unique(df$Team)[as.numeric(hdat$PANEL)]

p <- p + geom_hline(data=hdat, aes(yintercept = hline), col="red")
p

在此处输入图像描述

于 2012-12-04T16:17:55.043 回答