spss第十四章生存分析

第十四章存分LifeTables調(diào)用此過(guò)程時(shí)，系統(tǒng)將采用即表分析法，完成對(duì)病例隨訪資料在任意指定時(shí)點(diǎn)的生存狀況是否是否否2是2是否694

6 激活數(shù)據(jù)管理窗口，定義變量名：隨訪月數(shù)的變量名為T(mén)IME，是否的變量名為DEATH，1012。14.1time，點(diǎn)擊timeDisplayTimeIntervals欄中定義需要顯在by框中2death，點(diǎn)擊鈕使之Status點(diǎn)擊DefineEvent...鈕彈出LifeTables:DefineEventforStatusVariable框，在Singlevalue欄中輸入1，表明death=1為發(fā)生事件者；點(diǎn)擊Continue鈕返回LifeTables框。選group，點(diǎn)擊鈕使之進(jìn)入Factor框，點(diǎn)擊DefineRange...鈕，彈出LifeTables:DefineRangeforFactorVariable框，定義分組的范圍，在Mininum框中輸入1，在Maxinum框中輸入2，點(diǎn)擊Continue鈕返回LifeTables框。圖圖14.1生存資料 點(diǎn)擊Options...鈕彈出LifeTables:Options框，在Plot欄中選Survival項(xiàng)，要求繪制生存率曲線圖；在CompareLevelsofFactor欄中選Overall項(xiàng)，要求作組間生存狀況的比較。之后點(diǎn)擊Continue鈕返回LifeTables框，再點(diǎn)擊OK鈕即完成分析。0242個(gè)月（0—484204248個(gè)月的生存狀況不再顯示，分別顯示進(jìn)入該時(shí)點(diǎn)例數(shù)（NumberEntrngthisIntrvl、從IntrvlRiskEventserminatingSurvivingEndDenstyRateSurvivingDenstyRate例，用中藥+化療的方式治療白血病患者，至8個(gè)月時(shí)，率為17.39%，生存概率為82.61%，生50%19.44個(gè)月。1600216Thissubfile 26LifeSurvival NumberNumberNumber SE SEIntrvlEntrngWdrawnExposd PropnPropn SEof bility nating atEndDensty viving Themediansurvivaltimeforthesedatais LifeTableSurvival NumberNumberNumber SE SEIntrvlEntrngWdrawnExposd PropnPropn SEof bility atEnd atEnd Themediansurvivaltimeforthesedata ComparisonofsurvivalexperienceusingtheWilcoxon(Gehan)statisticSurvivalVariable groupedOverall PctMean188.-14.2再又下降。但在治療中加用中藥，對(duì)個(gè)別患者而言，20個(gè)月后依然有一定的生存率。14.2中藥組與對(duì)照組生存率曲線的比14.2中藥組與對(duì)照組生存率曲線的比0Cum號(hào)隨訪天數(shù)是否［例14-2］25例某在不同時(shí)期經(jīng)隨機(jī)化分配號(hào)隨訪天數(shù)是否是TIME，DEATH10，TREATA1B2。激活Statistics菜單選Survival中的Kaplan-Meier...項(xiàng)，彈出Kaplan-Meier框（圖14.3。從點(diǎn)擊DefineEvent...鈕彈出Kaplan-Meier:DefineEventforStatusVariable框，在Singlevalue欄中輸入1，表明death=1為發(fā)生事件者；點(diǎn)擊Continue鈕返回Kaplan-Meier框。選treat，點(diǎn)擊Factor框。14.3Kaplan-Meier框點(diǎn)擊Save...鈕彈出Kaplan-Meier:SaveNewVariables框，選Survival項(xiàng)，要求將各觀察樣例的生存率存入原始數(shù)據(jù)庫(kù)中。點(diǎn)擊Continue鈕返回Kaplan-Meier框。點(diǎn)擊Options...鈕彈出Kaplan-Meier:Options框，在Plot欄中選Survival項(xiàng)，要求繪制生存率曲線圖。之后點(diǎn)擊Continue鈕返回LifeTables框，再點(diǎn)擊OK鈕即完成分析。先對(duì)A2SurvivalErrorRemaining—1563天。B組共13人，12人，生存1人，存活率為7.69%；平均生存時(shí)間為607天，標(biāo)準(zhǔn)誤為226，95%區(qū)間為163—1051天。8118121391481571.5064063061060Numberof ( Events:SurvivalStandard95%Confidence 1563(Limited1976 .FactorTREAT=11121314915816717618519413120110(1051 352A66B17Numberof Numberof Events:SurvivalStandard 95%ConfidenceB兩組的生存率相近，而后，A50%的生存率，BSurvival-- 14.4兩種治療方式生存率曲線比14.5用戶(hù)從中可見(jiàn)，如A組治療8天者，其8Survival-- 14.4兩種治療方式生存率曲線比14.5CoxRegressionCox（亦稱(chēng)比例風(fēng)險(xiǎn)模型。方式）對(duì)生存時(shí)間的影響。而在Cox回歸模型中，某一時(shí)點(diǎn)t，除了有一個(gè)本底風(fēng)險(xiǎn)量h0(t)外， i個(gè)影響因素可使該本底風(fēng)險(xiǎn)量h0(t)增至 h(t)=h0(t)[14-3]33Ag陽(yáng)、分組（Ag陽(yáng)性組17例，Ag組16例，同時(shí)白細(xì)胞數(shù)的影響作用。試據(jù)下表資料作CoxAgAg白細(xì)胞數(shù)（×109白細(xì)胞（×1091是2是1是3是4是3是5是3是是4否是4是否4是是7是是8是否是是是是是是是是是否是是是是激活數(shù)據(jù)管理窗口，定義變量名：生存月數(shù)為T(mén)IME，是否為DEATH，白細(xì)胞數(shù)為WBC，Ag陽(yáng)性與否為AG。TIME按原數(shù)據(jù)輸入，DEATH是的輸入1、否的輸入0，WBC亦按原數(shù)據(jù)輸入，AG陽(yáng)性的輸入1、的輸入2。14.6ag，點(diǎn)擊Covariates框。圖圖14.6Cox回歸模型分 2、Forward:Conditional：以假定參數(shù)為基礎(chǔ)作似然比概率檢驗(yàn)，向前逐步選擇自變量；3、Forward:LR：以最大局部似然為基礎(chǔ)作似然比概率檢驗(yàn)，向前逐步選擇自變量；5、Backward:Conditional：以假定參數(shù)為基礎(chǔ)作似然比概率檢驗(yàn)，向后逐步選擇自變量；6、Backward:LR：以最大局部似然為基礎(chǔ)作似然比概率檢驗(yàn)，向后逐步選擇自變量；Enter法。點(diǎn)擊Plots...鈕彈出CoxRegression:Plots框，在PoltsType欄中選Survival項(xiàng)，要求繪制生存率曲線圖，同時(shí)選HazardContinueCoxRegression點(diǎn)擊Save...鈕彈出CoxRegression:SaveNewVariables框，在Survival欄中選Function項(xiàng)，DiagnosticsHazardfunction項(xiàng)，要求將風(fēng)險(xiǎn)函數(shù)計(jì)算結(jié)果存入原數(shù)據(jù)庫(kù)；再選X*Beta項(xiàng)，要求計(jì)算各自變量與其系數(shù)的乘積并存盤(pán)。完成選擇后點(diǎn)擊Continue鈕返回CoxRegression框。點(diǎn)擊Options...鈕彈出CoxRegression:Options框，在ModelStatistics欄中選Atlaststep項(xiàng)，Displaybaselinefunction項(xiàng)，要求顯示各樣本的本底風(fēng)險(xiǎn)量。之后點(diǎn)擊Continue鈕返回CoxRegression框，再點(diǎn)擊OK鈕即完成分析。系統(tǒng)顯示共有33例樣本進(jìn)入分析，其中29例在觀察期內(nèi)，4例仍存活，存活率為12.1%。Cox回歸方程擬合結(jié)果的χ2檢驗(yàn)，χ211.773，P=0.0028AG與WBC兩自變量對(duì)生存(0.0089×WBC-狀態(tài)均有作用。得到風(fēng)險(xiǎn)量增加倍數(shù)為 Totalcases0000CaseswithmissingValidcaseswithnon-positiveCensoredcasesbeforetheearliesteventinastratumTotalcasesdroppedCasesavailableforDependent 4BeginningBlockNumber InitialLogLikelihood-2LogOverall2Change(-2LL)Previous2Previous2BR-11BeginningBlockNumber Variable(s)EnteredatStepNumber-2LogVariablesintheBeginningBlockNumber Variable(s)EnteredatStepNumber-2LogVariablesinthe ----Atmeanofcovariates--CumSEofCum1234578...CovariateMeans ...CovariateMeans 3Newvariableshavebeenadded: SurvivalHAZ1 CumulativeHazardFunction SurvivalFunctionatmeanof-- 14.7Cox曲Cox模型的風(fēng)險(xiǎn)量曲線圖，其趨勢(shì)也十分明顯，即隨時(shí)間的延長(zhǎng)，患者在生存上所經(jīng)歷的風(fēng)險(xiǎn)愈來(lái)愈大，到140個(gè)月時(shí)，大約是起初（0個(gè)月）的5HazardFunctionatmeanof543210014.8Cox曲Cumsur_1haz_1xbe_114.9HazardFunctionatmeanof543210014.8Cox曲Cum14.9Cox第十五章計(jì)圖的繪資料所反映的變化趨勢(shì)、數(shù)量多少、分布狀態(tài)和相互關(guān)系等形象直觀地來(lái)，以便于讀者的閱正以保證準(zhǔn)確無(wú)誤等特點(diǎn)，故在、報(bào)告等寫(xiě)作中有著十分重要的應(yīng)用價(jià)值。調(diào)用Graphs菜單的Bar過(guò)程，可繪制直條圖。直條圖用直條的長(zhǎng)短來(lái)表示非連續(xù)性資料（該資態(tài)為BP。RATE按原數(shù)據(jù)輸入，DISEASE按冠狀動(dòng)脈機(jī)能不全=1、猝死=2、心絞痛=3、心肌梗塞=4輸入，BP=1、臨界=2、異常=3Graphs單的Bar...過(guò)程，彈出BarChart定義選項(xiàng)框（圖15.1。在定義選項(xiàng)框的下方有315.1Summariesforgroupsofcases：以組為單位體現(xiàn)數(shù)據(jù)；Summariesofseparatevariables：以變量為單位體現(xiàn)數(shù)據(jù)；Valuesofindividualcases：以觀察樣例為單位體現(xiàn)數(shù)據(jù)。圖可選：SimpleClusteredStacked15.2在左側(cè)的變量列表rate鈕使之進(jìn)入BarsRepresentOthersnmmaryfunction的Variable框，選disease點(diǎn)擊鈕使之進(jìn)入CategoryAxis框，選bp點(diǎn)擊鈕使之進(jìn)入DefineClustersby圖圖15.2直條圖繪 點(diǎn)擊Titles...鈕，彈出Titles框，在Title欄內(nèi)輸入“血壓狀態(tài)與冠心病各臨床型ContinueDefineClusteredChart:SummariesforGroupsofCases框，再點(diǎn)擊OK鈕即完成。的1、2、3……等數(shù)字表示，故體現(xiàn)在統(tǒng)計(jì)圖中的分組條目會(huì)讓讀者感到不理解。為此，用戶(hù)可點(diǎn)左鍵，系統(tǒng)即彈出相應(yīng)的編輯框。編輯過(guò)程簡(jiǎn)便易行，用戶(hù)不妨一試。本章對(duì)此內(nèi)容的介紹從下圖為經(jīng)編輯（主要是將分組的標(biāo)目改為中文）后血壓狀態(tài)與冠心病各臨床型標(biāo)化發(fā)生率0冠脈GraphsLine激活數(shù)據(jù)管理窗口，定義變量名：心理問(wèn)題檢出率為RATE，分組為AGE，為SEX，AGE65-、75-，SEXMF。Graphs單的Line...過(guò)程，彈出LineChart義選項(xiàng)框3種線圖可選：Simple為單一線圖、MultipleDrop-line15.3rateLinesRepresentOthersnmmaryfunction項(xiàng)的Variableage點(diǎn)擊鈕使之CategoryAxissex鈕使之進(jìn)入DefineLinesby圖圖15.3線圖繪 點(diǎn)擊Titles...鈕，彈出Titles框，在Title欄內(nèi)輸入“某地男女性別心理問(wèn)題檢出率比較”，點(diǎn)擊Continue鈕返回DefineMultipleLine:SummariesforGroupsofCases框，OK 8GraphsArea有無(wú)激活數(shù)據(jù)管理窗口，定義變量名：有無(wú)的人數(shù)為NUMBER，分組為AGE，現(xiàn)狀 ，AGEN。Graphs單的Area...過(guò)程，彈出AreaChart義選項(xiàng)框2種線圖可選：Simple為簡(jiǎn)單區(qū)域圖、Stacked15.4numberAreasRepresentOthersnmmaryfunction選項(xiàng)的Variableage點(diǎn)擊鈕使之進(jìn)入CategoryAxiscontra鈕使之進(jìn)入DefineAreasby圖圖15.4區(qū)域圖繪 點(diǎn)擊Titles...鈕彈出Titles框在Title欄內(nèi)輸入“某市已婚育齡婦女狀況分析”，點(diǎn)擊Continue鈕返回DefineStackedArea:SummariesforGroupsofCases框，再點(diǎn)擊OK絕大多數(shù)，而45歲以后，無(wú)人數(shù)又開(kāi)始增加。某市已某市已婚育齡婦狀況分現(xiàn)狀0分組（歲Sum調(diào)用Graphs菜單的Pie過(guò)程，可繪制構(gòu)成圖。構(gòu)成圖也稱(chēng)餡餅圖，用一個(gè)圓來(lái)表現(xiàn)百分構(gòu)成，百分構(gòu)成DATA，TEXT，TEXT為字符變量。DATATEXT1、2、3、4。選Graphs菜單的Pie...過(guò)PieChart定義選項(xiàng)框構(gòu)成圖僅有一故直接點(diǎn)擊Define15.5，dataSlicesRepresentOthersnmmaryfunctionVariabletext點(diǎn)擊鈕使之進(jìn)入DefineSlicesby框。點(diǎn)擊Titles...鈕，彈出Titles框，在欄內(nèi)輸入“中草藥治療慢性支氣管炎效果構(gòu)成圖”ContinueDefinePie:SummariesforGroupsofCases框，再點(diǎn)擊OK鈕即完成。圖圖15.5構(gòu)成圖繪 中草藥治療慢性支氣管炎效果構(gòu)成圖中草藥治療慢性支氣管炎效果構(gòu)成圖無(wú)效控制好轉(zhuǎn)33顯效區(qū)域，如一組測(cè)定值的范圍（最小值—最大值、95%區(qū)間值（低限、X±1.96·SD（低

-

Bq·L湖GROUP，1、2、3。SimpleHigh-Low-Close：簡(jiǎn)單線型高低區(qū)域圖；ClusteredHigh-Low-Close：復(fù)式線型高低區(qū)域圖；SimpleRangeBar：簡(jiǎn)單直條型高低區(qū)域圖；ClusteredRangeBar：復(fù)式直條型高低區(qū)域圖；DifferenceLine：差異線區(qū)域圖。Define鈕，彈出Define15.6，點(diǎn)擊BarsRepresentOthersnmmaryfunctionVariablecatCategoryAxisgroupDefineHigh-Low-Closeby點(diǎn)擊Continue鈕返回DefineSimpleHigh-Low-Close:SummariesforGroupsofCases框，OK圖圖15.6高低區(qū)域圖繪 某市測(cè)某市測(cè)定不同水體放射性元素鍶的含量圖例高值低值均值自來(lái) 湖 水庫(kù)水源點(diǎn)百分構(gòu)成春夏秋冬DATA，輸入具體數(shù)字；季節(jié)的變量GraphsPareto...ParetoChart2Simple為單一直條構(gòu)成線圖，Stacked為堆積式直條構(gòu)成線圖，本例選用單一直條構(gòu)成線圖。然后點(diǎn)擊Define鈕，彈出DefineSimplePareto:SummariesforGroupsofCases框（圖15.7dataSumsofvariableseason鈕使之進(jìn)入CategoryAxis框。點(diǎn)擊Titles...鈕，彈出Titles框，在Title欄內(nèi)輸入“1000名20-25歲男性各季節(jié)感冒發(fā)生人數(shù)分析”ContinueDefineSimplePareto:SummariesforGroupsofCases框，再點(diǎn)擊OK鈕即完成。圖圖15.7直條構(gòu)成線圖繪 10%。1000名1000名20-25各季節(jié)感冒發(fā)生人數(shù)分00季節(jié)調(diào)用Graphs菜單的ontrol1234515.815.8X-Bar,R,s：均數(shù)控制圖和極差（標(biāo)準(zhǔn)差）X圖，用于控制重復(fù)測(cè)Rs本例選用X-Bar,R,s選項(xiàng)框的下方為數(shù)據(jù)類(lèi)型選擇（DataOrganizationCasesareunits表示數(shù)據(jù)文件中各觀察樣例只是一個(gè)值，其分組需要再定義；Casesaresubgroups表示數(shù)據(jù)文件15.9，選data鈕使之進(jìn)入ProcessMeasurementgroup鈕使之進(jìn)入SubgroupsDefinedbyChartsX-Barandrange圖。點(diǎn)擊Titles...鈕，彈出Titles框，在Title欄內(nèi)輸入“樣品測(cè)定的質(zhì)量控制圖”，點(diǎn)擊Continue鈕返回X-Bar,R,s:CasesAreUnits框，再點(diǎn)擊OK鈕即完成。圖圖15.9質(zhì)量控制圖定 樣品測(cè)定的質(zhì)量控樣品測(cè)定的質(zhì)量控制測(cè)定值均數(shù)上控制限=平均中心線=下控制限= Sigmalevel:樣品測(cè)樣品測(cè)定的質(zhì)量控制極差值上控制限=平均中心線下控制限=. Sigmalevel:0[例15-8]研究甲基對(duì)肝臟脂質(zhì)過(guò)氧化的毒性作用，選用25只大白鼠，隨機(jī)分成五組，按不LPO（nmol/L，資料如下表，試?yán)L制箱圖。12345GROUP1、2、3、4、5。Graphs單的Boxplot...過(guò)程，彈出BoxplotChart定義選2圖可選：Simple為簡(jiǎn)單箱圖，Clustered為復(fù)式箱圖，本例選用簡(jiǎn)單箱圖。然后點(diǎn)擊Define鈕，彈出DefineSimple15.10，VariablegroupCategoryAxisOK圖圖15.10箱圖定 不同染毒劑量組大白鼠LPO含量的0-20染毒劑量組（mg/kg體重DOSE(n不同染毒劑量組大白鼠LPO含量的0-20染毒劑量組（mg/kg體重DOSE(n調(diào)用Graphs菜單的ErrorBar過(guò)程，可繪制均值相關(guān)區(qū)間圖。正態(tài)分布資料的描述性指標(biāo)：如均值、標(biāo)準(zhǔn)差、標(biāo)準(zhǔn)誤，并由此求得的參照值范圍、總體均值的區(qū)間等，都可用均值相關(guān)區(qū)間15-9]食品中微量砷（As）DDC-Ag比較兩種方法測(cè)定不同濃度As標(biāo)準(zhǔn)液（μg/5ml）的光密度值區(qū)間，試?yán)L制均值相關(guān)區(qū)間圖。DDC-Ag12345激活數(shù)據(jù)管理窗口，定義變量名：數(shù)據(jù)的變量名為DATA，將新銀鹽法和DDC-Ag法的測(cè)定所得CATE1、2。GraphsErrorBar...ErrorBar2然后點(diǎn)擊Define鈕，彈出DefineClusteredErrorBar:SummariesforGroupsofCases框15.11CategoryAxiscateDefineClustersby圖圖15.11均值相關(guān)區(qū)間圖定 BarRepresent3Confidenceintervalformean：繪出總體均值的區(qū)間，要求輸入?yún)^(qū)間的百分?jǐn)?shù)；Standarderrorofmean：繪出均值+/-兩倍標(biāo)準(zhǔn)誤的區(qū)間；本例選用Confidenceintervalformean。之后點(diǎn)擊Titles...鈕，彈出Titles框，在Title欄內(nèi)輸入“兩種方法測(cè)定微量砷的光密度95%區(qū)間比較”，點(diǎn)擊Continue鈕返回DefineClusteredErrorBar:SummariesforGroupsofCases框，再點(diǎn)擊OK鈕即完成。測(cè)定方法新銀鹽法-標(biāo)準(zhǔn)液As濃度兩種方法測(cè)定微量砷的光密度 區(qū)間比95%CIGraphsScatter123CATE1，123；再設(shè)一變CATE2，用于定義微量元素種類(lèi)，Zn1、Fe2、Mn3。15.12Simple為單層散點(diǎn)圖，Overlay為多層散點(diǎn)圖，Matrix為矩陣散點(diǎn)圖，3-D為散點(diǎn)圖，本例選用單層散點(diǎn)圖。然后點(diǎn)擊Define鈕，彈出SimpleScatterplot框（圖15.13，在左側(cè)的變量dataYAxiscate1XAxis（指定飲茶狀況為橫軸標(biāo)目）,cate2SetMarkersby（指定微量元素種類(lèi)為散點(diǎn)標(biāo)志。點(diǎn)擊Titles...鈕，彈出Titles框，在Title欄內(nèi)輸入“飲茶與微量元素的關(guān)系”，點(diǎn)擊Continue鈕返回SimpleScatterplot框，再點(diǎn)擊OK鈕即完成。15.12圖圖15.13散點(diǎn)圖繪 飲茶與微量元素的關(guān)系0微量元素- 飲茶狀況下圖橫軸標(biāo)目中1.0為多喝茶組、2.0飲茶與微量元素的關(guān)系0微量元素- 飲茶狀況調(diào)用Graphs菜單的Histogram過(guò)程，可繪制直方圖。直方圖是用直條的長(zhǎng)短來(lái)表示連續(xù)性的絕10-15-20-25-30-35-745-355-1NUMBER，將每歲患者人數(shù)資料輸入；設(shè)一變量為GROUP，用于定義組，將各分組的下限值輸入。為使頻數(shù)資料在作圖中生效，應(yīng)選Data菜單的WeightCases...命令項(xiàng)，在彈出的WeightCases框中取Weightcasesby項(xiàng)，NUMBERFrequencyVariableOK（圖15.14，在左側(cè)的變量列表中選group點(diǎn)擊鈕使之進(jìn)入Variable框；點(diǎn)擊Titles...鈕，彈出Titles框，在Title欄內(nèi)輸入“某地某年流行性乙型腦炎患者分布”，點(diǎn)擊Continue鈕返回Histogram框，再點(diǎn)擊OK鈕即完成。圖圖15.14直方圖繪 ChartCarousel（即每歲一組系統(tǒng)只按默認(rèn)的每5歲一組方式輸出圖形，所以需要安用戶(hù)的要求對(duì)統(tǒng)計(jì)圖進(jìn)行編輯。點(diǎn)擊ChartCarousel窗口上端中的Edit鈕，彈出!Chart1窗口，將鼠標(biāo)箭頭指向圖的橫軸并雙擊鼠標(biāo)左15.15，IntervalsCustomDefine..IntervalAxis:DefineCustomInterval框（15.16DefinitionIntervalwidth1，要求按每歲一組的方式作圖；在Range欄的Minimum處輸入0在um處輸入65要求橫軸從0歲開(kāi)始至65歲止點(diǎn)擊Continue鈕返回IntervalAxis框，再點(diǎn)擊OK鈕即可。圖圖15.15橫軸修 15.16某地某年流行性乙型某地某年流行性乙型腦炎患者分布Std.Dev=4.06Mean=5.6 N=分GraphsNormalP-P[例15-12]某醫(yī)師測(cè)得30名健康女大學(xué)生總蛋白含量（g/L）資料如下，試?yán)L制正態(tài)概率選Graphs菜單的NormalP-P...過(guò)程，彈出NormalP-PPlot框（圖15.17，在左側(cè)的量列表中選data點(diǎn)擊鈕使之進(jìn)入Variable框。在Transform欄中，系統(tǒng)