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

によって出力されたものと同じプロットを作成するにはどうすればよい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))