在 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请参阅此页面上的注释，了解为什么这对您不起作用。到目前为止，对于复杂模型，本机不支持此功能。但是，您可以通过手动而不是使用该bootstrap方法来解决此问题。你可以关注这个notebook来看看它是如何实现的。