文章预览
今天来学习使用mstate包绘制竞争风险曲线。 这个包可以借助ggplot2绘制图形。 1. 加载R包 library (mstate) library (dplyr) 2. 加载数据 本案例数据来自R自带数据集(survival包的mgus2)，共计1341名单克隆丙种球蛋白病患者，结局事件定义为"发生浆细胞恶性肿瘤"，而某些患者在"发生浆细胞恶性肿瘤"之前因为其他原因死亡，那这些发生其他死亡的患者就无法观察到"发生浆细胞恶性肿瘤"的终点，也就是说"其他死亡"与"发生浆细胞恶性肿瘤"存在竞争风险，故采用竞争风险模型分析。 加载数据，并定义结局为因子变量。 dt dt$etime 0, dt$futime, dt$ptime) dt$event 0, 2 *dt$death, 1 ) dt$event 0: 2 ) glimpse(dt) # 数据说明 age # 年龄 sex # 性别 dxyr # 疾病诊断的年份 hgb # 血红蛋白；\ creat # 肌酐；\ mspike # 单克隆血清样的大小 ptime # 直至发展为浆细胞
………………………………
