2

嗨，我目前正在使用 gam 模型估计幸福和年龄。我的回归命令是

    fit <- gam(happy ~ s(age) + s(age, by=nochild) ,data=happy)

我想可视化交互式变量“s(age):nochild”的估计效果visreg。我正在努力这样做，因为在使用时只出现了对年龄和 nochild 的估计

    visreg(fit)

我希望图表看起来像这样

  • Y轴：s（年龄）：nochild

  • X轴：年龄

但是，我不确定如何执行此操作，visreg或者是否有任何其他命令可以帮助执行此操作？

非常感谢你提前

4

1 回答 1

1

我得到了我现在需要的代码，所以我想我可以分享一下。

    plot(fit, select=2, main="No Children")
于 2016-03-15T13:41:59.793 回答