1

以下代码显示了使用该effects库的交互图:

model <- lm(mpg ~ hp + wt + hp:wt, data=mtcars)
library(effects)
plot(effect("hp:wt", model, list(wt=c(2.2,3.2,4.2))), multiline=TRUE)

在此处输入图像描述

我尝试制作相同的模型但保持wt在 2.2,但该模型无法计算 的系数wt

mtcars$wt_2.2 <- 2.2
model2.2 <- lm(mpg ~ hp + wt_2.2 + hp:wt, data=mtcars)
coef(model2.2)

如何使用 ggplot2 制作相同的情节?

4

1 回答 1

5
tmp <- as.data.frame(effect("hp:wt", model, list(wt=c(2.2,3.2,4.2))))
ggplot(data=tmp, aes(x=hp, y=fit, colour=as.factor(wt))) +
       geom_line() +
       labs(colour="wt")

阴谋

于 2013-11-14T09:53:00.847 回答