r - ggplot2 を使用した相互作用プロット

Question

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 を使用して同じプロットを作成するにはどうすればよいですか?

Answer 1

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")