r - 具有竞争风险的生存分析 R 中的 Fine 和 Gray 模型 - 错误

Question

我的论文是关于研究影响 COVID-19 死亡率的因素。在 3 月 2 日至 12 月 31 日期间，对所有病毒检测呈阳性的人进行了追踪，直至 12 月 31 日。对我的研究很重要的变量是性别、年龄组和地区，所有因素。

我的数据如下所示：

369689 obs. of  5 variables:
 $ status   : Factor w/ 4 levels "Ativo","Morte COVID19",..: 1 4 4 1 4 4 4 4 4 4 ...
 $ time     : int  12 14 14 3 14 14 14 31 11 14 ...
 $ age_group: Factor w/ 5 levels "]20-60]","]60-70]",..: 1 1 1 1 1 1 1 1 1 1 ...
 $ sex      : Factor w/ 2 levels "F","M": 1 2 2 2 1 2 2 1 2 1 ...
 $ region   : Factor w/ 7 levels "Açores","Alentejo",..: 7 7 7 7 4 4 4 5 4 1 ...

有了这些数据，我正在做一个存在竞争风险的生存分析。我正在关注 Zhang Z. 在存在竞争风险的情况下的生存分析。Ann Transl Med 2017；5(3):47。doi: 10.21037/atm.2016.08.62 文章做我的分析。

我的问题是当我在 R 中执行命令时：

SH <- FGR(Hist(time,status)~sex+age+invasion,data=Melanoma)

当我适应我的数据时，它根本不运行。我尝试使用 R 云，但它打破了......

然后我尝试使用替代命令

cov<-model.matrix(~sex+age+invasion,data=Melanoma)[,-1]
crr.model<-crr(Melanoma$time,Melanoma$status,cov1=cov)

我得到这个错误：

Error in solve.default(v[[1]]) : 
  Lapack routine dgesv: system is exactly singular: U[1,1] = 0

我真的需要做精细和灰色模型......你知道另一种获得模型的方法吗？

标签： rsurvival-analysis