周磊

摘要：在醫(yī)學(xué)科研中，生存分析是較為常用的統(tǒng)計分析方法，對患者生存期進(jìn)行研究，收集患者結(jié)局、發(fā)生時間資料，進(jìn)而全面評定患者。同時，在臨床研究中，患者結(jié)局多呈現(xiàn)多樣性，且不同結(jié)局之間具有一定聯(lián)系，因此需要利用多結(jié)局生存分析模型分析數(shù)據(jù)，提高分析準(zhǔn)確性。本文主要對的多結(jié)局風(fēng)險比例模型方法進(jìn)行概述，并在此基礎(chǔ)上探究該模型在醫(yī)學(xué)生存分析中的具體應(yīng)用。

關(guān)鍵詞：風(fēng)險比例模型 多結(jié)局 生存分析 應(yīng)用

中圖分類號：G642 文獻(xiàn)標(biāo)識碼： A 文章編號：1672-1578（2015）02-0082-01

對兩個變量對象的結(jié)局情況進(jìn)行觀察是生存分析資料的主要特征， 即刻畫觀察對象預(yù)期結(jié)局發(fā)生的時間、結(jié)果。目前，醫(yī)學(xué)研究對生存分析廣泛應(yīng)用，且常用方法為Cox比例風(fēng)險模型、壽命表法、Kaplan-Meier法等[1]。其中，Cox比例風(fēng)險模型最常用。但是這些的方法均為單結(jié)局分析法，僅能對一種結(jié)局的是否發(fā)生進(jìn)行分析，這樣，研究應(yīng)用范圍越大，方法使用不足性越凸顯，不利于多個失效時間情況的處理，需要充分考慮多結(jié)局現(xiàn)象，如某醫(yī)院研究異體骨髓移植患者手術(shù)后感染受排異性、遺傳性的影響，觀察300名患者，患者出現(xiàn)病毒感染、細(xì)菌感染、真菌感染三種不同結(jié)局，多結(jié)局風(fēng)險比例模型的重要性凸顯。

1 多結(jié)局風(fēng)險比例模型應(yīng)用方法

通常情況下，利用風(fēng)險函數(shù)描述、刻畫生存分析，應(yīng)用多結(jié)局生存函數(shù)描述多結(jié)局生存狀態(tài)。因此，生存函數(shù)、風(fēng)險函數(shù)相互對應(yīng)，即可通過生存函數(shù)確定風(fēng)險函數(shù)。例如，多結(jié)局生存模型有兩種，在這兩種基礎(chǔ)上對不同模型進(jìn)行定義、限制，得到不同函數(shù)關(guān)系。當(dāng)觀察對象可能存在s個結(jié)局總數(shù)時，用X表示協(xié)變量，則得到表達(dá)式為X=（xk1，xk2，……xkm），而h（t，Xk）表示風(fēng)險函數(shù)，表達(dá)式為h（t，Xk）=h0k（t）exp（β1xk1+β2xk2+……βmxkm），其中，k=1，2，……m，而歸為系數(shù)為β=（β1，β2，……βm）[2]。同時，當(dāng)準(zhǔn)基風(fēng)險函數(shù)h0（t）= h0k（t），k=1，2，……s。此時，h（t，Xk）函數(shù)表示為h（t，Xk）=h0（t）exp（β1xk1+β2xk2+……βmxkm）。這樣，當(dāng)某因素xki增加一個單位時，HR為不同結(jié)局的風(fēng)險函數(shù)比，且HR=eβi，但該模型存在一些不符合實際研究的問題，需要在實際應(yīng)用中驗證。此外，在實際研究中，根據(jù)不同情況，多結(jié)局生存時間具有不同定義，表現(xiàn)為有序多結(jié)局、無序多結(jié)局、有序多結(jié)局三種情況，其中，有序多結(jié)局指的是結(jié)局發(fā)生存在先后順序，只有發(fā)生第一種結(jié)局，才會接著產(chǎn)生第二種結(jié)局。而無序多結(jié)局即結(jié)局發(fā)生時間順序無法確定，有序多結(jié)局則為各結(jié)局均先后發(fā)生，與其他結(jié)局之間存在一定聯(lián)系[3]。

2 多結(jié)局風(fēng)險比例模型在生存分析中的具體應(yīng)用

2.1在有序多結(jié)局中的應(yīng)用

某醫(yī)院研究出治療乳腺癌的新藥，為了探究該藥的有效性，選取該院在2009年3月至2013年3月收治的200例乳腺癌術(shù)后患者，同時將所有患者隨機(jī)分為100例對照組和100例觀察組，給予對照組藥物治療，給予觀察組新藥治療，并回顧分析手術(shù)方式、臨床資料，統(tǒng)計記錄不同手術(shù)對乳腺癌死亡的影響。在此基礎(chǔ)上，對兩組患者進(jìn)行隨訪觀察，統(tǒng)計觀察時間所有患者出現(xiàn)的乳腺癌復(fù)發(fā)、死亡結(jié)局。同時，在此研究中，手術(shù)方式（o），治療藥物（g），放療（r）為影響研究結(jié)果的因素，乳房切除術(shù)（o=1）、腫瘤切除術(shù)（o=2）、新藥物（g=1）、常規(guī)藥物（g=2）、放療（r=1）、未放療（r=2）六種單獨因素造成不同結(jié)局[4]。

針對該種情況，如果乳腺癌復(fù)發(fā)為本次研究結(jié)局，則風(fēng)險函數(shù)為h1（t，X）= h01（t）exp（β11g+β21o+β31r），而如果死亡為本次研究結(jié)局，則h2（t，X）= h02（t）exp（β12g+β22o+β32r），兩種結(jié)局函數(shù)函數(shù)形式相同，回歸系數(shù)不同。同時，兩種結(jié)局為有序結(jié)局，結(jié)合兩個函數(shù)，得多多結(jié)局比例風(fēng)險模型hk（t，X）=h0k（t）exp（β11g1+β21o1+β31r1+β12g2+β22o2+β32r2），利用SAS、Stata軟件，以此建立合理模型，得到模型結(jié)果，結(jié)果如表1所示。

根據(jù)g1回歸系數(shù)，即新藥治療回歸系數(shù)可知，服用新藥患者的中路復(fù)發(fā)、死亡風(fēng)險明顯低于常規(guī)治療患者，驗證新藥的有效性。

2.2在無序多結(jié)局中的應(yīng)用

例如在慢性乙肝疾病結(jié)局因素分析中，因素多樣、無序。在具體研究中，選取我院在2008年1月至2013年1月確診的300例慢性乙肝患者，研究體內(nèi)白蛋白含量、肝纖維化程度對患者終末型肝病、代償性肝硬化的影響。由于兩種結(jié)局完全不同，風(fēng)險函數(shù)也不同。在時間變量方面，由于所有研究的慢性乙肝患者均具有發(fā)生這兩種不同結(jié)局的可能，因此時間變了按總時間定義。在此基礎(chǔ)上，收集分析患者臨床資料，統(tǒng)計分析患者不同性別（x1）、年齡（x2）、體內(nèi)白蛋白含量（x3）、肝纖維化（x4）比例，建立分析模型，獲取HR、P值，分析x3、x4風(fēng)險比，進(jìn)而得出體內(nèi)白蛋白含量、肝纖維化與兩種結(jié)局的關(guān)系，嚴(yán)重兩種因素具有造成兩種結(jié)局的較大可能性，風(fēng)險較高。

3 結(jié)語

在醫(yī)學(xué)生存分析中，合理利用多結(jié)局風(fēng)險比例模型，可有效提高分析結(jié)果準(zhǔn)確性，幫助研究者獲取準(zhǔn)確數(shù)據(jù)，得到科學(xué)結(jié)局，驗證研究，以此得到科學(xué)的治療方法、試驗新藥，為醫(yī)學(xué)水平提高提供保障。

參考文獻(xiàn)：

[1]張欣峰，婁清濤，陸建邦等.現(xiàn)時生存分析方法的應(yīng)用實踐與評價[J].中國衛(wèi)生統(tǒng)計，2011，28（1）：26-28.

[2]李文琦.神經(jīng)網(wǎng)絡(luò)模型和Cox比例風(fēng)險模型在生存分析中應(yīng)用的比較[D].徐州醫(yī)學(xué)院，2011.

[3]沈玲，朱海超，黃建鋒等.TACE治療中晚期原發(fā)性肝癌108例生存分析[J].腫瘤防治研究，2011，38（10）：1163-1166.

[4]姚珊珊，劉慧敏，劉艷民等.乙型肝炎相關(guān)原發(fā)性肝癌患者的生存分析[J].中華實驗和臨床感染病雜志（電子版），2013，7（1）：69-74.