β在具有相应置信区间的不同值处的边际效应

问题描述

基本上，我有一个混合效应模型，其中包含2个预测变量，一个交互项和一个随机效应，如下所示：

x1<-abs(rnorm(300)) #continuous variable
x2<-sample(c(1,2,3),300,replace = T) #discrete variable 
x3<-sample(c("A","B","C"),replace = T) #categorical variable 
y<-abs(rnorm(300)) #continuous variable
df<-as.data.frame(cbind(as.numeric(x1),as.numeric(x2),x3,as.numeric(y)))

m1=lmer(y~x1*x2+(1|x3),data=df)
summary(m1)

我用lme4创建了这样一个模型。 我想在x2的不同级别（这是3个级别的因子）上找到x1对y（β系数）的边际效应。

我还想获得每个级别的beta的置信区间，并绘制出这样的结果。

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）