HH软件包似乎提供了一种简单的方法来可视化ANCOVA。
library(HH)
data(hotdog)
ancova(Sodium ~ Calories * Type, data=hotdog)
有没有一种舒适的方式将其与类似于的
panel.ablineq
组合起来? (http://latticeextra.r-forge.r-project.org#panel.ablineq),以获取具体的函数(斜率,截距)?HH软件包似乎提供了一种简单的方法来可视化ANCOVA。
library(HH)
data(hotdog)
ancova(Sodium ~ Calories * Type, data=hotdog)
panel.ablineq
组合起来? (http://latticeextra.r-forge.r-project.org#panel.ablineq),以获取具体的函数(斜率,截距)?一般来说,直接使用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))