ggplot中的ANCOVA图

Question

下面的代码绘制ANCOVA的结果：

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

​

​

如果我尝试在ggplot中绘制相同的图，则最佳拟合线的斜率在dose的不同级别上不是均匀的。

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))