我正在关注R中brms
软件包上的this教程。有一个使用marginal_plot
来绘制交互的函数。代码如下:
plot(marginal_effects(fit1, effects = "log_Base4_c:Trt"))
我可以使用其他参数points = TRUE
向其添加数据点:
plot(marginal_effects(fit1, effects = "log_Base4_c:Trt"),points = TRUE)
我还可以在交互中拟合连续预测变量(年龄)并绘制交互图。
fit2 <- brm(count ~ log_Age_c + log_Base4_c * Age + (1|patient),
data = epilepsy, family = poisson())
plot(marginal_effects(fit2, effects = "log_Base4_c:Age"),points=TRUE)
我想知道如何控制边际图产生的范围?我得到的曲线分别为34.6、28.34和22.08。也可以对它们进行颜色编码,以使数据点与曲线匹配。
谢谢