当无法有效使用rms中的lrm函数的一个应用程序时，如何防止仿真崩溃？

Question

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

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

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

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

Answer 1

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

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

与lrm相关的lrm将在每次迭代中更改（例如x ）。