在 tidycmprsk package 中使用 cuminc 时估计累积发病率

Question

我正在使用新的 package tidycmprsk并尝试 plot 累积发病率 function 具有大约 800 个观察值的数据集。 但是，当我包含超过 320 个观察值时，我收到以下错误：

Error in if (any(is.na(df$time)) || any(is.na(df$time_chr)) || max(abs(df$time -  : 
  missing value where TRUE/FALSE needed
In addition: Warning message:
In cuminc_matrix_to_df(., name = "estimate", times = times) :
  NAs introduced by coercion

起初，我认为问题出在 nr 320 的观测值上，但包含了这个观测值并减少了包含的观测值的数量，我没有收到任何错误。

知道问题出在哪里吗？

感谢转发并感谢 Daniel 和 Teng 提供的出色的 package。

这是我的代码：

cif <- tidycmprsk::cuminc(Surv(time, outcome) ~ strata, data = x1)

这是我数据的一部分

x1 <- structure(list(strata = c("T2", "T1", "T3", "T2", "T4", "T1", 
"T2", "T4", "T2", "T4", "T4", "T3", "T1", "T1", "T2", "T3"), 
    outcome = structure(c(1L, 1L, 1L, 2L, 3L, 3L, 1L, 2L, 1L, 
    2L, 2L, 2L, 1L, 1L, 1L, 1L), levels = c("Alive or censored", 
    "Death from oral cancer", "Death from other causes"), class = "factor"), 
    time = c(4.96986301369863, 9.4027397260274, 6.15068493150685, 
    1.51232876712329, 1.84383561643836, 8.94520547945205, 5.23561643835616, 
    1.65753424657534, 7.44109589041096, 0.583561643835616, 2.80547945205479, 
    12.6958904109589, 9.36986301369863, 6.3041095890411, 10.7342465753425, 
    6.87945205479452)), row.names = 310:325, class = "data.frame")