使用R ggplot2绘制单个数据图，并添加2条水平线和正确的图例。

Question

在Axeman的帮助下，我已经能够绘制一个具有平均数/中位数hlines的数据。为了完成注释，我想把平均数/中位数值放入图例标签中。我的代码：

require(data.table)
require(ggplot2)
wf <- data.table::fread('Date Gals Days  GpD    GpM
                     2016-10-21  6.0    1  6.0  186.0
                     2016-10-22  6.0    1  6.0  186.0
                     2016-10-23 12.4    1 12.4  384.4
                     2016-10-24 26.8    1 26.8  830.8
                     2016-10-25 33.3    1 33.3 1032.3
                     2016-10-26 28.3    1 28.3  877.3')
nmu <- mean(wf$Gals)
textmu <- paste("Mean:", round(nmu, 2))
nmed <- median(wf$Gals)
textmed <- paste("Median:", round(nmed, 2))

p <- ggplot(wf, aes(Date, Gals)) +  
     geom_line(aes(lty = "Gals", color = "Gals", group = 1)) +
     geom_hline(aes(yintercept = nmu,  linetype = textmu, color = textmu), 
                color='red',size=1)+
     geom_hline(aes(yintercept = nmed, linetype = textmed, color = textmed), 
                color='orange',size=1)+
     #scale_linetype('Water Usage') +
     scale_color_manual('Water Usage', 
                        values=c('Gals'='black',textmu='red', textmed='orange')) +
     xlab("") + 
     ylab("Gallons per Day")

print(p)

产生的结果： 我该如何纠正这个问题？

- user50506

Answer 1

准备一个标签向量，然后使用我在上一个答案中给出的代码：

labels <- c(Gals = 'Gals',
            Mean = paste("Mean:", round(mean(wf$Gals), 2)),
            Median = paste("Median:", round(median(wf$Gals), 2)))

ggplot(wf, aes(Date, Gals)) +  
  geom_line(aes(lty = "Gals", color = "Gals", group = 1)) +
  geom_hline(aes(yintercept = mean(wf$Gals), linetype = "Mean", color = "Mean"), size=1) +
  geom_hline(aes(yintercept = median(wf$Gals), linetype = "Median", color = "Median"), size = 1) +
  scale_linetype('Water Usage', labels = labels) +
  scale_color_manual(
    'Water Usage', 
    values = c('Gals' = 'black', 'Mean' = 'red', 'Median' = 'orange'),
    labels = labels
  ) +
  xlab("") + ylab("Gallons per Day")

- Axeman

---

很好的结果。这正是我想要做的，尽可能保持数据框的简单性，并学习更多ggplot2的知识。感谢Axeman和社区！ - user50506

Answer 2

通过融合你的 data.table，你可以：

require(data.table)
require(ggplot2)
wf <- data.table::fread('Date Gals Days  GpD    GpM
                    2016-10-21  6.0    1  6.0  186.0
                    2016-10-22  6.0    1  6.0  186.0
                    2016-10-23 12.4    1 12.4  384.4
                    2016-10-24 26.8    1 26.8  830.8
                    2016-10-25 33.3    1 33.3 1032.3
                    2016-10-26 28.3    1 28.3  877.3')

wf[,`:=`(Mean=round(mean(Gals),2),Median=round(median(Gals),2))]

p <- ggplot(melt(wf,id.vars = c(1,3))
[variable%in%c("Gals","Mean","Median")], aes(Date, 
value,color=variable,group=variable)) +  
geom_line() +
scale_color_manual('Water Usage',values=c('Gals'='black',Mean='red', 
Median='orange')) +
xlab("") + 
ylab("Gallons per Day")

print(p)

- Erdem Akkas

---

这不会将平均值和中位数放在图例标签中，与问题描述不符。 - Axeman