Kaplan Meier我在R markdown公司上班
我想知道我计算的存活时间是否正确。我的数据是MOWHTO_复杂度
我的数据包含以下变量:D_SURGERY(手术日期)、REV_ARTHROPLASTY(翻修日期)、SENSOR_STATUS(0 =删失或1 =翻修)。
翻修日期仅包含翻修病例的日期。其余单元格为空。x1c 0d1x的数据
我用下面的代码计算了存活时间:
MOWHTO_COMPLICATIONS$SURVIVAL_TIME_YEARS = as.numeric(difftime(MOWHTO_COMPLICATIONS$REV_ARTHROPLASTY, MOWHTO_COMPLICATIONS$D_SURGERY, units = "weeks"))/52.25
字符串
然后使用以下代码创建曲线:
survfit2(Surv(SURVIVAL_TIME_YEARS, CENSOR_STATUS) ~ 1, data = MOWHTO_COMPLICATIONS) %>%
ggsurvfit() +
labs(
x = "Years",
y = "Overall Survival Probability"
)+
add_confidence_interval()+
add_risktable()
型
然后我想看看10年后的生存率,用了下面的代码:
summary(survfit(Surv(SURVIVAL_TIME_YEARS, CENSOR_STATUS) ~ 1, data = MOWHTO_COMPLICATIONS), times = 10)
型
我有以下结果调用:生存期(公式=生存期(生存时间年数,CENSOR状态)~ 1,数据= MOWHTO并发症)
由于缺失时间n.风险n.事件生存期标准误95% CI下限10 8 55 0.127 0.0419 0.0665 95% CI上限0.243,删除了678个观察结果
因此,10年生存率为0.127,即13%。它不可能是正确的。应该是八九成以上吧?
我做错了什么?难道是生死存亡的时刻?我是否应该为未修订的病例提供日期?那应该是什么日子呢?
如有任何帮助,我们将不胜感激。
1条答案
按热度按时间vq8itlhq1#
空列需要有一个日期在它或当你采取的difftime,结果将是空的。
举例来说:
字符串
返回NA
填写未发生关注事件时观察受试者的最后日期。这样,当你减去日期时,你会得到一个时间值。
在输出的这一部分中有一个提示:
由于缺失时间,删除了678个观察结果
他们删除了那些记录,因为它们没有时间价值。