ggplot2中的中值回归线

Question

我正在尝试使用 ggplot 来拟合中值回归线。这使用了 quantreg 包的 rq() 函数。数据是来自 HSAUR3 包的云数据。这是我的代码：

seeding.no <- predict(rq(rainfall ~ sne, data = clouds, subset = seeding == "no"))
seeding.yes <- predict(rq(rainfall ~ sne, data = clouds, subset = seeding == "yes"))
clouds.predcit <- data.frame(Rainfall=clouds$rainfall, No=seeding.no, Yes=seeding.yes)

ggplot(clouds, aes(x=sne, y=rainfall, colour=clouds$seeding, shape=clouds$seeding)) + 
  geom_point(size=2) +
  labs(title="Rainfall vs S-Ne Criterion", x="S-Ne Criterion", y="Rainfall") +
  scale_color_manual(values=c("blue", "red")) + 
  scale_shape_manual(values=c(1,2)) +
  geom_line(color="blue",data=clouds.predcit, aes(x=No, y=Rainfall)) +
  geom_line(color="red",data=clouds.predcit, aes(x=Yes, y=Rainfall)) +
  theme_minimal()

这是我得到的：

这就是我想要得到的：

我的代码有什么问题？

Answer 1

我用过：

geom_quantile(quantiles = 0.5)

并让它工作。毕竟不需要预测。

Answer 2

我认为你使用了错误的 geom，使用 geom_smooth 会更幸运。这有很多选项，但您可能需要 +geom_smooth(method="lm")。即使安装了 HSAUR3 包，我也无法让您的代码按原样运行。

试试这个

seeding.no <- predict(rq(rainfall ~ sne, data = clouds, subset = seeding == "no"))
seeding.yes <- predict(rq(rainfall ~ sne, data = clouds, subset = seeding == "yes"))
clouds.predcit <- data.frame(Rainfall=clouds$rainfall, No=seeding.no, Yes=seeding.yes)

ggplot(clouds, aes(x=sne, y=rainfall, colour=clouds$seeding, shape=clouds$seeding) + 
  geom_point(size=2) +
  labs(title="Rainfall vs S-Ne Criterion", x="S-Ne Criterion", y="Rainfall") +
  scale_color_manual(values=c("blue", "red")) + 
  scale_shape_manual(values=c(1,2)) +
  geom_smooth(method="lm") +
  theme_minimal()

如果您想轻松地重命名您的图例，只需将 colour = "your legend title" 添加到 labs() 参数中。