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

Question

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

模拟因错误消息而崩溃： 无法使用＆＃34; lrm.fit＆＃34;。

来拟合模型

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

有没有办法可以做到这一点？

Answer 1

考虑使用try报告错误但不退出整个循环或函数。

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

与lrm相关的内容在每次迭代时都会发生变化（例如x）。