多项混合效应逻辑回归

Question

是否有任何软件包可以在 Python 中执行多项逻辑混合效应回归？ pymer4 可以进行混合效应逻辑回归，但我找不到实现其多项式版本的包。

Answer 1

您可以通过将广义线性模型的 GLM 类与混合效应建模的 MixedLM 类相结合来组合“statsmodels”来拟合混合效应多项逻辑回归模型

import statsmodels.api as sm
import statsmodels.formula.api as smf
    
formula = 'Response ~ Predictor1 + Predictor2 + (1 | RandomEffect)'
    
model = smf.mnlogit(formula, data=df).fit()
    
print(model.summary())