统计模型中的自举混合效应回归系数

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笔记本以了解如何实现。