r - ggplot 中的 ANCOVA 图

Question

以下代码绘制了 ANCOVA 的结果:

library(HH)
ancova(weight ~ gesttime + dose, data=litter)

如果我尝试在 ggplot 中绘制相同的图，最佳拟合线的斜率在不同的剂量水平上是不均匀的。

library(ggplot2)
ggplot(litter, aes(gesttime, weight)) + geom_point() + facet_grid(. ~ dose) + stat_smooth(method="lm")

如何使用 ggplot 绘制与 ancova() 输出相同的图？

Answer 1

使用 ggplot，您需要使用 geom_line() 绘制模型的预测值，如下所示以生成您想要的结果。

让你的模型合身:

data(litter)
mod <- ancova(weight ~ gesttime + dose, data=litter)
pred <- predict(mod)

绘制它:

ggplot(data = cbind(litter, pred),
    aes(gesttime, weight, color=dose)) + geom_point() +
    facet_grid(. ~ dose) + geom_line(aes(y=pred))