周磊
摘要:在醫(yī)學科研中,生存分析是較為常用的統(tǒng)計分析方法,對患者生存期進行研究,收集患者結(jié)局、發(fā)生時間資料,進而全面評定患者。同時,在臨床研究中,患者結(jié)局多呈現(xiàn)多樣性,且不同結(jié)局之間具有一定聯(lián)系,因此需要利用多結(jié)局生存分析模型分析數(shù)據(jù),提高分析準確性。本文主要對的多結(jié)局風險比例模型方法進行概述,并在此基礎(chǔ)上探究該模型在醫(yī)學生存分析中的具體應用。
關(guān)鍵詞:風險比例模型 多結(jié)局 生存分析 應用
中圖分類號:G642 文獻標識碼: A 文章編號:1672-1578(2015)02-0082-01
對兩個變量對象的結(jié)局情況進行觀察是生存分析資料的主要特征, 即刻畫觀察對象預期結(jié)局發(fā)生的時間、結(jié)果。目前,醫(yī)學研究對生存分析廣泛應用,且常用方法為Cox比例風險模型、壽命表法、Kaplan-Meier法等[1]。其中,Cox比例風險模型最常用。但是這些的方法均為單結(jié)局分析法,僅能對一種結(jié)局的是否發(fā)生進行分析,這樣,研究應用范圍越大,方法使用不足性越凸顯,不利于多個失效時間情況的處理,需要充分考慮多結(jié)局現(xiàn)象,如某醫(yī)院研究異體骨髓移植患者手術(shù)后感染受排異性、遺傳性的影響,觀察300名患者,患者出現(xiàn)病毒感染、細菌感染、真菌感染三種不同結(jié)局,多結(jié)局風險比例模型的重要性凸顯。
1 多結(jié)局風險比例模型應用方法
通常情況下,利用風險函數(shù)描述、刻畫生存分析,應用多結(jié)局生存函數(shù)描述多結(jié)局生存狀態(tài)。因此,生存函數(shù)、風險函數(shù)相互對應,即可通過生存函數(shù)確定風險函數(shù)。例如,多結(jié)局生存模型有兩種,在這兩種基礎(chǔ)上對不同模型進行定義、限制,得到不同函數(shù)關(guān)系。當觀察對象可能存在s個結(jié)局總數(shù)時,用X表示協(xié)變量,則得到表達式為X=(xk1,xk2,……xkm),而h(t,Xk)表示風險函數(shù),表達式為h(t,Xk)=h0k(t)exp(β1xk1+β2xk2+……βmxkm),其中,k=1,2,……m,而歸為系數(shù)為β=(β1,β2,……βm)[2]。同時,當準基風險函數(shù)h0(t)= h0k(t),k=1,2,……s。此時,h(t,Xk)函數(shù)表示為h(t,Xk)=h0(t)exp(β1xk1+β2xk2+……βmxkm)。這樣,當某因素xki增加一個單位時,HR為不同結(jié)局的風險函數(shù)比,且HR=eβi,但該模型存在一些不符合實際研究的問題,需要在實際應用中驗證。此外,在實際研究中,根據(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é)局風險比例模型在生存分析中的具體應用
2.1在有序多結(jié)局中的應用
某醫(yī)院研究出治療乳腺癌的新藥,為了探究該藥的有效性,選取該院在2009年3月至2013年3月收治的200例乳腺癌術(shù)后患者,同時將所有患者隨機分為100例對照組和100例觀察組,給予對照組藥物治療,給予觀察組新藥治療,并回顧分析手術(shù)方式、臨床資料,統(tǒng)計記錄不同手術(shù)對乳腺癌死亡的影響。在此基礎(chǔ)上,對兩組患者進行隨訪觀察,統(tǒng)計觀察時間所有患者出現(xiàn)的乳腺癌復發(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ā)為本次研究結(jié)局,則風險函數(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é)局比例風險模型hk(t,X)=h0k(t)exp(β11g1+β21o1+β31r1+β12g2+β22o2+β32r2),利用SAS、Stata軟件,以此建立合理模型,得到模型結(jié)果,結(jié)果如表1所示。
根據(jù)g1回歸系數(shù),即新藥治療回歸系數(shù)可知,服用新藥患者的中路復發(fā)、死亡風險明顯低于常規(guī)治療患者,驗證新藥的有效性。
2.2在無序多結(jié)局中的應用
例如在慢性乙肝疾病結(jié)局因素分析中,因素多樣、無序。在具體研究中,選取我院在2008年1月至2013年1月確診的300例慢性乙肝患者,研究體內(nèi)白蛋白含量、肝纖維化程度對患者終末型肝病、代償性肝硬化的影響。由于兩種結(jié)局完全不同,風險函數(shù)也不同。在時間變量方面,由于所有研究的慢性乙肝患者均具有發(fā)生這兩種不同結(jié)局的可能,因此時間變了按總時間定義。在此基礎(chǔ)上,收集分析患者臨床資料,統(tǒng)計分析患者不同性別(x1)、年齡(x2)、體內(nèi)白蛋白含量(x3)、肝纖維化(x4)比例,建立分析模型,獲取HR、P值,分析x3、x4風險比,進而得出體內(nèi)白蛋白含量、肝纖維化與兩種結(jié)局的關(guān)系,嚴重兩種因素具有造成兩種結(jié)局的較大可能性,風險較高。
3 結(jié)語
在醫(yī)學生存分析中,合理利用多結(jié)局風險比例模型,可有效提高分析結(jié)果準確性,幫助研究者獲取準確數(shù)據(jù),得到科學結(jié)局,驗證研究,以此得到科學的治療方法、試驗新藥,為醫(yī)學水平提高提供保障。
參考文獻:
[1]張欣峰,婁清濤,陸建邦等.現(xiàn)時生存分析方法的應用實踐與評價[J].中國衛(wèi)生統(tǒng)計,2011,28(1):26-28.
[2]李文琦.神經(jīng)網(wǎng)絡(luò)模型和Cox比例風險模型在生存分析中應用的比較[D].徐州醫(yī)學院,2011.
[3]沈玲,朱海超,黃建鋒等.TACE治療中晚期原發(fā)性肝癌108例生存分析[J].腫瘤防治研究,2011,38(10):1163-1166.
[4]姚珊珊,劉慧敏,劉艷民等.乙型肝炎相關(guān)原發(fā)性肝癌患者的生存分析[J].中華實驗和臨床感染病雜志(電子版),2013,7(1):69-74.