我已经运行glm（用于逻辑回归），结果我得到算法发散。我想试试这个功能logistf来检查 Firth 的方法是否可以改进我的数据分析。让dat表示我的数据框并y ~ f1 + f2成为我的模型。试
fit <- logistf(y ~ f1 + f2, dat)
或者
fit <- logistf(y ~ f1 + f2, data=dat)
我总是得到错误
Error in logistf.fit(x = x, y = y, weight = weight, offset = offset, firth, :
no loop for break/next, jumping to top level
由于glm运行顺利，我看不出我的数据结构如何成为错误的根源（这就是我没有dat明确描述的原因）。你看到任何潜在的冲突吗？我非常感谢你，阿维图斯