0
我在R中有一个线性混合效应模型(lme4
),我想使用MixedModels
包在Julia中运行它。我的问题是,Julia的MixedModels
不允许在lme4
中使用偏移量变量。有什么建议么?如何在Julia中的线性混合模型中包含偏移量分量?
例子:Y ~ X1 + X2 + X3 + X2|Group + X3|Group + offset(X4)
我在R中有一个线性混合效应模型(lme4
),我想使用MixedModels
包在Julia中运行它。我的问题是,Julia的MixedModels
不允许在lme4
中使用偏移量变量。有什么建议么?如何在Julia中的线性混合模型中包含偏移量分量?
例子:Y ~ X1 + X2 + X3 + X2|Group + X3|Group + offset(X4)
如果你真的运行线性混合模式(这是,我相信,所有的朱莉娅MixedModels
目前的可用),那么你应该只能够减去从响应变量的偏移,即对于线性模型(或更一般的广义线性(混合)机型的身份链接),你可以只使用
(Y-X4)~X1+X2+X3+X2|Group+X3|Group)
(我不知道是否MixedModels
是足够聪明的评估Formu的LHS la,或者您是否需要预先计算Y2 = Y-X4
,然后使用Y2 ~ ...
)
您的建议(从因变量减去偏移量)对于_linear_模型是正确的。起初,我用R中的'glmer'在logit链接函数中尝试了它,并在使用偏移功能和减法版本之间获得了不同的结果。在纸上做数学后,减法只适用于_linear_模型是合理的。 – VassilisGeorgiou
是的,这就是为什么为GLMM实施偏移量非常重要(而LMMs只是便利/糖) –