广义线性混合效应模型

广义线性混合效应 (GLIMMIX) 模型是线性预测器中包含随机效应的广义线性模型。statsmodels 目前支持使用两种贝叶斯方法估计二项式和泊松 GLIMMIX 模型:后验的拉普拉斯近似和后验的变分贝叶斯近似。两种方法都提供点估计(后验均值)和对不确定性的评估(后验标准差)。

当前的实现仅支持独立随机效应。

技术文档

与 statsmodels 混合线性模型不同,GLIMMIX 实现不是基于组的。组是通过使所有随机效应与分类变量交互来创建的。请注意,这会创建大型稀疏随机效应设计矩阵 exog_vc。在内部,exog_vc 被转换为 scipy 稀疏矩阵。当直接将参数传递给类初始化器时,可以传递稀疏矩阵。当使用公式时,会创建一个密集矩阵,然后将其转换为稀疏矩阵。对于非常大的问题,由于此密集中间矩阵的大小,可能无法使用公式。

参考

Blei, Kucukelbir, McAuliffe (2017)。变分推断:统计学家的综述 https://arxiv.org/pdf/1601.00670.pdf

模块参考

模型类为

BinomialBayesMixedGLM(endog, exog, exog_vc, ...)

具有贝叶斯估计的广义线性混合模型

PoissonBayesMixedGLM(endog, exog, exog_vc, ident)

具有贝叶斯估计的广义线性混合模型

结果类为

BayesMixedGLMResults(model, params, cov_params)

用于保存混合 GLM 模型的贝叶斯估计结果的类。


最后更新:2024 年 10 月 3 日