以下代码显示了使用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创建相同的图?
答案 0 :(得分: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")