我已经运行了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
)。你看到任何潜在的冲突吗?我非常感谢你,
阿维图斯