当前位置: 首页>>代码示例 >>用法及示例精选 >>正文


R transplant 肝移植等候名单


R语言 transplant 位于 survival 包(package)。

说明

1990-1999 年肝移植等待名单上的受试者及其处置:接受移植、等待期间死亡、从名单中退出或被审查。

用法

transplant
data(transplant, package="survival")

格式

包含以下 6 个变量的 815 个(移植)观察值的 DataFrame 。

age

加入候补名单时的年龄

sex

mf

abo

血型:ABABO

year

他们进入等候名单的年份

futime

从进入到最终处理的时间

event

最终处置:censoreddeathltxwithdraw

细节

这代表了特定地区在一段时期内的移植经历,在这段时期内,肝移植被广泛认为是一种可行的治疗方式。在此期间,肝移植的数量有所增加,但加入肝移植等待名单的受试者数量增长得更快。这些数据解决的重要问题是等待时间的变化、谁在等待,以及名单上的死亡人数是否会随之增加。

血型是一个重要的考虑因子。 O 型血受试者的供体肝脏可供 A、B、AB 或 0 型血患者使用,而 AB 型肝脏只能供 AB 型受血者使用。因此,等待名单上的 O 型受试者处于不利地位,因为 O 型供体肝脏的竞争者群体更大。

这些数据具有历史意义,并提供了竞争风险的有用示例,但与当前实践关系不大。肝脏分配政策已经发展,现在直接取决于每个患者的风险和需求,当患者在等待名单上时,其评估会定期更新。然而,总体器官短缺仍然严重。

transplant数据集是分析早期使用的版本,transplant2有一些补充和更正,是最终的数据集,与论文相符。

例子

#since event is a factor, survfit creates competing risk curves
pfit <- survfit(Surv(futime, event) ~ abo, transplant)
pfit[,2]  #time to liver transplant, by blood type
plot(pfit[,2], mark.time=FALSE, col=1:4, lwd=2, xmax=735,
       xscale=30.5, xlab="Months", ylab="Fraction transplanted",
       xaxt = 'n')
temp <- c(0, 6, 12, 18, 24)
axis(1, temp*30.5, temp)
legend(450, .35, levels(transplant$abo), lty=1, col=1:4, lwd=2)

# competing risks for type O
plot(pfit[4,], xscale=30.5, xmax=735, col=1:3, lwd=2)
legend(450, .4, c("Death", "Transpant", "Withdrawal"), col=1:3, lwd=2)

参考

Kim WR, Therneau TM, Benson JT, Kremers WK, Rosen CB, Gores GJ, Dickson ER. Deaths on the liver transplant waiting list: An analysis of competing risks. Hepatology 2006 Feb; 43(2):345-51.

相关用法


注:本文由纯净天空筛选整理自R-devel大神的英文原创作品 Liver transplant waiting list。非经特殊声明,原始代码版权归原作者所有,本译文未经允许或授权,请勿转载或复制。