r - 当 rms 中的 lrm 函数的一个应用程序不适合时，如何防止模拟崩溃？

Question

我正在运行 1000 次迭代的蒙特卡洛模拟。在每次迭代中，我使用 Harrell 的 rms 包中的 lrm 函数拟合加权逻辑回归模型。使用以下代码拟合模型： lrm(y ~ x, weights=wt,x=T,y=T) 从拟合模型中，我提取了一些信息，例如回归系数和估计的标准误差。

模拟因错误消息而崩溃：无法使用“lrm.fit”拟合模型。

我想通过仅在安全的情况下评估功能来防止模拟崩溃。在绝大多数迭代中，没有问题。不知何故，在每次迭代中，我想告诉 R 只有在可以安全完成的情况下才适合该函数。

有没有办法做到这一点？

Answer 1

Consider using try which will report an error but not exit the whole loop or function.

for (i in 1:10){
    try(lrm(y ~ x, weights=wt,x=T,y=T) )
}

Where something relevant to lrm would be changing on each iteration (such as x for example).