взаимодействие участка с использованием 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")