【原】一文掌握Kaplan Meier生存分析

Mark_killua 2020-07-09

展開全文

Kaplan Meier，是一種單因素生存分析。它可用于研究1個(gè)因素對(duì)于生存時(shí)間的影響，在醫(yī)療領(lǐng)域中使用廣泛。

一、案例說明

當(dāng)前某研究人員擬觀察一種新型癌癥藥物的療效情況，首先將100名癌癥患者隨機(jī)分成兩組，對(duì)照組使用傳統(tǒng)治療方式，實(shí)驗(yàn)組使用新式藥物治療方式。并且隨訪時(shí)間為2年。并且以‘是否死亡’為作為結(jié)局。希望通過研究了解到新式藥物是否對(duì)于生存時(shí)間帶來(lái)影響。

藥物組別：0代表傳統(tǒng)治療組；1代表新藥組。

生存狀態(tài)：0代表存活或失訪，1代表死亡。

生存時(shí)間：定量數(shù)據(jù)。

二、操作步驟

將數(shù)據(jù)上傳至SPSSAU平臺(tái)，選擇【實(shí)驗(yàn)研究】--【Kaplan Meier】。

SPSSAU默認(rèn)生存狀態(tài)中，用數(shù)字1代表死亡，0代表生存。

如果數(shù)據(jù)不是這樣設(shè)置，可以通過【數(shù)據(jù)處理】--【數(shù)據(jù)編碼】修改。

[數(shù)據(jù)處理]--[數(shù)據(jù)編碼]

將[藥物組別]放入【X(定類)】，[生存時(shí)間]放入【Y1生存時(shí)間】，[生存狀態(tài)]放入【Y2生存狀態(tài)】。點(diǎn)擊開始分析。

Kaplan Meier

三、結(jié)果分析

1 模型基本描述

模型基本描述

上表是簡(jiǎn)單描述數(shù)據(jù)基本情況。其中新藥、傳統(tǒng)治療分別有50個(gè)樣本。傳統(tǒng)治療組死亡43人；新式藥物組死亡44人。傳統(tǒng)治療組中位數(shù)為9周，新藥組中位數(shù)為65周，整體生存時(shí)間中位數(shù)為24.5周。

特別說明一點(diǎn)，通常情況下感興趣的結(jié)局事情比例（事件占比）一般希望大于20%。

2 生存時(shí)間估計(jì)

上表格為Kaplan Meier模型得到的生存時(shí)間中位數(shù)估計(jì)值及其置信區(qū)間。由于生存時(shí)間的波動(dòng)性，一般使用中位數(shù)表示整體情況（而不是使用平均值）。

從上表可知，整體研究病例生存時(shí)間估計(jì)值是36周，傳統(tǒng)治療的生存時(shí)間中位數(shù)是9周，但 ‘新式藥物’時(shí)中位數(shù)生存時(shí)間為69周，明顯高于傳統(tǒng)治療時(shí)的生存時(shí)間。

注：中位數(shù)是與表中小于或等于 0.50 的第一個(gè)生存概率相關(guān)聯(lián)的時(shí)間。

3 假設(shè)檢驗(yàn)

此處模型檢驗(yàn)的原定假設(shè)為：是否放入藥物組別兩種情況時(shí)模型質(zhì)量均一樣。

從上表可知，模型拒絕原定假設(shè)(χ2=10.691,p =0.001 < 0.05)，即說明本次模型構(gòu)建有意義。

配對(duì)Log Rank檢驗(yàn)用于分析兩兩組別生存時(shí)間估計(jì)值差異性。

從上表可知，模型拒絕原假設(shè)，P<0.05，說明兩組之間生存時(shí)間存在顯著性差異。

4 生存曲線

上圖可以明顯的看出，‘新式藥物’的生存率會(huì)明顯的高于‘傳統(tǒng)治療’方案。而且‘新式藥物’時(shí)，0~80周時(shí)生存率下降的幅度相對(duì)較小，但在80周后生存率下降幅度加大?！畟鹘y(tǒng)治療’方案來(lái)看，0~20周生存率下降非常明顯，20周后生存率下降幅度相對(duì)較為穩(wěn)定。無(wú)論如何，新式藥物的生存率都會(huì)明顯的高于傳統(tǒng)治療方案，整體說明，‘新式藥物’起著明顯的作用。