具有连续变量的多级主持调解

Question

我是R的初学者，所以如果我的问题反映出背景不足，请原谅我。

我正在尝试使用mediation和lme4库来运行仲裁中介模型。

我所有的变量都是连续的。我的数据具有嵌套结构，其中个体嵌套在分支（Branch）中。

在我要测试的模型中，我的预测变量/独立变量（abranch）在分支级别。我的调解员（bmed）和结果（cout）处于个人级别。另一个单独的级别变量（dmod）调节了中介的作用。因此，在我的模型中，我abranch预测bmed，而bmed*dmod预测cout。

这是我使用的语法：

med.fit <- glmer(
    bmed ~ abranch + (1|Branch),
    family = binomial(link = "logit"),
    data = Dataset
)

out.fit <- glmer(
    cout ~ dmod*bmed + (1+bmed|Branch),
    family = binomial(link = "logit"),
    data = Dataset
)

当时我正在考虑使用：

med.out <- mediate(med.fit, out.fit, treat = "abranch", mediator = "bmed",
+ sims = 100)
summary(med.out)

但是即使在到达最后两行之前，我仍然遇到以下错误：

Error in eval(family$initialize, rho) : y values must be 0 <= y <= 1

我现在意识到这是因为我使用的是“二项式” / logit系列，而我的DV是连续的，并且不在0到1之间。考虑到变量的性质，我该怎么办？