加載中........
×

您的生存曲線領導還滿意嗎

2019-10-8 作者:王曉曉,趙一鳴   來源:臨床流行病學和循證醫學 我要評論2
Tags: 生存曲線  

小伙伴抱怨道:我覺得我畫的生存曲線很好啊,為什么領導說我畫的很丑呢?這個,可能是因為審美不同吧......那領導到底喜歡什么樣的呢?可能像下圖這樣。



上圖是采用R軟件繪制的,繪圖之前需要安裝兩個程序包survminer和survival。在開始繪圖之前,需要通過library函數加載上述程序包。

library(survminer);library(survival)

在讀取數據后,需要通過survfit確定生存分析的對象。因為小編的程序和數據在同一文件下,所以讀取數據的程序為 surv<-read.csv("surv.csv")。

接下來,通過survfit確定生存分析的對象,其中time為生存時間,status為生存狀態。fit<- survfit(Surv(time, status) ~ 1,data = surv)。萬一,領導想讓我分組繪制生存曲線呢,比如按性別,只需要改為fit<- survfit(Surv(time, status) ~ gender, data =surv)。



還沒結束哦,各位領導們會不會更喜歡下圖?相對上圖,該圖增加了1年、5年和10年的生存率,提供的信息也更多。

首先,我們需要知道1年、5年和10年的生存率,通過summary(fit)可得到不同時間點的生存率,但因為系統默認設置的問題,往往無法得到完整的結果。簡單設置一下就搞定了,options(max.print = 999999999),然后再調用summary(fit)。

結果大致如下所示,大多時候,我們會發現結果中并沒有1年、5年、10年這3個整點的時間點,那是因為在整點時并未出現死亡。這時候找小于上述時間點且離上述時間點最近的即可。本示例中,1年、5年和10年的生存率分別為88.5%、57.0%和39.0%。

在上面生存曲線的基礎上,我們標注3個坐標點(1,88.5)、(5,57)、(10,39),并通過文本說明其生存率即可。



小提示:78%用戶已下載梅斯醫學APP,更方便閱讀和交流,請掃描二維碼直接下載APP

只有APP中用戶,且經認證才能發表評論!馬上下載

177BCyz****(暫無昵稱)

漲知識了

(來自:梅斯醫學APP)

2019-10-10 18:53:13 回復

qingfengqishi5

學習了,學習了

(來自:梅斯醫學APP)

2019-10-8 23:24:46 回復

web對話
500彩票网开奖结果