对于混合效应线性模型，我有收敛性问题。我想了解一下优化过程中发生的情况。有没有办法获取迭代日志？我现在能做的最好的就是通过设置disp = True
来获得优化的摘要。mdf = md.fit（full_output = True，reml = True，method ='cg'，disp = True）
这给了我
警告：由于精度损失而未必达到所需的误差。 当前功能值：2.989371 迭代次数：5 功能评估：73 梯度评估：62
谢谢
更新：这没有回答我的问题，但是使用其他求解器，我设法获得了一些收敛。但是，这提出了另一个问题。我希望得分，即对数似然函数的梯度很小，对于我的示例来说并非如此。因此，另一个问题是：
什么可以被信任？凭直觉，我的答案就是分数