でプロットを作成し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
はリストです。2番目の要素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
プロットし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