r 中 logistf.fit 中的错误是什么意思“Fisher 信息矩阵的行列式在数值上为 0”

Error in logistf.fit in r what does it mean "determinant of Fisher information matrix was numerically 0"

提问人:Xiaoyong Fu 提问时间:11/4/2023 最后编辑:Xiaoyong Fu 更新时间:11/6/2023 访问量:40

问:

我正在尝试使用 r 中的 logistf 包进行 Firth 校正的逻辑建模。我有一个数据框 (df),其中包含 110 个预测协变量,第一列 (nonR) 中有一个结果。我粘贴到此数据集的子集下面,该子集重现了相同的错误:

Error in logistpl(x, y, beta, i, LL.0, firth, -1, offset, weight, plcontrol,  : 
  In iteration 63: Determinant of Fisher information matrix was numerically 0
In addition: Warning message:
In logistf(formula = as.formula(paste0("nonR ~", paste0(names(df)[2:10],  :
  logistf.fit: Maximum number of iterations for full model exceeded. Try to increase the number of iterations or alter step size by passing 'logistf.control(maxit=..., maxstep=...)' to parameter control

我跑完之后:

lf <- logistf(formula= as.formula(paste0("nonR ~", paste0(names(df)[2:10], collapse = " + "))), data = df)
     nonR   APC ARID1A ARID1B   ATM  BAP1 BCL11B BIRC3  BRAF CDKN2A
   <dbl> <dbl>  <dbl>  <dbl> <dbl> <dbl>  <dbl> <dbl> <dbl>  <dbl>
 1     1     0      0      0     1     0      0     0     0      0
 2     1     0      0      0     0     0      0     1     0      0
 3     0     0      0      0     0     0      0     0     1      0
 4     0     0      0      0     0     0      0     0     0      1
 5     1     1      0      0     0     0      0     0     0      0
 6     0     0      0      1     0     1      0     0     0      0
 7     1     0      0      0     0     0      0     0     0      1
 8     0     0      0      0     0     0      1     0     0      0
 9     1     0      1      0     0     0      0     0     0      0
10     1     0      0      0     0     0      0     0     1      0
11     1     0      0      0     0     0      0     0     1      0
12     0     0      0      0     0     0      0     0     0      1

有人可以向我解释这意味着什么以及如何解决这个问题吗?我在网上搜索了一下,但没有这样的信息来解决这个问题。先谢谢你。

R 逻辑回归

评论

0赞 Daniel 11/5/2023
我认为这里提出了一个类似的问题:stackoverflow.com/questions/74570346/......与这个问题类似,为可重现的例子提供一些(假)数据会有所帮助。

答: 暂无答案