我创建了一个与 iv*sex 的交互项,并将数据与小鼠进行插补。然后使用估算数据运行逻辑回归模型(glm):
model <- with(data=imp, glm(dv~control+iv+sex+iv*sex, family="binomial"))
以下是变量名的缩写: 因变量=dv,自变量=iv,调节器=性别,交互项= iv*性别
iv*sex 有显着的相互作用,我想绘制一个相互作用的图表,但我找不到如何做。如果提供任何解决方案,我们将不胜感激。谢谢!
我刚刚遇到了同样的问题,通过
effects
包我解决了它。
e <- effects::effect("iv*sex", model)
e <- as.data.frame(e)
ggplot2::ggplot(e, ggplot2::aes(iv, fit, color=sex, group = sex)) +
ggplot2::geom_point() +
ggplot2::geom_line() +
“fit”是你的因变量,在你的例子中是“dv”。