如何查看基于cuminc对象的ggcompetingrisks模型计算的95%置信区间？

• 手动计算置信区间：cuminc对象的核心数据存于列表中，你可以用估计值结合方差推导置信区间（默认95%置信区间用1.96倍标准误）：

  # 假设你的cuminc对象名为cr_obj
  # 提取基础数据
  time_points <- cr_obj$time
  estimates <- cr_obj$est
  variances <- cr_obj$var
  
  # 计算标准误和95%置信区间
  se_values <- sqrt(variances)
  lower_ci <- estimates - 1.96 * se_values
  upper_ci <- estimates + 1.96 * se_values
  
  # 整理为数据框查看
  ci_results <- data.frame(
    time = time_points,
    estimate = estimates,
    lower_95ci = lower_ci,
    upper_95ci = upper_ci
  )
  print(ci_results)
  

• 用broom包快速提取：借助broom的tidy()函数可以直接将cuminc对象转换为包含置信区间的整洁数据框，无需手动计算：

  library(broom)
  
  # 转换并获取置信区间
  tidy_cr_data <- tidy(cr_obj)
  # 输出包含time、estimate、conf.low、conf.high等列
  print(tidy_cr_data)
  

  其中conf.low和conf.high就是对应时间点的置信区间上下限，和ggcompetingrisks(conf.int = TRUE)绘制的区间完全一致。

内容的提问来源于stack exchange，提问作者Alex