analysis生存分析ppt课件.ppt

乙方有义务提供相应的执业资格证、执业印章等配合甲方工程招投标、承担的工程项目报建及建设施工过程中各类手续办理。 生存分析的SAS编程操作 薛 富 波 2006/05/21 生存分析的有关概念 事件时间（time-to-event，TTE）：又称为生存时间（survival time）或失效时间（failure time），是指从研究开始到相应事件发生的时间。 删失值（censored value）：是指TTE在某些观测上其确切的取值是未知的，所能知道的仅仅是其测量值肯定大于或小于某个特定值（即其取值的下限或上限），或者位于某个特定的取值范围内（取值区间）。 右侧删失：反映真实值下限的数据。 左侧删失：反映真实值上限的数据。 区间删失：反映真实值所在区间的数据。 生存分析的方法分类 非参数法：乘积极限法、寿命表法等； 参数法：指数分布法、威布尔分布（Weibull Distribution）法、对数正态回归分析法以及对数logistic回归分析法等； 半参数法：Cox模型分析方法。 生存分析方法的SAS过程 非参数法：lifetest过程； 参数法：lifereg过程； 半参数法：phreg过程。 生存分析的非参数方法 非参数法生存分析示例1 45例乳腺癌患者中，免疫过氧化物酶检测结果显示9例阳性、36例阴性，比较其生存时间分布的SAS程序如下。 proc lifetest plots=(s); time time*censor(1); strata immuno; run; lifetest过程可包含的语句 proc lifetest options ; time variable *censor(list) ; by variables ; freq variable ; id variables ; strata variable (list) ... variable (list) ; survival options ; test variables ; run ; proc lifetest语句 proc lifetest语句为调用lifetest过程的开始，其主要语句选项及其可执行的功能见下表。 proc lifetest语句的主要选项 选 项 功能和用法 method= 设置用以计算生存函数估计值的方法。此选项可设置的值及其含义如下：“pl”（或“km”）——乘积极限法（或Kaplan-Meier法）；“act”（或“life”或“lt”）——寿命表法。默认设置为“method=pl”。 intervals= 指定用以计算寿命表的区间端点。区间端点必须为非负数，无论指定的区间端点中有无0，寿命表的计算总是以0为第一个区间的起点，而每一个区间只具有下端点而无上端点（即为半开半闭区间）。 ninterval= 指定寿命表区间的个数，此选项可被“width=”和“intervals=”所覆盖。当设置“ninterval=”选项时，lifetest过程将根据所设置的区间个数划分寿命表的区间，但区间端点将会作适当的调整，使其以整数的形式出现。因此，最终划分的区间数不一定与设置的区间数完全相符。默认设置为“ninterval=10”。 width= 以寿命表法计算生存函数估计值时，指定寿命表的区间宽度。此选项可被“intervals=”所覆盖。 plots= 为生存函数估计值或删失值绘制图形。其完整设置形式应为“PLOTS= ( type (NAME=name) , ..., type (NAME=name) )”，其中“type”代表绘图的类型，“NAME=name”表示将所绘图形以指定的名称存储为SAS目录的条目（entry），其中“NAME”为关键字（可为小写），“name”为用户指定的条目名称。关于绘图类型详见“plots=”选项的绘图类型。 nocensplot 当设置了“plots=”选项时，禁止将删失值显示在图形中。当选择寿命表法进行计算时，无需设置此选项。 timelist= 指定一系列时间点，从而在结果中针对这些时间点显示相应的Kaplan-Meier估计值，该时间点在输出结果中所在的列以“_TIME_”为标识。 alphaqt= 指定生存时间四分位数间距可信区间的置信水平（须在0.0001~0.9999之间），默认设置为“alpha=0.05”。 intervals=选项的设置方法 端点值间以空格分隔，如“intervals=10 30 50 70”； 端点之间以逗号分隔，如“intervals=10,30,50,70”； “a to b”的形式，如“intervals=1 to 7”； “a to b by c”的形式，如“inter