R 可视化包括公式在内的ANCOVA（例如库HH）

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

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

---

是否有一种舒适的方式将其与

latticeExtra

中的

panel.ablineq

这样的东西结合起来？（）要获得具体的函数（斜率、截距）

一般的答案是直接使用

HH:：ancovaplot（）

（而不是通过

HH:：ancova（）

隐式地）以及latticeExtra的便捷

层（）

函数和重载的

+

操作符。我将把它留给你去解决需要的微妙细节，使它看起来正是你想要的

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