为什么 geeglm 不在这里收敛？

问：

我正在尝试使用 geeglm 在 GEE 模型中实现倾向得分加权。我还使用引导来估计标准误差。但是，我在 geeglm 中遇到了一个奇怪的收敛相关问题。

我相信我已经编写了正确实现这些方法所需的一切。但是，当使用 bootstrap 示例运行 geeglm 时，偶尔会出现一个特殊的示例，这似乎会导致不收敛。但是，该函数似乎在未收敛时卡住，但不会中断并给出收敛错误。因此，停止的唯一方法是按下 Rstudio 中的停止按钮，这需要我中止所有内容并重新启动 Rstudio 从头开始运行所有内容。

下面是我从控制设置中的跟踪选项获得的迭代序列示例，它允许我打印出拟合算法每次迭代的详细信息。

iter = 0 贝塔 = -6.129982 0.157005 0.069532 0.226057 -0.403785 -0.050105 -0.163753 -0.205799 0.617615 -0.331346 伽玛 = 653.576761 阿尔法 = 0.000000

ITER = 1 贝塔 = -6.129982 0.157005 0.069532 0.226057 -0.403785 -0.050105 -0.163753 -0.205799 0.617615 -0.331346 伽玛 = 0.932907 阿尔法 = 0.459056

ITER = 2 贝塔 = -6.145171 0.078880 -0.014521 0.315024 -0.338607 -0.046403 -0.150386 -0.185138 0.463438 -0.317506 伽玛 = 0.955960 阿尔法 = 1.002303

热膨胀系数 = 3 贝塔 = -11.616787 1.595249 -0.368715 2.298262 0.478158 -0.002019 -0.842909 -4.365958 -0.839819 -1.269375 伽玛 = 1065.562333 阿尔法 = 1.400400

ITER = 4 贝塔 = 396.013316 -72.850710 -52.460387 -240.090210 -34.443043 0.425689 -74.907351 414.566790 30.545857 45.713360 伽玛 = 南 阿尔法 = 南

当它进入迭代 4 时，它显然没有收敛，但它似乎只是停留在保持模式中，没有爆发并给出非收敛错误。唯一的方法是完全中止并重新启动 Rstudio。

我已经能够通过指定 scale.fix=TRUE 来避免非收敛，但这对我作为解决方案没有帮助。

我预计如果 R 卡住并且没有收敛，那么它就会爆发并给出非收敛错误。但是对于geeglm，情况似乎并非如此？关于如何获得非收敛错误的任何建议都会有所帮助。