使用plot_model()绘制原始比例的边际效应

时间:2019-02-19 11:33:24

标签: r lme4 sjplot marginal-effects

我正在尝试使用plot_model()绘制模型(用lme4完成)的边际效应

模型(相同模型的2个版本)如下

版本1:

m1 <- glmer(Y~Xs*Zs+(1|random.factor), family=binomial(), data=dat)

版本2:

m2 <- glmer(Y~scale(X)*scale(Z)+(1|random.factor), family=binomial(), data=dat)

在两个版本中,连续预测变量X和Z均使用scale()函数进行缩放,但是在版本1中,它是在创建模型m1之前完成的,并在{{1中保留为单独的列(X和Z)) }},而在版本2中则在模型m2中指定

我现在的问题是 如果我按原样绘制m1,我会得到

Plot without rescaled x-axis

我想绘制X的边际效应,但我想用原始值代替X轴上X的缩放比例值 我尝试过

dataframe=dat

a <- plot_model(m1, type="pred", terms="Xs [myfun]") b <- plot_model(m2, type="pred", terms="X [myfun]") 被定义为myfun的反向功能

Plot a

Plot b

所以,这给了我两个不同的情节。由于图b看起来与没有重新缩放x轴的m1图相似,我想说图b是正确的图。但是我不明白术语“ ...”中的“ [myfun]”实际上是如何工作的以及它的作用。 感谢您的帮助!

0 个答案:

没有答案