spss第十四章生存分析

第十四章存分LifeTables调用此过程时，系统将采用即表分析法，完成对病例随访资料在任意指定时点的生存状况是否是否否2是2是否694

6 激活数据管理窗口，定义变量名：随访月数的变量名为TIME，是否的变量名为DEATH，1012。14.1time，点击timeDisplayTimeIntervals栏中定义需要显在by框中2death，点击钮使之Status点击DefineEvent...钮弹出LifeTables:DefineEventforStatusVariable框，在Singlevalue栏中输入1，表明death=1为发生事件者；点击Continue钮返回LifeTables框。选group，点击钮使之进入Factor框，点击DefineRange...钮，弹出LifeTables:DefineRangeforFactorVariable框，定义分组的范围，在Mininum框中输入1，在Maxinum框中输入2，点击Continue钮返回LifeTables框。图图14.1生存资料 点击Options...钮弹出LifeTables:Options框，在Plot栏中选Survival项，要求绘制生存率曲线图；在CompareLevelsofFactor栏中选Overall项，要求作组间生存状况的比较。之后点击Continue钮返回LifeTables框，再点击OK钮即完成分析。0242个月（0—484204248个月的生存状况不再显示，分别显示进入该时点例数（NumberEntrngthisIntrvl、从IntrvlRiskEventserminatingSurvivingEndDenstyRateSurvivingDenstyRate例，用中药+化疗的方式治疗白血病患者，至8个月时，率为17.39%，生存概率为82.61%，生50%19.44个月。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再又下降。但在治疗中加用中药，对个别患者而言，20个月后依然有一定的生存率。14.2中药组与对照组生存率曲线的比14.2中药组与对照组生存率曲线的比0Cum号随访天数是否［例14-2］25例某在不同时期经随机化分配号随访天数是否是TIME，DEATH10，TREATA1B2。激活Statistics菜单选Survival中的Kaplan-Meier...项，弹出Kaplan-Meier框（图14.3。从点击DefineEvent...钮弹出Kaplan-Meier:DefineEventforStatusVariable框，在Singlevalue栏中输入1，表明death=1为发生事件者；点击Continue钮返回Kaplan-Meier框。选treat，点击Factor框。14.3Kaplan-Meier框点击Save...钮弹出Kaplan-Meier:SaveNewVariables框，选Survival项，要求将各观察样例的生存率存入原始数据库中。点击Continue钮返回Kaplan-Meier框。点击Options...钮弹出Kaplan-Meier:Options框，在Plot栏中选Survival项，要求绘制生存率曲线图。之后点击Continue钮返回LifeTables框，再点击OK钮即完成分析。先对A2SurvivalErrorRemaining—1563天。B组共13人，12人，生存1人，存活率为7.69%；平均生存时间为607天，标准误为226，95%区间为163—1051天。8118121391481571.5064063061060Numberof ( Events:SurvivalStandard95%Confidence 1563(Limited1976 .FactorTREAT=11121314915816717618519413120110(1051 352A66B17Numberof Numberof Events:SurvivalStandard 95%ConfidenceB两组的生存率相近，而后，A50%的生存率，BSurvival-- 14.4两种治疗方式生存率曲线比14.5用户从中可见，如A组治疗8天者，其8Survival-- 14.4两种治疗方式生存率曲线比14.5CoxRegressionCox（亦称比例风险模型。方式）对生存时间的影响。而在Cox回归模型中，某一时点t，除了有一个本底风险量h0(t)外， i个影响因素可使该本底风险量h0(t)增至 h(t)=h0(t)[14-3]33Ag阳、分组（Ag阳性组17例，Ag组16例，同时白细胞数的影响作用。试据下表资料作CoxAgAg白细胞数（×109白细胞（×1091是2是1是3是4是3是5是3是是4否是4是否4是是7是是8是否是是是是是是是是是否是是是是激活数据管理窗口，定义变量名：生存月数为TIME，是否为DEATH，白细胞数为WBC，Ag阳性与否为AG。TIME按原数据输入，DEATH是的输入1、否的输入0，WBC亦按原数据输入，AG阳性的输入1、的输入2。14.6ag，点击Covariates框。图图14.6Cox回归模型分 2、Forward:Conditional：以假定参数为基础作似然比概率检验，向前逐步选择自变量；3、Forward:LR：以最大局部似然为基础作似然比概率检验，向前逐步选择自变量；5、Backward:Conditional：以假定参数为基础作似然比概率检验，向后逐步选择自变量；6、Backward:LR：以最大局部似然为基础作似然比概率检验，向后逐步选择自变量；Enter法。点击Plots...钮弹出CoxRegression:Plots框，在PoltsType栏中选Survival项，要求绘制生存率曲线图，同时选HazardContinueCoxRegression点击Save...钮弹出CoxRegression:SaveNewVariables框，在Survival栏中选Function项，DiagnosticsHazardfunction项，要求将风险函数计算结果存入原数据库；再选X*Beta项，要求计算各自变量与其系数的乘积并存盘。完成选择后点击Continue钮返回CoxRegression框。点击Options...钮弹出CoxRegression:Options框，在ModelStatistics栏中选Atlaststep项，Displaybaselinefunction项，要求显示各样本的本底风险量。之后点击Continue钮返回CoxRegression框，再点击OK钮即完成分析。系统显示共有33例样本进入分析，其中29例在观察期内，4例仍存活，存活率为12.1%。Cox回归方程拟合结果的χ2检验，χ211.773，P=0.0028AG与WBC两自变量对生存(0.0089×WBC-状态均有作用。得到风险量增加倍数为 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模型的风险量曲线图，其趋势也十分明显，即随时间的延长，患者在生存上所经历的风险愈来愈大，到140个月时，大约是起初（0个月）的5HazardFunctionatmeanof543210014.8Cox曲Cumsur_1haz_1xbe_114.9HazardFunctionatmeanof543210014.8Cox曲Cum14.9Cox第十五章计图的绘资料所反映的变化趋势、数量多少、分布状态和相互关系等形象直观地来，以便于读者的阅正以保证准确无误等特点，故在、报告等写作中有着十分重要的应用价值。调用Graphs菜单的Bar过程，可绘制直条图。直条图用直条的长短来表示非连续性资料（该资态为BP。RATE按原数据输入，DISEASE按冠状动脉机能不全=1、猝死=2、心绞痛=3、心肌梗塞=4输入，BP=1、临界=2、异常=3Graphs单的Bar...过程，弹出BarChart定义选项框（图15.1。在定义选项框的下方有315.1Summariesforgroupsofcases：以组为单位体现数据；Summariesofseparatevariables：以变量为单位体现数据；Valuesofindividualcases：以观察样例为单位体现数据。图可选：SimpleClusteredStacked15.2在左侧的变量列表rate钮使之进入BarsRepresentOthersnmmaryfunction的Variable框，选disease点击钮使之进入CategoryAxis框，选bp点击钮使之进入DefineClustersby图图15.2直条图绘 点击Titles...钮，弹出Titles框，在Title栏内输入“血压状态与冠心病各临床型ContinueDefineClusteredChart:SummariesforGroupsofCases框，再点击OK钮即完成。的1、2、3……等数字表示，故体现在统计图中的分组条目会让读者感到不理解。为此，用户可点左键，系统即弹出相应的编辑框。编辑过程简便易行，用户不妨一试。本章对此内容的介绍从下图为经编辑（主要是将分组的标目改为中文）后血压状态与冠心病各临床型标化发生率0冠脉GraphsLine激活数据管理窗口，定义变量名：心理问题检出率为RATE，分组为AGE，为SEX，AGE65-、75-，SEXMF。Graphs单的Line...过程，弹出LineChart义选项框3种线图可选：Simple为单一线图、MultipleDrop-line15.3rateLinesRepresentOthersnmmaryfunction项的Variableage点击钮使之CategoryAxissex钮使之进入DefineLinesby图图15.3线图绘 点击Titles...钮，弹出Titles框，在Title栏内输入“某地男女性别心理问题检出率比较”，点击Continue钮返回DefineMultipleLine:SummariesforGroupsofCases框，OK 8GraphsArea有无激活数据管理窗口，定义变量名：有无的人数为NUMBER，分组为AGE，现状 ，AGEN。Graphs单的Area...过程，弹出AreaChart义选项框2种线图可选：Simple为简单区域图、Stacked15.4numberAreasRepresentOthersnmmaryfunction选项的Variableage点击钮使之进入CategoryAxiscontra钮使之进入DefineAreasby图图15.4区域图绘 点击Titles...钮弹出Titles框在Title栏内输入“某市已婚育龄妇女状况分析”，点击Continue钮返回DefineStackedArea:SummariesforGroupsofCases框，再点击OK绝大多数，而45岁以后，无人数又开始增加。某市已某市已婚育龄妇状况分现状0分组（岁Sum调用Graphs菜单的Pie过程，可绘制构成图。构成图也称馅饼图，用一个圆来表现百分构成，百分构成DATA，TEXT，TEXT为字符变量。DATATEXT1、2、3、4。选Graphs菜单的Pie...过PieChart定义选项框构成图仅有一故直接点击Define15.5，dataSlicesRepresentOthersnmmaryfunctionVariabletext点击钮使之进入DefineSlicesby框。点击Titles...钮，弹出Titles框，在栏内输入“中草药治疗慢性支气管炎效果构成图”ContinueDefinePie:SummariesforGroupsofCases框，再点击OK钮即完成。图图15.5构成图绘 中草药治疗慢性支气管炎效果构成图中草药治疗慢性支气管炎效果构成图无效控制好转33显效区域，如一组测定值的范围（最小值—最大值、95%区间值（低限、X±1.96·SD（低

-

Bq·L湖GROUP，1、2、3。SimpleHigh-Low-Close：简单线型高低区域图；ClusteredHigh-Low-Close：复式线型高低区域图；SimpleRangeBar：简单直条型高低区域图；ClusteredRangeBar：复式直条型高低区域图；DifferenceLine：差异线区域图。Define钮，弹出Define15.6，点击BarsRepresentOthersnmmaryfunctionVariablecatCategoryAxisgroupDefineHigh-Low-Closeby点击Continue钮返回DefineSimpleHigh-Low-Close:SummariesforGroupsofCases框，OK图图15.6高低区域图绘 某市测某市测定不同水体放射性元素锶的含量图例高值低值均值自来 湖 水库水源点百分构成春夏秋冬DATA，输入具体数字；季节的变量GraphsPareto...ParetoChart2Simple为单一直条构成线图，Stacked为堆积式直条构成线图，本例选用单一直条构成线图。然后点击Define钮，弹出DefineSimplePareto:SummariesforGroupsofCases框（图15.7dataSumsofvariableseason钮使之进入CategoryAxis框。点击Titles...钮，弹出Titles框，在Title栏内输入“1000名20-25岁男性各季节感冒发生人数分析”ContinueDefineSimplePareto:SummariesforGroupsofCases框，再点击OK钮即完成。图图15.7直条构成线图绘 10%。1000名1000名20-25各季节感冒发生人数分00季节调用Graphs菜单的ontrol1234515.815.8X-Bar,R,s：均数控制图和极差（标准差）X图，用于控制重复测Rs本例选用X-Bar,R,s选项框的下方为数据类型选择（DataOrganizationCasesareunits表示数据文件中各观察样例只是一个值，其分组需要再定义；Casesaresubgroups表示数据文件15.9，选data钮使之进入ProcessMeasurementgroup钮使之进入SubgroupsDefinedbyChartsX-Barandrange图。点击Titles...钮，弹出Titles框，在Title栏内输入“样品测定的质量控制图”，点击Continue钮返回X-Bar,R,s:CasesAreUnits框，再点击OK钮即完成。图图15.9质量控制图定 样品测定的质量控样品测定的质量控制测定值均数上控制限=平均中心线=下控制限= Sigmalevel:样品测样品测定的质量控制极差值上控制限=平均中心线下控制限=. Sigmalevel:0[例15-8]研究甲基对肝脏脂质过氧化的毒性作用，选用25只大白鼠，随机分成五组，按不LPO（nmol/L，资料如下表，试绘制箱图。12345GROUP1、2、3、4、5。Graphs单的Boxplot...过程，弹出BoxplotChart定义选2图可选：Simple为简单箱图，Clustered为复式箱图，本例选用简单箱图。然后点击Define钮，弹出DefineSimple15.10，VariablegroupCategoryAxisOK图图15.10箱图定 不同染毒剂量组大白鼠LPO含量的0-20染毒剂量组（mg/kg体重DOSE(n不同染毒剂量组大白鼠LPO含量的0-20染毒剂量组（mg/kg体重DOSE(n调用Graphs菜单的ErrorBar过程，可绘制均值相关区间图。正态分布资料的描述性指标：如均值、标准差、标准误，并由此求得的参照值范围、总体均值的区间等，都可用均值相关区间15-9]食品中微量砷（As）DDC-Ag比较两种方法测定不同浓度As标准液（μg/5ml）的光密度值区间，试绘制均值相关区间图。DDC-Ag12345激活数据管理窗口，定义变量名：数据的变量名为DATA，将新银盐法和DDC-Ag法的测定所得CATE1、2。GraphsErrorBar...ErrorBar2然后点击Define钮，弹出DefineClusteredErrorBar:SummariesforGroupsofCases框15.11CategoryAxiscateDefineClustersby图图15.11均值相关区间图定 BarRepresent3Confidenceintervalformean：绘出总体均值的区间，要求输入区间的百分数；Standarderrorofmean：绘出均值+/-两倍标准误的区间；本例选用Confidenceintervalformean。之后点击Titles...钮，弹出Titles框，在Title栏内输入“两种方法测定微量砷的光密度95%区间比较”，点击Continue钮返回DefineClusteredErrorBar:SummariesforGroupsofCases框，再点击OK钮即完成。测定方法新银盐法-标准液As浓度两种方法测定微量砷的光密度 区间比95%CIGraphsScatter123CATE1，123；再设一变CATE2，用于定义微量元素种类，Zn1、Fe2、Mn3。15.12Simple为单层散点图，Overlay为多层散点图，Matrix为矩阵散点图，3-D为散点图，本例选用单层散点图。然后点击Define钮，弹出SimpleScatterplot框（图15.13，在左侧的变量dataYAxiscate1XAxis（指定饮茶状况为横轴标目）,cate2SetMarkersby（指定微量元素种类为散点标志。点击Titles...钮，弹出Titles框，在Title栏内输入“饮茶与微量元素的关系”，点击Continue钮返回SimpleScatterplot框，再点击OK钮即完成。15.12图图15.13散点图绘 饮茶与微量元素的关系0微量元素- 饮茶状况下图横轴标目中1.0为多喝茶组、2.0饮茶与微量元素的关系0微量元素- 饮茶状况调用Graphs菜单的Histogram过程，可绘制直方图。直方图是用直条的长短来表示连续性的绝10-15-20-25-30-35-745-355-1NUMBER，将每岁患者人数资料输入；设一变量为GROUP，用于定义组，将各分组的下限值输入。为使频数资料在作图中生效，应选Data菜单的WeightCases...命令项，在弹出的WeightCases框中取Weightcasesby项，NUMBERFrequencyVariableOK（图15.14，在左侧的变量列表中选group点击钮使之进入Variable框；点击Titles...钮，弹出Titles框，在Title栏内输入“某地某年流行性乙型脑炎患者分布”，点击Continue钮返回Histogram框，再点击OK钮即完成。图图15.14直方图绘 ChartCarousel（即每岁一组系统只按默认的每5岁一组方式输出图形，所以需要安用户的要求对统计图进行编辑。点击ChartCarousel窗口上端中的Edit钮，弹出!Chart1窗口，将鼠标箭头指向图的横轴并双击鼠标左15.15，IntervalsCustomDefine..IntervalAxis:DefineCustomInterval框（15.16DefinitionIntervalwidth1，要求按每岁一组的方式作图；在Range栏的Minimum处输入0在um处输入65要求横轴从0岁开始至65岁止点击Continue钮返回IntervalAxis框，再点击OK钮即可。图图15.15横轴修 15.16某地某年流行性乙型某地某年流行性乙型脑炎患者分布Std.Dev=4.06Mean=5.6 N=分GraphsNormalP-P[例15-12]某医师测得30名健康女大学生总蛋白含量（g/L）资料如下，试绘制正态概率选Graphs菜单的NormalP-P...过程，弹出NormalP-PPlot框（图15.17，在左侧的量列表中选data点击钮使之进入Variable框。在Transform栏中，系统