statsmodels混合模型需要很长时间

Question

我一直在尝试使用statsmodels拟合混合模型，例如：

import statsmodels.formula.api as smf
model = smf.mixedlm('Y ~ cov_1 + cov_2 + cov_3', 
        data, 
        groups = data['grouping_variable']).fit()

其中data是熊猫DataFrame。我有大约10,000个群组。对于我的数据的大多数子集来说，这可以正常工作（大约需要2分钟的时间），但偶尔会挂起并运行数小时以上，但未产生结果。我不知道为什么会这样或如何调试它。