运行标准最小二乘混合模型时,为什么会出现“收敛有问题:检查迭代”?大家好,
我是一名生物医学博士生,并假设特定细胞类型随着药物治疗时间的不同而增加(例如,将 PBS 对照与治疗第 12 天和第 19 天进行比较)。 我使用了 JMP pro12 并运行了标准最小二乘模型。 我进行了两个实验来测试这一点,并将实验代码设置为随机效应。 另一个作用是治疗。
我正在查看固定效果测试 -> 效果详细信息 -> LS 表示差异 Tukeys HSD -> 有序差异报告,以查看处理是否显着不同。 固定效应检验在 prob>f 列中有一个点而不是 ap 值,并且 Lower CL/Upper CL 列中也有点。 因此,我不确定是否相信 p 值(<0.0001)。
我认为这可能与“收敛有问题:检查迭代”错误有关?
你知道如何解决这个问题吗? 我是否可以相信 p 值并说我的所有组之间存在显着差异,或者我的数据不适合此测试?
多谢!
艾玛


回复:为什么在运行标准最小二乘混合模型时会出现“收敛有问题:检查迭代”?这一切都有道理,艾玛。
回复:为什么在运行标准最小二乘混合模型时会出现“收敛有问题:检查迭代”?真的很有帮助,非常感谢。 我收到 JMP 技术支持的以下回复,我认为它已经回答了我的问题:
嗨艾玛,
您有关收敛失败的问题已转交给 JMP 技术支持。
在检查你的模型时,我注意到了一些事情:
在最后一次迭代中,Expt 的方差分量估计值为负,这表明 Expt 之间可能没有足够的变异证据来对随机效应进行建模。
事实上,在绘制数据时,似乎没有明显的证据表明 Expt 之间存在任何差异。 该图还强调了只有几行数据的治疗组合估计的不确定性。

我认为这就是不稳定和无法收敛到解决方案的原因。
请注意,如果您在模型对话框中取消选择“无界方差分量”(绑定 vc),则模型将收敛到 Expt 随机效应方差估计为零的解。
该模型相当于删除 Expt 随机效应并仅拟合固定治疗效应的模型。 我认为这是一个合理的方法。
我希望这些信息对你有用。 如果您有后续问题或需要澄清,请告诉我。
谢谢,
亚当·莫里斯
JMP 技术支持
回复:为什么在运行标准最小二乘混合模型时会出现“收敛有问题:检查迭代”?其他一些想法:
1. 在 JMP,我们总是说您应该在进行统计分析之前绘制数据。 我们为此设计了 JMP。 当您绘制您的反应与治疗的图表时,看起来是否存在很大差异? 数据看起来合理吗? 有什么奇怪的异常值吗?
2. 我想绘制响应的分布图。 从快速查看表格的屏幕截图来看,它似乎在数量级上变化。 因此响应的转换可能对建模有用。回复:为什么在运行标准最小二乘混合模型时会出现“收敛有问题:检查迭代”?嗨艾玛,
如果您能附上 JMP 数据表将会有所帮助。
我想知道 Treatment 和 Expt 之间的相关性,如果您可以提供 JMP 表,那么探索这一点会容易得多。
谢谢,
菲尔