安装.包( " tidycomprsk公司 " )
# install.packages(“devtools”)
开发工具 :: 安装github( " MSKCC-Epi-Bio/tidycmprsk公司 " )
图书馆( tidycomprsk公司 ) crr修改 <- crr(测量( tt死亡 , 死亡c r ) ~ 年龄 + 信托收据 , 试验 ) # >11例因缺失值而遗漏
crr修改
# >
# >——crr() ───────────────────────────────────────────────────────────────────────
# >•呼叫Surv(ttdeath,death_cr)~年龄+trt
# >•兴趣失败类型“癌症死亡”
# >
# >可变系数SE HR 95%CI p值
# >年龄0.006 0.010 1.01 0.99,1.03 0.56
# >trt药物B 0.417 0.279 1.52 0.88,2.62 0.13
tbl(待定) <-
crr修改 % > % 全球技术法规摘要 :: tbl回归( 指数化 = 真的 ) % > % 全球技术法规摘要 :: add_global_p()% > % 添加(_n)( 位置 = " 水平 " )
全球技术法规摘要 :: 内联文本( tbl(待定) , 变量 = 年龄 ) # >[1]“1.01(95%置信区间0.99,1.03;p=0.6)”
cuminc(苏尔夫( tt死亡 , 死亡cr ) ~ 1 , 试验 ) # >
# >——cuminc() ────────────────────────────────────────────────────────────────────
# >•失败类型“癌症死亡”
# >时间n.风险估计标准误差95%置信区间
# >5.00 199 0.000 0.000不适用,不适用
# > 10.0 189 0.030 0.012 0.012, 0.061
# > 15.0 158 0.120 0.023 0.079, 0.169
# > 20.0 116 0.215 0.029 0.161, 0.274
# >•故障类型“死亡-其他原因”
# >时间n.风险估计标准误差95%置信区间
# > 5.00 199 0.005 0.005 0.000, 0.026
# > 10.0 189 0.025 0.011 0.009, 0.054
# > 15.0 158 0.090 0.020 0.055, 0.135
# > 20.0 116 0.205 0.029 0.152, 0.264
图书馆( ggsurvfit公司 ) # >加载所需包:ggplot2 cuminc(苏尔夫( tt死亡 , 死亡c r ) ~ 信托收据 , 试验 ) % > % ggcuminc() + 添加确认间隔() + add_risktable() + 缩放ggsurvfit( x _刻度 = 列表 ( 打破 = 序列号( 0 , 24 , 通过 = 6 ))) # >绘制结果“癌症死亡”。
tbl(待定) <- cuminc(苏尔夫( tt死亡 , 死亡cr ) ~ 信托收据 , 试验 ) % > % tbl_cuminc公司( 次 = c(c)( 12 , 24 ), 标签_标题 = " **月份{time}** " ) % > % add_p()% > % 添加_()