第十四章生存分析的sas實(shí)現(xiàn)【整理版】(參考版)
2025-06-19 17:37本頁(yè)面
【正文】 事實(shí)上，Logrank統(tǒng)計(jì)檢驗(yàn)方法要求資料的兩條或多條生存曲線不交叉，所以出現(xiàn)上述統(tǒng)計(jì)結(jié)果也就不奇怪了。教材中的說(shuō)明用KaplanMeier法計(jì)算兩組的生存率并作生存率曲線圖（如圖14－6所示），從兩組的生存率曲線可以發(fā)現(xiàn)兩條樣本估計(jì)的生存率曲線大約在40天處交叉，在40天之前保守組的患者生存率高于手術(shù)組，而40天之后則手術(shù)組生存率高于保守組，并且兩條生存曲線的差異是比較明顯的，這兩條交叉生存曲線提示保守治療有可能不會(huì)很快死亡，但手術(shù)不太成功者很快死亡，而手術(shù)成功能夠存活35天以上的對(duì)象可能生存期會(huì)較長(zhǎng)。由結(jié)果可以發(fā)現(xiàn)檢驗(yàn)的統(tǒng)計(jì)量為 ，p值= ，所以在a=，不能拒絕H0，不能認(rèn)為兩層(c=1,2)之間生存率存在差異。 指定分層變量run。 指定時(shí)間變量和截尾變量并指出數(shù)據(jù)截尾時(shí)截尾變量的取值freq f。proc lifetest plots=(s)。cards。end。 輸入時(shí)間變量和頻數(shù)變量output。 定義截尾變量do i=1 to n。 定義分層變量input n。表14－6 204名晚期胰腺癌患者的生存資料（生存時(shí)間的單位為天）保守治療組生存時(shí)間23581013141521354045556175死亡人數(shù)1436044041268440失訪人數(shù)000230010028004手術(shù)治療組生存時(shí)間235678101213154663656775死亡人數(shù)11249444414348884失訪人數(shù)0000001300008012SAS分析程序data ch14_5。例14－5 某醫(yī)生研究手術(shù)是否可以延長(zhǎng)晚期胰腺癌患者的生存期。圖14－5 3期和4期淋巴瘤的病人生存曲線 Logrank檢驗(yàn)是一種非參數(shù)的檢驗(yàn)方法，對(duì)資料的分布基本沒(méi)有要求，但要求每組均含有失效事件發(fā)生的觀察對(duì)象資料，可以證明：Logrank檢驗(yàn)的檢驗(yàn)效能僅與發(fā)生失效事件的人數(shù)有關(guān)。教材中的說(shuō)明由于檢驗(yàn)統(tǒng)計(jì)量 , P , 所以在a=，拒絕H0，接受H1。可認(rèn)為兩個(gè)(c=1,2)總體生存率之間存在差異。SAS軟件輸出結(jié)果解釋The LIFETEST Procedure Testing Homogeneity of Survival Curves for time over Strata Rank Statistics c LogRank Wilcoxon 1 2 Covariance Matrix for the LogRank Statistics c 1 2 1 2 Covariance Matrix for the Wilcoxon Statistics c 1 2 1 2 Test of Equality over Strata Pr Test ChiSquare DF ChiSquare LogRank 1 Wilcoxon 1 2Log(LR) 1 SAS軟件輸出結(jié)果解釋這里只給出了兩個(gè)(c=1,2)總體生存率之間是否有差異的檢驗(yàn)結(jié)果，其他的輸出結(jié)果與前面類(lèi)似，不再列出了。 指定時(shí)間變量和截尾變量strata c。proc lifetest plots=(s)。cards。end。input time censor 。 定義分組變量input n。SAS分析程序data ch14_4。例14－4 在例14－3中介紹了McKelvey et al收集了19名3期淋巴瘤病人的生存資料，該研究者還同時(shí)收集了4期淋巴瘤病人的生存數(shù)據(jù)共61例：4，6，10，11，11，11，13，17，20，20，21，22，24，24，29，30，30，31，33，34，35，39，40，41+，43+，45，46，50，56，61+，61+，63，68，82，85，88，89，90，93，104，110，134，137，160+，169，171，173，175，184，201，222，235+，247+，260+，284+，290+，291+，302+，304+，341+，345+。③．壽命表法沒(méi)有確切的死亡時(shí)間和失訪時(shí)間，假定在每個(gè)時(shí)間段中的“死亡發(fā)生時(shí)間”和“失訪發(fā)生的時(shí)間”分別呈均勻分布，因此在每個(gè)時(shí)間段中的生存率呈線性變化，故各個(gè)組段間的生存率用直線相連接，Kaplan
高考資料相關(guān)推薦
鄂ICP備17016276號(hào)-1