Question

我正在使用cmprsk软件包中的crr（）函数来执行竞争风险回归。但是，它失败并显示警告“ crr convered：FALSE”。我已经进行了很多次竞争风险回归，这是我第一次遇到这个问题。

    x_total2 <- cbind(factor2ind(m1_sub3), factor2ind(age1, "1"), factor2ind(sex2, "1"), factor2ind(race2, "1"), 
              factor2ind(marital2, "1"), factor2ind(grade2, "1"), factor2ind(t_stage2, "1"), 
              factor2ind(n_stage2, "1"), factor2ind(ssct2, "0"))

    mod_nrt1 <- crr(surv_mo, sum_cod3, x_total2[,c(1:24)])

    summary.crr(mod_nrt1, digits = max(options()$digits - 4, 3))

这是警告。

    Competing Risks Regression

    Call:
    crr(ftime = surv_mo, fstatus = sum_cod3, cov1 = x_total2[, c(1:24)])

    crr converged: FALSE

此警告后，我使用了

    mod_nrt1[["converged"]] <- TRUE

通过此功能，我可以获得以下结果。

               exp(coef) exp(-coef)     2.5%    97.5%
    sex2:2      4.36e+00   2.30e-01 1.28e+00 1.48e+01
    race2:2     1.30e+00   7.67e-01 8.71e-01 1.95e+00
    race2:3     1.17e+00   8.58e-01 6.58e-01 2.07e+00
    race2:4     2.42e+00   4.13e-01 1.18e+00 4.98e+00
    race2:5     3.57e-05   2.80e+04 4.15e-06 3.07e-04
    t_stage2:2  2.60e-01   3.85e+00 1.10e-01 6.12e-01
    t_stage2:3  2.75e-01   3.64e+00 1.16e-01 6.51e-01
    t_stage2:4  4.08e-01   2.45e+00 1.53e-01 1.09e+00
    t_stage2:5  2.74e-01   3.65e+00 1.05e-01 7.14e-01
    t_stage2:6  8.58e-01   1.17e+00 2.80e-01 2.63e+00
    n_stage2:2  1.01e+00   9.91e-01 4.70e-01 2.16e+00
    n_stage2:3  6.97e-01   1.44e+00 2.74e-01 1.77e+00
    n_stage2:4  1.26e+00   7.95e-01 5.30e-01 2.98e+00
    n_stage2:5  3.22e+00   3.10e-01 1.01e+00 1.02e+01
    ssct2:0     3.74e+00   2.68e-01 1.75e+00 7.96e+00
    ssct2:1     2.92e+00   3.42e-01 1.36e+00 6.27e+00
    ssct2:3     9.93e-01   1.01e+00 6.49e-01 1.52e+00

我知道这是错误的。什么是完美的解决方案？原始数据有什么问题吗？

Answer 1

尝试增加 maxiter（默认为 10）

使用cmprsk软件包时竞争风险回归中的错误收敛

1 个答案: