在 statsmodels 中自举混合效应回归系数

Question

我有一个看起来像这样的混合效果模型：

import statsmodels.formula.api as smf

formula = "revised_error ~ C(condition, Treatment('solo_feedback'))*round_index"

model = smf.mixedlm(forumla, data=data, groups=data['player_id']).fit()

现在 model.conf_int() 得到了基于标准正态分布的拟合参数的置信区间。但是，我想要自举系数和置信度，但 model.boostrap() 抛出错误：

AttributeError: 'MixedLMResults' object has no attribute 'endog'

Answer 1

首先，您的代码有一个拼写错误，如果您没有更正它，可能会预先导致一些问题：formula

请参阅 this page 上的说明，了解为什么这对您不起作用。到目前为止，复杂模型还不支持此功能。但是，您可以通过手动而不是使用 bootstrap 方法来解决这个问题。

你可以关注 this notebook 看看它是如何实现的。