夜雨聆风
有人咨询，搜索了许多方法均绘制不出PCM CIF曲线
疑问：《R语言进行竞争风险模型分析》中只绘制了年龄分组的CIF曲线，没有绘制PCM CIF，怎么才能绘制出PCM CIF曲线
今天在此讲一下，如何用代码模版绘制自己需要的图形
原代码点击查看《R语言进行竞争风险模型分析》
1.加载必要的包
library(tidycmprsk)
library(ggsurvfit)
library(dplyr)
library(ggplot2)
2.加载数据
library(survival)
data(mgus2,package = "survival")
head(mgus2)
3.数据预处理
df <- mgus2 %>% mutate(
futime = ifelse(pstat == 1, ptime, futime),
event = case_when(
pstat == 1 ~ 1,
death == 1 ~ 2,
TRUE ~ 0 ),
event_fct = factor(event,
levels = c(0, 1, 2),
labels = c("censor", "pcm", "death")),
age_grp = ifelse(age > 65, ">65", "≤65"),
sex = factor(sex) )
table(df$event_fct)
4.累积发生率（CIF）
cuminc_all <-cuminc(Surv(futime,event_fct)~1,data =df)
print(cuminc_all)
5.绘制CIF曲线
ggcuminc(cuminc_all, outcome = "pcm") +
add_confidence_interval() +
add_risktable() +
labs(
title = "总体PCM累积发生率曲线",
x = "随访时间（月）",
y = "PCM累积发生率"
) +
theme_minimal()
6 如何替换模版？
原模版代码：
ggcuminc(cuminc_age) +
add_confidence_interval() +
add_risktable() +
labs(
title = "累积发生率曲线（按年龄分组）",
x = "随访时间（月）",
y = "累积发生率",
color = "年龄组"
) +
theme_minimal()
替换后代码：
ggcuminc(cuminc_all, outcome = "pcm") +
add_confidence_interval() +
add_risktable() +
labs(
title = "总体PCM累积发生率曲线",
x = "随访时间（月）",
y = "PCM累积发生率"
) +
theme_minimal()
替换内容：
配套数据点击《SPSS和R语言文章配套数据汇总（持续更新）》获取
不知道怎么选择统计方法点击《医学统计免费咨询》咨询
进入交流群点击《统计学习交流群 · 入群指南》