如何在ggplot2中添加均值和中位数的图例?

时间:2019-03-18 08:56:12

标签: r ggplot2

我正在绘制单个轨迹,并且还添加了均值和中位数的线。我不确定如何添加图例来标识这两行,非常感谢您提供使用ggplot2进行操作的建议。

我的数据:

> dput(junk)
structure(list(Patient = c("TLA001", "TLA001", "TLA001", "TLA001", 
"TLA001", "TLA001", "TLA002", "TLA002", "TLA002", "TLA002", "TLA002", 
"TLA003", "TLA003", "TLA003", "TLA003", "TLA003", "TLA003", "TLA004", 
"TLA004", "TLA004", "TLA004", "TLA004", "TLA004", "TLA005", "TLA005", 
"TLA005", "TLA005", "TLA005", "TLA005", "TLA006", "TLA006", "TLA006", 
"TLA006", "TLA006", "TLA006", "TLA006", "TLA008", "TLA008", "TLA008", 
"TLA008", "TLA008", "TLA008", "TLA009", "TLA009", "TLA009", "TLA009", 
"TLA009", "TLA009", "TLA010", "TLA010", "TLA010", "TLA010", "TLA010", 
"TLA010", "TLA011", "TLA011", "TLA011", "TLA011", "TLA011", "TLA011", 
"TLA012", "TLA012", "TLA012", "TLA012", "TLA012", "TLA012", "TLA012", 
"TLA013", "TLA013", "TLA013", "TLA013", "TLA013", "TLA013", "TLA013", 
"TLA014", "TLA014", "TLA014", "TLA014", "TLA015", "TLA015", "TLA015", 
"TLA015", "TLA015", "TLA015"), Night = c("TLA -3", "TLA -2", 
"TLA -1", "TLA 1", "TLA 2", "TLA 3", "TLA -3", "TLA -2", "TLA -1", 
"TLA 1", "TLA 2", "TLA -3", "TLA -2", "TLA -1", "TLA 1", "TLA 2", 
"TLA 3", "TLA -3", "TLA -2", "TLA -1", "TLA 1", "TLA 2", "TLA 3", 
"TLA -3", "TLA -2", "TLA -1", "TLA 1", "TLA 2", "TLA 3", "TLA -3", 
"TLA -2", "TLA -1", "TLA 1", "TLA 2", "TLA 3", "TLA 4", "TLA -3", 
"TLA -2", "TLA -1", "TLA 1", "TLA 2", "TLA 3", "TLA -3", "TLA -2", 
"TLA -1", "TLA 1", "TLA 2", "TLA 3", "TLA -3", "TLA -2", "TLA -1", 
"TLA 1", "TLA 2", "TLA 3", "TLA -3", "TLA -2", "TLA -1", "TLA 1", 
"TLA 2", "TLA 3", "TLA -3", "TLA -2", "TLA -1", "TLA 1", "TLA 2", 
"TLA 3", "TLA 4", "TLA -3", "TLA -2", "TLA -1", "TLA 1", "TLA 2", 
"TLA 3", "TLA 4", "TLA -3", "TLA -2", "TLA -1", "TLA 1", "TLA -3", 
"TLA -2", "TLA -1", "TLA 1", "TLA 2", "TLA 3"), var = c(0.1, 
0, 0, 0, 0, 0, 0.5, 0, 0, 2.8, 0.2, 0.3, 0, 0.4, 0.3, 0, 2.2, 
1.1, 0.6, 2.1, 1.4, 1.7, 0.9, 0.3, 0.3, 0, 1, 0, 0.3, 0, 0, 0, 
0, 0, 0, 0, 0.1, 0.2, 0.1, 0.1, 0.4, 0, 0.4, 0.1, 0, 0, 0, 0, 
1.1, 0.8, 0, 0, 0.6, 0, 0, 0, 0, 0, 0, 0, 2, 0.9, 0.2, 0.4, 0.3, 
1, 0.9, 0, 0, 0, 0.3, 0.3, 0, 0, 1.2, 0.9, 0.6, 0.7, 0.7, 0, 
0, 0, 0, 0), SF_Cluster = c("Responder", "Responder", "Responder", 
"Responder", "Responder", "Responder", "Non-responder", "Non-responder", 
"Non-responder", "Non-responder", "Non-responder", "Responder", 
"Responder", "Responder", "Responder", "Responder", "Responder", 
"Responder", "Responder", "Responder", "Responder", "Responder", 
"Responder", "Responder", "Responder", "Responder", "Responder", 
"Responder", "Responder", "Non-responder", "Non-responder", "Non-responder", 
"Non-responder", "Non-responder", "Non-responder", "Non-responder", 
"Responder", "Responder", "Responder", "Responder", "Responder", 
"Responder", "Non-responder", "Non-responder", "Non-responder", 
"Non-responder", "Non-responder", "Non-responder", "Responder", 
"Responder", "Responder", "Responder", "Responder", "Responder", 
"Responder", "Responder", "Responder", "Responder", "Responder", 
"Responder", "Non-responder", "Non-responder", "Non-responder", 
"Non-responder", "Non-responder", "Non-responder", "Non-responder", 
"Responder", "Responder", "Responder", "Responder", "Responder", 
"Responder", "Responder", "Non-responder", "Non-responder", "Non-responder", 
"Non-responder", "Responder", "Responder", "Responder", "Responder", 
"Responder", "Responder")), row.names = c(NA, -84L), class = c("tbl_df", 
"tbl", "data.frame"), spec = structure(list(cols = list(Patient = structure(list(), class = c("collector_character", 
"collector")), Night = structure(list(), class = c("collector_character", 
"collector")), var = structure(list(), class = c("collector_double", 
"collector")), SF_Cluster = structure(list(), class = c("collector_character", 
"collector"))), default = structure(list(), class = c("collector_guess", 
"collector"))), class = "col_spec"))

剧情代码:

ggplot(data = junk, aes(x = Night, y=var, group=Patient, col=SF_Cluster)) + 
  geom_line(alpha=0.5) + 
  geom_smooth(size=0.2, se=F)+
  ggtitle("Number of vardput(junk)") + 
  xlab("Night")+
  ylab("N") + 
  geom_vline(xintercept = 4, linetype = "dotted", color = "blue", size = 1.5) + 
  stat_summary(fun.y = median, geom = "line", lwd = 1, aes(group = 1)) + 
  stat_summary(fun.y = mean, geom = "line", lwd = 1, aes(group = 1), 
               col = "red", linetype = "dotted", size = 2.5) 

0 个答案:

没有答案