我是否可以信任仅使用bobyqa和对比度和编码收敛的完整glmer模型?

时间:2016-02-20 13:07:03

标签: r logistic-regression lme4

我正在使用R 3.2.0和lme4版本1.1.8。从心理语言学实验中对一些二项式数据(编码为0和1)运行混合效应逻辑回归模型。有2个分类预测因子(一个具有2个级别,一个具有3个级别)和两个随机术语(参与者和项目)。我使用和编码预测器(即contr.sum ..),它给了我感兴趣的效果和相互作用。

我发现完整的模型(具有固定的效果和交互,加上两个随机项的随机截距和斜率)仅在我指定(optimizer =“bobyqa”)时收敛。如果我没有指定优化器,模型只有在简化模型后才会收敛。当我使用默认的处理编码时,即使我指定optimizer =“bobyqa”,也会发生同样的事情。

我的第一个问题是为什么会发生这种情况,我能相信完整模型的输出吗?

我的第二个问题是,这是否可能是由于我的数据不完全平衡,因为我的条件没有完全相同的观察数。当数据不完全平衡时,是否必须采取特殊预防措施?可以建议对这个特例进行任何阅读吗?

非常感谢

1 个答案:

答案 0 :(得分:1)

您应该查看?convergence更新版本的lme4帮助页面(或者您可以阅读here)。如果两个使用不同的优化器拟合给出类似的估计参数(尽管一个给出收敛警告而另一个没有),并且不同对比的拟合给出相同的对数似然,那么你可能有合理的拟合。

一般来说,缺乏平衡会降低统计功效并使拟合变得更加困难,但温和至中度的不平衡数据应该不存在特殊问题。