當前位置: 首頁>>代碼示例 >>用法及示例精選 >>正文


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。非經特殊聲明,原始代碼版權歸原作者所有,本譯文未經允許或授權,請勿轉載或複製。