可视化 ANCOVA 包括公式（例如库 HH）

Question

HH 包似乎提供了一种可视化 ANCOVA 的简单方法

library(HH)
data(hotdog)
ancova(Sodium ~ Calories * Type, data=hotdog)

有没有一种方便的方法可以将它与 latticeExtra 中的 panel.ablineq 之类的东西结合起来？ (http://latticeextra.r-forge.r-project.org#panel.ablineq) 得到具体函数(斜率、截距)?

Answer 1

一般的答案是直接使用 HH::ancovaplot()（而不是通过 HH::ancova() 隐含地使用）以及 latticeExtra 的方便 layer()函数和重载 + 运算符。我将留给您解决使它看起来完全符合您的要求所需的繁琐细节。

ancovaplot(Sodium ~ Calories * Type, data=hotdog) +
layer(panel.ablineq(lm(y ~ x), rot = TRUE, at = 0.5, pos = 3))

Visualize ANCOVA incl formulas (e.g. library HH)