如何在Julia的线性混合模型中包含偏移分量？

Question

我在R（lme4）中有一个线性混合效果模型，我想使用MixedModels包在Julia中运行它。我的问题是Julia中的MixedModels不允许lme4中的偏移量变量。有什么建议吗？

示例：Y ~ X1 + X2 + X3 + X2|Group + X3|Group + offset(X4)

Answer 1

如果您真的在运行线性混合模型（我相信，目前Julia MixedModels中提供的所有内容都是如此），那么您应该能够减去从响应变量的偏移量，即对于线性模型（或更常见的具有身份链接的广义线性（混合）模型），您可以使用

(Y-X4)~X1+X2+X3+X2|Group+X3|Group)

（我不知道MixedModels是否足够聪明，可以为您评估公式的LHS，或者您是否必须预先计算Y2 = Y-X4然后使用Y2 ~ ...）