您好我正在使用gam模型估算幸福和年龄。 我对回归的命令是
fit <- gam(happy ~ s(age) + s(age, by=nochild) ,data=happy)
我希望可视化交互变量&#34; s(年龄)的估计效果:nochild&#34;使用visreg
。我正在努力这样做,因为只有在使用
visreg(fit)
我希望图表看起来像这样
Y轴:s(年龄):nochild
X轴:年龄
但是,我不知道如何使用visreg
执行此操作,还是有其他命令可以帮助执行此操作?
非常感谢你提前
答案 0 :(得分:1)
我得到了我现在需要的代码,所以我想我可能会分享这个。
plot(fit, select=2, main="No Children")