广义线性混合效应模型¶
广义线性混合效应 (GLIMMIX) 模型是线性预测器中包含随机效应的广义线性模型。statsmodels 目前支持使用两种贝叶斯方法估计二项式和泊松 GLIMMIX 模型:后验的拉普拉斯近似和后验的变分贝叶斯近似。两种方法都提供点估计(后验均值)和对不确定性的评估(后验标准差)。
当前的实现仅支持独立随机效应。
技术文档¶
与 statsmodels 混合线性模型不同,GLIMMIX 实现不是基于组的。组是通过使所有随机效应与分类变量交互来创建的。请注意,这会创建大型稀疏随机效应设计矩阵 exog_vc。在内部,exog_vc 被转换为 scipy 稀疏矩阵。当直接将参数传递给类初始化器时,可以传递稀疏矩阵。当使用公式时,会创建一个密集矩阵,然后将其转换为稀疏矩阵。对于非常大的问题,由于此密集中间矩阵的大小,可能无法使用公式。
参考¶
Blei, Kucukelbir, McAuliffe (2017)。变分推断:统计学家的综述 https://arxiv.org/pdf/1601.00670.pdf
模块参考¶
模型类为
|
具有贝叶斯估计的广义线性混合模型 |
|
具有贝叶斯估计的广义线性混合模型 |
结果类为
|
用于保存混合 GLM 模型的贝叶斯估计结果的类。 |
最后更新:2024 年 10 月 3 日