ggplot 中的 ANCOVA 图答案

【问题标题】：ANCOVA plot in ggplotggplot 中的 ANCOVA 图
【发布时间】：2013-05-30 08:04:58
【问题描述】：

以下代码绘制 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")

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

【问题讨论】：

嗨@luciano，你如何计算你的回归权重~ gesttime 它因不同的因素而不同？我的意思是，你如何证明这种相关性对于因子“dose = 0”在统计上是不同的

标签： r ggplot2

【解决方案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))

【讨论】：