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只是便利/糖) –