我正在尝试使用visreg
软件包进行三向交互。 (我可以使用我的R可复制代码来适应如下所示的双向交互。)
我想知道是否有可能使用visreg
:“ year”绘制三向交互关系; “年龄”; “教育”。使用ggeffects
包是可能的-但使用visreg也可以吗?
(我以前已经在一个星期前在r-help上发布了相同的请求,但没有得到任何答案,在CrossValidated上也没有得到任何答案(答案:“该网站的主题”))。
#install.packages("ISLR")
library(ISLR)
#install.packages("visreg")
library(visreg)
#install.packages("mgcv")
library(mgcv)
mod1 <- gam(wage ~ education+s(age,bs="ps") + year + te(age,year,bs="ps") +
s(year,bs="ps",by=education,m=1) + te(year,age,by=education,bs=rep("ps",2)),
data=Wage)
visreg(mod1, "year", by="education",partial=FALSE,overlay=TRUE)