HH包似乎提供了一种简单的方法来可视化ANCOVA
library(HH)
data(hotdog)
ancova(Sodium ~ Calories * Type, data=hotdog)

是否有一种舒适的方法将其与类似于panel.ablineq的latticeExtra?(http://latticeextra.r-forge.r-project.org#panel.ablineq)相结合,以获得具体的函数(斜率、拦截)?
发布于 2015-01-18 19:36:44
一般的答案将是直接使用HH::ancovaplot() (而不是通过HH::ancova()隐式地使用),同时使用HH::ancova()方便的layer()函数和重载的+操作符。我将把它留给您去解决所需的繁琐细节,以使它看起来与您想要的完全一致。
ancovaplot(Sodium ~ Calories * Type, data=hotdog) +
layer(panel.ablineq(lm(y ~ x), rot = TRUE, at = 0.5, pos = 3))

https://stackoverflow.com/questions/27988965
复制相似问题